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_0009132,713,5680TC96.8% 102.5 / NA 31W225R (TGG→CGG) srmBATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (18/12);  minor base A (0/1);  total (18/13)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CTGGCGTAAACAGACCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCG  >  NC_000913/2713432‑2713701
                                                                                                                                        |                                                                                                                                     
cTGGCGTAAACAGACCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCgg                                                                                                                                     >  2:230696/1‑139 (MQ=255)
        aaCAGACCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCTGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACcg                                                                                                                             >  1:301435/1‑139 (MQ=255)
           aGACCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCgcgc                                                                                                                          <  2:720511/139‑1 (MQ=255)
             aCCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCg                                                                                                                        >  1:212477/1‑139 (MQ=255)
             aCCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCg                                                                                                                        >  2:688919/1‑139 (MQ=255)
             aCCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCg                                                                                                                        >  1:392335/1‑139 (MQ=255)
             aCCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCg                                                                                                                        >  1:560184/1‑139 (MQ=255)
             aCCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTACCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCg                                                                                                                        >  2:113065/1‑139 (MQ=255)
             aCCCTGCTCTTTTCGGCAACGCTGAAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCg                                                                                                                        >  1:434836/1‑139 (MQ=255)
             aCCCTGCTCTTTTCGGAAACGCTGGAAGGCGATGCGATTCAGGATTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGAGGTATTACCGCGCCg                                                                                                                        <  2:564889/139‑1 (MQ=255)
                   ctctTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATc                                                                                                                  >  2:251894/1‑139 (MQ=255)
                     ctTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCtt                                                                                                                <  2:212477/139‑1 (MQ=255)
                     ctTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTACCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCtt                                                                                                                <  1:113065/139‑1 (MQ=255)
                             cAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATaa                                                                                                        >  2:207711/1‑139 (MQ=255)
                                               cgatTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCa                                                                                      <  1:577101/139‑1 (MQ=255)
                                                            tGCCGAGCGTCTGCTGGAAGGTCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAg                                                                         >  2:300477/1‑139 (MQ=255)
                                                            tGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGt                                                                                                <  2:561290/116‑1 (MQ=255)
                                                            tGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGt                                                                                                >  1:561290/1‑116 (MQ=255)
                                                            tGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAg                                                                         >  1:255655/1‑139 (MQ=255)
                                                                       tgctgGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGAc                                                              <  2:392335/139‑1 (MQ=255)
                                                                         ctgGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGAccc                                                            <  2:334203/139‑1 (MQ=255)
                                                                            gAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCt                                                         >  2:653443/1‑139 (MQ=255)
                                                                             aaGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTc                                                        <  1:693786/139‑1 (MQ=255)
                                                                              aGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCt                                                         >  1:404446/1‑137 (MQ=255)
                                                                              aGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCt                                                         <  2:404446/137‑1 (MQ=255)
                                                                                   cGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATtgt                                                  <  1:592902/139‑1 (MQ=255)
                                                                                   cGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATtgt                                                  >  1:712631/1‑139 (MQ=255)
                                                                                            tttCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCg                                         <  2:710673/139‑1 (MQ=255)
                                                                                                                         aaaaaaaTTCATCAGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCaa            >  1:454649/1‑139 (MQ=255)
                                                                                                                                  catcaGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTgc   >  2:740688/1‑139 (MQ=255)
                                                                                                                                   atcaGCGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTgcg  <  1:251894/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                     
CTGGCGTAAACAGACCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATCCGGTGGAAGTTTCTGCCAATCCCTCCACCCGTGAGCGCAAAAAAATTCATCAGTGGTATTACCGCGCCGATGATCTTGAGCATAAAACCGCGTTGCTGGTGCATCTGTTAAAACAGCCGGAAGCGACCCGCTCAATTGTGTTTGTGCGTAAGCGTGAGCGTGTGCATGAGCTGGCAAACTGGCTGCG  >  NC_000913/2713432‑2713701

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