Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,232,057 C→G 12.7% P213P (CCG→CCC)  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,0570CG12.7% 146.8 / 12.1 55P213P (CCG→CCC) PP_4610PepSY‑associated TM helix domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (21/27);  new base G (2/5);  total (23/32)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.18e-01

GCTGCAGCTCCTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGATGGCATTGTGCCCGTCCAGCCAGGAACGCTGCCCCTTGCCAGGGCGGAAGGTGAAGAACTCCTTGAA  >  NC_002947/5231912‑5232200
                                                                                                                                                 |                                                                                                                                               
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                                                                                                                               <  2:44346/149‑1 (MQ=255)
               ccTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATg                                                                                                                               <  7:132227/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                                                                                                <  1:146201/149‑1 (MQ=255)
                                              cggcgTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGAtcatca                                                                                                <  8:151151/149‑1 (MQ=255)
                                                       cACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCGGCGGGCATCat                                                                                                                                         <  2:222226/99‑2 (MQ=255)
                                                       cACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCGGCGGGCATCat                                                                                                                                         >  1:222226/1‑98 (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                                                                       <  2:259602/122‑1 (MQ=255)
                                                                                                  ttGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGcc                                                                       >  1:259602/1‑122 (MQ=255)
                                                                                                      aGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGATGGCATTGTGCCCGTCCAGCCAGGAACGc                                        <  7:181204/149‑1 (MQ=255)
                                                                                                       gTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATAATCAGGTGAAACGGCAGCACCAGCACGCCGATGGCATTGTGCCCGTCCAGCCAGGAACGCt                                       >  8:171747/1‑149 (MQ=255)
                                                                                                         aCTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGt                                                                                             <  3:177590/93‑1 (MQ=255)
                                                                                                         aCTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGt                                                                                             >  4:177590/1‑93 (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)
                                                                                                                                          ccagctacgGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGATGGCATTGTGCCCGTCCAGCCAGGAACGCTGCCCCTTGCCAGGGCGGAAGGTGAAGAACTCCTTg    <  6:161467/142‑1 (MQ=255)
                                                                                                                                            cccGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGATGGCATTGTGCCCGTCCAGCCAGGAACGCTGCCCCTTGCCAGGGCGGAAGGTGAAGAACTCCTTGaa  >  7:194161/1‑149 (MQ=255)
                                                                                                                                            cccGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGATGGCATTGTGCCCGTCCAGCCAGGAACGCTGCCCCTTGCCAGGGCGGAAGGTGAAGAACTCCTTGa   <  3:217821/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GCTGCAGCTCCTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGATGGCATTGTGCCCGTCCAGCCAGGAACGCTGCCCCTTGCCAGGGCGGAAGGTGAAGAACTCCTTGAA  >  NC_002947/5231912‑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: ‑