Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 4,031,017 Δ1 bp intergenic (+159/‑625) RPPX_17905 → / → RPPX_17915 23S rRNA methyltransferase/sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099744,031,0170A.93.5% 84.6 / 3.7 32intergenic (+159/‑625)RPPX_17905/RPPX_1791523S rRNA methyltransferase/sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base . (14/15);  total (15/17)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGG  >  CP009974/4030871‑4031167
                                                                                                                                                  |                                                                                                                                                       
gtTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGccc                                                                                                                                                    >  1:131616/1‑151 (MQ=255)
gtTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGccc                                                                                                                                                    >  1:126644/1‑151 (MQ=16)
 ttctgcttctgcttctgcttctgcttctgcttctgcTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGcccc                                                                                                                                                   >  2:3737682/1‑151 (MQ=16)
 ttctgcttctgcttctgcttctgcttctgcttctgcTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGcccc                                                                                                                                                   <  1:4167019/151‑1 (MQ=255)
  tctgcttctgcttctgcttctgcttctgcttctgctTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCa                                                                                                                                                  <  1:4923055/151‑1 (MQ=16)
  cttgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCTGGCGATGCAGATTGCGAGTTCAGGAGGCCGAGCGCAGGGGTTGCGGAGGGAGGGGACGGGCATGGATGCCCGGCAAGCGCTGgggctcg                                                                                                                                                   >  1:3419376‑M2/33‑148 (MQ=255)
   ctgcttctgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAg                                                                                                                                                 <  2:433220/151‑1 (MQ=255)
      cttctgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGAt                                                                                                                                              >  1:166489/1‑151 (MQ=255)
      cttctgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGAt                                                                                                                                              <  1:3848694/151‑1 (MQ=255)
       ttctgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATg                                                                                                                                             >  2:536949/1‑151 (MQ=16)
       ttctgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATg                                                                                                                                             >  2:1409882/1‑151 (MQ=16)
        tctgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATgg                                                                                                                                            <  1:3076321/151‑1 (MQ=255)
        tctgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATgg                                                                                                                                            <  1:4738945/151‑1 (MQ=255)
        tctgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATgg                                                                                                                                            <  1:3076308/151‑1 (MQ=255)
        tctgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATgg                                                                                                                                            <  2:5182366/151‑1 (MQ=16)
         ctgcgtctgcttctgctgctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATggg                                                                                                                                           <  1:3809995/151‑1 (MQ=255)
          tgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATgggg                                                                                                                                          >  2:3607263/1‑151 (MQ=19)
             ttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCt                                                                                                                                       >  1:4402058/1‑151 (MQ=19)
             ttctgcttctgctgctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCt                                                                                                                                       <  2:3313541/151‑1 (MQ=12)
              tctgcttctgctgctgcttctgcttTTGCTGCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCtg                                                                                                                                      <  2:3563485/151‑2 (MQ=255)
                 gcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCtgcag                                                                                                                                   >  2:1471906/1‑147 (MQ=9)
                 gcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCtgcag                                                                                                                                   >  2:776108/1‑147 (MQ=9)
                 gcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCtgcag                                                                                                                                   >  1:3774768/1‑147 (MQ=9)
                 gcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCtgcag                                                                                                                                   <  1:3440283/151‑5 (MQ=255)
                  cttctgcttctgcttgtggttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGGGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCtgcagc                                                                                                                                  <  1:3563508/151‑6 (MQ=11)
                  cttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGTGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCtgcagc                                                                                                                                  >  2:2445136/1‑146 (MQ=11)
                  cttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCtgcagc                                                                                                                                  >  2:4459029/1‑146 (MQ=9)
                  cttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCtgcagc                                                                                                                                  <  2:929555/151‑6 (MQ=0)
                  cttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCtgcagc                                                                                                                                  >  1:4979674/1‑146 (MQ=255)
                                         gCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGccc                                                                                                           <  2:247470/151‑1 (MQ=11)
                                                                  cGACGCAGATTGCGACTGCGGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGGGACGGGCAGGGATGCCCGTCAAGCGCGGGGCCCCAGGACGGGGCCTGCAGCGCGGTCCGCCCGGGAGCAAGCCCGGAGCCGAGGGGACCCCGgcgcgcc                                                                                  >  2:389343/1‑150 (MQ=255)
                                                                                                           gggggggAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTa                                         <  2:1102519/147‑1 (MQ=1)
                                                                                                                          ggtaTGGATGCCCGTCAAGCGGTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGttt                          <  1:2453876/148‑1 (MQ=1)
                                                                                                                                                  aGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAgg  <  1:3838757/151‑1 (MQ=31)
                                                                                                                                                  |                                                                                                                                                       
GTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGG  >  CP009974/4030871‑4031167

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