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%	31.6 / NA	13	L36Q (CTG→CAG)	folD	bifunctional methylenetetrahydrofolate dehydrogenase/methenyltetrahydrofolate cyclohydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (5/8);  total (5/8)

GTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCAGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGA  >  NC_000913/557438‑557841                                                                                                                                                                                                      |                                                                                                                                                                                                               gTAACGGCAGTTGAACCGGAATGCCGTCGTTTTTGTTGTCGGCCTTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTCGTTTCCGGGAGGTCATAATCGCGGGAGACCAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCTTCATAAATTTCCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTcc                                                                                                                                                                                                            <  1:75698/202‑1 (MQ=255)          gTTGAACCAGAATGCCATCGATGGTGTTGGCGGCATTCAGCGGATCGATAAGCTCCAGCAGCGCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACTAACCCGACTTCTGCACAAGCCTTGCGTTTGCTTGCGACATAAATTTTGGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTc                                                                                                                                                                                                             <  1:83982/192‑1 (MQ=255)                  aGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTGTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACctgctg                                                                                                                                          <  1:109493/251‑1 (MQ=255)                                                            aGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATCAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCACTAATCTTTGCTGCCCTCag                                                                                                   >  1:88356/1‑248 (MQ=255)                                                                    cagcTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCACTACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTc                                                                                        >  1:84049/1‑251 (MQ=255)                                                                     agcgccGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAGCTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTcn                                                                                       <  1:104786/247‑2 (MQ=255)                                                                                                           ncGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGt                                                 >  1:29685/2‑251 (MQ=255)                                                                                                                          naCTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGgcgcg                                  >  1:40712/2‑251 (MQ=255)                                                                                                                                                       cGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCGGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCg                                         <  1:156188/215‑1 (MQ=255)                                                                                                                                                           aTAAATTTGCGATGCAGGGTTACTACNCACCAGNNCANCGGCCTGTCCTGGTGCCCGNAGNNNGNNNGNCATACNCGNNTGCACTTTTNGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGtttt                                                                                                                              >  2:36341/1‑126 (MQ=255)                                                                                                                                                           aTAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCTCaagagtttt                                                                                                                              <  1:36341/126‑8 (MQ=255)                                                                                                                                                                             gTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTg   <  1:145984/231‑1 (MQ=255)                                                                                                                                                                                    cccACCAGCACAACGTCCTGTCCTGGTGCTCGCAGTCCGGCTGCAATACTCGCCTGAACTTTTTGAGCAACTTCAGGGCGCACCTTCTGCGCAATCGTTTTACCGTCAATAAGCTTTGCGGCCATCATAGAGAGGATTCTATCTATTACGTAGATCGGAGGTGTTGCGCCTTTATTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGa  <  1:167592/225‑1 (MQ=255)                                                                                                                                                                                                      |                                                                                                                                                                                                               GTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCAGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGA  >  NC_000913/557438‑557841

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