Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 4,165,118 Δ1 bp intergenic (‑312/‑163) RPPX_18525 ← / → RPPX_18530 transporter/16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099744,165,1160A.93.9% 54.7 / 4.2 33intergenic (‑310/‑165)RPPX_18525/RPPX_18530transporter/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base . (16/15);  total (17/16)
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.96e-01

ACGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGACGAAAGGTTCTTAACCAAACGCTCTTTAACAAATCGAATCAAGCAATTC  >  CP009974/4164967‑4165148
                                                                                                                                                     |                                
aCGCGGTAGCGAAGCAAAGCTAGCCACTGTATTTAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCaa                                 >  2:1248455/1‑151 (MQ=255)
aCGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCaa                                 <  2:2240074/151‑1 (MQ=255)
 cgcgGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCaac                                <  1:4411707/151‑2 (MQ=255)
 cgcgGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCaac                                <  1:3934365/151‑2 (MQ=255)
  gcgGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCaacg                               <  2:3534747/151‑3 (MQ=255)
   cgGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCaacgc                              >  1:1197275/1‑148 (MQ=255)
    ggTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCaacgct                             >  2:718664/1‑147 (MQ=255)
    ggTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCaacgct                             >  1:2020886/1‑147 (MQ=255)
    ggTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCaacgct                             <  1:2323275/151‑5 (MQ=17)
    ggTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGATTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCaacgct                             >  2:1107594/1‑147 (MQ=17)
     gTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGctc                           >  1:3103275/1‑151 (MQ=17)
       aGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCtt                         <  1:2919467/151‑1 (MQ=14)
         cGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGTATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTa                       >  2:4565996/1‑151 (MQ=14)
         cGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTa                       >  1:3099084/1‑151 (MQ=14)
         cGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTa                       <  1:4542838/151‑1 (MQ=14)
          gAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTaa                      >  2:2757666/1‑151 (MQ=14)
           aaGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAAc                     <  1:4871865/151‑1 (MQ=14)
            aGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACa                    <  2:1969329/151‑1 (MQ=14)
             gCAAAGGTAGCCACTGGATTGAGTGAAGTGTGTGAAGCTAAGCGAGTTTCTCGCAAAAAACTGCAAGATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACaa                   <  2:5211946/151‑1 (MQ=11)
             gCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACaa                   <  2:2010689/151‑1 (MQ=14)
              cAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACaaa                  <  1:1711507/151‑1 (MQ=14)
              cAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACaaa                  >  1:1754546/1‑151 (MQ=14)
               aaaGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAAt                 >  1:817184/1‑151 (MQ=14)
                aaGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATc                >  2:1029447/1‑151 (MQ=14)
                aaGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATc                <  2:5091906/151‑1 (MQ=14)
                 aGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCg               <  1:2123667/151‑1 (MQ=11)
                 aGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCg               <  2:2010591/151‑1 (MQ=11)
                 aGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCg               <  1:3403971/151‑1 (MQ=14)
                 aGCTAGCAACTGTATTGAGTAAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCg               <  2:4279832/151‑1 (MQ=11)
                  gCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGa              <  1:3726740/151‑1 (MQ=11)
                  gCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGa              >  2:714126/1‑151 (MQ=11)
                   cTAGCCACTGTATAGAGACAAGAGATTGAAGCTAAACGAGTTTCTAGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGaa             >  2:3577941/1‑151 (MQ=11)
                    tAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAAt            >  1:915195/1‑151 (MQ=11)
                    tAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAAt            >  1:636292/1‑151 (MQ=11)
                    tAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAAt            >  1:2612827/1‑151 (MQ=11)
                     aGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAATc           <  2:4453255/151‑1 (MQ=11)
                      gCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAATCa          <  2:3746211/151‑1 (MQ=11)
                         aCTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAATCAAGc       >  1:1199721/1‑151 (MQ=11)
                           tGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAATCAAGCaa     >  2:4758951/1‑151 (MQ=11)
                           tGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCGCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAATCAAGCaa     >  2:246453/1‑151 (MQ=11)
                            gaagcgAGTTAAGTGACTGAAGTTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGACGAAAGGTTCTTAACCAAACGCTCTTTAACAAATCGAATCAAGCaa     >  2:4206179/6‑151 (MQ=1)
                             tATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAATCAAGCAAtt   >  1:1367519/1‑151 (MQ=255)
                               gcgAGTTAAGTGACTGAAGTTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGACGAAAGGTTCTTAACCAAACGCTCTTTAACAAATCGAATCAAGCAATTc  <  1:3624635/149‑1 (MQ=11)
                                                                                                                                                     |                                
ACGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGACGAAAGGTTCTTAACCAAACGCTCTTTAACAAATCGAATCAAGCAATTC  >  CP009974/4164967‑4165148

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