Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,717,710 G→A 7.8% A137T (GCG→ACG)  rsmJ → ribosomal RNA small subunit methyltransferase J

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,717,7100GA 7.8% 147.8 / 2.4 51A137T (GCG→ACG) rsmJribosomal RNA small subunit methyltransferase J
Reads supporting (aligned to +/- strand):  ref base G (25/22);  new base A (2/2);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

GGCATTGCCCAAGGTATCCGGCCGCAGGTACTGGATGCCACGGCGGGGCTGGGCAAGGACGCGTTCGTGCTGGCCAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCATCGAGCGCATGCGCGCCTGGGAAGGGGAGGCGCCGCAGGTGATCTACCTCGACCCGATGTTTCCACACCGCGACAAGA  >  NC_002947/1717566‑1717857
                                                                                                                                                |                                                                                                                                                   
ggCATTGCCCAAGGTATCCGGCCGCAGGTACTGGATGCCACGGCGGGGCTGGGCAAGGACGCGTTCGTGCTGGCCAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTggcgc                                                                                                                                                  <  4:129309/148‑1 (MQ=255)
 gCATTGCCCAAGGTATCCGGCCGCAGGTACTGGATGCCACGGCGGGGCTGGGCAAGGACGCGTTCGTGCTGGCCAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTggcgcgg                                                                                                                                                <  6:244863/149‑1 (MQ=255)
   aTTGCCCAAGGTATCCGGCCGCAGGTACTGGATGCCACGGCGGGGCTGGGCAAGGACGCGTTCGTGCTGGCCAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTggcgcgggc                                                                                                                                              >  7:302911/1‑149 (MQ=255)
       cccAAGGTATCCGGCCGCAGGTACTGGATGCCACGGCGGGGCTGGGCAAGGACGCGTTCGTGCTGGCCAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGa                                                                                                                                          <  8:250707/149‑1 (MQ=255)
       cccAAGGTATCCGGCCGCAGGTACTGGATGCCACGGCGGGGCTGGGCAAGGACGCGTTCGTGCTGGCCAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGa                                                                                                                                          <  7:249912/149‑1 (MQ=255)
        ccAAGGTATCCGGCCGCAGGTACTGGATGCCACGGCGGGGCTGGGCAAGGACGCGTTCGTGCTGGCCAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAg                                                                                                                                         >  8:193206/1‑149 (MQ=255)
         cAAGGTATCCGGCCGCAGGTACTGGATGCCACGGCGGGGCTGGGCAAGGACGCGTTCGTGCTGGCCAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGc                                                                                                                                        <  2:292335/149‑1 (MQ=255)
                tCCGGCCGCAGGTACTGGATGCCACGGCGGGGCTGGGCAAGGACGCGTTCGTGCTGGCCAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATg                                                                                                                                   >  6:53602/1‑147 (MQ=255)
                    gCCGCAGGTACTGGATGCCACGGCGGGGCTGGGCAAGGACGCGTTCGTGCTGGCCAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGaagaag                                                                                                                             >  4:204754/1‑149 (MQ=255)
                          ggTACTGGATGCCACGGCGGGGCTGGGCAAGGACGCGTTCGTGCTGGCCAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGcc                                                                                                                       >  1:297762/1‑149 (MQ=255)
                              cTGGATGCCACGGCGGGGCTGGGCAAGGACGGGTTGGTGCTGGCCAGCCTTGGCTGCCAGATGACCCTGTTCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCTAGCGGATGAAGAAGTGGGCCCGAt                                                                                                                   <  1:118800/149‑1 (MQ=255)
                              cTGGATGCCACGGCGGGGCTGGGCAAGGACGCGTTCGTGCTGGCCAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGa                                                                                                                    <  4:178098/148‑1 (MQ=255)
                                ggATGCCACGGCGGGGCTGGGCAAGGACGCGTTCGTGCTGGCCAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTg                                                                                                                 >  7:217671/1‑149 (MQ=255)
                                  aTGCCACGGCGGGGCTGGGCAAGGACGCGTTCGTGCTGGCCAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTc                                                                                                               <  7:193206/149‑1 (MQ=255)
                                       aCGGCGGGGCTGGGCAAGGACGCGTTCGTGCTGGCCAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGgc                                                                                                           >  2:261314/1‑148 (MQ=255)
                                                 tGGGCAAGGACGCGTTCGTGCTGGCCAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCctg                                                                                                <  2:297762/149‑1 (MQ=255)
                                                   ggCAAGGACGCGTTCGTGCTGGCCAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCctgct                                                                                              >  5:217203/1‑149 (MQ=255)
                                                      aaGGACGCGTTCGTGCTGGCCAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGAc                                                                                           >  2:224327/1‑149 (MQ=255)
                                                         gACGCGTTCGTGCTGGCCAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCgg                                                                                        >  1:60046/1‑149 (MQ=255)
                                                              gTTCGTGCTGGCCAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACg                                                                                   >  4:296874/1‑149 (MQ=255)
                                                              gTTCGTGCTGGCCAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGACGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACg                                                                                   >  4:6602/1‑149 (MQ=255)
                                                              gTTCGTGCTGGCCAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGACGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACg                                                                                   >  4:2093/1‑149 (MQ=255)
                                                                  gTGCTGGCCAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCAt                                                                               >  1:115320/1‑149 (MQ=255)
                                                                    gCTGGCCAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCATc                                                                              >  5:159020/1‑148 (MQ=255)
                                                                     cTGGCCAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCATc                                                                              >  3:103204/1‑147 (MQ=255)
                                                                     cTGGCCAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCATc                                                                              >  4:306457/1‑147 (MQ=255)
                                                                        gCCAGCCTTGGCTGCCAGATGACCCTGAACGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGACGACGCCCTGGCCCGGGCGCGAGCGCATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCATCGAgcg                                                                         >  7:89758/1‑149 (MQ=255)
                                                                         ccAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGACGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCATCGAgcgc                                                                        <  3:6602/149‑1 (MQ=255)
                                                                         ccAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGACGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCATCGAgcgc                                                                        <  3:2093/149‑1 (MQ=255)
                                                                          cAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCctg                                                                                                <  2:8686/124‑1 (MQ=255)
                                                                          cAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCctg                                                                                                >  1:8686/1‑124 (MQ=255)
                                                                          cAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCATCGAGCGCa                                                                       >  3:6239/1‑149 (MQ=255)
                                                                                 ggCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCATCGAGCGCATgcgcg                                                                 <  6:159020/148‑1 (MQ=255)
                                                                                  gCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCATCGAGCGCATgcgcgc                                                                <  3:306457/148‑1 (MQ=255)
                                                                                      ccAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCATCGAGCGCATGCGCGCCTgg                                                            >  6:86708/1‑148 (MQ=255)
                                                                                             aCCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCATCGAGCGCATGCGCGCCTGGGAAggg                                                      >  1:113696/1‑147 (MQ=255)
                                                                                             aCCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCATCGAGCGCATGCGCGCCTGGGAAggg                                                      <  2:113696/147‑1 (MQ=255)
                                                                                             aCCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCATCGAGCGCATGCGCGCCTGGGAAGGGGa                                                    <  8:302911/149‑1 (MQ=255)
                                                                                             aCCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCATCGAGCGCATGCGCGCCTGGGAAGGGGa                                                    <  3:296874/149‑1 (MQ=255)
                                                                                               ccTGATCGAGCGCCAGTCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCATCGAGCGCATGCGCGCCTGGGAAGGGGAgg                                                  >  5:156948/1‑149 (MQ=255)
                                                                                               ccTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCATCGAGCGCATGCGCGCCTGGGAAGGGGAgg                                                  >  6:151524/1‑149 (MQ=255)
                                                                                               ccTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCATCGAGCGCATGCGCGCCTGGGAAGGGGAgg                                                  <  4:6239/149‑1 (MQ=255)
                                                                                                 tGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCATCGAGCGCATGCGCGCCTGGGAAGGGGAGgcg                                                >  8:129726/1‑149 (MQ=255)
                                                                                                        gcgcCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCATCGAGCGCATGCGCGCCTGGGAAGGGGAGGCGCCGCAgg                                         <  5:36299/149‑1 (MQ=255)
                                                                                                                cGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCATc                                                                              >  8:155166/1‑104 (MQ=255)
                                                                                                                cGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCATc                                                                              <  7:155166/104‑1 (MQ=255)
                                                                                                                cGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCATCGAGCGCATGCGCGCCTGGGAAGGGGAGGCGCCGCAGGTGATCTAc                                 <  2:115320/149‑1 (MQ=255)
                                                                                                                      ttGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACTCCATCGAGCGCATGCGCGCCTGGGAAGGGGAGGCGCCGCAGGTGATCTACCTCGAc                           <  5:53602/149‑1 (MQ=255)
                                                                                                                           gcgcTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCATCGAGCGCATGCGCGCCTGGGAAGGGGAGGCGCCGCAGGTGATCTACCTCGACCCGAt                      <  7:129726/149‑1 (MQ=255)
                                                                                                                           gcgcTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCATCGAGCGCATGCGCGCCTGGGAAGGGGAGGCGCCGCAGGTGATCTACCTCGACCCGAt                      <  3:204754/149‑1 (MQ=255)
                                                                                                                             gctgctGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCATCGAGCGCATGCGCGCCTGGGAAGGGGAGGCGCCGCAGGTGATCTACCTCGACCCGATGt                    <  6:156948/149‑1 (MQ=255)
                                                                                                                                       gaCGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCATCGAGCGCATGCGCGCCTGGGAAGGGGAGGCGCCGCAGGTGATCTACCTCGACCCGATGTTTCCACACcg          >  4:292529/1‑149 (MQ=255)
                                                                                                                                        aCGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCATCGAGCGCATGCGCGCCTGGGAAGGGGAGGCGCCGCAGGTGATCTACCTCGACCCGATGTTTCCACACcgc         >  6:24491/1‑149 (MQ=255)
                                                                                                                                          ggCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCATCGAGCGCATGCGCGCCTGGGAAGGGGAGGCGCCGCAGGTGATCTACCTCGACCCGATGTTTCCACACCGCGa       <  6:217203/149‑1 (MQ=255)
                                                                                                                                               ggcgcgggcgcgAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCATCGAGCGCATGCGCGCCTGGGAAGGGGAGGCGCCGCAGGTGATCTACCTCGACCCGATGTTTCCACACCGCGACAaga  <  5:151524/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GGCATTGCCCAAGGTATCCGGCCGCAGGTACTGGATGCCACGGCGGGGCTGGGCAAGGACGCGTTCGTGCTGGCCAGCCTTGGCTGCCAGATGACCCTGATCGAGCGCCAGCCGCTGATTGCCGCGCTGCTGGAGGACGGCCTGGCGCGGGCGCGAGCGGATGAAGAAGTGGGCCCGATTGTCGGGCGCATGCGCCTGCTGACCGGCAACGCCATCGAGCGCATGCGCGCCTGGGAAGGGGAGGCGCCGCAGGTGATCTACCTCGACCCGATGTTTCCACACCGCGACAAGA  >  NC_002947/1717566‑1717857

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