Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,630,852	G→A	S176N (AGC→AAC)	yfgI →	nalidixic acid resistance protein YfgI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,630,852	0	G	A	100.0%	87.0 / NA	25	S176N (AGC→AAC)	yfgI	nalidixic acid resistance protein YfgI
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (7/18);  total (7/18)

GACAAAAAACGCGATCAAAGTTATGAAGATGAGCTACGCAGCTTAGAGATTCAGAAACAAAAACTGGCGCTGGCGAAAGAAGAAGCCCGCGTTAAGCGAGAAAACGAATTTATCGATCAGGAACTGAAGCACAAAGCTGCGCAAACCGATGTGGTGCAATCTGAAGCTGACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAATCGTACGGCAGATGCGCCCAGTGTGCGGTC  >  NC_000913/2630605‑2631093                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  gACAAAAAACGCGATCAAAGTTATGAAGATGAGCTACGCAGCTTAGAGATTCAGAAACAAAAACTGGCGCTGGCGAAAGAAGAAGCCCGCGTTAAGCGAGAAAACGAATTTATCGATCAGGAACTGAAGCACAAAGCTGCGCAAACCGATGTGGTGCAATCTGAAGCTGACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAACTg                                                                                                                                                                                                                                                <  1:110223/251‑1 (MQ=255)                        tGAAGATGAGCTACGCAGCTTAGAGATTCAGAAACAAAAACTGGCGCTGGCGAAAGAAGAAGCCCGCGTTAAGCGAGAAAACGAATTTATCGATCAGGAACTGAAGCACAAAGCTGCGCAAACCGATGTGGTGCAATCTGAAGCTGACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAACTGGTTTAATTAATCGATGTTAGTa                                                                                                                                                                                                                          >  2:390422/1‑250 (MQ=255)                            gATGAGCTACGCAGCTTAGAGATTCAGAAACAAAAACTGGCGCTGGCGAAAGAAGAAGCCCGCGTTAAGCGAGAAAACGAATTTATCGATCAGGAACTGAAGCACAAAGCTGCGCAAACCGATGTGGTGAAATCTGAAGCTGACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGAGGGCAAAGCAGAAGAGAACAAATCGGACAACTGGTTTAATTAATCGATGTTAGTAACTTc                                                                                                                                                                                                                     >  2:205701/1‑251 (MQ=255)                                aGCTACGCAGCTTAGAGATTCAGAAACAAAAACTGGCGCTGGCGAAAGAAGAAGCCCGCGTTAAGCGAGAAAACGAATTTATCGATCAGGAACTGAAGCACAAAGCTGCGCAAACCGATGTGGTGCAATCTGAAGCTGACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAACTGGTTTAATTAATCGATGTTAGTAACTTCAATc                                                                                                                                                                                                                 <  1:34685/251‑1 (MQ=255)                                             agagATTCAGAAACAAAAACTGGCGCTGGCGAAAGAAGAAGCCCGCGTTAAGCGAGAAAACGAATTTATCGATCAGGAACTGAAGCACAAAGCTGCGCAAACCGATGTGGTGCAATCTGAAGCTGACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAACTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGa                                                                                                                                                                                                    <  2:113861/251‑1 (MQ=255)                                                       aaaCAAAAACTGGCGCTGGCGAAAGAAGAAGCCCGCGTTAAGCGAGAAAACGAATTTATCGATCAGGAACTGAAGCACAAAGCTGCGCAAACCGATGTGGTGCAATCTGAAGCTGACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAACTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATcaaaaa                                                                                                                                                                                          <  2:343871/251‑1 (MQ=255)                                                           aaaaaCTGGCGCTGGCGAAAGAAGAAGCCCGCGTTAAGCGAGAAAACGAATTTATCGATCAGGAACTGAAGCACAAAGCTGCGCAAACCGATGTGGTGCAATCTGAAGCTGACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAACTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCct                                                                                                                                                                                      <  1:44482/251‑1 (MQ=255)                                                                    cgctggcgAAAGAAGAAGCCCGCGTTAAGCGAGAAAACGAATTTATCGATCAGGAACTGAAGCACAAAGCTGCGCAAACCGATGTGGTGCAATCTGAAGCTGACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAACTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTaa                                                                                                                                                                             <  1:390422/251‑1 (MQ=255)                                                                        ggcgAAAGAAGAAGCCCGCGTTAAGCGAGAAAACGAATTGATCGATCAGGAACTGAAGCACAAAGCTGCGCAAACCGATGTGGTGCAATCTGAAGCTGACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAATAGAACAAATCGGACAACTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACaga                                                                                                                                                                         <  1:67178/251‑1 (MQ=255)                                                                               gaagaagCCCGCGTTAAGCGAGATAACGAATTTATCGATCAGGAACTGAAGCCCAAAGCTGCGCAAACCGATTTGGTGCAATCTGAAGCTGACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGTCAAAGCAGAAGAGAACATATCTGACAACTGGTTTAATTAATCGATGTTAGTAACTTCTATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGtttt                                                                                                                                                                  <  1:274528/251‑1 (MQ=255)                                                                                      cccGCGTTAAGCGAGAAAACGAATTTATCGATCAGGAACTGAAGCACAAAGCTGCGCAAACCGATGTGGTGCAATCTGAAGCTGACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAACTGGTTTAATTAATCGATGTTAGTAACTGCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGttcatt                                                                                                                                                           <  2:136349/251‑1 (MQ=255)                                                                                                                aTCGATCAGGAACTGAAGCACAAAGCTGCGCAAACCGATGTGGTGCAATCTGAAGCTGACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTTGGCAAAGCAGAAGAGAACAAATCGGACAACTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCa                                                                                                                                 <  1:2510/251‑1 (MQ=255)                                                                                                                                    caAAGCTGCGCAAACCGATGTGGGGCAATCTGAAGCTGACGCCAACAGAAATATGACGGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAACTGGTTTAATTAATCGATGTTAGTAACTTCAAGCCTATAATTCTTGAAGATAAAAAACCCTCTTTAGTAACAGAGGGTTTTGGTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCt                                                                                                             <  1:389928/251‑1 (MQ=255)                                                                                                                                                  ccGATGTGGTGCAATCTGAAGCTGACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAACTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCa                                                                                               <  1:393136/251‑1 (MQ=255)                                                                                                                                                  ccGATGTGGTGCAATCTGAAGCTGACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAACTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACGCAATGGTAGCTGGCGGCTTGCCGCTGATGTCa                                                                                               <  2:209608/251‑1 (MQ=255)                                                                                                                                                            gCAATCTGAAGCTGACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAACTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCTCTGATGTCATACACCAcg                                                                                      <  1:113221/250‑1 (MQ=255)                                                                                                                                                                                aCAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAACTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTc                                                                  >  1:290077/1‑250 (MQ=255)                                                                                                                                                                                     aaTATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAACTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATGCCGTTCACTTCATTGAt                                                            <  2:180055/251‑1 (MQ=255)                                                                                                                                                                                         tGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAACTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATg                                                        >  1:160501/1‑251 (MQ=255)                                                                                                                                                                                                 ggCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAACTGGTTTAATTAAGCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGTGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGa                                                <  1:218695/251‑1 (MQ=255)                                                                                                                                                                                                                   aaaaGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAACTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAATc                              <  1:108134/251‑1 (MQ=255)                                                                                                                                                                                                                    aaaGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAACTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAATCg                             >  2:166679/1‑251 (MQ=255)                                                                                                                                                                                                                           gggCAAAGCAGAAGAGAACAAATCGGACAACTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAATCGTACGGc                       <  2:352827/250‑1 (MQ=255)                                                                                                                                                                                                                           gggCAAAGCAGAAGAGAACAAATCGGACAACTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTATTGACGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAATCGTACGGc                       >  1:432643/1‑250 (MQ=255)                                                                                                                                                                                                                                                aTCGGACAACTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAATCGTACGGCAGATGCGCCCAGTGTGCGGTc  >  2:152336/1‑250 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  GACAAAAAACGCGATCAAAGTTATGAAGATGAGCTACGCAGCTTAGAGATTCAGAAACAAAAACTGGCGCTGGCGAAAGAAGAAGCCCGCGTTAAGCGAGAAAACGAATTTATCGATCAGGAACTGAAGCACAAAGCTGCGCAAACCGATGTGGTGCAATCTGAAGCTGACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAATCGTACGGCAGATGCGCCCAGTGTGCGGTC  >  NC_000913/2630605‑2631093

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