Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,513,956 T→C 7.5% L42P (CTG→CCG)  rsmH → ribosomal RNA small subunit methyltransferase H

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,513,9560TC 7.5% 148.9 / 2.2 54L42P (CTG→CCG) rsmHribosomal RNA small subunit methyltransferase H
Reads supporting (aligned to +/- strand):  ref base T (26/23);  new base C (2/2);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

CGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACAAAGATCCACAAGCGATTGCCACGGGGCAAGCGCTGGCGGCCGAAGACGGCCGCTTTGTCATTGTGCAGCGCAGTTTTGCCGAGCTGGGTGCGGAAGT  >  NC_002947/1513808‑1514100
                                                                                                                                                    |                                                                                                                                                
cgacgaATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCt                                                                                                                                                  >  3:88919/1‑149 (MQ=255)
    gaATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATa                                                                                                                                              <  5:102264/149‑1 (MQ=255)
    gaATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATa                                                                                                                                              >  1:9555/1‑149 (MQ=255)
     aaTTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGGTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATAc                                                                                                                                             <  4:88919/149‑1 (MQ=255)
         gCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCa                                                                                                                                          <  5:144125/148‑1 (MQ=255)
          cGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGc                                                                                                                                        <  4:245937/149‑1 (MQ=255)
          cGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGAGGCGGACACAGCCGTCTGATACTCAg                                                                                                                                         <  1:8921/148‑1 (MQ=255)
               ccTGAACCTGAGCCCATAGATAGCGCCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGACTAGAGCGCCGACGGTTGCAATCTTGACGGCACTTTCGGCCGGGGCGGCCACACCCGTCTGATACTCAGCAACCt                                                                                                                                   >  4:288404/1‑149 (MQ=255)
                cTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTc                                                                                                                                  >  5:21216/1‑149 (MQ=255)
                 tGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCg                                                                                                                                 >  5:286480/1‑149 (MQ=255)
                       tgtgACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCCGATACTCAGCAAGCTCgg                                                                                                                                >  1:105784/1‑144 (MQ=255)
                       tgtgACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCCGATACTCAGCAAGCTCgg                                                                                                                                >  1:101338/1‑144 (MQ=255)
                       tgtgACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCCGATACTCAGCAAGCTCgg                                                                                                                                <  2:101338/144‑1 (MQ=255)
                       tgtgACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCCGATACTCAGCAAGCTCgg                                                                                                                                <  2:105784/144‑1 (MQ=255)
                        gtgACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCCGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGc                                                                                                                           >  6:127805/1‑148 (MQ=255)
                                  atagCGGCTTCAACCACATTACCGTCCTTCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGctc                                                                                                                <  7:184658/149‑1 (MQ=255)
                                      cGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTgg                                                                                                            <  7:102500/149‑1 (MQ=255)
                                             aaCCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCTTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTgg                                                                                                            >  7:17963/1‑142 (MQ=255)
                                             aaCCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCTTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTgg                                                                                                            <  8:17963/142‑1 (MQ=255)
                                             aaCCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTGGGGCCGc                                                                                                                           >  5:323273/1‑127 (MQ=255)
                                             aaCCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGc                                                                                                                           <  6:323273/127‑1 (MQ=255)
                                                       ccGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCCGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACAAAGATCca                                                                                           <  5:127805/149‑1 (MQ=255)
                                                            cTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACAAAg                                                                                                >  6:14064/1‑139 (MQ=255)
                                                            cTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACAAAg                                                                                                <  5:14064/139‑1 (MQ=255)
                                                              gCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACAAAGATCCACAAGCGa                                                                                    >  3:257074/1‑149 (MQ=255)
                                                              gCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACAAAGATCCACAAGCGa                                                                                    >  4:151122/1‑149 (MQ=255)
                                                                  gacgaAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACAAAGATCCACAAGCGATTGc                                                                                <  4:257074/149‑1 (MQ=255)
                                                                       aGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACAAAGATCCACAAGCGATTGCCACgg                                                                           <  2:9555/149‑1 (MQ=255)
                                                                        gCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACAAAGATCCACAAGCGATTGCCACggg                                                                          <  3:151122/149‑1 (MQ=255)
                                                                                 gCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACAAAGATCCACAAGc                                                                                      <  8:289534/128‑1 (MQ=255)
                                                                                 gCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACAAAGATCCACAAGc                                                                                      >  7:289534/1‑128 (MQ=255)
                                                                                  cATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACAAAGATCCACAAGCGATTGCCACGGGGCAAGCGCTg                                                                <  1:210369/149‑1 (MQ=255)
                                                                                   aTTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTg                                                                                                             >  4:44587/1‑103 (MQ=255)
                                                                                   aTTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTg                                                                                                             <  3:44587/103‑1 (MQ=255)
                                                                                    ttGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACAAAGATCCACAAGCGATTGCCACGGGGCAAGCGCTggc                                                              >  5:27684/1‑149 (MQ=255)
                                                                                       gCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACaaa                                                                                                 <  1:51297/111‑1 (MQ=255)
                                                                                       gCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACaaa                                                                                                 >  2:51297/1‑111 (MQ=255)
                                                                                          ttGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACAAAGATCCACAAGCGATTGCCACGGGGCAAGCGCTGGCGGCCg                                                         >  1:150590/1‑148 (MQ=255)
                                                                                            gcgagcCGACGGGTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACATCCGTCGGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACAAAGATCCACAAGCGATTGCCACGGGGCAAGCGCTGGCGGCCGaa                                                       <  7:298743/148‑1 (MQ=255)
                                                                                                  cGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACAAAGATCCACAAGCGATTGCCACGGGGCAAGCGCTGGCGGCCGAAGACGGcc                                                <  6:21216/149‑1 (MQ=255)
                                                                                                                tGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACaaa                                                                                                 >  7:96131/1‑86 (MQ=255)
                                                                                                                tGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACaaa                                                                                                 <  8:96131/86‑1 (MQ=255)
                                                                                                                     ggCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACAAAGATCCACAAGCGATTGCCACGGGGCAAGCGCTGGCGGCCGAAGACGGCCGCTTTGTCAtt                                     <  1:266515/141‑1 (MQ=255)
                                                                                                                     ggCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACAAAGATCCACAAGCGATTGCCACGGGGCAAGCGCTGGCGGCCGAAGACGGCCGCTTTGTCAtt                                     >  2:266515/1‑141 (MQ=255)
                                                                                                                     ggCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACAAAGATCCACAAGCGATTGCCACGGGGCAAGCGCTGGCGGCCGAAGACGGCCGCTTTGTCATTGTGCAgcg                             >  7:292849/1‑149 (MQ=255)
                                                                                                                        aCTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACAAAGATCCACAAGCGATTGCCACGGGGCAAGCGCTGGCGGCCGAAGACGGCCGCTTTGTCATTGTGCAGCGCAg                          >  7:98295/1‑149 (MQ=255)
                                                                                                                               gCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACAAAGATCCACAAGCGATTGCCACGGGGCAAGCGCTGGCGGCCGAAGACGGCCGCTTTGTCATTGTGCAGCGCAGTTTTGcc                   >  3:101541/1‑149 (MQ=255)
                                                                                                                                ccGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACAAAGATCCACAAGCGATTGCCACGGGGCAAGCGCTGGCGGCCGAAGACGGCCGCTTTGTCATTGTGCAGCGCAGTTTTGcc                   >  5:86191/1‑148 (MQ=255)
                                                                                                                                      cGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACAAAGATCCACAAGCGATTGCCACGGGGCAAGCGCTGGCGGCCGAAGACGGCCGCTTTGTCATTGTGCAGCGCAGTTTTGCCGAGCTgg            >  1:76800/1‑149 (MQ=255)
                                                                                                                                        gACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACAAAGATCCACAAGCGATTGCCACGGGGCAAGCGCTGGCGGCCGAAGACGGCCGCTTTGTCATTGTGCAGCGCAGTTTTGCCGAGCTGGGt          <  4:101541/149‑1 (MQ=255)
                                                                                                                                           acaGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACAAAGATCCACAAGCGATTGCCACGGGGCAAGCGCTGGCGGCCGAAGACGGCCGCTTTGTCATTGTGCAGCGCAGTTTTGCCGAGCTGGGTGCg       >  7:267694/1‑149 (MQ=255)
                                                                                                                                           acaGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACAAAGATCCACAAGCGATTGCCACGGGGCAAGCGCTGGCGGCCGAAGACGGCCGCTTTGTCATTGTGCAGCGCAGTTTTGCCGAGCTGGGTGCg       >  2:316623/1‑149 (MQ=255)
                                                                                                                                             aGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACAAAGATCCACAAGCGATTGCCACGGGGCAAGCGCTGGCGGCCGAAGACGGCCGCTTTGTCATTGTGCAGCGCAGTTTTGCCGAGCTGGGTGCgg      >  5:117786/1‑148 (MQ=255)
                                                                                                                                               ccGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACAAAGATCCACAAGCGATTGCCACGGGGCAAGCGCTGGCGGCCGAAGACGGCCGCTTTGTCATTGTGCAGCGCAGTTTTGCCGAGCTGGGTGCGGAAg   >  6:51071/1‑149 (MQ=255)
                                                                                                                                                cGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACAAAGATCCACAAGCGATTGCCACGGGGCAAGCGCTGGCGGCCGAAGACGGCCGCTTTGTCATTGTGCAGCGCAGTTTTGCCGAGCTGGGTGCGGAAGt  <  2:76800/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CGACGAATTGCGTGACCTGATCCTGTGACCATAGATAGCGGCTTCAACCACATTACCGTCCTGCTCGACGAAGCTGTCGAGGCATTGGCCTTGCGCGCCGACGGTTGCTATCTGGACGGCACTTTCGGCCGGGGCGGACACAGCCGTCTGATACTCAGCAAGCTCGGGCCGCAAGGGCGGCTCTTGGGCTTCGACAAAGATCCACAAGCGATTGCCACGGGGCAAGCGCTGGCGGCCGAAGACGGCCGCTTTGTCATTGTGCAGCGCAGTTTTGCCGAGCTGGGTGCGGAAGT  >  NC_002947/1513808‑1514100

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: ‑