Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,414,745 Δ1 bp 10.5% intergenic (+287/‑58) PP_2114 → / → PP_2115 transposase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,414,7450T.10.5% 101.6 / 8.0 38intergenic (+287/‑58)PP_2114/PP_2115transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (15/19);  new base . (2/2);  total (17/21)
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.00e-01

TCCATTCGCGATGGCGACACTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAG  >  NC_002947/2414597‑2414892
                                                                                                                                                    |                                                                                                                                                   
tCCATTCGCGATGGCGACACTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCt                                                                                                                                                     >  8:94246/1‑149 (MQ=39)
tCCATTCGCGATGGCGACACTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCt                                                                                                                                                     >  6:122580/1‑149 (MQ=39)
tCCATTCGCGATGGCGACACTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCt                                                                                                                                                     >  5:163186/1‑149 (MQ=39)
   aTTCGCGATGGCGACACTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTccc                                                                                                                                                  >  7:73031/1‑149 (MQ=38)
         gATGGCGACACTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAg                                                                                                                                             <  1:140245/148‑1 (MQ=34)
         gATGGCGACACTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAg                                                                                                                                             >  3:215458/1‑148 (MQ=34)
         gATGGCGACACTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCAGCGCTGCGCACCGGGGTCCCCTCACAACCGCCTTGCTCCCGGCAg                                                                                                                                             >  2:140245/1‑148 (MQ=11)
             gCGACACTGACTCTCCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACcg                                                                                                                                        >  6:167268/1‑149 (MQ=21)
                      aCTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAgg                                                                                                                               >  8:181662/1‑149 (MQ=32)
                      aCTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAg                                                                                                                                <  1:300889/148‑1 (MQ=32)
                       cTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGTGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGc                                                                                                                              <  1:13074/149‑1 (MQ=17)
                          aCCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGcccc                                                                                                                           <  8:271023/149‑1 (MQ=32)
                          aCCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGcc                                                                                                                             <  7:44548/147‑1 (MQ=32)
                            cTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCa                                                                                                                                              >  6:119298/1‑128 (MQ=32)
                            cTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCa                                                                                                                                              <  5:119298/128‑1 (MQ=32)
                                tCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCt                                                                                                                      <  8:162605/148‑1 (MQ=31)
                                       ctattgctgtagctattgctgttgctgttgctgttgatcttgatcttgatcttgatcttgatcttgatcttgatcttcgcgacttcaggaggccgaacgcaggccttgcTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCt                                                                                                               >  7:303448‑M2/110‑148 (MQ=255)
                                                                     ctgttgatcttgatcttgatcttgatcttgatcttgatcttgatcttcgcgacttcaggaggccgaacgcaggccttgcTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCAcc                                                                                <  8:303448‑M2/70‑1 (MQ=255)
                                                                                  tcttgatcttgatcttgatcttgatcttgatcttcgcgacttcaggaggccgaacgcaggccttgcGAAGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCAc                                                                                 <  8:211186‑M2/68‑1 (MQ=255)
                                                                                  tcttgatcttgatcttgatcttgatcttgatcttcgcgacttcaggaggccgaacgcaggccttgcGAAGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCAc                                                                                 >  1:211419‑M2/67‑134 (MQ=255)
                                                                                  tcttgatcttgatcttgatcttgatcttgatcttcgcgacttcaggaggccgaacgcaggccttgcGAAGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCAc                                                                                 <  2:211419‑M2/68‑1 (MQ=255)
                                                                                  tcttgatcttgatcttgatcttgatcttgatcttcgcgacttcaggaggccgaacgcaggccttgcGAAGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCAc                                                                                 >  7:211186‑M2/67‑134 (MQ=255)
                                                                                   aaaCGCTGCGCTTCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGc                                                                  <  7:176980/149‑1 (MQ=17)
                                                                                                 catcatCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGcctc                                                    <  5:65522/149‑1 (MQ=32)
                                                                                                           ccccTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTcgc                                          <  8:73031/149‑1 (MQ=32)
                                                                                                            cccTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCGCAGCGCTTGACCGGCCTCCAGGCCCGTCCCCTGCCGTCGCCAGGCCGGCGTTCGGCCGCGGGAAGTcgcg                                         >  7:240339/1‑149 (MQ=11)
                                                                                                                        ctCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatc                             <  5:167268/149‑1 (MQ=32)
                                                                                                                           cGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaag                          <  5:163660/149‑1 (MQ=32)
                                                                                                                           cGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaa                           <  8:271203/148‑1 (MQ=32)
                                                                                                                             gggTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTg                                                 <  6:312479/124‑1 (MQ=32)
                                                                                                                             gggTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTg                                                 >  5:312479/1‑124 (MQ=32)
                                                                                                                                tCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGTGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaa                                 <  5:101962/137‑1 (MQ=17)
                                                                                                                                  cccTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatc                       >  5:137733/1‑145 (MQ=32)
                                                                                                                                  cccTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatc                       >  5:6410/1‑145 (MQ=32)
                                                                                                                                    cTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatc                 <  3:224241/149‑1 (MQ=32)
                                                                                                                                           cGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGcctcc                                                   >  7:21132/1‑108 (MQ=32)
                                                                                                                                           cGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGcctcc                                                   <  8:21132/108‑1 (MQ=32)
                                                                                                                                            cgttctgcTCCCGGCAGGACCGCGCCTAAGGCCCCATCCTGGGGCCTAAGCGCTTGACGGGCATCCATGCCCGTCACTTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTGGCTatgatcaagatcaagatcaagatcaagatcaa         <  1:19671‑M2/141‑1 (MQ=255)
                                                                                                                                             gccttgcTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCGGAAGTCGCGaagatcaagatcaagatcaagatcaagatc           >  1:131481‑M2/8‑146 (MQ=255)
                                                                                                                                                 tgcTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCAGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaagatcaagatc     >  6:94318‑M2/4‑148 (MQ=255)
                                                                                                                                                   cTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAACtcgcgaagatcaagatcaagatcaagatcaagatcaagagaacg  >  2:19671‑M2/2‑144 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCCATTCGCGATGGCGACACTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAG  >  NC_002947/2414597‑2414892

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: