Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,676,428 T→C 8.0% K140E (AAG→GAG)  PP_3241 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,676,4280TC 8.0% 139.3 / 2.7 50K140E (AAG→GAG) PP_3241hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (23/23);  new base C (2/2);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01

GGGGTGGGCATGCAGCAGGCGGGTAAGTACTTCGGTCTTCAATGCATCGTTACTCATGTGCGGCCTCACTTGCCGTGCTTGGACTGAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCACCGCTTTTTCGTGAGCGTCGACCTGGTTGTTGAGGTACGACTTGTCGAAGGACTCATCGCGCCACTCCAGGATCACTTTCTTGATCTTGTCGG  >  NC_002947/3676280‑3676568
                                                                                                                                                    |                                                                                                                                            
ggggTGGGCATGCAGCAGGCGGGTAAGTACTTCGGTCTTCAATGCATCGTTACTCATGTGCGGCCTCACTTGCCGTGCTTGGACTGAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCtt                                                                                                                                              <  2:125369/149‑1 (MQ=255)
ggggTGGGCATGCAGCAGGCGGGTAAGTACTTCGGTCTTCAATGCATCGTTACTCATGTGCGGCCTCACTTGCCGTGCTTGGACTGAAGGGCTTTGGCCTGTTCCAGGTGCTGGTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCtt                                                                                                                                              <  2:15986/149‑1 (MQ=255)
ggggTGGGCATGCAGCAGGCGGGTAAGTACTTCGGTCTTCAATGCATCGTTACTCATGTGCGGCCTCACTTGCCGTGCTTGGACTGAAGGGCTTTGGCCTGTTCCAGGTGCTGATCGAGTTTGGGCAGCGTTTCGATGGCAAAGGCCtt                                                                                                                                              >  7:202155/1‑149 (MQ=255)
  ggTGGGCATGCTGCTGGGGGGTATGTACTTGGGTCTTCTATGCATCGTTACTCTTGTGTGGCCTCACTTGCCGTGCTTGGACTGAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCa                                                                                                                                            <  8:202155/149‑1 (MQ=255)
       gCATGCAGCAGGCGGGTAAGTACTTCGGTCTTCAATGCATCGTTACTCATGTGCGGCCTCACTTGCCGTGCTTGGACTGAAGGGCTTTGGCCTGTTCCAGGTGCTGTTAGAGTTTGGGCAGCGTTTCGCTGGCAAAGTCCTTCCGCTCg                                                                                                                                       >  8:227654/1‑149 (MQ=255)
          tGCAGCAGGCGGGTAAGTACTTCGGTCTTCAATGCATCGTTACTCATGTGCGGCCTCACTTGCCGTGCTTGGACTGAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGcc                                                                                                                                    >  5:75259/1‑149 (MQ=255)
                agGCGGGTAAGTACTTCGGTCTTCAATGCATCGTTACTCATGTGCGGCCTCACTTGCCGTGCTTGGACTGAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCg                                                                                                                              <  8:68049/149‑1 (MQ=255)
                      gTAAGTACTTCGGTCTTCAATGCATCGTTACTCATGTGCGGCCTCACTTGCCGTGCTTGGACTGAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGAt                                                                                                                        >  3:157199/1‑149 (MQ=255)
                         aGTACTTCGGTCTTCAATGCATCGTTACTCATGTGCGGCCTCACTTGCCGTGCTTGGACTGAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGc                                                                                                                      <  5:43646/148‑1 (MQ=255)
                                cGGTCTTCAATGCATCGTTACTCATGTGCGGCCTCACTTGCCGTGCTTGGACTGAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGcc                                                                                                                                    >  5:151750/1‑127 (MQ=255)
                                cGGTCTTCAATGCATCGTTACTCATGTGCGGCCTCACTTGCCGTGCTTGGACTGAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGcc                                                                                                                                    <  6:151750/127‑1 (MQ=255)
                                   tCTTCAATGCATCGTTACTCATGTGCGGCCTCACTTGCCGTGCTTGGACTGAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCtctttct                                                                                                           <  5:251257/149‑1 (MQ=255)
                                        aaTGCATCGTTACTCATGTGCGGCCTCACTTGCCGTGCTTGGACTGAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGaa                                                                                                       >  4:129149/1‑148 (MQ=255)
                                        aaTGCATCGTTACTCATGTGCGGCCTCACTTGCCGTGCTTGGACTGAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGaa                                                                                                       <  5:281236/148‑1 (MQ=255)
                                        aaTGCATCGTTACTCATGTGCGGCCTCACTTGCCGTGCTTGGACTGAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGaa                                                                                                       <  7:227654/148‑1 (MQ=255)
                                         aTGCATCGTTACTCATGTGCGGCCTCACTTGCCGTGCTTGGACTGAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACa                                                                                                     <  4:257016/149‑1 (MQ=255)
                                                 ttACTCATGTGCGGCCTCACTTGCCGTGCTTGGACTGAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGTCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCAcc                                                                                             >  1:32159/1‑149 (MQ=255)
                                                     tCATGTGCGGCCTCACTTGCCGTGCTTGGACTGAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCACCGCtt                                                                                         >  5:174816/1‑149 (MQ=255)
                                                       aTGTGCGGCCTCACTTGCCGTGCTTGGACTGAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCACCGCtttt                                                                                       >  3:9315/1‑149 (MQ=255)
                                                         gtgCGGCCTCACTTGCCGTGCTTGGACTGAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCACCGCTTTTTc                                                                                     >  4:11800/1‑149 (MQ=255)
                                                                              ttGGACTGAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCACCGCTTTTTCGTGAGCGTCGACCTGgttgtt                                                                <  5:8846/149‑1 (MQ=255)
                                                                              ttGGACTGAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGACTCTTTCTTGAACAGCTCCACAGCTTTTTCGTGAGCGTCGACCTGgttgtt                                                                >  6:167400/1‑149 (MQ=255)
                                                                                ggACTGAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCACCGCTTTTTCGTGAGCGTCGACCTGgttgtt                                                                <  6:75259/147‑1 (MQ=255)
                                                                                 gACTGAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCACCGCTTTTTCGTGAGCGTCGACCTGGTTGTTGa                                                              <  6:174816/148‑1 (MQ=255)
                                                                                 gACTGAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCACCGCTTTTTCGTGAGCGTCGACCTGGTTGTTGa                                                              <  4:137506/148‑1 (MQ=255)
                                                                                  aaggAAGGTCTTTGGCCTGTTCAAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCATCCTCTTTCTTGAACAGCTCCTCCGCTTTTTCGTGAGCGTCGACATGTTTGTTGAgg                                                            <  2:32159/146‑1 (MQ=255)
                                                                                  aCTGAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCtctttctt                                                                                                          >  5:77274/1‑103 (MQ=255)
                                                                                  aCTGAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCtctttctt                                                                                                          <  6:77274/103‑1 (MQ=255)
                                                                                    tGAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCACCGCTTTTTc                                                                                     <  3:150343/122‑1 (MQ=255)
                                                                                    tGAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCACCGCTTTTTc                                                                                     >  4:150343/1‑122 (MQ=255)
                                                                                    tGAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCACCGCTTTTTCGTGAGCGTCGACCTGgttgt                                                                 >  5:46483/1‑142 (MQ=255)
                                                                                    tGAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCACCGCTTTTTCGTGAGCGTCGACCTGgttgt                                                                 <  6:46483/142‑1 (MQ=255)
                                                                                     gAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCACCGCtttt                                                                                       <  1:296978/119‑1 (MQ=255)
                                                                                     gAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCACCGCtttt                                                                                       >  2:296978/1‑119 (MQ=255)
                                                                                      aaGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCACCGCTTTTTCGt                                                                                   <  1:285224/122‑1 (MQ=255)
                                                                                      aaGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCACCGCTTTTTCGt                                                                                   >  2:285224/1‑122 (MQ=255)
                                                                                      aaGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCACCGCTTTTTCGTGAGCGTCGACCTGGTTGTTGAGGtcc                                                         >  7:141320/1‑146 (MQ=255)
                                                                                               ggCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCACCGCTTTTTCGTGAGCGTCGACCTGGTTGTTGAGGTACGACTTGTCg                                                >  2:162791/1‑148 (MQ=255)
                                                                                                  cTGATCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCACCGCTTTTTCGTGAGCGTCGACCTGGTTGTTGAGGTACGACTTGTCGAAg                                             >  4:247928/1‑148 (MQ=255)
                                                                                                     ttCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGTCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCACCGCTTTTTCGTGAGCGTCGACCTGGTTGTTGAGGTACGACTTGTCGAAGGAc                                          <  3:247928/148‑1 (MQ=255)
                                                                                                     ttCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCACCGCTTTTTCGTGAGCGTCGACCTGGTTGTTGAGGTACGACTTGTCGAAGGACt                                         >  2:298083/1‑149 (MQ=255)
                                                                                                          ggTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCACCGCTTTTTCGTGAGCGTCGACCTGGTTGTTGAGGTACGACTTGTCGAAGGACTCATcg                                    <  4:9315/149‑1 (MQ=255)
                                                                                                                gTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCACCGCTTTTTCGTGAGCGTCGACCTGGTTGTTGAGGTACGACTTATCGAAGGACTCATCGCGCCa                               >  7:155092/1‑148 (MQ=255)
                                                                                                                gTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTCCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCACCGCTTTTTc                                                                                     >  2:3886/1‑94 (MQ=255)
                                                                                                                gTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTCCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCACCGCTTTTTc                                                                                     >  2:9191/1‑94 (MQ=255)
                                                                                                                gTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTCCAGCTCGGCCTTGTCGGACGATGCAGCCTCCTTCTTGAACAGCTCCACCGCTTTTTc                                                                                     <  1:9191/94‑1 (MQ=255)
                                                                                                                gTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTCCAGCTCGGCCTTGTCGGACGATGCAGCCTCCTTCTTGAACAGCTCCACCGCTTTTTc                                                                                     <  1:3886/94‑1 (MQ=255)
                                                                                                                 ttCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCACCGCTTTTTCGTGAGCGTCGACCTGGTTGTTGAGGTACGACTTGTCGAAGGACTCATCGCGCCACt                             >  8:118013/1‑149 (MQ=255)
                                                                                                                 ttCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCACCGCTTTTTCGTGAGCGTCGACCTGGTTGTTGAGGTACGACTTGTCGAAGGACTCATCGCGCCACt                             >  5:147457/1‑149 (MQ=255)
                                                                                                                         tGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCACCGCTTTTTCGTGAGCGTCGACCTGGTTGTTGAGGTACGACTTGTCGAAGGACTCATCGCGCCACTCCAGGAt                      <  7:118013/148‑1 (MQ=255)
                                                                                                                           ggCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGTCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCACCGCTTTTTCGTGAGCGTCGACCTGGTTGTTGAGGTACGACTTTTCGAAGGACTCATCGCGCCACTCCAGGATCa                    <  8:155092/148‑1 (MQ=255)
                                                                                                                                     cGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCACCGCTTTTTCGTGAGCGTCGACCTGGTTGTTGAGGTACGACTTGTCGAAGGACTCATCGCGCCACTCCAGGATCACTTTCTTGAt          <  4:157199/148‑1 (MQ=255)
                                                                                                                                            aaaGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAACTCCACCGCTTTTTCGTGAGCGTCGACCTGGTTGTTGAGGTACGACTTGTCGAAGGACTCATCGCGCCACTCCAGGATCACTTTCTTGATCTTGTCgg  >  2:310021/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
GGGGTGGGCATGCAGCAGGCGGGTAAGTACTTCGGTCTTCAATGCATCGTTACTCATGTGCGGCCTCACTTGCCGTGCTTGGACTGAAGGGCTTTGGCCTGTTCCAGGTGCTGTTCGAGTTTGGGCAGCGTTTCGCTGGCAAAGGCCTTCAGCTCGGCCTTGTCGGACGATGCAGCCTCTTTCTTGAACAGCTCCACCGCTTTTTCGTGAGCGTCGACCTGGTTGTTGAGGTACGACTTGTCGAAGGACTCATCGCGCCACTCCAGGATCACTTTCTTGATCTTGTCGG  >  NC_002947/3676280‑3676568

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