Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 527,898 Δ1 bp 25.0% noncoding (1133/2902 nt) PP_23SC → 23S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947527,8980C.25.0% 109.5 / 41.2 48noncoding (1133/2902 nt)PP_23SC23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (17/19);  new base . (4/8);  total (21/27)
Fisher's exact test for biased strand distribution p-value = 5.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

TAGACAGCTAGGAGGTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGC  >  NC_002947/527801‑528033
                                                                                                 |                                                                                                                                          
tAGACAGCTAGGAGGTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑Cgatg                                                                                                   <  1:299042/136‑1 (MQ=32)
   aCAGCTAGGAGGTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGatgatg                                                                                                <  1:419982/136‑1 (MQ=32)
         acgatgTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCAACGGGTATCATCTTA‑‑CGATGATGCGGTag                                                                                          >  1:460286/6‑136 (MQ=12)
          ggatgTTGGCTTAGAAGCAGCCATCATTTAAAGAAAGCGTAATAGCTCACTGGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTaga                                                                                         >  1:642119/1‑136 (MQ=255)
          ggatgTTGGCTTAGAAGCAGCCATCATTTAAAGAAAGCGTAATAGCTCACTGGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTaga                                                                                         >  1:627790/1‑136 (MQ=255)
          ggatgTTGGCTTAGAAGCAGCCATCATTTAAAGAAAGCGTAATAGCTCACTGGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑Cgatg                                                                                                   <  2:160880/125‑1 (MQ=11)
          ggatgTTGGCTTAGAAGCAGCCATCATTTAAAGAAAGCGTAATAGCTCACTGGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑Cgatg                                                                                                   >  1:160880/1‑125 (MQ=11)
          ggatgTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTaga                                                                                         >  1:163685/1‑136 (MQ=17)
          ggatgTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTaga                                                                                         >  1:226471/1‑136 (MQ=17)
          ggaggTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTaga                                                                                         >  1:562375/1‑136 (MQ=32)
                tGGCTTAGAAGCAGCCATCATTTAAAGAAAGCGTAATAGCTCACTGGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCg                                                                                   >  1:117511/1‑136 (MQ=12)
                    ttAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCt                                                                               >  1:135122/1‑136 (MQ=32)
                       taaGCAGCCACCCTTTAAAGAAAGCGTATTAGCTCACTAGTCGAGTCGGGCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTa                                                                            <  1:79475/135‑1 (MQ=12)
                         agcagcCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACAGGTATCATCTTA‑‑CGATGATGCTGTACAGGATCGTTCTGTAAg                                                                          >  2:270513/1‑136 (MQ=255)
                            agcCATCATTTAAAGAAAGCGTAATAGCTCACTGGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCt                                                                       >  1:273469/1‑136 (MQ=12)
                                aCCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGa                                                                   >  2:460001/1‑136 (MQ=32)
                                aCCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGa                                                                   >  2:57000/1‑136 (MQ=32)
                                           aaagCGTAATAGCTCACTGGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGCGAAGGTgtgctg                                                         >  1:303388/1‑131 (MQ=0)
                                            aagCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAATAATGTGAAGGTGAGTTgag                                                       <  1:135315/136‑1 (MQ=12)
                                                gTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGc                                                   >  1:281630/1‑136 (MQ=32)
                                                     aGCTCACTGGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCt                                              <  2:117511/136‑1 (MQ=17)
                                                     aGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCt                                              <  2:579762/136‑1 (MQ=32)
                                                     aGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCt                                              >  1:391247/1‑136 (MQ=32)
                                                        tCACTGGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGCGAAGGTGTGCTGTGAGGCATGCTgg                                            >  1:551397/1‑136 (MQ=0)
                                                                    tcgGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGt                               >  1:460044/1‑136 (MQ=32)
                                                                    tcgGCCTGCGCGGAAGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGCGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAg                                >  2:650262/1‑136 (MQ=0)
                                                                         cTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGaa                          >  2:354611/1‑136 (MQ=32)
                                                                              cgGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTg                     <  1:140452/136‑1 (MQ=32)
                                                                                gAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGAc                   <  2:661182/136‑1 (MQ=32)
                                                                                gAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGAc                   <  2:454076/136‑1 (MQ=32)
                                                                                 aaGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGCGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGAc                   <  2:30331/136‑1 (MQ=0)
                                                                                  tgATGTAACGGGGCTCAAACCTTGCTCCGAAGGTACGGGTATGATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                 <  2:627790/135‑1 (MQ=255)
                                                                                  aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCTTCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                 <  2:358307/136‑1 (MQ=17)
                                                                                  aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                 <  2:511246/136‑1 (MQ=32)
                                                                                  aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                 <  2:163685/136‑1 (MQ=32)
                                                                                  aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                 <  2:546473/136‑1 (MQ=32)
                                                                                  aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                 <  2:416454/136‑1 (MQ=32)
                                                                                  aGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACa                  <  2:544602/136‑1 (MQ=0)
                                                                                   gATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGTATGCTGACATg                <  2:415240/136‑1 (MQ=17)
                                                                                   gATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATg                <  2:562375/136‑1 (MQ=32)
                                                                                   gATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATg                <  2:226471/136‑1 (MQ=32)
                                                                                   gATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                 <  2:623273/136‑1 (MQ=0)
                                                                                   gATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGCGAAGGTGTGCTGTGAGGCATGCTGGATGTATCAGAAGTGCGAATGCTGACAt                 <  2:24214/136‑1 (MQ=0)
                                                                                   gATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGCGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                 <  2:231097/136‑1 (MQ=0)
                                                                                   gATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGCGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                 <  2:21523/136‑1 (MQ=0)
                                                                                    aTGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGCGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACAta                <  2:233666/136‑2 (MQ=0)
                                                                                              gCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACaa     <  2:339759/136‑1 (MQ=32)
                                                                                                tCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATg   >  1:262908/1‑136 (MQ=32)
                                                                                                 cAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGc  >  2:262715/1‑136 (MQ=32)
                                                                                                 |                                                                                                                                          
TAGACAGCTAGGAGGTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCATCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGC  >  NC_002947/527801‑528033

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