Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,697,440 C→T 35.5% P262P (CCG→CCA mrcA ← fused penicillin‑binding protein 1a murein transglycosylase and murein transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,697,4400CT35.5% 29.0 / 23.0 45P262P (CCG→CCAmrcAfused penicillin‑binding protein 1a murein transglycosylase and murein transpeptidase
Reads supporting (aligned to +/- strand):  ref base C (20/9);  new base T (6/10);  total (26/19)
Fisher's exact test for biased strand distribution p-value = 6.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.42e-01

ACCGTCTTCATAGGCACTTTCGCCATAACGGTTATACATCTCCTGGCGCACCATTTCGCTCAGGTACGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGTTAATCGCCTCAGTGCGTGTCTGATCG  >  minE/2697374‑2697506
                                                                  |                                                                  
aCCGTCTTCATAGGCACTTTCGCCATAACGGTTATACATCTCCTGGCGCACCATTTCGCTCAGGTACgg                                                                  >  1:821746/1‑69 (MQ=255)
              cACTTTCGCCATAACGGTTATACATCTCCTGGCGCACCATTTCGCTCAGGTACGGCGCAGAGAAAGCa                                                     <  1:2336871/68‑1 (MQ=255)
               aCTTTCGCCATAACGGTTATACATCTCCTGGCGCACCATTTCGCTCAGGTACGGCGCATAGAAAGCAAt                                                   >  1:1131414/1‑69 (MQ=255)
               aCTTTCGCCATAACGGTTATACATCTCCTGGCGCACCATTTCGCTCAGGTACGGCGCAGAGAAAGCAAt                                                   >  1:3024178/1‑69 (MQ=255)
               aCTTTCGCCATAACGGTTATACATCTCCTGGCGCACCATTTCGCTCAGGTACGGCGCAGAGAAAGCAAt                                                   >  1:2310543/1‑69 (MQ=255)
                          aaCGGTTATACATCTCCTGGCGCACCATTTCGCTCAGGTACGGCGCAGAGAAAGCAATCTCCGGCGCGt                                        <  1:1164663/69‑1 (MQ=255)
                          aaCGGTTATACATCTCCTGGCGCACCATTTCGCTCAGGTACGGCGCAGAGAAAGCAATCTCCGGCGCGt                                        <  1:2051465/69‑1 (MQ=255)
                           aCGGTTATACATCTCCTGGCGCACCATTTCGCTCAGGTACGGCGCag                                                             >  1:2971202/1‑47 (MQ=255)
                           aCGGTTATACATCTCCTGGCGCACCATTTCGCTCAGGTACGGCGCAGAGAAAGCa                                                     >  1:1043457/1‑55 (MQ=255)
                           aCGGTTATACATCTCCTGGCGCACCATTTCGCTCAGGTACGGCGCAGAGAAAGCAATCTCCGGCGCGt                                        >  1:2794529/1‑68 (MQ=255)
                           aCGGTTATACATCTCCTGGCGCACCATTTCGCTCAGGTACGGCGCAGAGAAAGCAATCTCCGGCGCGt                                        >  1:75966/1‑68 (MQ=255)
                           aCGGTTATACATCTCCTGGCGCACCATTTCGCTCAGGTACGGCGCAGAGAAAGCAATCTCCGGCGCGt                                        >  1:3409289/1‑68 (MQ=255)
                           aCGGTTATACATCTCCTGGCGCACCATTTCGCTCAGGTACGGCGCAGAGAAAGCAATCTCCGGCGCGt                                        >  1:2395817/1‑68 (MQ=255)
                           aCGGTTATACATCTCCTGGCGCACCATTTCGCTCAGGTACGGCGCAGAGAAAGCAATCTCCGGCGCGTg                                       <  1:1056502/69‑1 (MQ=255)
                           aCGGTTATACATCTCCTGGCGCACCATTTCGCTCAGGTACGGCGCAGAGAAAGCAATCTCCGGCGCGTg                                       <  1:2406834/69‑1 (MQ=255)
                           aCGGTTATACATCTCCTGGCGCACCATTTCGCTCAGGTACGGCGCAGAGAAAGCAATCTCCGGCGCGTg                                       <  1:2402660/69‑1 (MQ=255)
                           aCGGTTATACATCTCCTGGCGCACCATTTCGCTCAGGTACGGCGCAGAGAAAGCAATCTCCGGCGCGTg                                       <  1:2891539/69‑1 (MQ=255)
                           aCGGTTATACATCTCCTGGCGCACCATTTCGCTCAGGTACGGCGCAGAGAAAGCAATCTCCGGCGCGTg                                       <  1:1955877/69‑1 (MQ=255)
                           aCGGTTATACATCTCATGGCGCACAATTTCGCTCAGGTACGGCGCAGAGAAAGCaa                                                    >  1:1727320/1‑56 (MQ=255)
                              gTTATACATCTCCTGGCGCACCATTTCGCTCAGGTACGGCGCAGAGAAAGCAATCTCCGGCGCGTGATa                                    <  1:2634306/69‑1 (MQ=255)
                                 ataCATCTCCTGGCGCACCATTTCACTCAGGTATGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGtt                                 >  1:926125/1‑69 (MQ=255)
                                 ataCATCTCCTGGCGCACCATTTCACTCAGGTATGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGtt                                 >  1:991127/1‑69 (MQ=255)
                                 ataCATCTCCTGGCGCACCATTTCACTCAGGTATGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGtt                                 >  1:2739004/1‑69 (MQ=255)
                                 ataCATCTCCTGGCGCACCATTTCACTCAGGTATGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGtt                                 >  1:1243432/1‑69 (MQ=255)
                                 ataCATCTCCTGGCGCACCATTTCACTCAGGTATGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGtt                                 >  1:2413633/1‑69 (MQ=255)
                                 ataCATCTCCTGGCGCACCATTTCACTCAGGTATGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGtt                                 >  1:2122312/1‑69 (MQ=255)
                                            ggCGCACCATTTCACTCAGGTATGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGGTAATc                       <  1:247004/68‑1 (MQ=255)
                                            ggCGCACCATTTCACTCAGGTATGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGGTAATc                       <  1:1489378/68‑1 (MQ=255)
                                            ggCGCACCATTTCACTCAGGTATGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGGTAATc                       <  1:1747853/68‑1 (MQ=255)
                                            ggCGCACCATTTCACTCAGGTATGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGGTAATc                       <  1:2393485/68‑1 (MQ=255)
                                            ggCGCACCATTTCACTCAGGTATGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGGTAATc                       <  1:2923474/68‑1 (MQ=255)
                                            ggCGCACCATTTCACTCAGGTATGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGGTAATc                       <  1:2140013/68‑1 (MQ=255)
                                            ggCGCACCATTTCACTCAGGTATGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGGTAATc                       <  1:3015040/68‑1 (MQ=255)
                                            ggCGCACCATTTCACTCAGGTATGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGGTAATc                       <  1:1793823/68‑1 (MQ=255)
                                            ggCGCACCATTTCACTCAGGTATGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGGTAATc                       <  1:3235904/68‑1 (MQ=255)
                                            ggCGCACCATTTCACTCAGGTATGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGGTAATc                       <  1:1825824/68‑1 (MQ=255)
                                                            cAGGTACGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGTTAATCGCCTCAGTGCGTGTCTg      >  1:2521129/1‑69 (MQ=255)
                                                            cAGGTACGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGTTAATCGCCTCAGTGCGTGTCTg      >  1:532265/1‑69 (MQ=255)
                                                            cAGGTACGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGTTAATCGCCTCAGTGCGTGTCTg      >  1:2762204/1‑69 (MQ=255)
                                                            cAGGTACGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGTTAATCGCCTCAGTGCGTGTCTg      >  1:2476704/1‑69 (MQ=255)
                                                                tACGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGTTAATCGCCTCa                 >  1:498299/1‑54 (MQ=255)
                                                                tACGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGTTAATCGCCTCAGTGCGTGTCTGATCg  >  1:3181443/1‑69 (MQ=255)
                                                                tACGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGTTAATCGCCTCAGTGCGTGTCTGATCg  >  1:310340/1‑69 (MQ=255)
                                                                tACGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGTTAATCGCCTCAGTGCGTGTCTGATCg  >  1:848410/1‑69 (MQ=255)
                                                                tACGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGTTAATCGCCTCAGTGCGTGTCTGATCg  >  1:2426338/1‑69 (MQ=255)
                                                                  |                                                                  
ACCGTCTTCATAGGCACTTTCGCCATAACGGTTATACATCTCCTGGCGCACCATTTCGCTCAGGTACGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGTTAATCGCCTCAGTGCGTGTCTGATCG  >  minE/2697374‑2697506

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