Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,860,630	(C)5→4	coding (348/894 nt)	yidL →	putative DNA‑binding transcriptional regulator YidL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,860,626	0	C	.	100.0%	78.0 / NA	19	coding (344/894 nt)	yidL	putative DNA‑binding transcriptional regulator YidL
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (10/9);  total (10/9)

AGGATGTGGGAAAATTCTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTACCCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCACC  >  NC_000913/3860489‑3860731                                                                                                                                          |                                                                                                          aGGATGTGGGAAAATTCTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTAccc                                                                                                         <  1:74539/140‑1 (MQ=255)       tgGGAAAATTCTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTAACACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAg                                                                                                  >  1:213228/1‑140 (MQ=255)             aaTTCTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTatgat                                                                                            >  1:313364/1‑140 (MQ=255)                  ttttAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGata                                                                                       <  1:15293/140‑1 (MQ=255)                                    aaaTAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTAACACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCaaa                                                                     <  2:213228/140‑1 (MQ=255)                                       tAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCg                                                                  <  2:313364/140‑1 (MQ=255)                                       tAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCg                                                                  <  2:222932/140‑1 (MQ=255)                                       tAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCg                                                                  <  1:249553/140‑1 (MQ=255)                                                     aaCTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGc                                                    <  1:164036/140‑1 (MQ=255)                                                         gTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGaaa                                                >  1:35360/1‑140 (MQ=255)                                                                      aaaGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑TCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGGAATTTATAATCAGGa                                   <  1:313627/140‑1 (MQ=255)                                                                      aaaGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑TCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGGAATTTATAATCAGGa                                   >  2:313627/1‑140 (MQ=255)                                                                                    ttCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTg                     <  2:35360/140‑1 (MQ=255)                                                                                      cACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCt                   <  1:193769/140‑1 (MQ=255)                                                                                         tctTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAg                >  2:334706/1‑140 (MQ=255)                                                                                         tctTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAg                >  2:384139/1‑140 (MQ=255)                                                                                                   gtgAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACtt      >  1:177453/1‑140 (MQ=255)                                                                                                     ggaGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGGGTTAATAACTTCa    >  1:266789/3‑140 (MQ=255)                                                                                                      aaGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑‑CCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCAcc  >  2:130103/1‑140 (MQ=255)                                                                                                       aGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCAcc  >  2:169345/1‑140 (MQ=255)                                                                                                                                          |                                                                                                          AGGATGTGGGAAAATTCTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTACCCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCACC  >  NC_000913/3860489‑3860731

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