Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,406,852 G→T 24.2% G58C (GGC→TGC)  dbpA → ATP‑dependent RNA helicase, specific for 23S rRNA

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,406,8520G→T24.2% 4.1 ‑4.4 0.033172dbpAATP‑dependent RNA helicase, specific for 23S rRNA
Reads supporting (aligned to +/- strand):  new base (9/0):  ref base (12/12):  total (21/12)
Fisher's exact test strand distribution p-value = 1.22e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GGCGAAAACCGGCAGCGGCAAAACGGCGGCTTTTGGCCTCGGCTTGTTACAGCAAATTGATGCGTCGCTA  >  REL606/1406818‑1406887
                                  |                                   
ggCGAAAACCGGCAGCGTCAAAACGGCGGCTTTTgg                                    <  1:582232/36‑1 (MQ=255)
ggCGAAAACCGGCAGCGGCAAAACGGCGGCTTTTgg                                    >  1:7436617/1‑36 (MQ=255)
ggCGAAAACCGGCAGCGGCAAAACGGCGGCCTTTgg                                    >  1:1617598/1‑36 (MQ=255)
 gCGAAAACCGGCAGCGGCAAAACGGCGGCTTTTGGc                                   >  1:1393337/1‑36 (MQ=255)
   gAAAACCGGCAGCGGCAAAACGGCGGGTTTTTGccc                                 >  1:5351558/1‑35 (MQ=21)
    aagaCCGGCAGCGGCAAAACGGCGGCTTTTGGCCTc                                <  1:4864917/33‑1 (MQ=255)
    aaaaCCGGCAGCGGCAAAACGGCGGTTTTATGCCTc                                >  1:149931/1‑36 (MQ=18)
       aCCGGCAGCGGCAAAACGGCGGCTTTTTGCCTCGGc                             >  1:5003047/1‑36 (MQ=255)
       aCCGGCAGCGGCAAAACGGCGGCTTCTGGCCTCGGc                             <  1:2499715/36‑1 (MQ=255)
         cGGCAGCGGCAAAACGGCGGCTTTTGGCCTCGGCtt                           <  1:6460559/36‑1 (MQ=255)
          ggCAGCGGCAAAACGGCGGCTTTTTGCCGCGGATTg                          >  1:7905400/1‑36 (MQ=21)
          ggCAGCGGCAAAACGGCGGCTTTTGGCCTCGGCTTg                          <  1:2432026/36‑1 (MQ=255)
          ggCAGCGGCAAAACGGCGGCTTTTGGCCTCGGCTTg                          <  1:2129906/36‑1 (MQ=255)
           gCAGCGGCAAAACGGCGGCTTTTTGCCTCGGCtttt                         >  1:2724529/1‑34 (MQ=37)
           gCAGCGGCAAAACGGCGGCTTTTTGCCTCGGCTTGt                         >  1:5844583/1‑36 (MQ=255)
             aGCGGCAAAACGGCNNNTTTTTGCCTCGGCTNNTTa                       >  1:7333375/1‑36 (MQ=21)
             aGCGGCAAAACGGCGGCTTTTTGCCTCGGCTTTttc                       >  1:6795989/1‑35 (MQ=37)
             aGCGGCAAAACGGCGGCTTTTTGCCTCGGCTTGttt                       >  1:2309685/1‑35 (MQ=255)
             aGCGGCAAAACGGCGGCTTTTGGCCTCGGCTTTttt                       >  1:6873722/1‑35 (MQ=255)
              gCGGCAAAACGGCGGCTTTTGGCCTCGGCTTGTTAc                      <  1:3934382/36‑1 (MQ=255)
              gCGGCAAAACGGCGGCTTTTGGCCTCGGCTTGTTAc                      <  1:4022656/36‑1 (MQ=255)
                ggCAAAACGGCGGCTTTTGGCCTCGGCTTGTTACat                    >  1:919553/1‑35 (MQ=255)
                 gCAAAACGGCGGCTTTTGGCCTCGGCTTTTTACAGc                   >  1:908059/1‑36 (MQ=255)
                 gCAAAACGGCGGCTTGTGGCCTCGGCTTGTTTCAGc                   >  1:318212/1‑36 (MQ=37)
                       cggcggCTTTTGGCCTCGGCTTGTTACAGCAAATTg             <  1:7470432/36‑1 (MQ=255)
                         gcggcTTTTGGCCTCGGCTTGTTACAGCAAATTGAt           >  1:1081861/1‑36 (MQ=255)
                          cggcTTTTGGCCTCGGCTTGTTACAGCAAATTGATg          <  1:6843233/36‑1 (MQ=255)
                            gcTTTTGGCCTCGGCTTGTTACAGCAAATTGATGCg        >  1:326376/1‑36 (MQ=255)
                             cTTTTGGCCTCGGCTTGTTACAGCAAATTGATGCGt       <  1:513510/36‑1 (MQ=255)
                             cTTTTGGCCTCGGCTTGTTACAGCAAATTGATGCGt       <  1:4415163/36‑1 (MQ=255)
                             cTTTTGGCCTCGGCTTGTTACAGCAAATTGATGCGt       >  1:7076380/1‑36 (MQ=255)
                             cTTTTGGCCTCGGCTTGTTACAGCAAATTGATGCGt       >  1:3562299/1‑36 (MQ=255)
                              ttttGGCCTCGGCTTGTTACAGCAAATTGATGCGTc      >  1:544319/1‑36 (MQ=255)
                              ttttGGCCTCGGCTTGTTACAGCAAATTGATGCGTc      >  1:4632122/1‑36 (MQ=255)
                                ttGGCCTCGGCTTGTTACAGCAAATTGATGCGTCGc    >  1:4699558/1‑36 (MQ=255)
                                 tGGCCTCGGCTTGTTACAGCAAATTGATGCGTCGCt   <  1:7884072/36‑1 (MQ=255)
                                  ggCCTCGGCTTGTTACAGCAAATTGATGCGTCGCTa  <  1:3404505/36‑1 (MQ=255)
                                  ggCCTCGGCTTGTTACAGCAAATTGATGCGTCGCTa  <  1:2371049/36‑1 (MQ=255)
                                  ggCCTCGGCTTGTTACAGCAAATTGATGCGTCGCTa  >  1:1235222/1‑36 (MQ=255)
                                  |                                   
GGCGAAAACCGGCAGCGGCAAAACGGCGGCTTTTGGCCTCGGCTTGTTACAGCAAATTGATGCGTCGCTA  >  REL606/1406818‑1406887

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 19 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG