Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 846,580 C→T 8.1% G207D (GGC→GAC)  PP_0729 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947846,5800CT 8.1% 135.7 / 3.2 49G207D (GGC→GAC) PP_0729hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (24/21);  new base T (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.13e-01

GCAGATGGTGTTGCTTGCAGATCAGCTCTTCGAGCTGGGAAACCTCGTGCTCCAAAACTTTGCGTAGTTGGGCGTCGGCAAGTGCGCTGGCATTACGGGGCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCG  >  NC_002947/846440‑846726
                                                                                                                                            |                                                                                                                                                  
gCAGATGGTGTTGCTTGCAGATCAGCTCTTCGAGCTGGGAAACCTCGTGCTCCAAAACTTTGCGTAGTTGGGCGTCGGCAAGTGCGCTGGCATTACGGGGCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACa                                                                                                                                            <  1:282049/149‑1 (MQ=255)
gCAGATGGTGTTGCTTGCAGATCAGCTCTTCGAGCTGGGAAACCTCGTGCTCCAAAACTTTGCGTAGTTGGGCGTCGGCAAGTGCGCTGGCATTACGGGGCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACa                                                                                                                                            <  1:286744/149‑1 (MQ=255)
     tGGTGTTGCTTGCAGATCAGCTCTTCGAGCTGGGAAACCTCGTGCTCCAAAACTTTGCGTAGTTGGGCGTCGGCAAGTGCGCTGGCATTACGGGGCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACAGTGAt                                                                                                                                       >  4:143467/1‑149 (MQ=255)
        tgtTGCTTGCAGATCAGCTCTTCGAGCTGGGAAACCTCGTGCTCCAAAACTTTGCGTAGTTGGGCGTCGGCAAGTGCGCTGGCATTACGGGGCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATcgc                                                                                                                                    <  4:51692/149‑1 (MQ=255)
         gtTGCTTGCAGATCAGCTCTTCGAGCTGGGAAACCTCGTGCTCCAAAACTTTGCGTAGTTGGGCGTCGGCAAGTGCGCTGGCATTACGGGGCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATcgcg                                                                                                                                   >  4:216598/1‑149 (MQ=255)
         gctgcatgcaGATCAGCTCCTCGAGCTGGGAAACCTCGTGCTCCAAAACTTTGCGTAGTTGGGCGTCGGCAAGTGCGCTGGCATTACGGGGCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATcgcg                                                                                                                                   <  6:86541/147‑1 (MQ=255)
              ttgcAGATCAGCTCTTCGAGCTGGGAAACCTCGTGCTCCAAAACTTTGCGTAGTTGGGCGTCGGCAAGTGCGCTGGCATTACGGGGCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGTCATCGACAGTGATCGCGAActc                                                                                                                              >  3:19415/1‑149 (MQ=255)
                gcAGATCAGCTCTTCGAGCTGGGAAACCTCGTGCTCCAAAACTTTGCGTAGTTGGGCGTCGGCAAGTGCGCTGGCATTACGGGGCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTa                                                                                                                            >  1:157277/1‑149 (MQ=255)
                 cAGATCAGCTCTTCGAGCTGGGAAACCTCGTGCTCCAAAACTTTGCGTAGTTGGGCGTCGGCAAGTGCGCTGGCATTACGGGGCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTAt                                                                                                                           <  2:157277/149‑1 (MQ=255)
                   gATCAGCTCTTCGAGCTGGGAAACCTCGTGCTCCAAAACTTTGCGTAGTTGGGCGTCGGCAAGTGCGCTGGCATTACGGGGCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGt                                                                                                                         >  5:197596/1‑149 (MQ=255)
                    aTCAGCTCTTCGAGCTGGGAAACCTCGTGCTCCAAAACTTTGCGTAGTTGGGCGTCGGCAAGTGCGCTGGCATTACGGGGCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGTCATCGACAGTGATCGCGAACTCTATGTc                                                                                                                        >  5:69231/1‑149 (MQ=255)
                    aTCAGCTCTTCGAGCTGGGAAACCTCGTGCTCCAAAACTTTGCGTAGTTGGGCGTCGGCAAGTGCGCTGGCATTACGGGGCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTc                                                                                                                        >  3:96103/1‑149 (MQ=255)
                              cGAGCTGGGAAACCTCGTGCTCCAAAACTTTGCGTAGTTGGGCGTCGGCAAGTGCGCTGGCATTACGGGGCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATc                                                                                                               <  4:90336/148‑1 (MQ=255)
                                    gggAAACCTCGTGCTCCAAAACTTTGCGTAGTTGGGCGTCGGCAAGTGCGCTGGCATTACGGGGCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCg                                                                                                        >  3:201918/1‑149 (MQ=255)
                                          ccTCGTGCTCCAAAACTTTGCGTAGTTGGGCGTCGGCAAGTGCGCTGGCATTACGGGGCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAg                                                                                                   >  1:288130/1‑148 (MQ=255)
                                                    cAAAACTTTGCGTAGTTGGGCGTCGGCAAGTGCGCTGGCATTACGGGGCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTAc                                                                                          <  4:96103/147‑1 (MQ=255)
                                                    cAAAACTTTGCGTAGTTGGGCGTCGGCAAGTGCGCTGGCATTACGGGGCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACct                                                                                        <  7:29132/149‑1 (MQ=255)
                                                    cAAAACTTTGCGTAGTTGGGCGTCGGCAAGTGCGCTGGCATTACGGGGCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACAGTCATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGATAGAAATTGTAcc                                                                                         >  1:274570/1‑148 (MQ=255)
                                                       aaCTTTGCGTAGTTGGGCGTCGGCAAGTGCGCTGGCATTACGGGGCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGTCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCtt                                                                                     <  6:69231/149‑1 (MQ=255)
                                                         cTTTGCGTAGTTGGGCGTCGGCAAGTGCGCTGGCATTACGGGGCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCg                                                                                   <  5:132941/149‑1 (MQ=255)
                                                                     gggCGTCGGCAAGTGCGCTGGCATTACGGGGCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCTCGCCCGGTTAGAACTTGTACCTCTTCGa                                                                                  <  4:71285/138‑1 (MQ=255)
                                                                     gggCGTCGGCAAGTGCGCTGGCATTACGGGGCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGa                                                                                  >  3:71285/1‑138 (MQ=255)
                                                                         gTCGGCAAGTGCGCTGGCATTACGGGGCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAAcc                                                                   >  2:148836/1‑149 (MQ=255)
                                                                          tCGGCAAGTGCGCTGGCATTACGGGGCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAAccc                                                                  <  3:143467/149‑1 (MQ=255)
                                                                             gCAAGTGCGCTGGCATTACGGGGCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTgg                                                               <  4:201918/149‑1 (MQ=255)
                                                                                     gcTGGCATTACGGGGCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCg                                                       >  2:56380/1‑149 (MQ=255)
                                                                                      cTGGCATTACGGGGCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGc                                                      <  1:148836/149‑1 (MQ=255)
                                                                                         gCATTACGGGGCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAgg                                                   >  3:252655/1‑149 (MQ=255)
                                                                                           aTTACGGGGCCCATCTGCGGTCATGGGGCGTACTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCt                                                 >  1:15834/1‑149 (MQ=255)
                                                                                            ttACGGGGCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTc                                                >  2:248021/1‑149 (MQ=255)
                                                                                                ggggCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCtt                                            >  2:287258/1‑149 (MQ=255)
                                                                                                       aTCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCtt                                     >  3:273889/1‑149 (MQ=255)
                                                                                                            cGGTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTg                                >  8:79299/1‑149 (MQ=255)
                                                                                                             ggTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGc                               <  6:197596/149‑1 (MQ=255)
                                                                                                                 aTGGGACGTACTTGGAATAGCAGAGGGTCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCaa                            <  4:19415/148‑1 (MQ=255)
                                                                                                                  tGGGACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAgg                          <  3:216598/149‑1 (MQ=255)
                                                                                                                   gggACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGc                                          >  7:61023/1‑132 (MQ=255)
                                                                                                                   gggACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGc                                          <  8:61023/132‑1 (MQ=255)
                                                                                                                   gggACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAgg                          >  4:16754/1‑148 (MQ=255)
                                                                                                                        gTACTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGAc                     >  5:12751/1‑148 (MQ=255)
                                                                                                                           cTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGaaa                 <  3:16754/149‑1 (MQ=255)
                                                                                                                            ttGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAg                >  8:29510/1‑149 (MQ=255)
                                                                                                                               gAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCAt             <  4:252655/149‑1 (MQ=255)
                                                                                                                               gAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCAt             <  5:258651/149‑1 (MQ=255)
                                                                                                                               gAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCAt             <  6:12751/149‑1 (MQ=255)
                                                                                                                               gAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCAt             >  7:101050/1‑149 (MQ=255)
                                                                                                                                       gagGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGtc     >  4:276134/1‑149 (MQ=255)
                                                                                                                                       gagGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGtc     <  1:56380/149‑1 (MQ=255)
                                                                                                                                          ggCCATCGACAGTGATCGCGAACTCTATGTCCAAGACATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCg  >  1:49281/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
GCAGATGGTGTTGCTTGCAGATCAGCTCTTCGAGCTGGGAAACCTCGTGCTCCAAAACTTTGCGTAGTTGGGCGTCGGCAAGTGCGCTGGCATTACGGGGCCCATCTGCGGTCATGGGACGTACTTGGAATAGCAGAGGGCCATCGACAGTGATCGCGAACTCTATGTCCAAGGCATCGCGCCCGGTTAGAACTTGTACCTCTTCGAGCAAAGCAAACAACCCTGGTAATTCCGCAGGCTCGCTTGCCGCTTCGCTGCCAAGGGAGACGAAAGCATGCACAGTCTCG  >  NC_002947/846440‑846726

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