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%	81.3 / NA	20	coding (344/894 nt)	yidL	putative DNA‑binding transcriptional regulator YidL
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (14/6);  total (14/6)

AAGGATGTGGGAAAATTCTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTACCCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCACCAGAAGCAATAAAAAATAATCTGGCAGTCG  >  NC_000913/3860488‑3860760                                                                                                                                           |                                                                                                                                       aaGGATGTGGGAAAATTCTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTAc                                                                                                                                        >  2:317724/1‑139 (MQ=255) aaGAATGTGGGAAAATTCTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTAc                                                                                                                                        <  1:317724/139‑1 (MQ=255)        tgGGAAAATTCTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAg                                                                                                                               >  1:198967/1‑140 (MQ=255)                                        tAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCg                                                                                               <  2:401886/140‑1 (MQ=255)                                              tGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGttattt                                                                                         >  1:366669/1‑140 (MQ=255)                                                cATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGttatttat                                                                                       >  1:275376/1‑140 (MQ=255)                                                     cAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATgg                                                                                  >  1:314589/1‑140 (MQ=255)                                                      aaCTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGc                                                                                 <  2:41419/140‑1 (MQ=255)                                                       aCTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCg                                                                                <  1:41326/140‑1 (MQ=255)                                                          gTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGaaa                                                                             >  1:24987/1‑140 (MQ=255)                                                              tataTTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTta                                                                         <  2:394236/140‑1 (MQ=255)                                                                   tCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGTAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATc                                                                    >  1:366703/1‑140 (MQ=255)                                                                                      tCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGc                                                 <  2:314589/140‑1 (MQ=255)                                                                                          tctTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAg                                             >  2:43584/1‑140 (MQ=255)                                                                                                   tgtgAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACt                                    <  2:366669/140‑1 (MQ=255)                                                                                                    gtgAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACtt                                   >  2:351940/1‑140 (MQ=255)                                                                                                    gtgAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACtt                                   >  1:284605/1‑140 (MQ=255)                                                                                                      gAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCa                                 >  2:336985/1‑140 (MQ=255)                                                                                                      gAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCa                                 >  2:197395/1‑140 (MQ=255)                                                                                                             tAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCACCAGAAg                          >  1:340841/1‑140 (MQ=255)                                                                                                             tAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCACCAGAAg                          >  1:277402/1‑140 (MQ=255)                                                                                                                                     aaTTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCACCAGAAGCAATAAAAAATAATCTGGCAGTCg  >  2:244754/1‑140 (MQ=255)                                                                                                                                           |                                                                                                                                       AAGGATGTGGGAAAATTCTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTACCCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCACCAGAAGCAATAAAAAATAATCTGGCAGTCG  >  NC_000913/3860488‑3860760

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