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%	177.4 / NA	51	W225R (TGG→CGG)	srmB	ATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/24);  total (27/24)

ACAGACCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCG  >  NC_000913/2713441‑2713703                                                                                                                                |                                                                                                                                        aCAGACCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACcgc                                                                                                                              <  1:1418090/139‑1 (MQ=255)     aCCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCg                                                                                                                          >  1:438173/1‑139 (MQ=255)     aCCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCg                                                                                                                          >  1:818139/1‑139 (MQ=255)                          cTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACcg                                                                                                     >  1:1509469/1‑139 (MQ=255)                            ggAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACcgcg                                                                                                   <  2:1509469/139‑1 (MQ=255)                              aaGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGtt                                                                                                 <  2:484834/139‑1 (MQ=255)                                       cgatTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCa                                                                                        <  2:304021/139‑1 (MQ=255)                                               gACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCga                                                                                                                         <  2:891005/98‑1 (MQ=255)                                               gACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCga                                                                                                                         >  1:891005/1‑98 (MQ=255)                                               gACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTaa                                                                                >  2:296372/1‑139 (MQ=255)                                                    tGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACcgcg                                                                                                   >  2:4309/1‑115 (MQ=255)                                                    tGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACcgcg                                                                                                   <  1:4309/115‑1 (MQ=255)                                                    tGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACcg                                                                                                     >  2:286919/1‑113 (MQ=255)                                                    tGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACcg                                                                                                     >  1:1369151/1‑113 (MQ=255)                                                    tGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACcg                                                                                                     <  2:1369151/113‑1 (MQ=255)                                                    tGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACcg                                                                                                     <  1:286919/113‑1 (MQ=255)                                                    tGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGtt                                                                                                 <  1:219936/117‑1 (MQ=255)                                                    tGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGtt                                                                                                 >  2:219936/1‑117 (MQ=255)                                                    tGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAg                                                                           >  1:1048126/1‑139 (MQ=255)                                                    tGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAg                                                                           >  2:342829/1‑139 (MQ=255)                                                       cGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCg                                                                        <  1:342829/139‑1 (MQ=255)                                                              ctgctgGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGa                                                                 >  1:679657/1‑139 (MQ=255)                                                                gctgGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCgcgc                                                                                                                            >  1:1430506/1‑78 (MQ=255)                                                                gctgGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCgcgc                                                                                                                            <  2:1430506/78‑1 (MQ=255)                                                                      aGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATc                                                                                                                    >  1:753182/1‑80 (MQ=255)                                                                      aGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATc                                                                                                                    <  2:753182/80‑1 (MQ=255)                                                                         tCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGtt                                                                                                 <  2:1172894/96‑1 (MQ=255)                                                                         tCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGtt                                                                                                 >  1:1172894/1‑96 (MQ=255)                                                                         tCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAAtt                                                      >  1:1196837/1‑139 (MQ=255)                                                                           cGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGc                                                            <  1:265555/131‑1 (MQ=255)                                                                           cGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGc                                                            >  2:265555/1‑131 (MQ=255)                                                                             gTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTgtgt                                                  <  2:1196837/139‑1 (MQ=255)                                                                                     ttCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGt                                          <  2:421556/139‑1 (MQ=255)                                                                                           cAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAagcgt                                    <  1:296372/139‑1 (MQ=255)                                                                                             aTCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAagcgtga                                  >  1:557543/1‑139 (MQ=255)                                                                                                  tCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAagcgtgagcgtg                             >  2:254458/1‑139 (MQ=255)                                                                                                         gTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAagcgtgagc                                >  2:899372/1‑129 (MQ=255)                                                                                                         gTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAagcgtgagc                                <  1:899372/129‑1 (MQ=255)                                                                                                          tGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCgatg                                                                                                                       >  2:3469/1‑41 (MQ=255)                                                                                                          tGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCgatg                                                                                                                       <  1:3469/41‑1 (MQ=255)                                                                                                          tGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGTATAAAACCGCGTTGCTGGTGCa                                                                                        >  2:1132371/1‑72 (MQ=255)                                                                                                          tGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGTATAAAACCGCGTTGCTGGTGCa                                                                                        <  1:1132371/72‑1 (MQ=255)                                                                                                          tGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAg                     >  1:151444/1‑139 (MQ=255)                                                                                                            aGCGCAAAAAATTTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCa                                                                                                             >  1:522288/1‑49 (MQ=38)                                                                                                            aGCGCAAAAAATTTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCa                                                                                                             <  2:522288/49‑1 (MQ=38)                                                                                                            aGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGc                                                                          <  2:1245516/84‑1 (MQ=255)                                                                                                            aGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGc                                                                          >  1:1245516/1‑84 (MQ=255)                                                                                                            aGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCt                   >  1:1491771/1‑139 (MQ=255)                                                                                                                cAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGc                                                            >  2:1260586/1‑94 (MQ=255)                                                                                                                cAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGc                                                            <  1:1260586/94‑1 (MQ=255)                                                                                                                             caGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGcgcg  <  2:557543/139‑1 (MQ=255)                                                                                                                                |                                                                                                                                        ACAGACCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCGCG  >  NC_000913/2713441‑2713703

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