Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	557,635	A→T	L36Q (CTG→CAG)	folD ←	bifunctional methylenetetrahydrofolate dehydrogenase/methenyltetrahydrofolate cyclohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	557,635	0	A	T	100.0%	36.3 / NA	14	L36Q (CTG→CAG)	folD	bifunctional methylenetetrahydrofolate dehydrogenase/methenyltetrahydrofolate cyclohydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (7/7);  total (7/7)

GACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCAGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGT  >  NC_000913/557418‑557823                                                                                                                                                                                                                          |                                                                                                                                                                                             gACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATAcgcg                                                                                                                                                             >  1:111602/1‑251 (MQ=255)     ttactaatccccGACGGTAACGGCAGTTGATCCAAAATGCCATCGATGGTTTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTGCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCt                                                                                                                                                          <  2:167575/241‑1 (MQ=255)      tatatatcccctccgGTAACGTCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCCGCATCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGTGAGACGATCCCGACTTCTTCACCAGCCTTGCGTTTGCTTGCTACATAAATTTGCGATGCAGGGTTACTCCCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTg                                                                                                                                                         <  1:1326/238‑1 (MQ=255)             aCCCGCCGGTAACGGCAGTTGAACCAGACTGCCATCGATGGTGTTTTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTAGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCCACATAAATTTGCGATGCAGGGTTACTACCCACCAGAACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACttttt                                                                                                                                                 >  2:162715/1‑251 (MQ=255)                     gTAACGGCAGTTGAACCAGAATTCCATCGATGGTGTTGTCGGAATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGCGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAAACCGACTTCTTCACAAACCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCCCCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCACTAAGCACCTTAACCTTTTGAGCaaat                                                                                                                                         >  2:165335/1‑249 (MQ=255)                         cGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATA‑GCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTc                                                                                                                                       >  2:21353/1‑248 (MQ=255)                          ggCAGTTGAACCAGAATGCAATCGATGGAGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCAACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCATGAACTTTTTGAGCAACTTc                                                                                                                                       >  2:9327/1‑248 (MQ=255)                                                cGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGTGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGt                                                                                                               <  2:138497/250‑1 (MQ=255)                                                                                aGCTCCCGCAGCTGCGCTTCGCCGGTGGTTTCCGGGCGGTCCTAAGAGCGGGAGAAGCACCCGACTTCTTCCCAAGCATTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCCCCTGATGCGCAATCGTTTTACCGCCACTAATCTTTGCTGCCATCagtgg                                                                              >  2:95933/1‑248 (MQ=255)                                                                                             ccGCTTCGCGTGTGGTGTCCGGGAGGTCATAAGATCGGGAGACGAACCCGACTTCTTCACAAGCCTTCCGTTTGCTTGCGACATAAATTTGCTATGCAGGGTTACTACCCTCCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCACCTGACCTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGTATTCCATCt                                                                 <  2:127640/251‑1 (MQ=255)                                                                                                                                   gagaCGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGTTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACACGCCTGCACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGCTCGAAGGGGATGCGCCTATTTTGTc                              >  2:168174/1‑248 (MQ=255)                                                                                                                                   gagaCGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGa                           <  2:111602/251‑1 (MQ=255)                                                                                                                                           cccGACTTCTTCACAAGCCTTGCGTTTTCTTGCGACATAAATTTGCGATGCAGTTTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGATCTTTTTGAGCAACTTCAGATCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTTCCATCAGAGAGATGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGgc                   <  2:73805/251‑1 (MQ=255)                                                                                                                                                            ccTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGgtttcgt  <  1:63014/251‑1 (MQ=255)                                                                                                                                                                                                                          |                                                                                                                                                                                             GACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCAGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGT  >  NC_000913/557418‑557823

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