Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,697,431 G→A 43.2% S265S (AGC→AGT 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,4310GA43.2% 7.8 / 25.1 37S265S (AGC→AGTmrcAfused penicillin‑binding protein 1a murein transglycosylase and murein transpeptidase
Reads supporting (aligned to +/- strand):  ref base G (12/9);  new base A (6/10);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 3.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

ACCGTCTTCATAGGCACTTTCGCCATAACGGTTATACATCTCCTGGCGCACCATTTCGCTCAGGTACGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGTTAATC  >  minE/2697374‑2697485
                                                         |                                                      
aCCGTCTTCATAGGCACTTTCGCCATAACGGTTATACATCTCCTGGCGCACCATTTCGCTCAGGTACgg                                             >  1:821746/1‑69 (MQ=255)
              cACTTTCGCCATAACGGTTATACATCTCCTGGCGCACCATTTCGCTCAGGTACGGCGCAGAGAAAGCa                                <  1:2336871/68‑1 (MQ=255)
               aCTTTCGCCATAACGGTTATACATCTCCTGGCGCACCATTTCGCTCAGGt                                                 >  1:2356158/1‑50 (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:2395817/1‑68 (MQ=255)
                           aCGGTTATACATCTCCTGGCGCACCATTTCGCTCAGGTACGGCGCAGAGAAAGCAATCTCCGGCGCGt                   >  1:2794529/1‑68 (MQ=255)
                           aCGGTTATACATCTCCTGGCGCACCATTTCGCTCAGGTACGGCGCAGAGAAAGCAATCTCCGGCGCGt                   >  1:3409289/1‑68 (MQ=255)
                           aCGGTTATACATCTCCTGGCGCACCATTTCGCTCAGGTACGGCGCAGAGAAAGCAATCTCCGGCGCGt                   >  1:75966/1‑68 (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:1056502/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:991127/1‑69 (MQ=255)
                                 ataCATCTCCTGGCGCACCATTTCACTCAGGTATGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGtt            >  1:2122312/1‑69 (MQ=255)
                                 ataCATCTCCTGGCGCACCATTTCACTCAGGTATGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGtt            >  1:2413633/1‑69 (MQ=255)
                                 ataCATCTCCTGGCGCACCATTTCACTCAGGTATGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGtt            >  1:926125/1‑69 (MQ=255)
                                 ataCATCTCCTGGCGCACCATTTCACTCAGGTATGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGtt            >  1:2739004/1‑69 (MQ=255)
                                 ataCATCTCCTGGCGCACCATTTCACTCAGGTATGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGtt            >  1:1243432/1‑69 (MQ=255)
                                            ggCGCACCATTTCACTCAGGTATGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGGTAATc  <  1:2923474/68‑1 (MQ=255)
                                            ggCGCACCATTTCACTCAGGTATGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGGTAATc  <  1:1489378/68‑1 (MQ=255)
                                            ggCGCACCATTTCACTCAGGTATGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGGTAATc  <  1:3235904/68‑1 (MQ=255)
                                            ggCGCACCATTTCACTCAGGTATGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGGTAATc  <  1:1747853/68‑1 (MQ=255)
                                            ggCGCACCATTTCACTCAGGTATGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGGTAATc  <  1:3015040/68‑1 (MQ=255)
                                            ggCGCACCATTTCACTCAGGTATGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGGTAATc  <  1:1793823/68‑1 (MQ=255)
                                            ggCGCACCATTTCACTCAGGTATGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGGTAATc  <  1:2393485/68‑1 (MQ=255)
                                            ggCGCACCATTTCACTCAGGTATGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGGTAATc  <  1:1825824/68‑1 (MQ=255)
                                            ggCGCACCATTTCACTCAGGTATGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGGTAATc  <  1:247004/68‑1 (MQ=255)
                                            ggCGCACCATTTCACTCAGGTATGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGGTAATc  <  1:2140013/68‑1 (MQ=255)
                                                         |                                                      
ACCGTCTTCATAGGCACTTTCGCCATAACGGTTATACATCTCCTGGCGCACCATTTCGCTCAGGTACGGCGCAGAGAAAGCAATCTCCGGCGCGTGATAGTTAGCGTTAATC  >  minE/2697374‑2697485

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