Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,251,032 C→A 20.0% G26V (GGT→GTT)  cidA ← murein hydrolases holin regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,251,0320CA20.0% 37.4 / 4.9 20G26V (GGT→GTT) cidAmurein hydrolases holin regulator
Reads supporting (aligned to +/- strand):  ref base C (6/10);  new base A (2/2);  total (8/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.35e-01

AGCCATCAGCCAGCCGGCCCCCAATTGCAGCATGGCCGGTTTGAGCACACCCGAAGCGAACAGTATCAACAGCAGTGCCAGGCCCATCACACCACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCACCAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATAATGCAGGTTTCATGGCGGGGTTCCTCGACAGGCCCTCATTTTAAGATTTGGCTGCCTAGGCCAGAAGCGAATTGTTAGACTGCTC  >  NC_002947/5250899‑5251180
                                                                                                                                     |                                                                                                                                                    
agccatcagccaGCCGGCCCCCAATTGCAGCATGGCCGGTTTGAGCACACCCGAAGCGAACAGTATCAACAGCAGTGCCAGGCCCATCACACCACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCACCAAGCAGGAACAgcgc                                                                                                                                      >  2:53796/1‑150 (MQ=255)
   catcagccaGCCGGCCCCCAATTGCAGCATGGCCGGTTTGAGCACACCCGAAGCGAACAGTATCAACAGCAGTGCCAGGCCCATCACACAACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCACCAAGCAGGAACAGCGcaca                                                                                                                                   >  1:83008/1‑150 (MQ=255)
    atcagccaGCCGGCCCCCAATTGCTGCATGGCCGGTTTGAGCACACCCGAAGCGAACAGTATCAACAGCAGTGCCAGGCCCATCACACCACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCACCAAGCAGGAACAGCGCACag                                                                                                                                  >  1:176200/1‑150 (MQ=255)
            gccGGCCCCCAATTGCAGCATGGCCGGTTTGAGCACACCCGAAGCGAACAGTATGAACAGCAGTGCCACGCTCATCACTGCACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCACCAAGCAGGAACAGCGCACAGAGAACCGc                                                                                                                          <  2:147928/150‑1 (MQ=255)
                         tGCAGCATGGCCGGTTTGAGCACACCCGAAGCGAACAGTATCAACAGCAGTGCCAGGCCCATCACACAACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCACCAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACaaac                                                                                                             <  2:83008/150‑3 (MQ=255)
                                           aGCACACCCGAAGCGAACAGTATCAACAGCAGTGCCAGGCCCATCACACCACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCAACAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTc                                                                                           >  1:295509/1‑150 (MQ=255)
                                           aGCACACCCGAAGCGAACAGTATCAACAGCAGTGCCAGGCCCATCACACCACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCAACAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTc                                                                                           >  1:295473/1‑150 (MQ=255)
                                                    gAAGCGAACAGTATCAACAGCAGTGCCAGGCCCATCACACCACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCACCAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATAATGCa                                                                                  <  1:129213/150‑1 (MQ=255)
                                                       gCGAACAGTATCAACAGCAGTGCCAGGCCCATCACACCACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCACCAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATAATGCAGGt                                                                               >  2:57388/1‑150 (MQ=255)
                                                                 tCAACAGCAGTGCCAGGCCCATCACACCACCGGGGATCGGCCAGCCCAGCCAACTGGCGATCTGGCCACCAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATAATGCAGGTTTCATGGCgg                                                                     <  2:176200/150‑1 (MQ=255)
                                                                  cAACAGCAGTGCCAGGCCCATCACACCACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCAACAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATAATGCAGGTTTCATGGCggg                                                                    <  2:295473/150‑1 (MQ=255)
                                                                  cAACAGCAGTGCCAGGCCCATCACACCACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCAACAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATAATGCAGGTTTCATGGCggg                                                                    <  2:295509/150‑1 (MQ=255)
                                                                   aaCAGCAGTGCCAGGCCCATCACACCACCGGGGATCGGCCAGCCCAGCCAACGTGCAATCTGGCCACCAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATAATGCAGGTTTCATGGCgggg                                                                   <  1:442247/150‑1 (MQ=255)
                                                                              cAGGCCCATCACACCACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCACCAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATAATGCAGGTTTCATGGCGGGGTTCCTCGACAg                                                        >  1:583246/1‑150 (MQ=255)
                                                                              cAGGCCCATCACACCACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCACCAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATAATGCAGGTTTCATGGCGGGATTCCTCGACAg                                                        <  2:104286/150‑1 (MQ=255)
                                                                                       cacaCCACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCACCAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATAATGCAGGTTTCATGGCGGGGTTCCTCGACAGGCCCTCAtt                                               <  2:557840/150‑1 (MQ=255)
                                                                                                                 gccAACTGGCAATCTGGCCACCAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATAATGCAGGTTTCATGGCGGGGTTCCTCGACAGGCCCTCATTTTAAGATTTGGCTGCCTAGGCCAGaa                     <  1:57388/150‑1 (MQ=255)
                                                                                                                   cAACTGGCAATCTGGCCACCAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATAATGCAGGTTTCATGGCGGGGTTCCTCGACAGGCCCTCATTTTAAGATTTGGCTGCCTAGGc                          <  2:211673/143‑1 (MQ=255)
                                                                                                                   cAACTGGCAATCTGGCCACCAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATAATGCAGGTTTCATGGCGGGGTTCCTCGACAGGCCCTCATTTTAAGATTTGGCTGCCTAGGc                          >  1:211673/1‑143 (MQ=255)
                                                                                                                       tGGCAATCTGGCCACCAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATACTGCAGGTTTCATGGCGGGGTTCCTCGACAGGCCCTCATTTTAAGATTTGGCTGCCTAGGCCAGAAGCGAAt               <  1:53796/150‑1 (MQ=255)
                                                                                                                                  ccaacaAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTCTTTTCAATCATGCAGGTTTCATGGCGGGGTTCCTCGACAGGCCCTCATTTTAAGATTTGGCTGCCTAGGCCAGAAGCGAATTGTTAGACTGc    <  2:201416/150‑1 (MQ=255)
                                                                                                                                    accaAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATAATGCAGGTTTCATGGCGGGGTTCCTCGACAGGCCCTCATTTTAAGATTTGGCTGCCTAGGCCAGAAGCGAATTGTTAGACTGCTc  <  2:56142/150‑1 (MQ=255)
                                                                                                                                    accaAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGAtttttt                                                                                            <  2:27270/60‑1 (MQ=255)
                                                                                                                                    accaAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGAtttttt                                                                                            >  1:27270/1‑60 (MQ=255)
                                                                                                                                     |                                                                                                                                                    
AGCCATCAGCCAGCCGGCCCCCAATTGCAGCATGGCCGGTTTGAGCACACCCGAAGCGAACAGTATCAACAGCAGTGCCAGGCCCATCACACCACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCACCAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATAATGCAGGTTTCATGGCGGGGTTCCTCGACAGGCCCTCATTTTAAGATTTGGCTGCCTAGGCCAGAAGCGAATTGTTAGACTGCTC  >  NC_002947/5250899‑5251180

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