Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,368,475 A→C 10.5% R580R (AGG→CGG)  PP_2964 → transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,368,4750AC10.5% 99.6 / 3.6 38R580R (AGG→CGG) PP_2964transposase
Reads supporting (aligned to +/- strand):  ref base A (16/18);  new base C (2/2);  total (18/20)
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

TCCCATCATGCCGAGGGGACGTATCGATTAACGAAGCTGGAGAACCGGTGTTGCGCCGAGTAGCGGCGCGGGACATCCCCCCTTCGGCCATCTCGCTGCAGACGGCGCTTATGCAGCGCATGCCGGCCAGGCATGTGCTCGATATCATGGCTAACATCGAGCATTGGATTCAGTTCACGCGGCATTTCGGGCCCATGTCCGGTAACGAGCCGAAGCTCAAAGAGCCTGCCGAGCGCTACCTGATGACGATCTTCGCCATGGGC  >  NC_002947/3368347‑3368609
                                                                                                                                |                                                                                                                                      
tCCCATCATGCCGAGGGGACGTATCGATTAACGAAGCTGGAGAACCGGTGTTGCGCCGAGTAGCGGCGCGGGACATCCCCCCTTCGGCCATCTCGCTGCAGACGGCGCTTATGCAGCGCATGCCGGCCAGGCATGTGCTCGATATCATgg                                                                                                                   <  1:1200445/150‑1 (MQ=255)
  ccATCATGCCGAGGGGACGTATCGATTAACGAAGCTGGAGAACCGGTGTTGCGCCGAGTAGCGGCGCGGGACATCCCCCCTTCGGCCATCTCGCTGCAGACGGCGCTTATGCAGCACATGCCTGGCCGGCATGCGCtgcata                                                                                                                         >  2:793559/1‑137 (MQ=255)
  ccATCATGCCGAGGGGACGTATCGATTAACGAAGCTGGAGAACCGGTGTTGCGCCGAGTAGCGGCGCGGGACATCCCCCCTTCGGCCATCTCGCTGCAGACGGCGCTTATGCAGCACATGCCTGGCCGGCATGCGCtgcata                                                                                                                         <  1:793559/142‑6 (MQ=255)
     tcatGCCGAGGGGACGTATCGATTAACGAAGCTGGAGAACCGGTGTTGCGCCGAGTAGCGGCGCGGGACATCCCCCCTTCGGCCATCTCGCTGCAGACGGCGCTTATGCAGCGCATGCCTGGCCGGCAtgcgc                                                                                                                               >  2:217433/1‑130 (MQ=255)
     tcatGCCGAGGGGACGTATCGATTAACGAAGCTGGAGAACCGGTGTTGCGCCGAGTAGCGGCGCGGGACATCCCCCCTTCGGCCATCTCGCTGCAGACGGCGCTTATGCAGCGCATGCCTGGCCGGCATgcgc                                                                                                                               <  1:217433/133‑1 (MQ=255)
     tcatGCCGAGGGGACGTATCGATTAACGAAGCTGGAGAACCGGTGTTGCGCCGAGTAGCGGCGCGGGACATCCCCCCTTCGGCCATCTCGCTGCAGACGGCGCTTATGCAGCGCATGCCGGCCAGGCATGTGCTAGATATCATGGCTAAc                                                                                                              <  1:240896/150‑1 (MQ=255)
            gAGGGGACGTATCGATTAACGAAGCTGGAGAACCGGTGTTGCGCCGAGTAGCGGCGCGGGACATCCCCCCTTCGGCCATCTCGCTGCAGACGGCGCTTATGCAGCGCATGCCGGCCAGGCATGTGCTCGATATCATGGCTAACATCGAGc                                                                                                       >  1:112278/1‑150 (MQ=255)
                    gTATCGATTAACGAAGCTGGAGAACCGGTGTTGCGCCGAGTAGCGGCGCGGGACATCCCCCCTTCGGCCATCTCGCTGCAGACGGCGCTTATGCAGCGCATGCCGGCCAGGCATGTGCTCGATATCATGGCTAACATCGAGCATTGGAtt                                                                                               <  2:512806/150‑1 (MQ=255)
                    gTATCGATTAACGAAGCTGGAGAACCGGTGTTGCGCCGAGTAGCGGCGCGGGACATCCCCCCTTCGGCCATCTCGCTGCAGACGGCGCGTATGCAGCGCATGCCGGCCAGGCATGTGCTCGATATCATGGCTAACATCGAGCATTGGAtt                                                                                               <  1:531242/150‑1 (MQ=255)
                      aTCGATTAACGAAGCTGGAGAACCGGTGTTGCGCCGAGTAGCGGCGCGGGACATCCCCCCTTCGGCCATCTCGCTGCAGACGGCGCTTATGCAGCGCATGCCGGCCAGGCATGTGCTCGATATCATGGCTAACATCGAGCATTGGATTCa                                                                                             >  2:1001489/1‑150 (MQ=255)
                         gATTAACGAAGCTGGAGAACCGGTGTTGCGCCGAGTAGCGGCGCGGGACATCCCCCCTTCGGCCATCTCGCTGCAGACGGCGCTTATGCAGCGCATGCCGGCCAGGCATGTGCTCGa                                                                                                                           <  2:47222/117‑1 (MQ=255)
                         gATTAACGAAGCTGGAGAACCGGTGTTGCGCCGAGTAGCGGCGCGGGACATCCCCCCTTCGGCCATCTCGCTGCAGACGGCGCTTATGCAGCGCATGCCGGCCAGGCATGTGCTCGa                                                                                                                           >  1:47222/1‑117 (MQ=255)
                          aTTAACGAAGCTGGAGAACCGGTGTTGCGCCGAGTAGCGGCGCGGGACATCCCCCCTTCGGCCATCTCGCTGCAGACGGCGCTTATGCAGCGCATGCCGGCCAGGCAt                                                                                                                                   <  1:859426/108‑1 (MQ=255)
                          aTTAACGAAGCTGGAGAACCGGTGTTGCGCCGAGTAGCGGCGCGGGACATCCCCCCTTCGGCCATCTCGCTGCAGACGGCGCTTATGCAGCGCATGCCGGCCAGGCAt                                                                                                                                   >  2:859426/1‑108 (MQ=255)
                               cGAAGCTGGAGAACCGGTGTTGCGCCGAGTAGCGGCGCGGGACATCCCCCCTTCGGCCATCTCGCTGCAGACGGCGCTTATGCAGCGCATGCCGGCCAGGCATGTGCTCGATATCAt                                                                                                                     <  2:593947/117‑1 (MQ=255)
                               cGAAGCTGGAGAACCGGTGTTGCGCCGAGTAGCGGCGCGGGACATCCCCCCTTCGGCCATCTCGCTGCAGACGGCGCTTATGCAGCGCATGCCGGCCAGGCATGTGCTCGATATCAt                                                                                                                     >  1:593947/1‑117 (MQ=255)
                                         gaACCGGTGTTGCGCCGAGTAGCGGCGCGGGACATCCCCCCTTCGGCCATCTCGCTGCAGACGGCGCTTATGCAGCGCATGCCGGCCAGGCATGTGCTCGATATCATGGc                                                                                                                  >  1:1000777/1‑110 (MQ=255)
                                         gaACCGGTGTTGCGCCGAGTAGCGGCGCGGGACATCCCCCCTTCGGCCATCTCGCTGCAGACGGCGCTTATGCAGCGCATGCCGGCCAGGCATGTGCTCGATATCATGGc                                                                                                                  <  2:1000777/110‑1 (MQ=255)
                                                    gcgcCGAGTAGCGGCGCGGGACATCCCCCCTTCGGCCATCTCGCTGCAGACGGCGCTTATGCAGCGCATGCCGGCCAGGCATGTGCTCGATATCATGGCTAACATCGAGCATTGGATTCAGTTCACGCGGCATTTCGGGCCCATGTCCgg                                                               <  2:311271/150‑1 (MQ=255)
                                                              gcggcgCGGGACATCCCCCCTTCGGCCATCTCGCTGCAGACGGCGCTTATGCAGCGCATGCCGGCCAGGCATGTGCTCGATATCATGGCTAACATCGAGCATTGGATTCAGTTCACGCGGCATTTCGGGCCCATGTCCg                                                                >  1:504147/1‑139 (MQ=255)
                                                              gcggcgCGGGACATCCCCCCTTCGGCCATCTCGCTGCAGACGGCGCTTATGCAGCGCATGCCGGCCAGGCATGTGCTCGATATCATGGCTAACATCGAGCATTGGATTCAGTTCACGCGGCATTTCGGGCCCATGTCCg                                                                <  2:504147/139‑1 (MQ=255)
                                                                ggcgCGGGACATCCCCCCTTCGGCCATCTCGCTGCAGACGGCGCTTATGCAGCGCATGCCGGCCAGGCATGTGCGTATTATCAGGGGTAACATCGAGCATTGGATTCAGTTCACGCGGCATTTCGGGCCCATGTCCGGTAACGCGCCGaa                                                   >  1:749097/1‑150 (MQ=255)
                                                                            ccccccTTCGGCCATCTCGCTGCAGACGGCGCTTATGCAGCGCATGCCGGCCAGGCATGTGCTCGATATCATGGCTAACATCGAg                                                                                                        >  1:1100053/1‑85 (MQ=255)
                                                                            ccccccTTCGGCCATCTCGCTGCAGACGGCGCTTATGCAGCGCATGCCGGCCAGGCATGTGCTCGATATCATGGCTAACATCGAg                                                                                                        <  2:1100053/85‑1 (MQ=255)
                                                                              ccccTTCGGCCATCTCGCTGCAGACGGCGCTTATGCAGCGCATGCCGGCCAGGCATGTGCTCGATATCATGGCTAACATCGAGCATTGGAtt                                                                                               <  2:239679/92‑1 (MQ=255)
                                                                              ccccTTCGGCCATCTCGCTGCAGACGGCGCTTATGCAGCGCATGCCGGCCAGGCATGTGCTCGATATCATGGCTAACATCGAGCATTGGAtt                                                                                               >  1:239679/1‑92 (MQ=255)
                                                                                ccTTCGGCCATCTCGCTGCAGACGGCGCTTATGCAGCGCATGCCGGCCAGGCATGTGCTCGATATCATGGCTAACATCGAGCATTGGATTCAGTTCACGCGGCATTTCGGGCCAATTTCCGGTAACGAGCCGAAGCTCAAAGAGCCTGcc                                   <  1:1001489/150‑1 (MQ=255)
                                                                                 cTTCGGCCATCTAGCTGCAGACGGCGCTTATGCAACGCATGCCGGCCAGGCATGTGCTCGATATCATGGCTAACATCGAGCATTGGATTCAGTTCACGCGGCATTTCGGGCCCATGTCCGGTAACGAGCCGAAGCTCAAAGAGCCTGCCg                                  >  1:833907/1‑150 (MQ=255)
                                                                                   tCGGCCATCTCGCTGCAGACGGCGCTTATGCAGCGCATGCCGGCCAGGCATGTGCTCGATATCATGGCTAACATCGAGCATTGGATTCAGTTCACGCGGCATTTCGGGCCCATGTCCGGTAACGAGCCGAAGCTCAAAGAGCCTGCCGAg                                <  2:112278/150‑1 (MQ=255)
                                                                                        cATCTCGCTGCAGACGGCGCTTATGCAGCGCATGCCGGCCAGGCATGTGCTCGATATCATGGCTAACATCGAGCATTGGATTCAGTTCACGCGGCAttt                                                                              <  2:600135/99‑1 (MQ=255)
                                                                                        cATCTCGCTGCAGACGGCGCTTATGCAGCGCATGCCGGCCAGGCATGTGCTCGATATCATGGCTAACATCGAGCATTGGATTCAGTTCACGCGGCAttt                                                                              >  1:600135/1‑99 (MQ=255)
                                                                                         aTCTCGCTGCAGACGGCGCTTATGCAGCGCATGCCGGCCAGGCATGTGCTCGATATCATGGCTAACATCGAGCATTGGATTCAGTTCACGCGGCATTTCGGGCCCATGTCCGGTAACGAGCCGAAGCTCAAAGAGCCTGCCGAGCGCTAc                          >  2:1094532/1‑150 (MQ=255)
                                                                                                     aCGGCGCTTATGCAGCGCATGCCGGCCAGGCATGTGCTCGATATCATGGCTAACATCGAGCATTGGATTCAGTTCACGCGGCATTTCGGGCCCATGTCCGGTAACGAGCCGAAGCTCAAAGAGCCTGCCGAGCGCTACCTGATGACGATc              <  1:1094532/150‑1 (MQ=255)
                                                                                                       ggCGCTTATGCAGCGCATGCCGGCCAGGCATGTGCTCGATATCATGGCTAACATCGAGCATTGGATTCAGTTCACGCGGCATTTCGGGCCCATGTCCGGTAACGAGCCGAAGCTCAAAGAGCCTGCCGAGCGCTACCTGATGACGATCtt            >  2:788450/1‑150 (MQ=255)
                                                                                                        gcgcTTATGCAGCGCATGCCGGCCAGGCATGTGCTCGATATCATGGCTAACATCGAGCATTGGATTCAGTTCACGCGGCATTTCGGGCCCATGTCCgg                                                               <  2:135506/98‑1 (MQ=255)
                                                                                                        gcgcTTATGCAGCGCATGCCGGCCAGGCATGTGCTCGATATCATGGCTAACATCGAGCATTGGATTCAGTTCACGCGGCATTTCGGGCCCATGTCCgg                                                               >  1:135506/1‑98 (MQ=255)
                                                                                                             tATGCAGCGCATGCCGGCCAGGCATGTGCTCGATATCATGGCTAACATCGAGCATTGGATTCAGTTCACGCGGCATTTCGGGCCCATGTCCGGTAACGAGCCGAAGCTCAAAGAGCCTGCCGAGCGCTACCTGATGACGATCTTCGCCAt      <  1:788450/150‑1 (MQ=255)
                                                                                                                 gaGCGCATGCCGGCCAGGCATGTGCTCGATATCATGGCTAACATCGAGCATTGGATTCAGTTCACGCGGCATTTCGGGCCCATGTCCGGTAACGAGCCGAAGCTCAAAGAGCCTGCCGAGCGCTACCTGATGACGATCTTCGCCATGGGc  >  2:172381/2‑150 (MQ=255)
                                                                                                                                |                                                                                                                                      
TCCCATCATGCCGAGGGGACGTATCGATTAACGAAGCTGGAGAACCGGTGTTGCGCCGAGTAGCGGCGCGGGACATCCCCCCTTCGGCCATCTCGCTGCAGACGGCGCTTATGCAGCGCATGCCGGCCAGGCATGTGCTCGATATCATGGCTAACATCGAGCATTGGATTCAGTTCACGCGGCATTTCGGGCCCATGTCCGGTAACGAGCCGAAGCTCAAAGAGCCTGCCGAGCGCTACCTGATGACGATCTTCGCCATGGGC  >  NC_002947/3368347‑3368609

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