Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,892,543 C→G 8.0% A509A (GCC→GCG)  PP_3435 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,892,5430CG 8.0% 150.9 / 4.5 50A509A (GCC→GCG) PP_3435hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (24/22);  new base G (2/2);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.49e-01

GGGCAAGCTGGACCTGGGCATCGTTGCCGAGGGTGTGGAAACCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATGGCGTCCTGCCTGTTCGCTGTGAT  >  NC_002947/3892399‑3892689
                                                                                                                                                |                                                                                                                                                  
gggCAAGCTGGACCTGGGCATCGTTGCCGAGGGTGTGGAAACCGACATGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGCCCCATGCCGGCCGCCCAGTTCCTCGAAGCCATGCCCACTc                                                                                                                                                >  2:239259/1‑149 (MQ=255)
   cAAGCTGGACCTGGGCATCGTTGCCGAGGGTGTGGAAACCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGcc                                                                                                                                             >  7:94245/1‑149 (MQ=255)
          gACCTGGGCATCGTTGCCGAGGGTGTGGAAACCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGgc                                                                                                                                        <  4:79204/147‑1 (MQ=255)
             cTGGGCATCGTTGCCGAGGGTGTGGAAACCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCGGGCGAGTGGCCAGGgcttcga                                                                                                                                   <  6:352953/149‑6 (MQ=255)
                ggCATCGTTGCCGAGGGTGTGGAAACCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCa                                                                                                                                <  6:103548/149‑1 (MQ=255)
                ggCATCGTTGCCGAGGGGGTGGAAACCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCGGGCGAGTGGCCAGGgcttcgag                                                                                                                                  <  5:371687/147‑7 (MQ=255)
                  cATCGTTGCCGAGGGTGTGGAAACCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGt                                                                                                                              <  1:68665/149‑1 (MQ=255)
                         gCCGAGGGTGTGGAAACCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGc                                                                                                                        >  3:282420/1‑148 (MQ=255)
                           cGAGGGTGTGGAAACCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGc                                                                                                                        >  2:47731/1‑146 (MQ=255)
                               ggTGTGGAAACCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATgc                                                                                                                 >  6:320894/1‑149 (MQ=255)
                                     gAAACCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGccccc                                                                                                           >  2:1648/1‑149 (MQ=255)
                                     gAAACCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGATCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGccccc                                                                                                           >  1:305247/1‑149 (MQ=255)
                                      aaaCCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCt                                                                                                          <  1:47731/149‑1 (MQ=255)
                                          cGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCt                                                                                                          >  5:177930/1‑145 (MQ=255)
                                             ccTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCa                                                                                                   <  2:352386/149‑1 (MQ=255)
                                                gCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATgc                                                                                                >  3:192488/1‑149 (MQ=255)
                                                 cAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATg                                                                                                 <  8:94245/147‑1 (MQ=255)
                                                          tACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGAt                                                                                       >  4:15875/1‑148 (MQ=255)
                                                            ccTGGCCGCGCTCGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCTTGCCGGCCGCCCAGTTCCTCGAAGACCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCTCAAGATGCGCCCCCTGAGATCTTGCGCGGCTgtt                                                                                       <  4:260672/146‑3 (MQ=255)
                                                                  cgcgCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGc                                                                               >  3:15303/1‑148 (MQ=255)
                                                                          gggTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTccc                                                                      <  7:288178/149‑1 (MQ=255)
                                                                            gTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTACCg                                                                     >  7:228543/1‑148 (MQ=255)
                                                                                      agcagGGTTACCTGTTCGCCAGACCCATGCCCGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGACCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCa                                                          >  6:58940/1‑149 (MQ=255)
                                                                                                  tGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCaa                                                >  3:282475/1‑147 (MQ=255)
                                                                                                    ttCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAgg                                            >  7:294476/1‑149 (MQ=255)
                                                                                                        ccAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCAtt                                        <  1:239259/149‑1 (MQ=255)
                                                                                                         cAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTg                                       >  2:315539/1‑149 (MQ=255)
                                                                                                            aCCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGtt                                     >  8:78252/1‑148 (MQ=255)
                                                                                                              ccATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCa                                  >  4:338916/1‑149 (MQ=255)
                                                                                                               cATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCAt                                 >  1:69244/1‑149 (MQ=255)
                                                                                                                aTGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGt                                      >  1:342669/1‑143 (MQ=255)
                                                                                                                aTGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGt                                      <  2:342669/143‑1 (MQ=255)
                                                                                                                aTGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATc                                >  5:275785/1‑149 (MQ=255)
                                                                                                                   ccggccgCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCgt                                                     <  5:54936/125‑1 (MQ=255)
                                                                                                                   ccggccgCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCgt                                                     >  6:54936/1‑125 (MQ=255)
                                                                                                                       ccgccCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATg                         <  5:105966/149‑1 (MQ=255)
                                                                                                                       ccgccCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATg                         <  3:15875/149‑1 (MQ=255)
                                                                                                                            cAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATGGCGTc                    <  1:1648/149‑1 (MQ=255)
                                                                                                                              gTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATGGCGTcct                  <  6:177930/149‑1 (MQ=255)
                                                                                                                               ttCCTCGCAGCCCTGGCCACTCGCCCGGGCGCAACGCAGGTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAACGGCATTGTTTCATCGGGCATGGCGTcctg                 <  2:305247/149‑1 (MQ=255)
                                                                                                                               ttCCTCGAAGCCCTGGCCACGCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATGGCGTcctg                 >  3:36383/1‑149 (MQ=255)
                                                                                                                                 gcgttgcgCCCGGGCGAGTGGCCAGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCACGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATGGCGTcctgcc               >  6:203564/8‑149 (MQ=255)
                                                                                                                                  cgttgcgCCCGGGCGAGTGGCCAGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATGGCGTcctgcct              >  4:350476/7‑149 (MQ=255)
                                                                                                                                  cTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATGGCGTcctgcct              >  2:36762/1‑149 (MQ=255)
                                                                                                                                         cccTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATGGCGTCCTGCCTGTTCGCt       <  4:158279/149‑1 (MQ=255)
                                                                                                                                          ccTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATGGCGTCCTGCCTGTTCGCtg      <  4:282420/149‑1 (MQ=255)
                                                                                                                                              gCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATGGCGTCCTGCCTGTTCGCTGTGa   <  1:36762/148‑1 (MQ=255)
                                                                                                                                              gCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATGGCGTCCTGCCTGTTCGCTGTGAt  <  4:192488/149‑1 (MQ=255)
                                                                                                                                              gCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATGGCGTCCTGCCTGTTCGCTGTGAt  <  2:62989/149‑1 (MQ=255)
                                                                                                                                              gCCACGCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATGGCGTCCTGCCTGTTCGCtgtg    <  4:36383/147‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GGGCAAGCTGGACCTGGGCATCGTTGCCGAGGGTGTGGAAACCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATGGCGTCCTGCCTGTTCGCTGTGAT  >  NC_002947/3892399‑3892689

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