Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,232,045 G→C 9.9% G233G (GGC→GGG)  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,0450GC 9.9% 106.9 / 2.9 38G233G (GGC→GGG) PP_4610PepSY‑associated TM helix domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (15/19);  new base C (2/2);  total (17/21)
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.62e-01

AAGGTAACGCGGGCATTGCTGTCGCCCGGGTTCTGCACCTGGATGTTGCCGACACGCGCCCCCGGCTGCAGCTCCTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGATGGCA  >  NC_002947/5231896‑5232186
                                                                                                                                                     |                                                                                                                                             
aaGGTAACGCGGGCATTGCTGTCGCCCGGGTTCTGCACCTGGATGTTGCCGACACGCGCCCCCGGCTGCAGCTCCTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCgg                                                                                                                                               <  2:625078/150‑1 (MQ=255)
     accgcgGGCATTGCTGTCGCCCGGGTTCTGCACCTGGATGTTGCCGACACGCGCCCCCGGCTGCAGCTCCTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGCCGCAACGATGCGCCCaag                                                                                                                                          <  1:451914/148‑2 (MQ=255)
              aTTGCTGTCGCCCGGGTTCTGCACCTGGATGTTGCCGACACGCGCCCCCGGCTGCAGCTCCTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGtt                                                                                                                                 >  2:667170/1‑150 (MQ=255)
               ttGCTGTCGCCCGGGTTCTGCACCTGGATGTTGCCGACACGCGCCCCCGGCTGCAGCTCCTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTg                                                                                                                                <  2:130203/150‑1 (MQ=255)
                 gCTGTCGCCCTGGTTCTGCACCTGGATGTTGCCGACACGCGCCCCACCCTGCAGCTCCTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGaa                                                                                                                              <  1:520340/150‑1 (MQ=255)
                       gCCCGGCTTCTGCACCTGGATGTTGCCGACACGCGCCCCCGGCTGCAGCTCCTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACtt                                                                                                                        <  2:472990/150‑1 (MQ=255)
                         ccGGGTTCTGCACCTGGATGTTGCCGACACGCGCCCCCGGCTGCAGCTCCTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGt                                                                                                                      >  2:515051/1‑150 (MQ=255)
                                                 cGACACGCGCCCCCGGCTGCAGCTCCTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCACCCTGACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCa                                                                                              >  1:52975/1‑150 (MQ=255)
                                                    cacGCGCCCCCGGCTGCAGCTCCTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCatga                                                                                           <  1:348762/150‑1 (MQ=255)
                                                    cacGCGCCCCCGGCTGCAGCTCCTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCatga                                                                                           <  1:837941/150‑1 (MQ=255)
                                                          cccccGGCTGCAGCTCCTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGCCGCAACGATGCGCCCaagg                                                                                                                                         <  1:194244/98‑3 (MQ=255)
                                                          cccccGGCTGCAGCTCCTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGCCGCAACGATGCGCCCaagg                                                                                                                                         >  2:194244/1‑96 (MQ=255)
                                                             ccGGCTGCAGCTCCTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCccgccg                                                                                  >  1:600266/1‑150 (MQ=255)
                                                                   gCAGCTCCTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGTCATGATGCCGGCGGGCATCa                                                                            <  1:667170/150‑1 (MQ=255)
                                                                     aGCTCCTGCACCTCGGCGTACAGGCTGGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCGGCGGGCATCatg                                                                          <  2:51529/150‑3 (MQ=255)
                                                                      gCTCCTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGAAACATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCa                                                                         >  1:437233/1‑150 (MQ=255)
                                                                       cTCCTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCAt                                                                        <  2:723628/150‑1 (MQ=255)
                                                                           tGCACCTTGGCCTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTAc                                                                    >  1:843618/1‑150 (MQ=255)
                                                                            gCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTGATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACa                                                                   >  2:615616/1‑150 (MQ=255)
                                                                            gCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACa                                                                   >  2:296342/1‑150 (MQ=255)
                                                                             cACCTTGGCGCACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCATTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACAt                                                                  <  1:337188/150‑1 (MQ=255)
                                                                                          aGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAgg                                                     <  2:600266/150‑1 (MQ=255)
                                                                                          aGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGGACATGAACAGCACCAgg                                                     >  2:105051/1‑150 (MQ=255)
                                                                                              tCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCa                                                                                              <  1:241074/105‑1 (MQ=255)
                                                                                              tCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCa                                                                                              >  2:241074/1‑105 (MQ=255)
                                                                                                        tACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGat                                       >  2:709425/1‑150 (MQ=255)
                                                                                                           cAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGt                                                                                                                      >  2:686496/1‑68 (MQ=255)
                                                                                                           cAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGt                                                                                                                      <  1:686496/68‑1 (MQ=255)
                                                                                                           cAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCatga                                                                                           <  1:602798/95‑1 (MQ=255)
                                                                                                           cAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCatga                                                                                           >  2:602798/1‑95 (MQ=255)
                                                                                                              cggcgTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAg                                 <  2:843618/150‑1 (MQ=255)
                                                                                                                         ccTGGGCCGCCTTGGGCGCATCGTTGCGGCCACACACGTGGGAGTAGAACGTCTCGGTTTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGgca                      <  1:192478/150‑1 (MQ=255)
                                                                                                                                 gccTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGcagcaccagc              >  2:45431/1‑150 (MQ=255)
                                                                                                                                       ggCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGt                               >  2:637181/1‑127 (MQ=255)
                                                                                                                                       ggCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGt                               <  1:637181/127‑1 (MQ=255)
                                                                                                                                       ggCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGgc                       <  1:196997/135‑1 (MQ=255)
                                                                                                                                       ggCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGgc                       >  2:196997/1‑135 (MQ=255)
                                                                                                                                        gcgcATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGa       <  2:197763/150‑1 (MQ=255)
                                                                                                                                             tCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGATGGCa  <  1:45431/150‑1 (MQ=255)
                                                                                                                                                     gCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGt                               >  2:63788/1‑113 (MQ=255)
                                                                                                                                                     gCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGt                               <  1:63788/113‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
AAGGTAACGCGGGCATTGCTGTCGCCCGGGTTCTGCACCTGGATGTTGCCGACACGCGCCCCCGGCTGCAGCTCCTGCACCTTGGCGTACAGGCTCGGCAGCGCTACCAGCGGCGTGGCCACCTGGGCCGCCTTGGGCGCATCGTTGCGGCCAAACAGGTCGTTGAAGTACTTGTCGGTGTCATTGCCGTAGCTGGCCATGATGCCCGCCGGCATCACCATGTACATGAACAGCACCAGGCTGCTGTAGCTGATCATCAGGTGAAACGGCAGCACCAGCACGCCGATGGCA  >  NC_002947/5231896‑5232186

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑