Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,449,172 C→G 10.2% P14R (CCC→CGC) ‡ ybeZ → nucleoside triphosphate hydrolase domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,449,1720CG10.2% 137.2 / 7.7 49P14R (CCC→CGC) ‡ybeZnucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (23/21);  new base G (2/3);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.13e-01

GCCTGTACATACCGTTGAACGTTTCAGTGCTTTCGCCACTGCCATTCAAGGGTTATTCTTGAGTTCAACTCCATCGCCGTCGAGCGGCTAAAAACGACCTTGAACGCACCCATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGCCATCGAGATCCGCAACCGCGGCAATCAGTTCGAACTGATCGGCGAACCCAGGACCACCTCCGCTGCAGAA  >  NC_002947/5449033‑5449317
                                                                                                                                           |                                                                                                                                                 
gCCTTTACATACCGTTGCACGTTTCAGTGCTTTCGCCACTGCCATTCAAGGGTTATTCTTGAGTTCAACTCCATCGCCGTCGAGCGGCTAAAAACGACCTTGAACGCACCCATAGAACCCCATCGTTTCATCCTCGAAGGGTTCGAgga                                                                                                                                          <  6:85290/149‑2 (MQ=255)
gCCTGTACATACCGTTGAACGTTTCAGTGCTTTCGCCACTGCCATTCAAGGGTTATTCTTGAGTTCAACTCCATCGCCGTCGAGCGGCTAAAAACGACCTTGAACGCACCCATAGAACCCCATCGTTTCATCCTCGAAGGGTTCGAgga                                                                                                                                          <  6:90367/149‑2 (MQ=255)
 ccTGTAGATAACGTTGAACGTTTCAGTGCTTTCGCCACTGCCATTCAAGGGTTATTCTTGAGTTCAACTCCATCGCCGTCGAGCGGCTAAAAACGACCTTGAACGCACCCATAGAACCCCATCGTTTCATCCTCGAAGGGTTCGAggat                                                                                                                                         <  3:189232/149‑3 (MQ=255)
                     tttCAGTGCTTTCGCCACTGCCATTCAAGGGTTATTCTTGAGTTCAACTCCATCGCCGTCGAGCGGCTAAAAACGACCTTGAACGCACCCATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTtgtg                                                                                                                     <  5:285358/149‑1 (MQ=255)
                              tttCGCCACTGCCATTCAAGGGTTATTCTTGAGTTCAACTCCATCGCCGTCGAGCGGCTAAAAACGACCTTGAACGCACCCATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGtt                                                                                                            >  7:251248/1‑149 (MQ=255)
                                  gCCACTGCCATTCAAGGGTTATTCTTGAGTTCAACTCCATCGCCGTCGAGCGGCTAAAAACGACCTTGAACGCACCCATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCAAACTTGTGCGGCCAGATagac                                                                                                        >  4:230544/1‑149 (MQ=255)
                                                aGGGTTATTCTTGAGTTCAACTCCATCGCCGTCGAGCGGCTAAAAACGACCTTGAACGCACCCATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCt                                                                                          >  1:88396/1‑149 (MQ=255)
                                                    ttATTCTTGAGTTCAACTCCATCGCCGTCGAGCGGCTAAAAACGACCTTGAACGCACCCATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGa                                                                                        >  6:204964/1‑147 (MQ=255)
                                                    ttATTCTTGAGTTCAACTCCATCGCCGTCGAGCGGCTAAAAACGACCTTGAACGCACCCATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGa                                                                                        <  5:204964/147‑1 (MQ=255)
                                                          ttGAGTTCAACTCCATCGCCGTCGAGCGGCTAAAAACGACCTTGAACGCACCCATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAg                                                                                                                                            >  5:54861/1‑89 (MQ=255)
                                                          ttGAGTTCAACTCCATCGCCGTCGAGCGGCTAAAAACGACCTTGAACGCACCCATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAg                                                                                                                                            <  6:54861/89‑1 (MQ=255)
                                                          ttGAGTTCAACTCCATCGCCGTCGAGCGGCTAAAAACGACCTTGAACGCACCCATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAg                                                                                >  7:99075/1‑149 (MQ=255)
                                                             aGTTCAACTCCATCGCCGTCGAGCGGCTAAAAACGACCTTGAACGCACCCATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGaa                                                                                   <  3:56966/143‑1 (MQ=255)
                                                             aGTTCAACTCCATCGCCGTCGAGCGGCTAAAAACGACCTTGAACGCACCCATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGaa                                                                                   >  4:56966/1‑143 (MQ=255)
                                                               ttCAACTCCATCGCCGTCGAGCGGCTAAAAACGACCTTGACCGCACCCATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCt                                                                           >  7:43875/1‑149 (MQ=255)
                                                               ttCAACTCCATCGCCGTCGAGCGGCTAAAAACGACCTTGAACGCACCCATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCt                                                                           >  5:151898/1‑149 (MQ=255)
                                                               ttCAACTCCATCGCCGTCGAGCGGCTAAAAACGACCTTGAACGCACCCATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCt                                                                           >  2:184698/1‑149 (MQ=255)
                                                                 cAACTCCATCGCCGTCGAGCGGCTAAAAACGACCTTGAACGCACCCATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCTAACAGCGCCTg                                                                          <  8:43875/148‑1 (MQ=255)
                                                                  aaCTCCATCGCCGTCGAGCGGCTAAAAACGACCTTGAACGCACCCATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTgg                                                                         <  2:88396/148‑1 (MQ=255)
                                                                  aaCTCCATCGCCGTCGAGCGGCTAAAAACGACCTTGAACGCACCCATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGc                                                                        >  3:319010/1‑149 (MQ=255)
                                                                    cTCCATCGCCGTCGAGCGGCTAAAAACGACCTTGAACGCACCCATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGc                                                                        >  4:303535/1‑147 (MQ=255)
                                                                    cTCCATCGCCGTCGAGCGGCTAAAAACGACCTTGAACGCACCCATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGCCa                                                                      >  1:88747/1‑149 (MQ=255)
                                                                       cATCGCCGTCGAGCGGCTAAAAACGACCTTGAACGCACCCATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGCCATc                                                                    <  3:230544/148‑1 (MQ=255)
                                                                        aTCGCCGTCGAGCGGCTAAAAACGACCTTGAACGCACCCATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGCCATc                                                                    >  7:57580/1‑147 (MQ=255)
                                                                           gCCGTCGAGCGGCTAAAAACGACCTTGAACGCACCCATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGCCATCgaga                                                                <  4:319010/148‑1 (MQ=255)
                                                                              gTCGAGCGGCTAAAAACGACCTTGAACGCACCCATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGCCATCGAGATCCg                                                            >  4:189527/1‑149 (MQ=255)
                                                                                       cTAAAAACGACCTTGAACGCACCCATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGCCATCGAGATCCGCAACCGCgg                                                   <  8:99075/149‑1 (MQ=255)
                                                                                            aaCGACCTTGAACGCACCCATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGCCATCGAGATCCGCAACCGCGGCAAt                                               >  7:69083/1‑148 (MQ=255)
                                                                                                     gAACGCACCCATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGCCATCGAGATCCGCAACCGCGGCAATCAGTTCGaa                                      <  1:184698/148‑1 (MQ=255)
                                                                                                          cACCCATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGCCATCGAGATCCGCAACCGCGGCAATCAGTTCGAACTGATc                                <  6:151898/149‑1 (MQ=255)
                                                                                                          cACCCATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGCCATCGAGATCCGCAACCGCGGCAATCAGTTCGAACTGATc                                <  8:69083/149‑1 (MQ=255)
                                                                                                           aCCCATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGCCATCGAGATCCGCAACCGCGGCAATCAGTTCGAACTGATCg                               <  8:251248/149‑1 (MQ=255)
                                                                                                             ccATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGCCATCGAGATCCGCAACCGCGGCAAt                                               >  8:208716/1‑131 (MQ=255)
                                                                                                             ccATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGCCATCGAGATCCGCAACCGCGGCAAt                                               <  7:208716/131‑1 (MQ=255)
                                                                                                              cATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGCCATc                                                                    <  4:327004/109‑1 (MQ=255)
                                                                                                              cATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGCCATc                                                                    >  3:327004/1‑109 (MQ=255)
                                                                                                                  gAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGCCATCGAGATcc                                                             >  1:168522/1‑112 (MQ=255)
                                                                                                                  gAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGCCATCGAGATcc                                                             <  2:168522/112‑1 (MQ=255)
                                                                                                                   aaCCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGCCATCGAGATCCGCAACCGCGGCAATCAGTTCGAACTGATCGGCGAAccc                       <  2:88747/149‑1 (MQ=255)
                                                                                                                     ccccATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGCCATCGAGATCCGCAACCGCGGCAATCAGTTCGAACTGATCGGCGAACCCAg                     <  8:57580/149‑1 (MQ=255)
                                                                                                                              gtgggccTCGAAGGGTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGCCATCGAGATCCGCAACCGCGGCAATCAGTTCGAACTGATCGGCGAACCCAGGACCACCTc            >  6:322244/6‑149 (MQ=255)
                                                                                                                                cATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGCCATCGAGATCCGCAACCGCGGCAATCAGTTCGAACTGATCg                               <  5:170169/128‑1 (MQ=255)
                                                                                                                                cATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGCCATCGAGATCCGCAACCGCGGCAATCAGTTCGAACTGATCg                               >  6:170169/1‑128 (MQ=255)
                                                                                                                                 ggccTCGAAGGGTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGCCATCGAGATCCGCAACCGCGGCAATCAGTTCGAACTGATCGGCGAACCCAGGACCACCTCCGc         >  4:54700/3‑149 (MQ=255)
                                                                                                                                 aTCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGCCATCGAGATCCGCAACCGCGGCAATCAGTTCGAACTGATCGGCGAACCCAGGACCACCTCCg          <  3:189527/148‑1 (MQ=255)
                                                                                                                                     tCGAACCCTTCGAGGCCGACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGCCATCGAGATCCGCAACCGCGGCAATCAGTTCGAACTGATCGGCGAACCCAGGACCACCTc            <  8:103232/142‑1 (MQ=255)
                                                                                                                                     tCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGCCATCGAGATCCGCAACCGCGGCAATCAGTTCGAACTGATCGGCGAACCCAGGACCACCTc            >  7:103232/1‑142 (MQ=255)
                                                                                                                                       gAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGCCATCGAGATCCGCAACCGCGGCAATCAGTTCGAACTGATCGGCGAACCCAGGACCACCTCCGCTGCAg    >  6:4835/1‑148 (MQ=255)
                                                                                                                                        aaCCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGCCATCGAGATCCGCAACCGCGGCAATCAGTTCGAACTGATCGGCGAACCCAGGACCACCTCCGCTGCAGaa  >  6:38607/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                 
GCCTGTACATACCGTTGAACGTTTCAGTGCTTTCGCCACTGCCATTCAAGGGTTATTCTTGAGTTCAACTCCATCGCCGTCGAGCGGCTAAAAACGACCTTGAACGCACCCATAGAACCCCATCGTTTCATCCTCGAACCCTTCGAGGCCCACCGTTTCGCCAACTTGTGCGGCCAGTTCGACGAGCACCTGCGCCTGATCGAACAGCGCCTGGCCATCGAGATCCGCAACCGCGGCAATCAGTTCGAACTGATCGGCGAACCCAGGACCACCTCCGCTGCAGAA  >  NC_002947/5449033‑5449317

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