Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,545,072 T→C 10.8% D344G (GAC→GGC)  ampG ← muropeptide permease AmpG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,545,0720TC10.8% 93.5 / 3.4 37D344G (GAC→GGC) ampGmuropeptide permease AmpG
Reads supporting (aligned to +/- strand):  ref base T (17/16);  new base C (2/2);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

CCGGCAGAGGCTACGCCGCCGATGAACAGGATGGGCATGATGCCGAAGCGTACGATCAGCAGGCCGCCAACCCCGGCGCCGACCAGGGTCATGACCAGGCCGAAGATCTTGCTGACGCTGGCAATCTGGTCCTTGGTGAAGCCCATGTCGATATAGAACACGTTGGCCATCACCCCCATCACCGTGTCCGACATGCGATAGGTGGCGATCAGGCCCAGCAGCAGCAGAGCCTGCCAGCGGTAGCGGGTGATGAAGTCGTTGATCGGGGTCAGCACGGGGGCCAGGCCGCGA  >  NC_002947/1544925‑1545215
                                                                                                                                                   |                                                                                                                                               
ccGGCAGAGGCTACGCCGCCGATGAACAGGATGGGCATGATGCCGAAGCGTACGATCAGCAGGCCGCCAACCCCGGCGCCGACCAGGGTCATGACCAGGCCGAAGATCTTGCTGACGCTGGCAATCTGGTCCTTGGTGAAGCCCATGTc                                                                                                                                                >  6:100591/1‑149 (MQ=255)
         gCTACGCCGCCGATGAACAGGATGGGCATGATGCCGAAGCGTACGATCAGCAGGCCGCCAACCCCGGCGCCGACCAGGGTCATGACCAGGCCGAAGATCTTGCTGACGCTGGCAATCTGGTCCTTGGTGAAGCCCATGTCGATATAGaa                                                                                                                                       <  3:188133/149‑1 (MQ=255)
             cgccgcCGATGAACAGGATGGGCATGATGCCGAAGCGTACGATCAGCAGGCCGCCAACCCCGGCGCCGACCAGGGTCATGACCAGGCCGAAGATCTTGCTGACGCTGGCAATCTGGTCCTTGGTGAAGCCCATGTCGATATAGAaca                                                                                                                                     <  5:100591/147‑1 (MQ=255)
                                    atgatgCCGAAGCGTACGATCAGCAGGCCGCCAACCCCGGCGCCGACCAGGGTCATGACCAGGCCGAAGATCTTGCTGACGCTGGCAATCTGGTCCTTGGTGAAGCCCATGTCGATATAGTACACGTTGGCCATCACCCCCATCACCgt                                                                                                            <  1:188481/149‑1 (MQ=255)
                                      gatgCCGAAGCGTACGATCAGCAGGCCGCCAACCCCGGCGCCGACCAGGGTCATGACCAGGCCGAAGATCTTGCTGACGCTGGCAATCTGGTCCTTGGTGAAGCCCATGTCGATATAGAACACGTTGGCCATCACCCCCATCACCgtgt                                                                                                          >  5:156065/1‑149 (MQ=255)
                                        tgCCTAAGCGTACGATCAGCAGGCCGCCAACCCCGGCGCCGACCAGTGTCATGACCAGGCCTAAGATCTTGCTGACGCTGGCAATCTGTTCCTTGGTGAAGCCCATGTCGatat                                                                                                                                           <  2:199573/114‑1 (MQ=255)
                                        tgCCGAAGCGTACGATCAGCAGGCCGCAAACCCCGGCGCCGACCAGGGTCATGACCAGGCCGAAGATCTTGCTGACGCTGGCAATCTGGTCCTTGGTGAAGCCCATGTCGatat                                                                                                                                           >  1:199573/1‑114 (MQ=255)
                                              aGCGTACGATCAGCAGGCCGCCAACCCCGGCGCCGACCAGGGTCATGACCAGGCCGAAGATCTTGCTGACGCTGGCAATCTGGTCCTTGGTGAAGCCCATGTCGATATAGAACACGTTGGCCATCACCCCCATCACCGTGTCCGACATg                                                                                                  <  3:104754/149‑1 (MQ=255)
                                               gCGTACGATCAGCAGGCCGCCAACCCCGGCGCCGACCAGGGTCATGACCAGGCCGAAGATCTTGCTGACGCTGGCAATCTGGTCCTTGGTGAAGCCCATGTCGATATAGAACACGTTGGCCATCACCCCCATCACCGTGTCCGACATGc                                                                                                 >  2:315875/1‑149 (MQ=255)
                                                 gTACGATCAGCAGGCCGCCAACCCCGGCGCCGACCAGGGTCATGACCAGGCCGAAGATCTTGCTGACGCTGGCAATCTGGTCCTTGGTGAAGCCCATGTCGATATAGAACACGTTGGCCATCACCCCCATCACCGTGTCCGACATGCGa                                                                                               >  5:290074/1‑149 (MQ=255)
                                                     gATCAGCAGGCCGCCAACCCCGGCGCCGACCAGGGTCATGACCAGGCCGAAGATCTTGCTGACGCTGGCAATCTGGTCCTTGGTGAAGCCCATGTCGATATAGAACACGTTGGCCATCACCCCCATCACCGTGTCCGACATGCGATAgg                                                                                           >  6:252397/1‑149 (MQ=255)
                                                          gcagGCCGCCAACCCCGGCGCCGACCAGGGTCATGACCAGGCCGAAGATCTTGCTGACGCTGGCAATCTGGTCCTTGGTGAAGCCCATGTCGATATAGAACACGTTGGCCATCACCCCCATCACCgtg                                                                                                           <  1:322528/128‑1 (MQ=255)
                                                          gcagGCCGCCAACCCCGGCGCCGACCAGGGTCATGACCAGGCCGAAGATCTTGCTGACGCTGGCAATCTGGTCCTTGGTGAAGCCCATGTCGATATAGAACACGTTGGCCATCACCCCCATCACCgtg                                                                                                           >  2:322528/1‑128 (MQ=255)
                                                            agGCCGCCAACCCCGGCGCCGACCAGGGTCATGACCAGGCCGAAGATCTTGCTGACGCTGGCAATCTGGTCCTTGGTGAAGCCCATGCCGATATAGAACACGTTGGCCATCACCCCCATCACCGTGTCCGACATGCGATAGGTGGCGAt                                                                                    >  4:24511/1‑149 (MQ=255)
                                                            agGCCGCCAACCCCGGCGCCGACCAGGGTCATGACCAGGCAGAAGATCTTGCTGACGCTGGCAATCTGGTCCTAGGTGAAGCCCATGCAGATAAAGAACACGATGGCCATCACCCCCATCACCGATACCGACATGAGATAGGTg                                                                                         >  4:218583/1‑144 (MQ=255)
                                                              gccgccAACCCCGGCGCCGACCAGGGTCATGACCAGGCCGAAGATCTTGCTGACGCTGGCAATCTGGTCCTTGGTGAAGCCCATGTCGATATAGAACACGTTGGCCATCACCCCCATCACCGTGTCCGACATGCGATAGGTGGCGATCa                                                                                  <  8:285258/149‑1 (MQ=255)
                                                                 gccAACCCCGGCGCCGACCAGGGTCATGACCAGGCCGAAGATCTTGCTGACGCTGGCAATCTGGTCCTTGGTGAAGCCCATGTCGATATAGAACACGTTGGCCATCACCCCCATc                                                                                                                 >  7:350/1‑115 (MQ=255)
                                                                 gccAACCCCGGCGCCGACCAGGGTCATGACCAGGCCGAAGATCTTGCTGACGCTGGCAATCTGGTCCTTGGTGAAGCCCATGTCGATATAGAACACGTTGGCCATCACCCCCATc                                                                                                                 <  8:350/115‑1 (MQ=255)
                                                                           gcgcCGACCAGGGTCATGACCAGGCCGAAGATCTTGCTGACGCTGGCAATCTGGTCCTTGGTGAAGCCCATGTCGATATAGAACACGTTGGCCATCACCCCCATCACCGTGTCCGACATGCGATAGGTGGCGATCAGGCCcagcagcag                                                                     <  1:104748/149‑1 (MQ=255)
                                                                             gcCGACCAGGGTCATGACCAGGCCGAAGATCTTGCTGACGCTGGCAATCTGGTCCTTGGTGAAGCCCATGTCGATATAGAACACGTTGGCCATCACCCCCATCACCGTGTCCGACATGCGATAGGTGGCGATCAGGCCcagcagcag                                                                     >  4:117450/1‑147 (MQ=255)
                                                                               cGACCAGGGTCATGACCAGGCCGAAGATCTTGCTGACGCTGGCAATCTGGTCCTTGGTGAAGCCCATGCCTATATAGAACACGTTGGCCATCACCCCCATCACCGTGTCCGACATGCGATAGGTGGCGATCAGGCCCAGCAGCAGCaga                                                                 <  3:218583/149‑1 (MQ=255)
                                                                               cGACCAGGGTCATGACCAGGCCGAAGATCTTGCTGACGCTGGCAATCTGGTCCTTGGTGAAGCCCATGCCGATATAGAACACGTTGGCCATCACCCCCATCACCGTGTCCGACATGCGATAGGTGGCGATCAGGCCCAGCAGCAGCaga                                                                 <  3:24511/149‑1 (MQ=255)
                                                                                    aGGGTCATGACCAGGCCGAAGATCTTGCTGACGCTGGCAATCTGGTCCTTGGTGAAGCCCATGTCGATATAGAACACGTTGGCCATCACCCCCATCACCGTGTCCGACATGCGATAGGTGGCGATCAGGCCCAGCAGCAGCAGAGCCTg                                                            <  7:184525/149‑1 (MQ=255)
                                                                                     gggTCATGACCAGGCCGAAGATCTTGCTGACGCTGGCAATCTGGTCCTTGGTGAAGCCCATGTCGATATAGAACACGTTGGCCATCACCCCCATCACCGTGTCCGACATCCGATAGGTGGCGATCAGGCCCAGCAGCAGCAGAGCCTg                                                            >  1:47134/1‑148 (MQ=255)
                                                                                             aCCAGGCCGAAGATCTTGCTGACGCTGGCAATCTGGTCCTTGGTGAAGCCCATGTCGATATAGAACACGTTGGCCATCACCCCCATCACCGTGTCCGACATGCGATAGGTGGCGATCAGGCCCAGCAGCAGCAGAGCCTGCCAGCGGt                                                    <  1:315875/148‑1 (MQ=255)
                                                                                               cAGGCCGAAGATCTTGCTGACGCTGGCAATCTGGTCCTTGGTGAAGCCCATGTCGATATAGAACACGTTGGCCATCACCCCCATCACCGTGTCCGACATGCGATAGGTGGCGATCAGGCCCAGCAGCAGCAGAGCCTGCCAGCGGTAGc                                                 <  6:156065/149‑1 (MQ=255)
                                                                                                                  aCGCTGGCAATCTGGTCCTTGGTGAAGCCCATGTCGATATAGAACACGTTGGCCATCACCCCCCTCACCGTGTCCGACATGCGATAGGTGGCGATCAGGCCCAGCAGCAGCAGAGCCTGCCAGCGGTAGCGGGtg                                            >  7:211622/1‑135 (MQ=255)
                                                                                                                  aCGCTGGCAATCTGGTCCTTGGTGAAGCCCATGTCGATATAGAACACGTTGGCCATCACCCCCATCACCGTGTCCGACATGCGATAGGTGGCGATCAGGCCCAGCAGCAGCAGAGCCTGCCAGCGGTAGCGGGtg                                            <  8:211622/135‑1 (MQ=255)
                                                                                                                    gCTGGCAATCTGGTCCTTGGTGAAGCCCATGTCGATATAGAACACGTTGGCCATCACCCCCATCACCGTGTCCGACATGCGATAGGTGGCGATCAgg                                                                                <  4:123310/97‑1 (MQ=255)
                                                                                                                    gCTGGCAATCTGGTCCTTGGTGAAGCCCATGTCGATATAGAACACGTTGGCCATCACCCCCATCACCGTGTCCGACATGCGATAGGTGGCGATCAgg                                                                                >  3:123310/1‑97 (MQ=255)
                                                                                                                       ggCAATCTGGTCCTTGGTGAAGCCCATGTCGATATAGAACACGTTGGCCATCACCCCCATCACCGTCTCCGACATGCGATAGGTGGCGATCAGGCCCAGCAGCAGCAGAGCCTGCCAGCGGTAGCGGGTGATGAAGTCGTAGAGCgggg                         >  1:302719/1‑149 (MQ=255)
                                                                                                                        gCAATCTGGTCCTTGGTGAAGCCCATGTCGATATAGAACACGTTGGCCATCACCCCCATCACCGTGTCCGACATGCGATAGGTGGCGATCAGGCCCAGCAGCAGCAGAGCCTGCCAGCGGTAGCGGGTGATGAAGTCGTTGATCggg                          <  5:252397/147‑1 (MQ=255)
                                                                                                                                    ttGGTGAAGCCCATGTCGATATAGAACACGTTGGCCATCACCCCCATCACCGTGTCCGACATGCGATAGGTGGCGATCAGGCCCAGCAGCAGCAGAGCCTGCCAGCGGTAGCGGGTGATGAAGTCGTTGATCGGGGTCAGCACGGGGGc            >  6:204850/1‑149 (MQ=255)
                                                                                                                                          aaGCCCATGTCGATATAGAACACGTTGGCCATCACCCCCATCACCGTGTCCGACATGCGATAGGTGGCGATCAGGCCCAGCAGCAGCAGAGCCTGCCAGCGGTAGCGGGTGATGAAGTCGTTGATCGGGGTCAGCACGGGGGCCAGGcc      >  1:35202/1‑149 (MQ=255)
                                                                                                                                            gCCCATGTCGATATAGAACACGTTGGCCATCACCCCCATCACCGTGTCCGACATGCGATAGGTGGCGATCAGGCCCAGCAGCAGCAGAGCCTGCCAGCGGTAGCGGGtga                                           >  2:164554/1‑110 (MQ=255)
                                                                                                                                            gCCCATGTCGATATAGAACACGTTGGCCATCACCCCCATCACCGTGTCCGACATGCGATAGGTGGCGATCAGGCCCAGCAGCAGCAGAGCCTGCCAGCGGTAGCGGGtga                                           <  1:164554/110‑1 (MQ=255)
                                                                                                                                              ccATGTCGATATAGAACACGTTGGCCATCACCCCCATCACCGTGTCCGACATGCCATACGTGGCGATCAGGCCCAGCACAACCAGAGCCTGCCAGCGGAAGCGGGTCATGAAGTCGTTGATCGGGGTCAGCACGGGGGCCAGGCCGCGa  >  5:206261/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CCGGCAGAGGCTACGCCGCCGATGAACAGGATGGGCATGATGCCGAAGCGTACGATCAGCAGGCCGCCAACCCCGGCGCCGACCAGGGTCATGACCAGGCCGAAGATCTTGCTGACGCTGGCAATCTGGTCCTTGGTGAAGCCCATGTCGATATAGAACACGTTGGCCATCACCCCCATCACCGTGTCCGACATGCGATAGGTGGCGATCAGGCCCAGCAGCAGCAGAGCCTGCCAGCGGTAGCGGGTGATGAAGTCGTTGATCGGGGTCAGCACGGGGGCCAGGCCGCGA  >  NC_002947/1544925‑1545215

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