Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,971,753 C→T 7.8% intergenic (+136/‑244) PP_3501 → / → PP_3502 transposase/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,971,7530CT 7.8% 144.2 / 2.3 51intergenic (+136/‑244)PP_3501/PP_3502transposase/transposase
Reads supporting (aligned to +/- strand):  ref base C (23/24);  new base T (2/2);  total (25/26)
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.87e-01

CCTGCCTAACTTACGCAAGGTGAGTTGGGCGGACGCTTACTTCCTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTCGCGCAGTACAGACCGTAGGAGAAATCACCCCTCCTCCTTTCAACACCC  >  NC_002947/3971609‑3971895
                                                                                                                                                 |                                                                                                                                              
cctgcctAACTTACGCAAGGTGAGTTGGGCGGACGCTTACTTCCTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTg                                                                                                                                            <  2:38273/149‑1 (MQ=255)
    cctAACTTACGCAAGGTGAGTTGGGCGGACGCTTACTTCCTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCt                                                                                                                                         >  2:246954/1‑148 (MQ=255)
    cctAACTTACGCAAGGTGAGTTGGGCGGACGCTTACTTCCTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCt                                                                                                                                         <  1:246954/148‑1 (MQ=255)
            aCGCAAGGTGAGTTGGGCGGACGCTTACTTCCTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGc                                                                                                                                  <  6:35581/147‑1 (MQ=255)
              gCAAGGTGAGTTGGGCGGACGCTTACTTCCTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCgg                                                                                                                              >  2:52355/1‑149 (MQ=255)
                aaGGTGAGTTGGGCGGACGCTTACTTCCTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGATATTGGCTAGTAGGCCCGGct                                                                                                                            >  2:119227/1‑149 (MQ=255)
                aaGGTGAGTTGGGCGGACGCTTACTTCCTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGATATTGGCTAGTAGGCCCGGct                                                                                                                            >  2:115137/1‑149 (MQ=255)
                          gggCGGACGCTTACTTCCTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCgg                                                                                                                   <  2:94670/148‑1 (MQ=255)
                              ggACGCTTACTTCCTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCgg                                                                                                                   <  1:3016/144‑1 (MQ=255)
                              ggACGCTTACTTCCTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCgg                                                                                                                   >  2:3016/1‑144 (MQ=255)
                               caCGCCTACTTCCTATGTTCAAAAACGCAGCTAGCCAGTGTGCGCGAGGCAATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGGGGACATTGGCTAGTAGGCCCGGCTCTTTCTCGGGTAAAc                                                                                                             <  2:224252/148‑1 (MQ=25)
                                  gCTTACTTCCTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCg                                                                                                          >  5:296503/1‑149 (MQ=255)
                                   cTTACTTCCTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGc                                                                                                         <  6:36506/149‑1 (MQ=255)
                                    ttACTTCCTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCt                                                                                                        >  6:78463/1‑149 (MQ=255)
                                    ttACTTCCTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGGAGCCGCGCCAAGGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCt                                                                                                        >  1:297107/1‑149 (MQ=255)
                                             aTGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGATTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATgg                                                                                                 >  1:270317/1‑148 (MQ=255)
                                              tGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGt                                                                                                <  6:296503/147‑1 (MQ=255)
                                              tGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAg                                                                                              <  5:78463/149‑1 (MQ=255)
                                                   aaaaaaGCAGCTAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACtt                                                                                         >  7:272171/1‑149 (MQ=255)
                                                   aaaaaaGCAGCTAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACtt                                                                                         <  5:101735/149‑1 (MQ=255)
                                                     aaaaGCAGCTAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGcc                                                                                                                                 <  4:46317/107‑1 (MQ=255)
                                                     aaaaGCAGCTAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGcc                                                                                                                                 >  3:46317/1‑107 (MQ=255)
                                                          cagcTAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGgga                                                                                  <  6:10161/149‑1 (MQ=255)
                                                           agcTAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGggag                                                                                 <  8:272171/149‑1 (MQ=255)
                                                                                 catacataGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGATATTGGCTNGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTc                                                           <  1:115137/149‑1 (MQ=255)
                                                                                 catacataGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGATATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTc                                                           <  1:119227/149‑1 (MQ=255)
                                                                                 catacataGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCgg                                                                                                                   >  1:307301/1‑93 (MQ=255)
                                                                                 catacataGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCgg                                                                                                                   <  2:307301/93‑1 (MQ=255)
                                                                                  atacataGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGcc                                                                                                                                 >  8:185586/1‑78 (MQ=255)
                                                                                  atacataGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGcc                                                                                                                                 <  7:185586/78‑1 (MQ=255)
                                                                                     cataGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGa                                                        <  2:297107/148‑1 (MQ=255)
                                                                                          gAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTTGCTAGTAGGCCAGGATCTTTAGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGAAGGAGGGCCATAACTCGCTTCTGaatttta                                                  >  3:48/1‑143 (MQ=25)
                                                                                          gAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACt                                                                 >  2:101009/1‑134 (MQ=255)
                                                                                          gAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACt                                                                 <  1:101009/134‑1 (MQ=255)
                                                                                               ggTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCAt                                                                     <  5:191921/125‑1 (MQ=255)
                                                                                               ggTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCAt                                                                     >  6:191921/1‑125 (MQ=255)
                                                                                               ggTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTCGCGCa                                             <  1:52355/149‑1 (MQ=255)
                                                                                                       ccGCGAATACGGTAGCCGCGACAATGGTGAACGGCTGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTCGCGCAGTACAGa                                      >  5:29537/1‑149 (MQ=255)
                                                                                                       ccGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTc                                                           >  1:315450/1‑128 (MQ=255)
                                                                                                       ccGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTc                                                           <  2:315450/128‑1 (MQ=255)
                                                                                                       ccGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTCGCGCAGTACAGa                                      >  8:267849/1‑149 (MQ=255)
                                                                                                       ccGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTCGCGCAGTACAGa                                      >  6:303513/1‑149 (MQ=255)
                                                                                                                      ccGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGa                                                        <  3:269905/116‑1 (MQ=255)
                                                                                                                      ccGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGa                                                        >  4:269905/1‑116 (MQ=255)
                                                                                                                             aaTGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTCGCGCAGTACAGACCGTAGGAGAAATCACCcctcc                >  5:111072/1‑149 (MQ=255)
                                                                                                                                  tGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTCGCGCa                                             <  6:224161/115‑1 (MQ=255)
                                                                                                                                  tGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTCGCGCa                                             >  5:224161/1‑115 (MQ=255)
                                                                                                                                  tGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTCGCGCAGTACAGACCGTAGGAGAAATCACCCCTCCTCCtt           >  3:110993/1‑149 (MQ=255)
                                                                                                                                    aaCGGCAGTTGTACATTTGCTAGTAGGCCCGGCTCTTTAGAGGGTAAACCCGATTGTATGGTAGGACTTGAAGGGAGGATGAGGGCCATAACTCGCTTCTGACTGTTCGCGCAGTACAGACCGTAGGATAAATCACCCCTCCTCCttt          <  4:48/148‑1 (MQ=17)
                                                                                                                                      cGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTCGCGCAGTACAGACCGTAGGAGAAATCACCCCTCCTCCTTTCaa       <  7:267849/149‑1 (MQ=255)
                                                                                                                                            gTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTCGCGCAGTACAGACCGTAGGAGAAATCACCCCTCCTCCTTTCAACAccc  >  3:30674/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CCTGCCTAACTTACGCAAGGTGAGTTGGGCGGACGCTTACTTCCTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTCGCGCAGTACAGACCGTAGGAGAAATCACCCCTCCTCCTTTCAACACCC  >  NC_002947/3971609‑3971895

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: