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.91.7% 109.3 / 8.2 36intergenic (+159/‑625)RPPX_17905/RPPX_1791523S rRNA methyltransferase/sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base . (17/16);  total (18/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGG  >  CP009974/4030871‑4031167
                                                                                                                                                  |                                                                                                                                                       
gtTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGccc                                                                                                                                                    >  1:3106084/1‑151 (MQ=255)
gtTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGccc                                                                                                                                                    >  1:3681641/1‑151 (MQ=16)
gtTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGccc                                                                                                                                                    >  2:270408/1‑151 (MQ=16)
 ttctgcttctgcttctgcttctgcttctgcttctgcTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGcccc                                                                                                                                                   <  1:697406/151‑1 (MQ=255)
 ttctgcttctgcttctgcttctgcttctgcttctgcTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGcccc                                                                                                                                                   <  1:852286/151‑1 (MQ=255)
  tctgcttctgcttctgcttctgcttctgcttctgctTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCa                                                                                                                                                  <  2:5059353/151‑1 (MQ=16)
   ttgcttcggattctgcttctgcttctgcttctgcttTTGCTTCTCAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAg                                                                                                                                                <  1:2073309/150‑1 (MQ=11)
    tgcttctgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGCCGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAgg                                                                                                                                               >  2:2142805/1‑151 (MQ=255)
     gcttctgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGa                                                                                                                                              >  1:1763211/1‑151 (MQ=16)
     gctgctgcttctgcttctgcttctgcttctgcttTTGCTTCTCAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGACGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGa                                                                                                                                              <  2:3481242/151‑1 (MQ=11)
       ttctgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATg                                                                                                                                            >  2:1052902/1‑151 (MQ=16)
       ttctgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATg                                                                                                                                            <  2:3049066/151‑1 (MQ=16)
       ttctgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATg                                                                                                                                            <  1:5125066/151‑1 (MQ=16)
         ctgcttctgcttctgcttctgcttctgcttTTGCTTCTAGGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATggg                                                                                                                                          <  1:5295473/151‑1 (MQ=255)
         ctgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATggg                                                                                                                                          >  1:323683/1‑151 (MQ=255)
         ctgcttctgcttcggcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATggg                                                                                                                                          <  1:3858899/151‑1 (MQ=255)
           gcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGc                                                                                                                                        <  2:1763211/151‑1 (MQ=16)
            cttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGcc                                                                                                                                       >  2:317539/1‑151 (MQ=255)
            cttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGcc                                                                                                                                       >  2:587619/1‑151 (MQ=12)
             ttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCt                                                                                                                                      <  2:714098/151‑1 (MQ=16)
             ttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCt                                                                                                                                      >  1:1211761/1‑151 (MQ=255)
             ttcggcttctgcttcggcttctgcttTTGCGTCAACGCGCGCGATGGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGCGCGCAGGGCTGGCGAAGGGGGGTGACGGGCTTTGATGCCCGTCACGCGCTGAGGACCCATGGATGGGGc                                                                                                                                        <  2:2858868/151‑1 (MQ=11)
              tctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCtg                                                                                                                                     >  2:3273285/1‑150 (MQ=16)
              tctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCtg                                                                                                                                     >  1:1182011/1‑150 (MQ=255)
                tgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGCGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCtgca                                                                                                                                   <  1:4093380/151‑4 (MQ=255)
                tgcgtctgcttctggttctgcttGTGCTTCTAGGGACGCGATAGGTGGGGCGACTCAGTTTGCTGCTTGAGGAGGCGGAGTGCAGGGCTTGTGAAGGGAGGTGACGTGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCtgca                                                                                                                                   <  2:4349027/151‑4 (MQ=2)
                  cttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCtgcagc                                                                                                                                 <  2:3074794/151‑6 (MQ=2)
                  cttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCtgcagc                                                                                                                                 >  1:3338203/1‑146 (MQ=255)
                   ttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCtgcagcg                                                                                                                                <  1:3939573/151‑7 (MQ=255)
                   ttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCtgcagcg                                                                                                                                >  2:346277/1‑145 (MQ=255)
                   ttctgctgcggcttcggcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCtgcagcg                                                                                                                                <  2:559177/151‑7 (MQ=1)
                      tgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCgg                                                                                                                             >  1:4216857/1‑151 (MQ=255)
                       gcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGt                                                                                                                            >  1:1951132/1‑151 (MQ=255)
                          tctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCt                                                                                                                         >  1:399776/1‑151 (MQ=255)
                                                                                                                                                 gAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAg   <  1:850631/151‑1 (MQ=11)
                                                                                                                                                 gAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAg   >  2:4130074/1‑151 (MQ=11)
                                                                                                                                                  aGGCCGGAGGATGGGGCGTGCGGGGCGGTCCTGCCGGGAGGAAGGCCGGAGTGAGGGGACGCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAgg  <  2:3162067/151‑1 (MQ=11)
                                                                                                                                                  |                                                                                                                                                       
GTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGG  >  CP009974/4030871‑4031167

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