Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,909,921 C→G 9.1% intergenic (+314/+165) PP_4318 → / ← PP_4319 transposase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,909,9210CG 9.1% 131.5 / 7.6 55intergenic (+314/+165)PP_4318/PP_4319transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (24/25);  new base G (3/2);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.71e-01

CCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCA  >  NC_002947/4909778‑4910064
                                                                                                                                               |                                                                                                                                                
ccGCCCTTACGGCGGGTCACTTTTTGTTAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATgg                                                                                                                                             >  2:176822/1‑149 (MQ=18)
 cGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATggg                                                                                                                                            <  3:180369/149‑1 (MQ=32)
 cGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATgg                                                                                                                                             <  7:155951/148‑1 (MQ=32)
  gCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTGCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATgggg                                                                                                                                           >  1:252254/1‑149 (MQ=17)
  gCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATgggg                                                                                                                                           >  8:200806/1‑149 (MQ=32)
     cTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCt                                                                                                                                        >  5:79958/1‑149 (MQ=32)
        aCGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCa                                                                                                                                     <  5:116718/148‑1 (MQ=32)
         cggcggGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAgc                                                                                                                                   <  5:108966/149‑1 (MQ=32)
         cggcggGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAgc                                                                                                                                   >  3:218835/1‑149 (MQ=32)
          ggcggGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAgcg                                                                                                                                  >  2:294815/1‑149 (MQ=32)
           gcggGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAgcg                                                                                                                                  >  4:49957/1‑148 (MQ=32)
                   cTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACggg                                                                                                                         <  1:169910/149‑1 (MQ=32)
                    tttttGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACggg                                                                                                                         >  5:151014/1‑148 (MQ=32)
                          tCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATcc                                                                                                                    <  2:56525/147‑1 (MQ=32)
                          tCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATcc                                                                                                                    >  1:56525/1‑147 (MQ=32)
                             aaggcgACAAAAAGTAACCAAAAAACGCTGGCTCCTTTCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACgg                                                                                                                          <  8:37602/135‑1 (MQ=17)
                             aaCGCGACAAAAAGTAACAAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTAGATACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACgg                                                                                                                          >  7:37602/1‑138 (MQ=12)
                              aCGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGccc                                                                                                              <  4:218835/149‑1 (MQ=32)
                                 cgACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATGGATGCCCGt                                                                                                            >  4:261634/1‑148 (MQ=17)
                                  gACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATGGATGCCCGTCa                                                                                                          <  3:261634/149‑1 (MQ=17)
                                  gACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCa                                                                                                          <  3:49957/149‑1 (MQ=32)
                                   aCAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTTCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCAc                                                                                                         <  5:156082/149‑1 (MQ=17)
                                   aCAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATGGATGCCCGTCaa                                                                                                         <  1:294815/149‑2 (MQ=17)
                                      aaaaGTAACCAAAACACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCTGGGTTCCCTCGCTACGGTCTTTCTCCCGGTAGTACGCGCTGCAGGCCCCATCCATGGGGCCtgcagc                                                                                                                                   <  7:156598/121‑6 (MQ=11)
                                       aaaGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTcc                                                                                                     <  1:176822/149‑1 (MQ=32)
                                                              ccTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCAAGGGGCCTCAGCGCTTGACGGGCATCCATGCCAGTCACCTCCCTTCGCAAGCCCTGAGCTCGGcc                                                                              >  8:210190/1‑149 (MQ=255)
                                                                  tcatcCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACTCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTc                                                                           <  7:83025/148‑1 (MQ=25)
                                                                    atcCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTg                                                                         <  6:151014/148‑1 (MQ=37)
                                                                     tcCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCCCTAAGCTCATGCTCCCGGGAGTACGCCCTGCAGCCCCCAAACAAGCGGCCTCAGCGCTTGAACGGCATCCATCCCAGTCCCCTAACTTCGAAAGCCCTGCGCACGGCCATCTg                                                                         >  6:180690/1‑147 (MQ=2)
                                                                              tGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATc                                                              >  2:103112/1‑149 (MQ=37)
                                                                                 gcTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCCGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGc                                                           >  5:23243/1‑149 (MQ=37)
                                                                                 gcTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGc                                                           >  5:27715/1‑149 (MQ=37)
                                                                                        tCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGTTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCt                                                    >  2:192581/1‑149 (MQ=37)
                                                                                          cGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCt                                                    <  6:79958/147‑1 (MQ=37)
                                                                                            gggTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAAc                                                <  7:60590/149‑1 (MQ=37)
                                                                                                            gTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACaa                                >  5:113910/1‑149 (MQ=40)
                                                                                                                 gCTCCCGGGAGTACGCGCTGCAGGCCGCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGc                              <  5:180690/146‑1 (MQ=255)
                                                                                                                 gCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCaa                            >  2:203935/1‑148 (MQ=39)
                                                                                                                 gCTCCCGGGAGTACGCGCTGCAGGCCCCAGCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCaa                            <  4:165988/148‑1 (MQ=255)
                                                                                                                  cTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCaaa                           <  7:200806/148‑1 (MQ=40)
                                                                                                                     ccGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCg                                                                                 <  4:236260/91‑1 (MQ=33)
                                                                                                                     ccGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCg                                                                                 >  3:236260/1‑91 (MQ=33)
                                                                                                                     ccGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaag                         >  8:102265/1‑147 (MQ=255)
                                                                                                                     ccGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaag                         <  7:134266/147‑1 (MQ=255)
                                                                                                                           tcaagcgcTGAGGCCCCATGGATGGGGCCTGCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGaa                                                 >  8:8100/5‑117 (MQ=11)
                                                                                                                           tcaagcgcTGAGGCCCCATGGATGGGGCCTGCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGaa                                                 <  7:8100/113‑1 (MQ=11)
                                                                                                                               gcgcTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaag             >  8:132560/1‑149 (MQ=255)
                                                                                                                                cgcTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaagc            <  7:102265/149‑1 (MQ=255)
                                                                                                                                     gaGGCCCCATGGATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGc                                                      <  3:280091/101‑1 (MQ=255)
                                                                                                                                     gaGGCCCCATGGATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGc                                                      >  4:280091/2‑102 (MQ=255)
                                                                                                                                     gaGGCCCCATGGATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaagcagaag       >  6:38709/2‑149 (MQ=255)
                                                                                                                                     cAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCgcgcgc                                      >  6:267964/1‑118 (MQ=39)
                                                                                                                                     cAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaagcaga         <  2:252254/147‑1 (MQ=255)
                                                                                                                                          cccATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCCCCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTGGCAATCTGCGTCGCCTGAACTACCGCGCGCTTCCAAGCAAAAGCAGAAGCAGAagcagaagccgcc  >  8:221417/1‑148 (MQ=255)
                                                                                                                                          cccATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaagcagaagcag    >  8:16877/1‑147 (MQ=255)
                                                                                                                                               |                                                                                                                                                
CCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCA  >  NC_002947/4909778‑4910064

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: