Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,232,054 C→G 13.2% A214A (GCG→GCC PP_4610 ← PepSY‑associated TM helix domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,232,0540CG13.2% 139.3 / 12.5 53A214A (GCG→GCCPP_4610PepSY‑associated TM helix domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (20/26);  new base G (2/5);  total (22/31)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.26e-01

CCCCCGGCTGCAGCTCCTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGATGGCATTGTGCCCGTCCAGCCAGGAACGCTGCCCCTTGCCAGGGCGGAAGGTGAAGAACTCCTTGAA  >  NC_002947/5231906‑5232200
                                                                                                                                                    |                                                                                                                                                  
cccccGGCTGCAGCTCCTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGccc                                                                                                                                                    <  6:72339/149‑1 (MQ=255)
 ccccGGCTGCAGCTCCTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCccg                                                                                                                                                   >  2:86554/1‑149 (MQ=255)
  cccGGCTGCAGCTCCTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCccgc                                                                                                                                                  >  2:186111/1‑149 (MQ=255)
      gCTGCAGCTCCTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCgg                                                                                                                                               >  7:148147/1‑148 (MQ=255)
      gCTGCAGCTCCTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGc                                                                                                                                              <  8:101544/149‑1 (MQ=255)
           aGCTCCTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCAc                                                                                                                                         <  1:186111/149‑1 (MQ=255)
            gCTCCTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCGGCGGGCATCatg                                                                                                                                        <  4:38127/149‑3 (MQ=255)
             cTCCTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCGGCGGGCATCatgg                                                                                                                                       <  1:175565/149‑4 (MQ=255)
             cTCCTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCGGCGGGCATCatgg                                                                                                                                       <  2:224978/149‑4 (MQ=255)
                cTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCGGCGGGCATCatggcca                                                                                                                                    <  8:76656/149‑7 (MQ=255)
                   cACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTAc                                                                                                                                  <  8:243060/148‑1 (MQ=255)
                    aCCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTTCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACAt                                                                                                                                <  7:47419/149‑1 (MQ=255)
                     ccTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATg                                                                                                                               <  7:132227/149‑1 (MQ=255)
                     ccTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATg                                                                                                                               <  2:44346/149‑1 (MQ=255)
                      cTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATg                                                                                                                               >  2:224645/1‑148 (MQ=255)
                          gCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACa                                                                                                                           <  1:133634/148‑1 (MQ=255)
                            gTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCa                                                                                                                        >  2:146201/1‑149 (MQ=255)
                               cAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCAc                                                                                                                       >  6:180304/1‑147 (MQ=255)
                                    tCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGgctg                                                                                                                >  8:181204/1‑149 (MQ=255)
                                      ggCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGgctg                                                                                                                <  5:180304/147‑1 (MQ=255)
                                         aGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGctgctg                                                                                                             >  5:111166/1‑147 (MQ=255)
                                                ccAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGatc                                                                                                    <  4:78414/149‑1 (MQ=255)
                                                    cggcgTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGAtcatca                                                                                                <  1:86554/149‑1 (MQ=255)
                                                    cggcgTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGAtcatca                                                                                                <  8:151151/149‑1 (MQ=255)
                                                    cggcgTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGAtcatca                                                                                                <  1:146201/149‑1 (MQ=255)
                                                             cACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCGGCGGGCATCat                                                                                                                                         >  1:222226/1‑98 (MQ=255)
                                                             cACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCGGCGGGCATCat                                                                                                                                         <  2:222226/99‑2 (MQ=255)
                                                             cACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACg                                                                                       >  3:30623/1‑149 (MQ=255)
                                                               ccTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGgc                                                                                     <  4:287020/149‑1 (MQ=255)
                                                                            gggCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGc                                                                        >  6:166535/1‑149 (MQ=255)
                                                                                 cATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGATg                                                                   <  6:111166/149‑1 (MQ=255)
                                                                                  aTCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCt                                                                                                        <  6:274299/111‑1 (MQ=255)
                                                                                  aTCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCt                                                                                                        >  5:274299/1‑111 (MQ=255)
                                                                                  aTCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGATg                                                                   >  8:257615/1‑148 (MQ=255)
                                                                                  aTCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGATg                                                                   <  7:257615/148‑1 (MQ=255)
                                                                                         ctgCCAAACAGGTCGTTGAAGTACTTGTAGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCTGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGATGGCATtgtg                                                           <  3:198976/147‑1 (MQ=255)
                                                                                          ggCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGATGGCATTGTGc                                                          >  4:293851/1‑149 (MQ=255)
                                                                                                    gTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGATGGCATTGTGCCCGTccagcc                                                <  3:293851/149‑1 (MQ=255)
                                                                                                        ttGATGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCGCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGATGGCATTGTGCCCGTCCAGCCAgg                                             >  1:203796/1‑148 (MQ=255)
                                                                                                        ttGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGcc                                                                       >  1:259602/1‑122 (MQ=255)
                                                                                                        ttGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGcc                                                                       <  2:259602/122‑1 (MQ=255)
                                                                                                            aGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGATGGCATTGTGCCCGTCCAGCCAGGAACGc                                        <  7:181204/149‑1 (MQ=255)
                                                                                                             gTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATAATCAGGTGAAACGGCAGCACCAGCACGCCGATGGCATTGTGCCCGTCCAGCCAGGAACGCt                                       >  8:171747/1‑149 (MQ=255)
                                                                                                               aCTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGt                                                                                             >  4:177590/1‑93 (MQ=255)
                                                                                                               aCTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGt                                                                                             <  3:177590/93‑1 (MQ=255)
                                                                                                                   gTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGATGGCATTGTGCCCGTCCAGCCAGGAACGCTGCCCCt                                 >  4:260999/1‑149 (MQ=255)
                                                                                                                    tCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGATGGCATTGTGCCCGTCCAGCCAGGAACGCTGCCCCtt                                <  5:166535/149‑1 (MQ=255)
                                                                                                                     cGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGATGGCATTGTGCCCGTCCAGCCAGGAACGCTGCCCCTTg                               <  4:30623/149‑1 (MQ=255)
                                                                                                                      ggTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGATGGCATTGTGCCCGTCCAGCCAGGAACGCTGCCCCTTGc                              >  5:275546/1‑149 (MQ=255)
                                                                                                                      ggTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGATGGCATTGTGCCCGTCCAGCCAGGAACGCTGCCCCTTGc                              <  1:224645/149‑1 (MQ=255)
                                                                                                                            aTTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGATGGCATTGTGCCCGTCCAGCCAGGAACGCTGCCCCTTGCCAggg                         <  3:260999/148‑1 (MQ=255)
                                                                                                                            aTTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGATGGCATTGTGCCCGTCCAGCCAGGAACGCTGCCCCTTGCCAGGGc                        >  5:48624/1‑149 (MQ=255)
                                                                                                                                ccGTAGCTGGCCATGATGCCCTCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCAACAGCACGCCGAAGGCATTGTGCCCGTCCAGCCAGGACCGCTGACCATTGCCAGGGCgg                      >  1:137876/1‑147 (MQ=255)
                                                                                                                                   tAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGATGGCATTGTGCCCGTCCAGCCAGGAACGCTGCCCCTTGCCAGGGCGGAAGGt                 <  6:48624/149‑1 (MQ=255)
                                                                                                                                         gCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGATGGCATTGTGCCCGTCCAGCCAGGAACGCTGCCCCTTGCCAGGGCGGAAGGTgaagaa           >  4:188677/1‑149 (MQ=255)
                                                                                                                                             tgatgCCGGCGGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGATGGCATTGTGCCCGTCCAGCCAGGAACGCTGCCCCTTGCCAGGGCGGAAGGTGAAGAACTcc       >  1:240276/1‑149 (MQ=255)
                                                                                                                                                  cccGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGATGGCATTGTGCCCGTCCAGCCAGGAACGCTGCCCCTTGCCAGGGCGGAAGGTGAAGAACTCCTTGaa  >  7:194161/1‑149 (MQ=255)
                                                                                                                                                  cccGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGATGGCATTGTGCCCGTCCAGCCAGGAACGCTGCCCCTTGCCAGGGCGGAAGGTGAAGAACTCCTTGa   <  3:217821/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCCCCGGCTGCAGCTCCTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGATGGCATTGTGCCCGTCCAGCCAGGAACGCTGCCCCTTGCCAGGGCGGAAGGTGAAGAACTCCTTGAA  >  NC_002947/5231906‑5232200

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: