Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,488,992 C→G 8.8% intergenic (+345/+226) PP_3981 → / ← PP_3982 transposase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,488,9920CG 8.8% 167.4 / 7.1 57intergenic (+345/+226)PP_3981/PP_3982transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (24/28);  new base G (2/3);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

TGCCCAAGGCTCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGTTTCACCACC  >  NC_002947/4488848‑4489135
                                                                                                                                                |                                                                                                                                               
tGCCCAAGGCTCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCTAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGTAGCTGAAGGATTAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTc                                                                                                                                              <  1:72356/148‑1 (MQ=255)
  cccAAGGCTCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTccc                                                                                                                                            <  3:193909/148‑1 (MQ=255)
   ccAAGGCTCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGc                                                                                                                                          >  3:342062/1‑149 (MQ=255)
    cAAGGCTCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGcc                                                                                                                                         <  4:56008/149‑1 (MQ=255)
     aaGGCTCGACATTGCGCCATGACGACTGGAAACAGCAGGTTGTCGAAAATCCACCGATCGCCAACATCAAGCATCGCGCCCACCGTGCAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAGGTTCCCGccc                                                                                                                                        <  2:153266/149‑1 (MQ=255)
       ggCTCGACATTGCGCCATGACGACTGGAAAATGCAGGTAGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAGGTTCCCGcccgt                                                                                                                                      <  3:329330/149‑3 (MQ=255)
         cTCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCatat                                                                                                                                    >  6:173689/1‑149 (MQ=255)
            gACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATgc                                                                                                                                  <  5:173689/148‑1 (MQ=255)
             aCATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATgcgc                                                                                                                                <  4:169658/149‑1 (MQ=255)
                ttGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTa                                                                                                                             <  8:40032/149‑1 (MQ=255)
                ttGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTa                                                                                                                             <  5:293397/149‑1 (MQ=255)
                  gcgcCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTAc                                                                                                                            <  8:12949/148‑1 (MQ=255)
                       aTGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGTGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATa                                                                                                                      <  7:342755/149‑1 (MQ=255)
                           cgacTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATACAGGGCGTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAgcgc                                                                                                                  <  2:151558/149‑1 (MQ=255)
                            gacTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGcgcg                                                                                                                 >  7:290405/1‑149 (MQ=255)
                                 gAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAAt                                                                                                            >  7:60331/1‑149 (MQ=255)
                                  aaaaTGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATc                                                                                                           <  2:184614/149‑1 (MQ=255)
                                  aaaaTGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATc                                                                                                           <  4:342062/149‑1 (MQ=255)
                                        cAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTc                                                                                                      <  1:130705/148‑1 (MQ=255)
                                               tCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGAc                                                                                                <  7:109609/147‑1 (MQ=255)
                                               tCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGAc                                                                                                >  8:109609/1‑147 (MQ=255)
                                                 gAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTca                                                                                            <  8:60331/149‑1 (MQ=255)
                                                  aaaaTCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAAGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTACCGCCCATATgct                                                                                                                                 >  1:60720/1‑110 (MQ=255)
                                                            gATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTg                                                                                 <  1:319188/149‑1 (MQ=255)
                                                             aTCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTgg                                                                                <  7:126684/149‑1 (MQ=255)
                                                              tCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTg                                                                                 <  1:314351/147‑1 (MQ=255)
                                                              tCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGCGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGc                                                                               <  5:206327/149‑1 (MQ=255)
                                                                 ccAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACATGGGCGGGAGGTTCCCGCCCGTGTgcaa                                                                                                                                <  6:249589/97‑3 (MQ=37)
                                                                 ccAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACATGGGCGGGAGGTTCCCGCCCGTGTgcaa                                                                                                                                >  5:249589/1‑95 (MQ=37)
                                                                           gCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTgg                                                                                >  5:140777/1‑135 (MQ=255)
                                                                           gCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTgg                                                                                <  6:140777/135‑1 (MQ=255)
                                                                                  gcCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAActca                                                           >  2:248451/1‑148 (MQ=255)
                                                                                    cctgCGTGGAGCTGTAGGATAAAGGGCTTATCTGAGGGCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAActc                                                            <  8:21246/143‑1 (MQ=255)
                                                                                    ccAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAActc                                                            >  7:21246/1‑146 (MQ=255)
                                                                                        cGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCACCATAg                                                     <  1:248451/149‑1 (MQ=255)
                                                                                           ggAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGt                                                   >  5:336032/1‑148 (MQ=255)
                                                                                                 gAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCt                                            >  6:220882/1‑149 (MQ=255)
                                                                                                  acggATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTa                                                                                                                             <  2:199241/65‑1 (MQ=255)
                                                                                                   aGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGa                                          >  7:151738/1‑149 (MQ=255)
                                                                                                     gATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACg                                        >  5:45919/1‑149 (MQ=255)
                                                                                                     gATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACg                                        >  7:168575/1‑149 (MQ=255)
                                                                                                         aaGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCt                                     <  5:220882/148‑1 (MQ=255)
                                                                                                            ggCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCca                                 >  5:212634/1‑149 (MQ=255)
                                                                                                                    tGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCa                            <  2:162403/146‑1 (MQ=255)
                                                                                                                    tGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCa                            >  1:162403/1‑146 (MQ=255)
                                                                                                                    tGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATa                         >  5:338740/1‑149 (MQ=255)
                                                                                                                    tGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCACACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATa                         >  6:143283/1‑149 (MQ=255)
                                                                                                                      aGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCaa                           <  8:290405/145‑1 (MQ=255)
                                                                                                                        gCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGa                     <  5:143283/149‑1 (MQ=255)
                                                                                                                          cttgcttgCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAgg                   >  4:360994/1‑149 (MQ=255)
                                                                                                                           ttgcttgcACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAgg                   >  3:82056/1‑148 (MQ=255)
                                                                                                                               ttgcACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAAc               >  5:345108/1‑148 (MQ=255)
                                                                                                                                   atatgggcgggAGGTTCCCGCCCGTATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGCCGCCTGCCACATCAATACCGAGGTAACGAGtt          >  5:273303/5‑149 (MQ=255)
                                                                                                                                      cgggcgggAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAACCTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGTTTca       >  7:342909/1‑149 (MQ=255)
                                                                                                                                         gcgggAACCTCCCGCCCATATGCGCGTAATACATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTCATGTCTGACGCCTGCCACATCAATACCGAGGTAACCAGTTTcacca    >  1:88179/1‑149 (MQ=255)
                                                                                                                                          cgggAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGTTTcaccac   <  6:336032/149‑1 (MQ=255)
                                                                                                                                           gggAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGTTTCaccacc  >  2:221787/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                               
TGCCCAAGGCTCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGTTTCACCACC  >  NC_002947/4488848‑4489135

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: ‑