Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,700,170 C→A 16.0% G106C (GGC→TGC)  ybhP‑I ← phosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,700,1700CA16.0% 59.4 / 2.8 25G106C (GGC→TGC) ybhP‑Iphosphohydrolase
Reads supporting (aligned to +/- strand):  ref base C (10/11);  new base A (2/2);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.84e-01

CCGCATGCACCTCGTCATGGCCCGGCACCTCCAGCTGGCAATGCAGCAGGCCACGCTGTTCATTGCCCTGCACCGAAACATCCAGGTTGTTGAAGCGCCTGATGGGGAATTTCGACAACAGCGCGTTGCCGTGATCGCCATGTGGATACACAGCATTGCGGCCATAGGCAAACTGCGGCCACATGCTGTCGGCCAGAAACTCGTACTGCGGCGTTTCGGGCCAGGCAGGGTGGCGCAAGGCGTGCTGCTGGTGACTGCCGTGCACCTCTTGCAAGAA  >  NC_002947/3700041‑3700317
                                                                                                                                 |                                                                                                                                                   
ccGCATGCACCTCGTCATGGCCCGGCACCTCCAGCTGGCAATGCAGCAGGCCACGCTGTTCATTGCCCTGCACCGAAACATCCAGGTTGTTGAAGCGCCTGATGGGGAATTTCGACAACAGCGCGTTGCCGTGATCGCCATGTGGATaca                                                                                                                                 >  1:480415/1‑150 (MQ=255)
          cTCGTCATGGCCCGGCACCTCCAGCTGGCAATGCAGCAGGCCACGCTGTTCATTGCCCTGCACCGAAACATCCAGGTTGTTGAAGCGCCTGATGGGGAATTTCGACAACAGCGCGTTGCCGTGATCGCCATGTGGATACACAGCATTGCg                                                                                                                       <  1:7787/150‑1 (MQ=255)
                     ccGGCACCTCCAGCTGGCAATGCAGCAGGCCACGCTGTTCATTGCCCTGCACCGAAACATCCAGGTTGTTGAAGCGCCTGATGGGGAATTTCGACAACAGCGCGTTGCCGTGATCGCCATGTGGATACACAGCATTGCGGCCATAGGCaa                                                                                                            >  1:422765/1‑150 (MQ=255)
                                aGCTGGCAATGCAGCAGGCCACGCTGTTCATTGCCCTGCACCGAAACATCCAGGTTGTTGAAGCGCCTGATGGGGAATTTCGACAACAGCGCGTTGCCGTGATCGCCATGTGGATa                                                                                                                                   <  2:56555/116‑1 (MQ=255)
                                aGCTGGCAATGCAGCAGGCCACGCTGTTCATTGCCCTGCACCGAAACATCCAGGTTGTTGAAGCGCCTGATGGGGAATTTCGACAACAGCGCGTTGCCGTGATCGCCATGTGGATa                                                                                                                                   >  1:56555/1‑116 (MQ=255)
                                aGCTGGCAATGCAGCAGGCCACGCTGTTCATTGCCCTGCACCGAAACATCCAGGTTGTTGAAGCGCCTGATGGGGAATTTCGACAACAGCGCGTTGCCGTGATCGCCATGTGGATACACAGCATTGCGGCCATAGGCAAACTGCGGCcac                                                                                                 <  1:441412/150‑1 (MQ=255)
                                         tGCAGCAGGCCACGCTGTTCATTGCCCTGCACCGAAACATCCAGGTTGTTGAAGCGCCTGATGGGGAATTTCGACAACAGCGCGTTGCCGTGATCGCCATGTGGATACACAGCATTGCGGCCATAGGCAAACTGCGGCCACATGCTGTCg                                                                                        <  2:23222/150‑1 (MQ=255)
                                           cagcagGCCACGCTGTTCATTGCCCTGCACCGAAACATCCAGGTTGTTGAAGCGCCTGATGGGGAATTTCGACAACAGCGCGTTGCCGTGATCGCCATGTGGATACACAGCATTGCGGCCATAGGCAAACTGCGGCCACATGCTGTCGGc                                                                                      <  1:352777/150‑1 (MQ=255)
                                                   cACGCTGTTCATTGCCCTGCACCGAAACATCCAGGTTGTTGAAGCGCCTGATGGGGAATTTCGACAACAGCGCGTTCCCGTGCTCCCCACGTGGCTCGACAGCATTGGGGCCAGAGGCAAACCGCGGCCACCTGGTGTCGGCCAGaaccg                                                                              >  2:30629/1‑147 (MQ=255)
                                                          ttCATTGCCCTGCACCGAAACATCCAGGTTGTTGAAGCGCCTGATGGGGAATTTCGACAACAGCGCGTTGCCGTGATCGCCATGTGGATACACAGCATTGCGGCCATAGGCAAACTGCGGCCACATGCTGTCGGCCAGAAACTCGTACTg                                                                       >  1:250026/1‑150 (MQ=255)
                                                             aTTGCCCTGCACCGAAACATCCAGGTTGTTGAAGCGCCTGATGGGGAATTTCGACAACAGCGCGTTGCCGTGATCGCCATGTGGATACACAGCATTGCGGCCATAGGCAAACTGCGGCCACATGCTGTCGGCCAGAAACTCGTACTgcgg                                                                    >  1:67783/1‑150 (MQ=255)
                                                                                                       ggggAATTTCGACAACAGCGCGTTGCAGTGATCGCCATGTGGATACACAGCATTGCGGCCATAGGCAAACTGCGGCCACATGCTGTCGGCCAGAAACTCGTACTGCGGCGTTTCGGGCCAGGCAGGGTGGCGCAAGGCGTGCTGCTgg                            <  1:46370/148‑1 (MQ=255)
                                                                                                       ggggAATTTCGACAACAGCGCGTTGCAGTGATCGCCATGTGGATACACAGCATTGCGGCCATAGGCAAACTGCGGCCACATGCTGTCGGCCAGAAACTCGTACTGCGGCGTTTCGGGCCAGGCAGGGTGGCGCAAGGCGTGCTGCTgg                            <  1:46489/148‑1 (MQ=255)
                                                                                                       ggggAATTTCGACAACAGCGCGTTGCAGTGATCGCCATGTGGATACACAGCATTGCGGCCATAGGCAAACTGCGGCCACATGCTGTCGGCCAGAAACTCGTACTGCGGCGTTTCGGGCCAGGCAGGGTGGCGCAAGGCGTGCTGCTgg                            >  2:46489/1‑148 (MQ=255)
                                                                                                       ggggAATTTCGACAACAGCGCGTTGCAGTGATCGCCATGTGGATACACAGCATTGCGGCCATAGGCAAACTGCGGCCACATGCTGTCGGCCAGAAACTCGTACTGCGGCGTTTCGGGCCAGGCAGGGTGGCGCAAGGCGTGCTGCTgg                            >  2:46370/1‑148 (MQ=255)
                                                                                                         ggAATTTCGACAACAGCGCGTTGCCGTGATCGCCATGTGGATACACAGCATTGCGGCCATAGGCAAACTGCGGCCACATGCTGTCGGCCAGAAACTCGTACTGCGGCGTTTCGGGCCAGGCAGGGTGGCGCAAGGCGtgc                                  >  1:224861/1‑140 (MQ=255)
                                                                                                         ggAATTTCGACAACAGCGCGTTGCCGTGATCGCCATGTGGATACACAGCATTGCGGCCATAGGCAAACTGCGGCCACATGCTGTCGGCCAGAAACTCGTACTGCGGCGTTTCGGGCCAGGCAGGGTGGCGCAAGGCGtgc                                  <  2:224861/140‑1 (MQ=255)
                                                                                                                 gACAACAGCGCGTTGCCGTGATCGCCATGTGGATACACAGCATTGCGGCCATAGGCAAACTGCGGCCACATGCTGTCGGCCAGAAACTCGTACTGCGGCGTTTCGGGCCAGGCAGGGTGGCGCAAGGCGTGCTGCTGGTGACTGCCGTGc                >  2:406898/1‑150 (MQ=255)
                                                                                                                  acaacaGCGCGTTGCCGTGATCGCCATGTGGATACACAGCATTGCGGCCATAGGCAAACTGCGGCCACATGCTGTCGGCCAGAAACTCGTACTGCGGCGTTTCGGGCCAGGCAGGGTGGCGCAAGGCGTGCTGCTGGTGACTGCCGTGCa               <  1:374299/150‑1 (MQ=255)
                                                                                                                       aGCGCGTTGCCGTGATCGCCATGTGGATACACAGCAtt                                                                                                                          <  1:451117/38‑1 (MQ=255)
                                                                                                                       aGCGCGTTGCCGTGATCGCCATGTGGATACACAGCAtt                                                                                                                          >  2:451117/1‑38 (MQ=255)
                                                                                                                          gcgTTGCCGTGATCGCCATGTGGATACACAGCATTGCGGCCATAGGCAAACTGCGGCCACATGCTGTCGGCCAGAAACTCGTACTGCGGCGTTTCGGGCCAGGCAGGGTGGCGCAAGGCGTGCTGCTGGTGACTGCCGTGCACCTCTTGc       <  1:406898/150‑1 (MQ=255)
                                                                                                                             ttGCCGTGATCGCCATGTGGATACACAGCATTGCGGCCATAGGCAAACTGCGGCCACATGCTGTCGGCCAGAAACTCGTACTGCGGCGTTTCGGGCCAGGCAGGGTGGCGCAAGGCGTGCTGCTGGt                           >  2:131263/1‑127 (MQ=255)
                                                                                                                             ttGCCGTGATCGCCATGTGGATACACAGCATTGCGGCCATAGGCAAACTGCGGCCACATGCTGTCGGCCAGAAACTCGTACTGCGGCGTTTCGGGCCAGGCAGGGTGGCGCAAGGCGTGCTGCTGGt                           <  1:131263/127‑1 (MQ=255)
                                                                                                                               gCCGTGATCGCCATGTGGATACACAGCATTGCGGCCATAGGCAAACTGCGGCCACATGCTGTCGGCCAGAAACTCGTACTGCGGCGTTTCGGGCCAGGCAGGGTGGCGCAAGGCGTGCTGCTGGTGACTGCCGTGCACCTCTTGCAAGaa  <  1:30629/150‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                                   
CCGCATGCACCTCGTCATGGCCCGGCACCTCCAGCTGGCAATGCAGCAGGCCACGCTGTTCATTGCCCTGCACCGAAACATCCAGGTTGTTGAAGCGCCTGATGGGGAATTTCGACAACAGCGCGTTGCCGTGATCGCCATGTGGATACACAGCATTGCGGCCATAGGCAAACTGCGGCCACATGCTGTCGGCCAGAAACTCGTACTGCGGCGTTTCGGGCCAGGCAGGGTGGCGCAAGGCGTGCTGCTGGTGACTGCCGTGCACCTCTTGCAAGAA  >  NC_002947/3700041‑3700317

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