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

TGGGAAAATTCTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTACCCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCACCAGAAGCAATAAAAAATA  >  NC_000913/3860495‑3860748                                                                                                                                    |                                                                                                                           tgGGAAAATTCTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTAcc                                                                                                                           >  1:479318/1‑133 (MQ=255) tgGGAAAATTCTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTAcc                                                                                                                           >  2:408524/1‑133 (MQ=255)          tCTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTa                                                                                                                 <  2:110744/133‑1 (MQ=255)            ttttAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTatg                                                                                                               <  2:479318/133‑1 (MQ=255)                              aaaTAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGt                                                                                             >  1:7421/1‑133 (MQ=255)                                  aaCACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGc                                                                                         >  2:295740/1‑133 (MQ=255)                                        gCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTTTGATGGATATTACCGTTGGTCAGCAAAGCg                                                                                   >  1:178312/1‑133 (MQ=255)                                        gCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCg                                                                                   <  1:349264/133‑1 (MQ=255)                                        gCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCg                                                                                   <  1:408524/133‑1 (MQ=255)                                        gCATGGCAACTGTATCATATTCTTAAAGCCCATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGTTGGATATTCCCGTTGGTCAGCAAAGCg                                                                                   <  1:142373/133‑1 (MQ=255)                                          aTGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CTCCACCAGTATGTCGGAAATTCCCGTTGGTCTGCAAAGCGtt                                                                                 >  2:142373/1‑133 (MQ=255)                                                            tCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATGTTCCCGTTGGTCAGCAAAGCGTTATGTATAGTGGCGAAAtt                                                               >  2:134613/1‑133 (MQ=255)                                                            tCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAAtt                                                               >  2:487837/1‑133 (MQ=255)                                                            tCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAAtt                                                               >  1:42967/1‑133 (MQ=255)                                                              ttAAAGCCAATGGATTTTCACTCTTATCCCTGTGAAGGTTTCGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTta                                                             <  2:131681/133‑1 (MQ=255)                                                              ttAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTta                                                             >  1:296940/1‑133 (MQ=255)                                                              ttAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTta                                                             >  1:297922/1‑133 (MQ=255)                                                              ttAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTta                                                             >  1:8652/1‑133 (MQ=255)                                                              ttAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTta                                                             >  1:121571/1‑133 (MQ=255)                                                                aaaGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTtata                                                           >  2:298968/1‑133 (MQ=255)                                                                aaaGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTCATTAATTATGGCGTAATTtata                                                           >  1:175499/1‑133 (MQ=255)                                                                        tGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGaa                                                   >  2:108040/1‑133 (MQ=255)                                                                               tCACTCTTATCACTGTGAAGGTTTCGTCAGGGTACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCg                                            <  2:178312/133‑1 (MQ=255)                                                                                     ttATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTg                                      <  2:121571/133‑1 (MQ=255)                                                                                     ttATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTg                                      <  2:297922/133‑1 (MQ=255)                                                                                       aTCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCAGCGAAGTTGCt                                    >  2:151766/1‑133 (MQ=255)                                                                                       aTCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCt                                    >  2:504813/1‑133 (MQ=255)                                                                                       aTCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCt                                    >  1:25474/1‑133 (MQ=255)                                                                                       aTCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCt                                    <  1:295740/133‑1 (MQ=255)                                                                                           cTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCt                                    <  1:42540/129‑1 (MQ=255)                                                                                           cTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCt                                    >  2:42540/1‑129 (MQ=255)                                                                                            tgtgAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGtt                               >  1:458457/1‑133 (MQ=255)                                                                                             gtgAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTta                              >  2:53168/1‑133 (MQ=255)                                                                                             gtgAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTta                              >  2:116049/1‑133 (MQ=255)                                                                                             gtgAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTACCGTTGGTCAACAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTta                              >  2:74910/1‑133 (MQ=255)                                                                                              tgAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTtaa                             <  1:116049/133‑1 (MQ=255)                                                                                               gAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTtaat                            >  2:401968/1‑133 (MQ=255)                                                                                               gAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTtaat                            <  1:298968/133‑1 (MQ=255)                                                                                                      tAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCa                     >  1:229547/1‑133 (MQ=255)                                                                                                           tGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCACCAGa                <  2:458457/133‑1 (MQ=255)                                                                                                            gggAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCACCAGaa               >  2:314221/1‑133 (MQ=255)                                                                                                             ggAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCACCAGAAg              >  1:34021/1‑133 (MQ=255)                                                                                                               aaCAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCACCAGAAGCa            <  1:487837/133‑1 (MQ=255)                                                                                                                     aCTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCACCAGAAGCAATaaaa      <  1:504813/133‑1 (MQ=255)                                                                                                                         gatggaATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAAGTCACCGAAGTTGCTGAGTTAATAACTTCACCAGAAGCAATAAAAaata  >  2:357444/1‑133 (MQ=255)                                                                                                                         gatggaATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCACCAGAAGCAAt          >  2:386290/1‑125 (MQ=255)                                                                                                                         gatggaATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCACCAGAAGCAAt          <  1:386290/125‑1 (MQ=255)                                                                                                                                    |                                                                                                                           TGGGAAAATTCTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTACCCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCACCAGAAGCAATAAAAAATA  >  NC_000913/3860495‑3860748

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