Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,085,591 A→G 11.3% intergenic (+300/‑295) PP_1864 → / → PP_1865 hypothetical protein/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,085,5910AG11.3% 112.5 / 5.9 44intergenic (+300/‑295)PP_1864/PP_1865hypothetical protein/transposase
Reads supporting (aligned to +/- strand):  ref base A (23/16);  new base G (3/2);  total (26/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 = 4.81e-01

GGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTT  >  NC_002947/2085449‑2085738
                                                                                                                                               |                                                                                                                                                       
gggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGAcc                                                                                                                                                    >  6:272110/1‑148 (MQ=34)
gggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCa                                                                                                                                                   <  5:239485/149‑1 (MQ=34)
      ggATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGCCCAAAGTaa                                                                                                                                             >  8:88720/1‑149 (MQ=21)
      ggATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGAGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTaa                                                                                                                                             >  6:54473/1‑149 (MQ=21)
       gATGCCCGTCAAGCGGTGGGCCCCTGGAGGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTaa                                                                                                                                             <  2:159237/148‑1 (MQ=255)
       gATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGGGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAAc                                                                                                                                            <  2:216530/149‑1 (MQ=21)
       gATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGTCGTGGCCAAAGTaa                                                                                                                                             <  5:98845/148‑1 (MQ=255)
       gATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAAc                                                                                                                                            >  8:141138/1‑149 (MQ=34)
           cccGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACc‑‑cgc                                                                                                                                      <  1:56044/149‑1 (MQ=34)
            ccGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACc‑‑cgcc                                                                                                                                     >  6:132445/1‑149 (MQ=34)
                           ccccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAgg                                                                                                                      >  2:140390/1‑149 (MQ=34)
                               aGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGaa                                                                                                                  >  8:223550/1‑149 (MQ=34)
                                          tGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTg                                                                                                        >  6:80471/1‑148 (MQ=35)
                                              gcgcGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCg                                                                                                   <  5:272110/149‑1 (MQ=37)
                                                       tccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGgcgc                                                                                                           <  4:212351/132‑1 (MQ=34)
                                                       tccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGgcgc                                                                                                           >  3:212351/1‑132 (MQ=34)
                                                                 aaGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑GGCCGGAAGGGCGGAAAGGTGAATAGGCGGCTGCCTCGTGAATGAAGGTCTACATGa                                                                                >  6:224594/1‑149 (MQ=11)
                                                                 aaGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCAGGATGAGGGGAGCCGTGCGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACtt                                                                                  >  1:208986/1‑148 (MQ=255)
                                                                             aGGGGACCCCGGAGCGTAGCTTAGGGGCCGGATGTGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAgcgc                                                                    <  7:141138/148‑1 (MQ=255)
                                                                                 gACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAAc                                                               >  4:225755/1‑149 (MQ=39)
                                                                                      cGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCaaa                                                          <  6:77113/149‑1 (MQ=39)
                                                                                      cGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCa                                                            >  6:112156/1‑147 (MQ=39)
                                                                                       ggAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAg                                                         <  5:132445/149‑1 (MQ=40)
                                                                                        gAGCGTAGCGTTGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAgg                                                        <  8:240081/149‑1 (MQ=255)
                                                                                                       gCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAGCCCAAAGGCCcaaaagctcaa                                          >  5:44351/1‑147 (MQ=255)
                                                                                                              gAGGGGACCCGACGTGTTTGGTTACTTTGCCCACGACCAAATTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTa                                  <  7:88720/148‑1 (MQ=14)
                                                                                                                      ccGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGtt                         >  5:150776/1‑149 (MQ=255)
                                                                                                                       cGACGGCTTTGGTTACTTTGGCCACGACCAAATTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTc                        <  2:208986/149‑1 (MQ=255)
                                                                                                                        gACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCt                       >  4:139788/1‑149 (MQ=255)
                                                                                                                            ggCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCCTCAAAAGCTTTAGCGCCGGTTCTGa                     <  1:140390/149‑1 (MQ=255)
                                                                                                                              cTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCcaaaagctcaaaagct                                     <  5:142768/130‑1 (MQ=255)
                                                                                                                              cTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCcaaaagctcaaaagct                                     >  6:142768/1‑130 (MQ=255)
                                                                                                                               gcgggTTACTTTGGTCGTGGCCAAAGTAACCAAAGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATc                   >  4:141510/4‑149 (MQ=17)
                                                                                                                                 gggTTACTTTGGTCGTGGCCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACtt                                                                                  >  1:150933/2‑84 (MQ=11)
                                                                                                                                    ttACTTTGGTCGTGGCCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTCGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAgcgc                              <  1:191641/132‑1 (MQ=14)
                                                                                                                                    ttACTTTGGTCGTGGCCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAgcgc                              >  2:191641/1‑131 (MQ=255)
                                                                                                                                    ttACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGGCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCCCCGGTTCTGATCGGGCGTa            >  2:27257/1‑149 (MQ=255)
                                                                                                                                       cTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCgg                 >  4:152065/1‑141 (MQ=255)
                                                                                                                                       cTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCgg                 <  3:152065/141‑1 (MQ=255)
                                                                                                                                          tGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTcc      >  3:101552/1‑149 (MQ=255)
                                                                                                                                           gcccACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTcc      <  7:223550/146‑1 (MQ=255)
                                                                                                                                            gCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGt    >  2:4575/1‑149 (MQ=255)
                                                                                                                                             ccACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGtt   >  3:77705/1‑149 (MQ=255)
                                                                                                                                              cACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAACGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGttt  >  8:101717/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                       
GGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTT  >  NC_002947/2085449‑2085738

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