Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_010473	142,348	N→A	?254* (TNG→TAG)	fhuA →	ferrichrome porin FhuA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_010473	142,348	0	N	A	93.8%	95.2 / ‑6.4	33	?254* (TNG→TAG)	fhuA	ferrichrome porin FhuA
Reads supporting (aligned to +/- strand):  ref base N (0/0);  major base A (14/16);  minor base C (2/0);  total (17/16)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTT‑‑CCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGTTCTGCCAATGCCCAGCAGAAAGGGTCAGAA‑GAGCAGCGTTATGCTATTGCACCGGCGTTCACCTNGCGTCC‑GGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCC‑GGAAACCGG‑TTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCC‑GCTGCCGAACGGTAAGCGTCTGCCGACAGACTTTAATGAAGGGGCGAAGAACAACACCTATTCTCGTAATGAGAAGATGGTCGGCTACAGCTTCGAT  >  NC_010473/142157‑142547                                                                                                                                                                                                   |                                                                                                                                                                                                            caccGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGTTCTGCCAATGCCCAGCAGAAAGGGTCAGAA‑GAGCAGCGTTATGCTATTGCACCGGCGTTCACCTAGCGTcc                                                                                                                                                                                                       >  1:21391S12/1‑200 (MQ=255)  accGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTT‑‑CCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGTTCTGCCAATGCCCAGCAGAAAGGGTCAGAA‑GAGCAGCGTTATGCTATTGCACCGGCGTTCACCTAGCGTCC‑Gga                                                                                                                                                                                                   <  1:16722S32/200‑1 (MQ=255)          gCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTT‑‑CCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGTTCTGCCAATGCCCAGCAGAAAGGGTTA‑AA‑AAGCAGCGTTATGCTATTGCAC‑GGCATTAGCCTCACGTCC‑GGATGATAAAAcc                                                                                                                                                                                         >  1:29723S105/1‑200 (MQ=255)                                 ccGGTACTGACAGCCTGTT‑‑CCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGTTCTGCCAATGCCCAGCAGAAAGGGTCAGAAAGAGCAGCGTTATGCTATTGCACCGGCGTTCACCTAGCGTCC‑GGATGATAAAACCAATTTTACCTTCCTTTCTTa                                                                                                                                                                     <  1:2345S49/200‑1 (MQ=255)                                 ccGGTACTGACAGCCTGTT‑‑CCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGTTCTGCCAATGCCCAGCAGAAAGGGTCAGAA‑GAGCAGCGTTATGCTATTGCACCGGCGTTCACCTAGCGTCC‑GGATGATAAAACCAATTTTACCTTCCTTTCTTAc                                                                                                                                                                    >  1:26280S5/1‑200 (MQ=255)                                       tcgACAGCCTGTT‑‑CCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTAC‑‑‑‑‑TGGCCTGAC‑GGTCTTGCGCGTTCTGCCAATGCCCAGCAGAAAGGGTCAGAA‑GAGCAGCGTTATGCTATTGCACCGGCGTTCACCTAGCGTCC‑GGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAg                                                                                                                                                        >  1:18042S25/3‑200 (MQ=255)                                        tGACAGCCTGTT‑‑CCAGACTGGTTTTGACTTTAGCGATCTGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGTTCTGCCAATGCCCAGCAGAAAGGGTCAGAA‑GAACAGCGTTATGCTATTGCACCGGCGTTCACCTAGCGTCC‑GGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGa                                                                                                                                                             <  1:27811S86/200‑1 (MQ=255)                                                      ccAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCT‑ATCGCCTGACCGGTCTTGCGCGTTCTGCCAATGCCCAGCAGAAAGGGTCAGAA‑GAGCAGCGTTATGCTATTGCACCGGCGTTCACCTAGCGTCC‑GGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCC‑GGAAAc                                                                                                                                               <  1:27047S2/200‑1 (MQ=255)                                                             ggTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGTTCTGCCAATGCCCAGCAGAAAGGGTCAGAA‑GAGCAGCGTTATGCTATTGCACCGGCGTTCACCGAGCGTCC‑GGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCC‑GGAAACCGG‑tta                                                                                                                                        <  1:1415S7/200‑1 (MQ=255)                                                                 ttGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGTTCTGCCAATGCCCAGCAGAAAGGGTCAGAA‑GAGCAGCGTTATGCTATTGCACCGGCGTTCACCTAGCGTCC‑GGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCC‑GGAAACCGG‑TTATTAc                                                                                                                                    >  1:5649S2/1‑200 (MQ=255)                                                                    aCTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGTTCTGCCAATGCCCAGCAGAAAGGGTCGAAA‑GAGCAGCGTTATGCTATTGCACCGGCGTTCACCTAGCGTCC‑GGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCC‑GGAAACCGG‑ttctgcggct                                                                                                                                 >  1:22434S127/1‑192 (MQ=255)                                                                                 gTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGTTCTGCCAATGCCCAGCAGAAAGGGTCAGAA‑GAGCAGCGTTATGCTATTGCACCGGCGTTCACCTAGCGTCC‑GGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCC‑GGAAACCGG‑TTATTACGGCTGGTTGCCGAAag                                                                                                                    <  1:139S38/200‑1 (MQ=255)                                                                                           gaCGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGTTCTGCCAATGCCCAGCAGAAAGGGTCAGAA‑GAGCAGCGTTATGCTATTGCACCGGCGTTCACCTAGCGTCC‑GGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCC‑GGAAACCGG‑TTATTACGGCTGGTTGCCGAAAGAggtaccgtt                                                                                                          <  1:25484S32/200‑8 (MQ=255)                                                                                             cGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGTTCTGCCAATGCCCAGCAGAAAGGGTCAGAA‑GAGCAGCGTTATGCTATTGCACCGGCGTTCACCTAGCGTCC‑GGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCC‑GGAACCCGGGTTATTACGGCTGGTTGCCGAAAGAGGGAACCGtt                                                                                                         >  1:28558S124/1‑200 (MQ=255)                                                                                                  tttACTCTTATCGCCTGACCGGTCTTGCGCGTTCTGCCAATGCCCAGCAGAAAGGGTCAGAA‑GAGCAGCGTTATGCTATTGCACCGGCGTTCACCTAGCGTCC‑GGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCC‑GGAAACCGG‑TTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCC‑g                                                                                                  >  1:2358S15/1‑200 (MQ=255)                                                                                                   ttACTCTTATCGCCTGACCGGTCTTGCGCGTTCTGCCAATGCCCAGCAGAAAGGGTCAGAA‑GAGCAGCGTTATGCTATTGCACCGGCGTTCACCTAGCGTCC‑GGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCC‑GGAAACCGG‑TTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCC‑Gc                                                                                                 >  1:29680S5/1‑200 (MQ=255)                                                                                                         ttATCGCCTGACCGGTCTTGCGCGTTCTGCCAATGCCCAGCAGAAAGGGTCAGAA‑GAGCAGCGTTATGCTATTGCACCGGCGTTCACCTAGCGTCC‑GGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCC‑GGAAACCGG‑TTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCC‑GCTGCCGa                                                                                           <  1:20162S18/200‑1 (MQ=255)                                                                                                              gCCTGACCGGTCTTGCGCGTTCTGCCAATGCCCAGCAGAAAGGGTCAGAA‑GAGCAGCGTTATGCTATTGCACCGGCGTTCACCTAGCGTCC‑GGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCC‑GGAAACCGG‑TTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCC‑GCTGCCGAACGGt                                                                                      >  1:5859S26/1‑200 (MQ=255)                                                                                                                     cGGTCTTGCGCGTTCTGCCAATGCCCAGCAGAAAGGGTCAGAA‑GAGCAGCGTTATGCTATTGCACCGGCGTTCACCTAGCGTCC‑GGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCC‑GGAAACCGG‑TTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCC‑GCTGCCGAACGGTAAGCGTc                                                                               <  1:5822S24/200‑1 (MQ=255)                                                                                                                          ttGCGCGTTCTGCCAATGCCCAGCAGAAAGGGTCAGAA‑GAGCAGCGTTATGCTATTGCGCCGGCATTAGCCTCACGTCCCGGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCC‑GGAAACCGG‑TTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCC‑GCTGCCGAACGGTAAGCGTCTGcc                                                                           >  1:14972S11/1‑200 (MQ=255)                                                                                                                             cgcgTTCTGCCAATGCCCAGCAGAAAGGGTCAGAA‑GAGCAGCGTTATGCTATTGCACCGGCGTTCACCTAGCGTCC‑GGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCC‑GGAAACCGG‑TTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCC‑GCTGCCGAACGGTAAGCGTCTGCCGACa                                                                       <  1:8378S87/200‑1 (MQ=255)                                                                                                                                gTTCTGCCAATGCCCAGCAGAAAGGGTCAGAA‑GAGCAGCGTTATGCTATTGCACCGGCGTTCACCTAGCGTCC‑GGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCC‑GGAAACCGG‑TTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCC‑GCTGCCGAACGGTAAGCGTCTGCCGACAGAc                                                                    <  1:9493S28/200‑1 (MQ=255)                                                                                                                                gTTCTGCCAATGCCCAGCAGAAAGGGTCAGAA‑GAGCAGCGTTATGCTATTGCACCGGCGTTCACCTAGCGTCC‑GGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCC‑GGAAACCGG‑TTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCC‑GCTGCCGAACGGTAAGCGTCTGCCGACAGAc                                                                    >  1:17833S82/1‑200 (MQ=255)                                                                                                                                       cAATGCCCAGCAGAAGGGGTCAGAA‑GAGCAGCGTTATGCTATTGCACCGGCGTTCACCTAGCGTCC‑GGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCC‑GGAAACCGG‑TTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCC‑GCTGCCGAACGGTAAGCGTCTGCCGACAGACTTTAATg                                                             <  1:4388S9/200‑1 (MQ=255)                                                                                                                                           gCCCAGCAGAAAGGGTCAGAA‑GAGCAGCGTTATGCTATTGCACCGGCGTTCACCTAGCGTCC‑GGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCC‑GGAAACCGG‑TTATTTAGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCC‑GCTGCCGAACGGTAAGCGTCTGCCGACAGACTTTAATGAAgg                                                         >  1:12397S56/1‑200 (MQ=255)                                                                                                                                            cccAGCAGAAAGGGTCAGAA‑GAGCAGCGTTATGCTGT‑‑CACCGGCGTTCGCCTGGCGTCC‑GGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCC‑GGAAACCGG‑TTATTACGGCTGGTTGCC‑AAAGAGGGAACCGTTGAGCC‑GCTGCCGAACGGTAAGCGTCTGCCGACAGACTTTAATGAAGGGGCg                                                     >  1:25665S21/1‑200 (MQ=255)                                                                                                                                                                     agcGTTATGCTATTGCACCGGCGTTCACCTAGCGTCC‑GGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCC‑GGAAACCGG‑TTATTACGGCTGGTTGCCGAAAGAGGGAACCGTT‑CGC‑‑GCTGCCGAACGGTAAGCGTCTGCCGACAGACTTTAATGAAGGGGCGAAGAACAACACCTATTCTCGTaa                              >  1:10334S30/1‑200 (MQ=255)                                                                                                                                                                            tGCTATTGCACCGGCGTTCACCTAGCGTCC‑GGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACAAGCCAGAAAACCGG‑TTACTGCGGCTGGTTGCCGAA‑G‑GGGAACCGTTGAGCC‑GCTGCCGAACGGTAAGCGTCTGCCGACAGACTTTAATGAAGGGGCGAAGAACAACACCTATTCTCGTAATGagaa                        >  1:13174S4/1‑200 (MQ=255)                                                                                                                                                                              cTATTGCACCGGCGTTCACCTAGCGTCC‑GGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCC‑GGAAACCGG‑TTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCC‑GCTGCCGAACGGTAAGCGTCTGCCGACAGACTTTAATGAAGGGGCGAAGAACAACACC‑AT‑‑TCGTAATGAGAAGATg                    <  1:18986S15/200‑1 (MQ=255)                                                                                                                                                                                 ttGCACCGGCGTTCACCTAGCGTCC‑GGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCC‑GGAAACCGG‑TTATTACGGCTGGTTGC‑AAAAGAGG‑ATCTGGTCAGTTTGCTGCCGAACGGTAAGCGTCTGCCGACAGACTTTAATGAAGGGGCGAAGAACAACACCTATTCTCGTAATGAGAAGATgg                   >  1:8048S39/1‑200 (MQ=255)                                                                                                                                                                                         gCGTTCACCTAGCGTCC‑GGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCC‑GGAAACCGG‑TTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCC‑GCTGCCGAACGGTAAGCGTCTGCCGACAGACTTTAATGAAGGGGCGAAGAACAACACCTATTCTCGTAATGAGAAGATGGTCGGCTa            <  1:24769S5/200‑1 (MQ=255)                                                                                                                                                                                          cGTTCACCTAGCGTCC‑GGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCC‑GGAAACCGG‑TTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCC‑GCTGCCGAACGGTAAGCGTCTGCCGACAGACTTTAATGAAGGGGCGAAGAACAACACCTATTCTCGTAATGAGAAGATGGTCGGCTAc           <  1:20677S9/200‑1 (MQ=255)                                                                                                                                                                                                ccTAGCGTCC‑GGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCC‑GGAAACCGG‑TTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCC‑GCTGCCGAACGGTAAGCGTCTGCCGACAGACTTTAATGAAGGGGCGAAGAACAACACCTA‑‑CTC‑CAATGAGAAGATGGTCGGCTACAGCTTCGAt  <  1:16576S36/200‑1 (MQ=255)                                                                                                                                                                                                   |                                                                                                                                                                                                            CACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTT‑‑CCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGTTCTGCCAATGCCCAGCAGAAAGGGTCAGAA‑GAGCAGCGTTATGCTATTGCACCGGCGTTCACCTNGCGTCC‑GGATGATAAAACCAATTTTACCTTCCTTTCTTACTTCCAGAACGAGCC‑GGAAACCGG‑TTATTACGGCTGGTTGCCGAAAGAGGGAACCGTTGAGCC‑GCTGCCGAACGGTAAGCGTCTGCCGACAGACTTTAATGAAGGGGCGAAGAACAACACCTATTCTCGTAATGAGAAGATGGTCGGCTACAGCTTCGAT  >  NC_010473/142157‑142547

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 50 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑