Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 217,959 A→G 23.9% intergenic (‑51/‑55) yafU ← / → yafF predicted inner membrane protein/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE217,9590AG23.9% 79.7 / 24.4 46intergenic (‑51/‑55)yafU/yafFpredicted inner membrane protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (20/15);  new base G (7/4);  total (27/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

AGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAAT  >  minE/217889‑218019
                                                                      |                                                            
aGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCa                                                              <  1:3579398/71‑1 (MQ=255)
aGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCa                                                              <  1:4895552/71‑1 (MQ=255)
      cGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCg                                                         >  1:4188791/1‑70 (MQ=255)
       gCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATcc                                                          <  1:109112/68‑1 (MQ=255)
       gCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGa                                                        <  1:4997074/70‑1 (MQ=255)
           tGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCa                                                   >  1:3517786/1‑71 (MQ=255)
           tGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCa                                                   >  1:3764827/1‑71 (MQ=255)
             aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                                  <  1:3644044/70‑1 (MQ=255)
             aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                                  <  1:4325023/70‑1 (MQ=255)
             aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                                  <  1:2679150/70‑1 (MQ=255)
             aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                                  <  1:3821054/70‑1 (MQ=255)
             aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACt                                                 >  1:837663/1‑71 (MQ=255)
             aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACt                                                 >  1:1708473/1‑71 (MQ=255)
             aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACt                                                 >  1:1639437/1‑71 (MQ=255)
             aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACt                                                 >  1:4511118/1‑71 (MQ=255)
             aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACt                                                 >  1:1398361/1‑71 (MQ=255)
             aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACt                                                 >  1:259609/1‑71 (MQ=255)
             aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTAGCCACAGCAATCCGAAACCACt                                                 >  1:2279135/1‑71 (MQ=255)
              aCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTg                                                >  1:919237/1‑71 (MQ=255)
              aCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTg                                                >  1:1790189/1‑71 (MQ=255)
              aCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTg                                                >  1:2108639/1‑71 (MQ=255)
                   taataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                                  <  1:3882919/64‑1 (MQ=255)
                   taataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACg                                           <  1:1283294/71‑1 (MQ=255)
                   aaataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACg                                           <  1:127667/70‑1 (MQ=255)
                    aataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAAc                                                     <  1:4770672/60‑1 (MQ=255)
                    aataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAAc                                                     <  1:2468057/60‑1 (MQ=255)
                     ataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAAc                                                     >  1:2325212/1‑59 (MQ=255)
                     ataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTg                                         >  1:4016527/1‑71 (MQ=255)
                               tGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATa                                  <  1:3891386/68‑1 (MQ=255)
                                ggTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCt                              >  1:197515/1‑71 (MQ=255)
                                       ctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAgg                       >  1:2937070/1‑71 (MQ=255)
                                        tgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGc                      >  1:2751679/1‑71 (MQ=255)
                                        tgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGc                      >  1:3919022/1‑71 (MQ=255)
                                          aTTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGt                    <  1:3954901/71‑1 (MQ=255)
                                                 cGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACg             <  1:2151764/71‑1 (MQ=255)
                                                           ttCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaa   >  1:3499276/1‑71 (MQ=255)
                                                           ttCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaa   >  1:975920/1‑71 (MQ=255)
                                                           ttCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaa   >  1:907351/1‑71 (MQ=255)
                                                           ttCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaa   >  1:4029888/1‑71 (MQ=255)
                                                           ttCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaa   >  1:3326733/1‑71 (MQ=255)
                                                           ttCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGGaatgaa   >  1:3503472/1‑71 (MQ=255)
                                                             cGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaatgaat  >  1:4132721/1‑70 (MQ=255)
                                                               ccACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaatga    <  1:3809680/66‑1 (MQ=255)
                                                                  caGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGc                      >  1:1450720/1‑45 (MQ=255)
                                                                   aGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTg                 <  1:43725/49‑1 (MQ=255)
                                                                   aGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTg                 <  1:2109642/49‑1 (MQ=255)
                                                                   aGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTg                 <  1:2076933/49‑1 (MQ=255)
                                                                   aGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTg                 <  1:3250859/49‑1 (MQ=255)
                                                                      |                                                            
AGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAAT  >  minE/217889‑218019

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: