Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,530,398 G→A 11.1% A26T (GCA→ACA)  mutM → formamidopyrimidine/5‑formyluracil/ 5‑hydroxymethyluracil DNA glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,530,3980GA11.1% 98.2 / 4.3 36A26T (GCA→ACA) mutMformamidopyrimidine/5‑formyluracil/ 5‑hydroxymethyluracil DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base G (23/9);  new base A (2/2);  total (25/11)
Fisher's exact test for biased strand distribution p-value = 5.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

ACCCGAAGTTGAAACCAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAGCGACC  >  minE/2530334‑2530461
                                                                |                                                               
aCCCGAAGTTGAAACCAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAgtg                                                            <  1:1715374/70‑1 (MQ=255)
 cccGAAGTTGAAACCAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATg                                                                 >  1:2180103/1‑64 (MQ=255)
    gAAGTTGAAACCAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTgcg                                                       >  1:835482/1‑71 (MQ=255)
       gTTGAAACCAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCaa                                                    >  1:388565/1‑71 (MQ=255)
       gTTGAAACCAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCaa                                                    >  1:2007260/1‑71 (MQ=255)
        ttGAAACCAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATg                                                                 <  1:1430973/57‑1 (MQ=255)
          gAAACCAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAAc                                                   >  1:1652547/1‑69 (MQ=255)
              ccAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGc                                             >  1:1068827/1‑71 (MQ=255)
                aGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAgtggtg                                                         >  1:171172/1‑57 (MQ=255)
                aGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCtt                                           >  1:1144615/1‑71 (MQ=255)
                  ccgccgCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTgc                                         <  1:1975675/71‑1 (MQ=255)
                   cgccgcGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCtt                                           >  1:687685/1‑68 (MQ=255)
                   cgccgcGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTgcg                                        <  1:230877/71‑1 (MQ=255)
                     ccgcGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATg                                                                 <  1:2011383/44‑1 (MQ=255)
                     ccgcGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATg                                                                 <  1:498137/44‑1 (MQ=255)
                     ccgcGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATg                                                                 <  1:82248/44‑1 (MQ=255)
                          gCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTgcgc                                       <  1:1015314/65‑1 (MQ=255)
                               gAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAgtgg                                                           >  1:1946299/1‑40 (MQ=255)
                                 aCCGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCag                          <  1:1078744/71‑1 (MQ=255)
                                  ccGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCag                          >  1:1949386/1‑70 (MQ=255)
                                   cGCATCTCGTTGGTGCAACCATTCTTCATACAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCagaa                        <  1:1997424/71‑1 (MQ=255)
                                     cATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCagaaga                      <  1:1759835/71‑1 (MQ=255)
                                     cATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCagaaga                      <  1:544075/71‑1 (MQ=255)
                                         tcGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCagaa                        >  1:2050360/1‑65 (MQ=255)
                                         tcGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCagaa                        >  1:400398/1‑65 (MQ=255)
                                           gTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTg                                     <  1:495068/50‑1 (MQ=255)
                                              ggTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCg             >  1:394937/1‑71 (MQ=255)
                                                 gCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGtt           >  1:2107962/1‑70 (MQ=255)
                                                 gCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGtt           >  1:356215/1‑70 (MQ=255)
                                                    aCCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAg       >  1:1901802/1‑71 (MQ=255)
                                                     ccATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAGc      >  1:541717/1‑71 (MQ=255)
                                                     ccATTCTTCATACAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAGc      <  1:1407520/71‑1 (MQ=255)
                                                      cATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAgaga                    >  1:1169710/1‑56 (MQ=255)
                                                      cATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAgaga                    >  1:1702270/1‑56 (MQ=255)
                                                      cATTCTTCATACAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAGc      >  1:870683/1‑70 (MQ=255)
                                                      cATTCTTCATACAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAGc      >  1:2172257/1‑70 (MQ=255)
                                                        ttcttcATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAagag                     <  1:1777135/53‑1 (MQ=255)
                                                        ttcttcATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAGCGac   >  1:329472/1‑71 (MQ=255)
                                                        ttcttcATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAGCGac   >  1:1811304/1‑71 (MQ=255)
                                                          cttcATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCagaa                        >  1:1553621/1‑48 (MQ=255)
                                                          cttcATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAGCGacc  >  1:444454/1‑70 (MQ=255)
                                                                |                                                               
ACCCGAAGTTGAAACCAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAGCGACC  >  minE/2530334‑2530461

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