Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 698,297 A→G 7.2% noncoding (474/1522 nt) PP_16SD → 16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947698,2970AG 7.2% 146.3 / 2.4 55noncoding (474/1522 nt)PP_16SD16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (27/24);  new base G (2/2);  total (29/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.75e-01

GAGGCAGCAGTGGGGAATATTGGACAATGGGCGAAAGCCTGATCCAGCCATGCCGCGTGTGTGAAGAAGGTCTTCGGATTGTAAAGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAGGTGGTTTGTTAAGTTGGATGTGAAAGCCCCGGGCTCAAC  >  NC_002947/698150‑698425
                                                                                                                                                   |                                                                                                                                
gAGGCAGCAGTGGGGAATATTGGACAATGGGCGAAAGCCTGATCCAGCCCTGCCGCGTGTGTGAAGAAGGTCTTCGGATTGTAAAGCACTTTATGTTGGGAGGAAGGGGAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGaa                                                                                                                                 >  5:112301/1‑149 (MQ=12)
  ggCAGCAGTGGGGAATATTGGACAATGGGCGAAAGCCTGATCCAGCCATGCCGCGTGTGTGAAGAAGGTCTTCGGATTGTAAAGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAAt                                                                                                                                >  4:213362/1‑148 (MQ=32)
    cagcagTGGGGAATATTGGACAATGGGCGAAAGCCTGATCCAGCCATGCCGCGTGTGTGAAGAAGGTCTTCGGATTGTAAAGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAg                                                                                                                             >  8:217075/1‑149 (MQ=32)
        agTGGGGAATATTGGACAATGGGCGAAAGCCTGATCCAGCCATGCCGCGTGTGTGAAGAAGGTCTTCGGATTGTAAAGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCAcc                                                                                                                         <  4:166124/149‑1 (MQ=32)
        agTGGGGAATATAGGACAATGGTCGAAAGCCTGATCCAGCCATGCCGCGTGTGTGAAGAAGGTCTTCGGATTGTAAAGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCAc                                                                                                                          <  4:311449/148‑1 (MQ=17)
              gAATACTGGACAAAGGGCGAAAGCCTGAACCAGCCAGGCCGCGTGTGTGAAGCAGGTCTTCGGATTGTACAGCACTTTCAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTaa                                                                                                                   <  8:339898/149‑1 (MQ=11)
                                  aaGCCTGATCCAGCCATGCCGCGTGTGTGAAGAAGGTCTTCGGATTGTAAAGCACTATAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTTACGTTACCGACAGAATAAGCAACGGCTAACTCTGTGCCAGCAGCCGCgg                                                                                               >  3:296882/1‑149 (MQ=12)
                                    gCCTGATCCAGCCATGCCGCGTGTGTGAAGAAGGTCTTCGGATTGTAAAGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGt                                                                                              >  5:220360/1‑148 (MQ=32)
                                         aTCCAGCCATGCCGCGTGTGTGAAGAAGGTCTTCGGATTGTAAAGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATa                                                                                          <  7:230175/147‑1 (MQ=32)
                                              gCCATGCCGCGTGTGTGAAGAAGGTCTTCGGATTGTAAAGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAggg                                                                                   <  6:220360/149‑1 (MQ=32)
                                                    ccGCGTGTGTGAAGAAGGTCTTCGGATTGTAAAGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAg                                                                             >  6:107541/1‑149 (MQ=32)
                                                      gcgTGTGTGAAGAAGGTCGTCGGATTGTAAAGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCg                                                                           >  2:43524/1‑149 (MQ=17)
                                                       cgTGTGTGAAGAAGGTCTTCGGATTGTAAAGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGt                                                                          <  3:213362/149‑1 (MQ=32)
                                                       cgTGTGTGAAGAAGGTCTTCGGATTGTAAAGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGt                                                                          <  1:43524/149‑1 (MQ=32)
                                                            gtgAAGAAGGTCTTCGGATTGTAAAGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATc                                                                     >  4:302434/1‑149 (MQ=32)
                                                            gtgAAGAAGGTCTTCGGATTGTAAAGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATc                                                                     >  8:138646/1‑149 (MQ=32)
                                                             tgAAGAAGGTCTTCGGATTGTAAAGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAActct                                                                                                               >  8:295218/1‑106 (MQ=32)
                                                             tgAAGAAGGTCTTCGGATTGTAAAGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAActct                                                                                                               <  7:295218/106‑1 (MQ=32)
                                                                 gaagGTCTTCGGATTGTAAAGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAAt                                                                <  3:104642/149‑1 (MQ=32)
                                                                 gaagGTCTTCGGATTGTAAAGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAAt                                                                <  6:38816/149‑1 (MQ=32)
                                                                     gTCTTCGGATTGTAAAGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACt                                                            >  8:116770/1‑149 (MQ=32)
                                                                       cTTCGGATTGTAAAGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTgg                                                          <  3:79744/149‑1 (MQ=32)
                                                                          cGGATTGTAAAGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTggg                                                         <  2:351510/147‑1 (MQ=32)
                                                                                 tAAAGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAgcg                                                >  5:360141/1‑149 (MQ=32)
                                                                                   aaGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAgcgcg                                              >  1:154453/1‑149 (MQ=32)
                                                                                   aaGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAgcgcg                                              >  2:174106/1‑149 (MQ=32)
                                                                                    aGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAgcgcgc                                             >  5:254430/1‑149 (MQ=32)
                                                                                     gCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGcgcgcg                                            >  2:193854/1‑149 (MQ=32)
                                                                                        cTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGc                                                                                <  4:180820/110‑1 (MQ=32)
                                                                                        cTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGc                                                                                >  3:180820/1‑110 (MQ=32)
                                                                                         tttAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAgg                                        >  5:51308/1‑149 (MQ=32)
                                                                                          ttAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAggt                                       >  2:81242/1‑149 (MQ=32)
                                                                                           tAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAggt                                       <  3:185231/148‑1 (MQ=32)
                                                                                            aaGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAggtg                                      <  3:302434/148‑1 (MQ=32)
                                                                                             aGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAggtggt                                    <  6:360141/149‑1 (MQ=32)
                                                                                               ttGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGGATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAgcgcgc                                             <  7:61185/138‑1 (MQ=17)
                                                                                               ttGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGGATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAgcgcgc                                             <  7:66381/138‑1 (MQ=17)
                                                                                               ttGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGGATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAgcgcgc                                             >  8:61185/1‑138 (MQ=17)
                                                                                               ttGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGGATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAgcgcgc                                             >  8:66381/1‑138 (MQ=17)
                                                                                               ttGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAGGTGGttt                                  <  7:217075/149‑1 (MQ=32)
                                                                                                  ggaggaAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAGGTGGTTTGtt                               <  7:116770/149‑1 (MQ=32)
                                                                                                  ggaggaAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAGGTGGTTTGt                                <  3:347875/148‑1 (MQ=32)
                                                                                                  ggaggaAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAGGTGGTTTGt                                >  7:23556/1‑148 (MQ=32)
                                                                                                    aggaAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAGGTGGTTTGTTaa                             <  2:154453/149‑1 (MQ=32)
                                                                                                     ggaAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAGGTGGTTTGTTAAg                            >  3:99062/1‑149 (MQ=32)
                                                                                                      gaACGGCAGTAAGTTAATACCTTGCGGTCTTGACGTTACCGACAGAATAAGCTCCGGCTAACTCTGTGCAAGCAGCCGCTGTTATACAGAGTGTTCAAGCGTTAATCGGAATTACTGGCCGTAAAGCGCGCGTAGGTGGTTTGTTAAGt                           <  6:112301/149‑1 (MQ=11)
                                                                                                          ggCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAGGTGGTTTGTTAAGTTGGa                       <  1:174106/149‑1 (MQ=32)
                                                                                                           gCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAGGTGGTTTGTTAAGTTGGAt                      <  7:138646/149‑1 (MQ=32)
                                                                                                           gCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAGGTGGTTTGTTAAGTTGGAt                      >  6:78847/1‑149 (MQ=32)
                                                                                                                   ttAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAGGTGGTTTGTTAAGTTGGATGTGAAAGc              >  4:328955/1‑149 (MQ=32)
                                                                                                                    tAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAGGTGGCTTGTTAAGTTGGATGTGAAAGcc             >  4:231668/1‑149 (MQ=17)
                                                                                                                       tACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAGGTGGTTTGTTAAGTTGGATGTGAAAGCCCCg          >  3:367676/1‑149 (MQ=32)
                                                                                                                        aCCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGTGTAGGTGGTTTGTTAAGTTGGATGTGAAAGCCCCgg         >  1:232064/1‑149 (MQ=17)
                                                                                                                        aCCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAGGTGGTTTGTTAAGTTGGATGTGaaa                <  2:254146/142‑1 (MQ=32)
                                                                                                                        aCCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAGGTGGTTTGTTAAGTTGGATGTGaaa                >  1:254146/1‑142 (MQ=32)
                                                                                                                               tGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAGGTGGTTTGTTAAGTTGGATGTGAAAGCCCCGGGCTCAAc  >  1:219263/1‑149 (MQ=32)
                                                                                                                                                   |                                                                                                                                
GAGGCAGCAGTGGGGAATATTGGACAATGGGCGAAAGCCTGATCCAGCCATGCCGCGTGTGTGAAGAAGGTCTTCGGATTGTAAAGCACTTTAAGTTGGGAGGAAGGGCAGTAAGTTAATACCTTGCTGTCTTGACGTTACCGACAGAATAAGCACCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAGGTGGTTTGTTAAGTTGGATGTGAAAGCCCCGGGCTCAAC  >  NC_002947/698150‑698425

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