Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,892,547 C→G 7.9% R511G (CGC→GGC)  PP_3435 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,892,5470CG 7.9% 146.5 / 4.5 51R511G (CGC→GGC) PP_3435hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (23/24);  new base G (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 = 2.63e-01

GGGCAAGCTGGACCTGGGCATCGTTGCCGAGGGTGTGGAAACCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATGGCGTCCTGCCTGTTCGCTGTGATGTAC  >  NC_002947/3892399‑3892693
                                                                                                                                                    |                                                                                                                                                  
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                                          <  2:342669/143‑1 (MQ=255)
                                                                                                                aTGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGt                                          >  1:342669/1‑143 (MQ=255)
                                                                                                                aTGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATc                                    >  5:275785/1‑149 (MQ=255)
                                                                                                                   ccggccgCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCgt                                                         >  6:54936/1‑125 (MQ=255)
                                                                                                                   ccggccgCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCgt                                                         <  5:54936/125‑1 (MQ=255)
                                                                                                                       ccgccCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATg                             <  3:15875/149‑1 (MQ=255)
                                                                                                                       ccgccCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATg                             <  5:105966/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      <  2:62989/149‑1 (MQ=255)
                                                                                                                                              gCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATGGCGTCCTGCCTGTTCGCTGTGAt      <  4:192488/149‑1 (MQ=255)
                                                                                                                                              gCCACGCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATGGCGTCCTGCCTGTTCGCtgtg        <  4:36383/147‑1 (MQ=255)
                                                                                                                                                  cTCGCCCTGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATGGCGTCCTGCCTGTTCGCTGTGATGTAc  <  4:15303/149‑1 (MQ=255)
                                                                                                                                                  cTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTCAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATGGCGTCCTGCCTGTTCGCTGTGATGTAc  <  5:320894/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGGCAAGCTGGACCTGGGCATCGTTGCCGAGGGTGTGGAAACCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATGGCGTCCTGCCTGTTCGCTGTGATGTAC  >  NC_002947/3892399‑3892693

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: ‑