Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,411,659 T→A 8.1% V58D (GTC→GAC)  PP_5551 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,411,6590TA 8.1% 138.9 / 2.7 49V58D (GTC→GAC) PP_5551hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (24/21);  new base A (2/2);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

TGGATTGCACTTGTCCTGTCTGTTGTAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGACCTTTCGGCAAGATGT  >  NC_002947/3411511‑3411807
                                                                                                                                                     |                                                                                                                                                    
tGGATTGCACTTGTCCTGTCTGTTGTAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGt                                                                                                                                                      >  5:77146/1‑149 (MQ=255)
 ggATTGCACTTGTCCTGTCTGTTGTAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGt                                                                                                                                                      <  1:138165/148‑1 (MQ=255)
  gATTGCTCTTGTCCTGTCTGTTGTAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTc                                                                                                                                                     <  5:38042/148‑1 (MQ=255)
   aTTGCACTTGTCCTGTCTGTTGTAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCaa                                                                                                                                                   >  7:202764/1‑149 (MQ=255)
         cTTGTCCTGTCTGTTGTAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCaa                                                                                                                                             <  7:334453/149‑1 (MQ=255)
                 gtctgtTGTAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGATTGACCTTCGCAAGCtc                                                                                                                                         <  5:204159/145‑2 (MQ=255)
                 gtctgtTGTAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGATTGACCTTCGCAAGCtc                                                                                                                                         >  6:204159/1‑144 (MQ=255)
                     gttgtAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGATTGACCTTCGCAAGCt                                                                                                                                          >  4:226462/1‑140 (MQ=255)
                     gttgtAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGATTGACCTTCGCAAGCt                                                                                                                                          <  3:226462/140‑1 (MQ=255)
                      ttgtAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGa                                                                                                                                 >  5:186481/1‑148 (MQ=255)
                        gtAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGaaa                                                                                                                               <  6:77146/148‑1 (MQ=255)
                         tAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCtt                                                                                                                                    >  7:157696/1‑142 (MQ=255)
                         tAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCtt                                                                                                                                    <  8:157696/142‑1 (MQ=255)
                           aCGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGa                                                                                                                            <  4:313077/148‑1 (MQ=255)
                            cGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGAtt                                                                                                                          <  7:308259/149‑1 (MQ=255)
                             gCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTa                                                                                                                         <  3:331609/149‑1 (MQ=255)
                                 gTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTc                                                                                                                      >  1:34737/1‑148 (MQ=255)
                                     cTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACt                                                                                                                       >  3:305924/1‑143 (MQ=255)
                                     cTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACt                                                                                                                       <  4:305924/143‑1 (MQ=255)
                                                         cGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGt                                                                                             <  5:226146/149‑1 (MQ=255)
                                                            tCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCAc                                                                                          >  5:28316/1‑149 (MQ=255)
                                                                cGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCActct                                                                                       <  6:186481/148‑1 (MQ=255)
                                                                     gATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGAGTCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGa                                                                                                                            >  1:248280/1‑107 (MQ=255)
                                                                     gATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGa                                                                                                                            <  2:248280/106‑1 (MQ=255)
                                                                           gCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCATAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGAc                                                                           >  8:79828/1‑149 (MQ=255)
                                                                           gCTAATTTAGCAAGTCAGGCCGCGAGGGATCATGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGAc                                                                           >  6:324957/1‑149 (MQ=255)
                                                                             tAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGa                                                                         >  2:309250/1‑149 (MQ=255)
                                                                                tttAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGt                                                                       >  4:659/1‑148 (MQ=255)
                                                                                        gTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGt                                                              <  6:28316/149‑1 (MQ=255)
                                                                                         tCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTg                                                             <  3:659/149‑1 (MQ=255)
                                                                                            ggCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAgg                                                          <  2:34737/149‑1 (MQ=255)
                                                                                                       aTCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTTCTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTcc                                               <  5:359228/149‑1 (MQ=255)
                                                                                                             gAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCtt                                          >  4:227788/1‑149 (MQ=255)
                                                                                                             gAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCtt                                          >  1:190929/1‑149 (MQ=255)
                                                                                                                gATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGAt                                       <  7:29566/149‑1 (MQ=255)
                                                                                                                     aaaGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAg                                                           >  2:18283/1‑124 (MQ=255)
                                                                                                                     aaaGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAg                                                           <  1:18283/124‑1 (MQ=255)
                                                                                                                     aaaGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTa                                   >  2:324644/1‑148 (MQ=255)
                                                                                                                      aaGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGa                                        >  6:271591/1‑142 (MQ=255)
                                                                                                                      aaGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGa                                        <  5:271591/142‑1 (MQ=255)
                                                                                                                                   aagaagAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAg                    >  7:124022/1‑149 (MQ=255)
                                                                                                                                      aagaGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGAc                 >  1:151152/1‑149 (MQ=255)
                                                                                                                                       agagCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGAc                 <  1:324644/148‑1 (MQ=255)
                                                                                                                                        gagCTTGCGAAGGTCAGTCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGACCt               >  8:338375/1‑149 (MQ=255)
                                                                                                                                         agCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGACCtt              >  4:367157/1‑149 (MQ=255)
                                                                                                                                          gCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGACCttt             >  6:219450/1‑149 (MQ=255)
                                                                                                                                           cTTGCGAAGGTCAGTCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGACCTTTc            <  7:338375/149‑1 (MQ=255)
                                                                                                                                            ttGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGACCTTTCg           >  2:340183/1‑149 (MQ=255)
                                                                                                                                              gCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGACCTTTCgg          >  6:278826/1‑148 (MQ=255)
                                                                                                                                                gAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCGGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGAACCTTTCGGc         >  1:319208/1‑147 (MQ=255)
                                                                                                                                                   ggTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGACCTTTCGGCAAGAt    <  8:202764/149‑1 (MQ=255)
                                                                                                                                                     tCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGACCTTTCGGCAAGATGt  <  1:340183/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGGATTGCACTTGTCCTGTCTGTTGTAACGCTGGTCGCTCAGAATAAGCAAAGCAAGCGATCTGCGGCTGATCAGGCTAATTTAGCAAGTCAGGCCGCGAGGGATCGTGAAGATCTAAAGAATTCTGCAAAAGAAGAGCTTGCGAAGGTCAATCGCAAGCTTCGCTTGGAAAAGATTACTCAAGATTTAGACGAGCTTTTGGCTGTCACTCTTCAGTATTGGACGAGTGACGGCCGTGAGGGCGCTACTTCCGCCTTGATGCTAAAAGTCAAGGTAAAGGACCTTTCGGCAAGATGT  >  NC_002947/3411511‑3411807

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: