Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	169,074	G→T	E531* (GAA→TAA)	fhuA →	ferrichrome outer membrane transporter/phage receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	169,074	0	G	T	100.0%	98.0 / NA	29	E531* (GAA→TAA)	fhuA	ferrichrome outer membrane transporter/phage receptor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (15/14);  total (15/14)

GGCAAACTCCGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCT  >  NC_000913/168833‑169315                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  ggCAAACTCCGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTAAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTg                                                                                                                                                                                                                                          <  1:154348/251‑1 (MQ=255)  gCAAACTCCGGCCCTTACCGCATTCTGAATAAAAAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTg                                                                                                                                                                                                                                          >  1:270414/1‑250 (MQ=255)          cGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACcttctt                                                                                                                                                                                                                                 <  1:111038/251‑1 (MQ=255)                    cGCATTCTGAATAAACAGAAACAAACGGACGTTTATGTTCAGGATCAGGCGCAGAGGGATAAAGAGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCCGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGCCAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTGCGCAAGTTgg                                                                                                                                                                                                                       >  2:273493/1‑251 (MQ=255)                     gCATTCTGAATAAACAGAAACAAACGGGCGTTTATGGTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAAGCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGTCGTTGTGGTGTTAACTACCTGTTGGATAATGGTGTAACACCTTACTTCAGCTAGAGCTAATCGGTTGAACCTTCTTCGctttttggg                                                                                                                                                                                                                      <  2:426829/251‑9 (MQ=255)                            gAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGAt                                                                                                                                                                                                                <  2:25109/250‑1 (MQ=255)                                 aaCAGAAACAAACGGGCGTTTATGTTCAGGAACAGGCGCAGTGGGATAAAGTGCTGGTCCCACTAGGCGGTCGTGATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACTGATAAACGTGATGACAAACAGTTTACCTGTCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACAACTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAat                                                                                                                                                                                                          <  1:273493/251‑1 (MQ=255)                                  aCAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAat                                                                                                                                                                                                          >  2:309205/1‑250 (MQ=255)                                                           cAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGc                                                                                                                                                                                 >  2:12861/1‑250 (MQ=255)                                                                      gcAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCg                                                                                                                                                                     <  1:340601/251‑1 (MQ=255)                                                                             gATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTTTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGa                                                                                                                                                               >  2:17449/1‑250 (MQ=255)                                                                                                  ggCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGt                                                                                                                                                                                                                          >  1:19709/1‑170 (MQ=255)                                                                                                  ggCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTGCGCAAGt                                                                                                                                                                                                                          <  2:19709/170‑1 (MQ=255)                                                                                                              gACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTAc                                                                                                                             <  1:162685/251‑1 (MQ=255)                                                                                                                          cAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAATTCGGCGGGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCgtgt                                                                                                                  <  1:305107/250‑1 (MQ=255)                                                                                                                                tctcTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAAtct                                                                                                           >  1:345901/1‑251 (MQ=255)                                                                                                                                  tctTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAAtctc                                                                                                          >  1:212201/1‑250 (MQ=255)                                                                                                                                     tAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCAc                                                                                                        <  2:341129/249‑1 (MQ=255)                                                                                                                                         cgcgTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAAGCTCACTAAAAc                                                                                                  <  1:276887/251‑1 (MQ=255)                                                                                                                                              tGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACcaac                                                                                              >  2:348028/1‑250 (MQ=255)                                                                                                                                                       gacCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATgg                                                                                    >  1:306677/1‑251 (MQ=255)                                                                                                                                                                                                         aCCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCcgcg                                    <  1:30348/249‑1 (MQ=255)                                                                                                                                                                                                               ttGATAAAGGTGTAACACCTTCCCTCAGCTATAGCTCATCGTTTGAACCTTCTTCGCCAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTACGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGCTTGTAGTTACTGGTGCCGTGTATAATCTCACTAAACCCAACAACCTGATGGCGGACCCTGAGGGTTCATTCTTCTCGGTTGCAGGTGGCGAGATCCGCGCACGTGGc                            >  2:457454/1‑251 (MQ=255)                                                                                                                                                                                                                  ataatggtgtcaccttttACTTCAGCTATAGCTAAGTGTTTGACCTTTCCTCGCATGTTGGGAAGGATTGTAATATTCTCGCCCTGTCTAAAGGTAAGCAGTATGAAGGCGGTGTGAAATATGTCCCGGAAGATCGTCCTATTGTAGTTTCTGGTGCTGTGTATAATCTCACTAAACCCAACACCCTGATCGCGCCCTCTGAGTGTCCCTTC‑‑CTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAg                        <  2:394888/234‑1 (MQ=255)                                                                                                                                                                                                                         tgtAACACCTTACTTAAGCTATAGCTAATCGTTTGACCCTTCTTCGCAAGTTGGGACGGATGGTAATATTTTCGCACCGTCTACAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCg                  >  2:231184/1‑251 (MQ=255)                                                                                                                                                                                                                             acacCTTACTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTCCCTTCTTCTCTGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGc              <  1:166263/251‑1 (MQ=255)                                                                                                                                                                                                                                    aCTTCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAgcg       >  1:145757/1‑251 (MQ=255)                                                                                                                                                                                                                                      ttCAGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAgcggc     >  1:207731/1‑251 (MQ=255)                                                                                                                                                                                                                                         aGCTATAGCTAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCt  >  1:226696/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  GGCAAACTCCGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTACCTGTTTGATAATGGTGTAACACCTTACTTCAGCTATAGCGAATCGTTTGAACCTTCTTCGCAAGTTGGGAAGGATGGTAATATTTTCGCACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCT  >  NC_000913/168833‑169315

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