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

AAGGATGTGGGAAAATTCTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTACCCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCACCAGAAGCAATAAAAAATAATCTGGCAGTCG  >  NC_000913/3860488‑3860760                                                                                                                                           |                                                                                                                                       aaGGATGTGGGAAAATTCTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTAcc                                                                                                                                       >  2:160462/1‑140 (MQ=255) aaGGATGTGGGAAAATTCTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTAcc                                                                                                                                       >  1:526964/1‑140 (MQ=255)  aGGATGTGGGAAAATTCTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTAccc                                                                                                                                      <  1:243465/140‑1 (MQ=255)       gtgGGAAAATTCTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTAccccacca                                                                                                                                 <  2:526964/140‑5 (MQ=255)        tgGGAAAATTCTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAg                                                                                                                               >  2:141137/1‑140 (MQ=255)            aaaaTTCTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTatg                                                                                                                           <  2:125438/140‑1 (MQ=255)             aaaTTCTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTatga                                                                                                                          >  2:159991/1‑140 (MQ=255)               aTTCTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTatgatg                                                                                                                        >  1:2893/1‑140 (MQ=255)               aTTCTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTatgatg                                                                                                                        >  1:314573/1‑140 (MQ=255)                     ttAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATAtt                                                                                                                  >  1:475007/1‑140 (MQ=255)                       aaaGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTcc                                                                                                                <  1:160462/140‑1 (MQ=255)                             tGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTgg                                                                                                          >  2:150019/1‑140 (MQ=255)                                   aCAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCa                                                                                                    <  1:439347/140‑1 (MQ=255)                                        tAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCg                                                                                               <  2:387477/140‑1 (MQ=255)                                                       aCTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCg                                                                                <  1:538683/140‑1 (MQ=255)                                                       aCTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCg                                                                                <  2:2893/140‑1 (MQ=255)                                                              tataTTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTta                                                                         <  2:314573/140‑1 (MQ=255)                                                                               tGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCg                                                        >  1:160494/1‑140 (MQ=255)                                                                               tGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCg                                                        <  2:160494/140‑1 (MQ=255)                                                                                     ttCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTg                                                  >  2:469988/1‑140 (MQ=255)                                                                                              aTCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGGATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTtaa                                         >  1:84184/1‑140 (MQ=255)                                                                                                    gtgAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACtt                                   >  2:455528/1‑140 (MQ=255)                                                                                                      gAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCa                                 >  2:409498/1‑140 (MQ=255)                                                                                                       aaGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCAc                                >  1:263176/1‑140 (MQ=255)                                                                                                        aGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCAcc                               >  1:385042/1‑140 (MQ=255)                                                                                                             tAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCACCAGAAg                          <  1:141137/140‑1 (MQ=255)                                                                                                                              tggatggaATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCACCAGAAGCAATAAAAAATAATCTg         >  1:81672/1‑140 (MQ=255)                                                                                                                                gatggaATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCACCAGAAGCAATAAAAAATAATCTGGc       <  1:409498/140‑1 (MQ=255)                                                                                                                                     aaTTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCACCAGAAGCAATAAAAAATAATCTGGCAGTCg  >  1:488468/1‑140 (MQ=255)                                                                                                                                     aaTTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCACCAGAAGCAATAAAAAATAATCTGGCAGTCg  >  2:447934/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: ‑