Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,713,568	T→C	W225R (TGG→CGG)	srmB →	ATP‑dependent RNA helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,713,568	0	T	C	100.0%	121.7 / NA	35	W225R (TGG→CGG)	srmB	ATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/17);  total (18/17)

ACCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACT  >  NC_000913/2713445‑2713694                                                                                                                            |                                                                                                                               aCCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCgcgc                                                                                                                   >  2:421611/1‑137 (MQ=255) aCCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCgcgc                                                                                                                   <  1:421611/137‑1 (MQ=255)       ctctTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATc                                                                                                           >  1:1161724/1‑139 (MQ=255)       ctctTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATc                                                                                                           <  1:22439/139‑1 (MQ=255)                        ggAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACcgcg                                                                                          <  1:368274/139‑1 (MQ=255)                          aaGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGtt                                                                                        <  1:1092880/139‑1 (MQ=255)                                 tgcgatTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTg                                                                                 >  1:745484/1‑139 (MQ=255)                                  gcgatTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGc                                                                                <  2:1173909/139‑1 (MQ=255)                                   cgatTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTg                                                                                       <  1:276834/131‑1 (MQ=255)                                   cgatTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTg                                                                                       >  2:276834/1‑131 (MQ=255)                                   cgatTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCa                                                                               <  2:745484/139‑1 (MQ=255)                                            aCTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTNNNCCAATNCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTaaa                                                                      <  1:81079/139‑1 (MQ=255)                                                tGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACcgcg                                                                                          >  1:819028/1‑115 (MQ=255)                                                tGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACcgcg                                                                                          <  2:819028/115‑1 (MQ=255)                                                tGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAg                                                                  >  2:855843/1‑139 (MQ=255)                                                tGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAg                                                                  >  2:769876/1‑139 (MQ=255)                                                tGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAg                                                                  >  1:1298845/1‑139 (MQ=255)                                                tGCCGAGCGCCTGCTGGAAAATCGGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAg                                                                  >  2:1228833/1‑139 (MQ=255)                                                               ggAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGc                                                   <  2:1093883/139‑1 (MQ=255)                                                               ggAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGc                                                   <  2:1161724/139‑1 (MQ=255)                                                               ggAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGc                                                   <  1:1228833/139‑1 (MQ=255)                                                                     tCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAAtt                                             >  2:11930/1‑139 (MQ=255)                                                                                      ccAATCCCTCCACCCGTGAGCGCAAAAAAATTCTTCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAagcg                            <  1:364186/139‑1 (MQ=255)                                                                                        aaTCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAagcgtg                          <  1:769876/139‑1 (MQ=255)                                                                                          tCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTc                                                 >  1:577200/1‑114 (MQ=255)                                                                                          tCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTc                                                 <  2:577200/114‑1 (MQ=255)                                                                                                   ccGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGt                                                                          <  2:631064/80‑1 (MQ=255)                                                                                                   ccGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGt                                                                          >  1:631064/1‑80 (MQ=255)                                                                                                     gTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATtttgtttgagtgt                                 <  2:763134/119‑6 (MQ=255)                                                                                                     gTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGt                                 >  1:763134/1‑119 (MQ=255)                                                                                                     gTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGa             >  1:34240/1‑139 (MQ=255)                                                                                                      tGAGCGCAAAAAAATTCATCAGCGGTATGACCGCGCCGATGATCTTGAGAATAAAACCGCGGTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAg            >  2:849276/1‑139 (MQ=255)                                                                                                        aGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCt          >  1:709827/1‑139 (MQ=255)                                                                                                        aGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCt          >  1:869137/1‑139 (MQ=255)                                                                                                             aaaaaaaTTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAA‑‑‑AGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACt  >  2:1093103/1‑139 (MQ=255)                                                                                                                            |                                                                                                                               ACCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACT  >  NC_000913/2713445‑2713694

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