Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,622,247 A→C 9.3% intergenic (‑174/+119) ypfE ← / ← maeB predicted carboxysome structural protein with predicted role in ethanol utilization/fused malic enzyme predicted oxidoreductase and predicted phosphotransacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,622,2470AC 9.3% 100.4 / 5.7 43intergenic (‑174/+119)ypfE/maeBpredicted carboxysome structural protein with predicted role in ethanol utilization/fused malic enzyme predicted oxidoreductase and predicted phosphotransacetylase
Reads supporting (aligned to +/- strand):  ref base A (19/20);  new base C (2/2);  total (21/22)
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.95e-01

AATGGAGAGATATTCGCTGTGGTGCATAAACTTTCATATGCAACGCATGAATATAAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCAAACTCAGTAATAAGCTTTACTAATATATTGCAGC  >  minE/1622177‑1622316
                                                                      |                                                                     
aaTGGAGAGATATTCGCTGTGGTGCATAAACTTTCATATGCAACGCATGAATATAAATAAAATAATATCga                                                                       <  1:3435478/71‑1 (MQ=255)
 aTGGAGAGATATTCGCTGTGGTGCATAAACTTTCATATGCAACGCATGAATATAAATAAAATAATATCga                                                                       <  1:2996037/70‑1 (MQ=255)
    gagagaTATTCGCTGTGGTGCATAAACTTTCATATGCAACGCATGAATATAAATAAAATAATATCGAGata                                                                   <  1:1654539/71‑1 (MQ=255)
     agagaTATTCGCTGTGGTGCATAAACTTTCATATGCAACGCATGAATATAAATAAAATAATATCTCGatat                                                                  >  1:434617/1‑71 (MQ=255)
       agaTATTCGCTGTGGTGCATAAACTTTCATATGCAACGCATGAATATAAATAAAATAATATCGAGAtatta                                                                >  1:673467/1‑71 (MQ=255)
         atatTCGCTGTGGTGCATAAACTTTCATATGCAACGCATGAATATAAATAAAATAATATCGAGATATtata                                                              <  1:3727289/71‑1 (MQ=255)
          tatTCGCTGTGGTGCATAAACTTTCATATGCAACGCATGAATATAAATAAAATAATATCGAGATATtata                                                              <  1:1856598/70‑1 (MQ=255)
            ttCGCTGTGGTGCATAAACTTTCATATGCAACGCATGAATATAAATAAAATAATATCGAGATATTATATAg                                                           <  1:2533865/71‑1 (MQ=255)
                 tgtgGTGCATAAACTTTCATATGCAACGCATGAATATAAATAAAATAATATCGAGATATTATATAGACTGc                                                      >  1:2132595/1‑71 (MQ=255)
                       gCATAAACTTTCATATGCAACGCATGAATATAAATAAAATAATATCGAGAtatt                                                                 <  1:3774380/54‑1 (MQ=255)
                             aCTTTCATATGCAACGCATGAATATAAATAAAATAATATCGAGATATtatata                                                            >  1:3799056/1‑53 (MQ=255)
                              cTTTCATATGCAACGCATGAATATAAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTTtattat                                         <  1:31134/71‑1 (MQ=255)
                              cTTTCATATGCAACGCATGAATATAAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTTtattat                                         >  1:817336/1‑71 (MQ=255)
                               tttCATATGCAACGCATGAATATAAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTttatta                                          <  1:1217132/69‑1 (MQ=255)
                                 tCATATGCAACGCATGAATATAAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTTTATTATTc                                       >  1:3300664/1‑70 (MQ=255)
                                   atatGCAACGCATGAATATAAATAAAATAATATCTCGAtatt                                                                 <  1:10914/42‑1 (MQ=38)
                                   atatGCAACGCATGAATATAAATAAAATAATATCTCGAtatt                                                                 <  1:2618428/42‑1 (MQ=38)
                                   atatGCAACGCATGAATATAAATAAAATAATATCGAGATATTATATAGACTGCtttttttt                                              <  1:2699563/61‑1 (MQ=255)
                                       gCAACGCATGAATATAAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCAAAc                                 >  1:1563975/1‑70 (MQ=255)
                                        cAACGCATGAATATAAATAAAATAATATCTCGAtattat                                                               >  1:1633906/1‑39 (MQ=25)
                                         aaCGCATGAATATAAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCAAACt                                >  1:3567042/1‑69 (MQ=255)
                                         aaCGCATGAATATAAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCAAACt                                >  1:3562468/1‑69 (MQ=255)
                                           cGCATGAATATAAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCAAACTCAg                             >  1:2269243/1‑70 (MQ=255)
                                              aTGAATATAAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCAAACTCAGtaa                          >  1:722552/1‑70 (MQ=255)
                                               tGAATATAAATAAAATAATATCGAGATATTATATAGACTGCtttt                                                  >  1:2812495/1‑45 (MQ=255)
                                                 aaTATAAATAAAATAATATCGAGATATTATATAGACTGCttttttt                                               <  1:3358594/46‑1 (MQ=255)
                                                 aaTATAAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCAAACTCAGtaataa                       <  1:2384585/70‑1 (MQ=255)
                                                 aaTATAAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCAAACTCAGTAATAAg                      <  1:779428/71‑1 (MQ=255)
                                                  atatAAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTttatta                                          <  1:564881/50‑1 (MQ=255)
                                                  atatAAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCAAACTCAGTAATAAGc                     <  1:2657812/71‑1 (MQ=255)
                                                  atatAAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCAAACTCAGTAATAAGc                     >  1:1185765/1‑71 (MQ=255)
                                                   tataAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTTTattatt                                        <  1:1125773/51‑1 (MQ=255)
                                                   tataAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCAAACTCAGTAATAAGCt                    >  1:2507465/1‑71 (MQ=255)
                                                   tataAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCAAACTCAGTAATAAGCt                    >  1:2491885/1‑71 (MQ=255)
                                                     taaataaaATAATATCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCAAACTCAGTAATAAGCttt                  <  1:3073492/71‑1 (MQ=255)
                                                        ataaaataaTATCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCAAACTCAGTAATAAGCTTTa                 <  1:2768502/69‑1 (MQ=255)
                                                           aaataaTATCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCAAACTCAGTAATAAGCTTTACTAat            <  1:2287831/71‑1 (MQ=255)
                                                            aataatATCGAGATATTATATAGACTGCTTTTTTTTattatt                                        >  1:2104773/1‑42 (MQ=255)
                                                            aataatATCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCAAACTCAGTAATAAGCTTTACTAat            <  1:117301/70‑1 (MQ=255)
                                                              taataTCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCaaa                                  >  1:3393806/1‑46 (MQ=255)
                                                              taataTCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCAAACTCAGTAATAAGCTTTACTAatata         >  1:630646/1‑71 (MQ=255)
                                                                atatCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCAAACTCAGTAATAAGCTTTACTAATATAtt       >  1:1486207/1‑71 (MQ=255)
                                                                atatCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCAAACTCAGTAATAAGCTTTACTAATATAtt       >  1:855072/1‑71 (MQ=255)
                                                                 tatCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCAAACt                                <  1:938218/45‑1 (MQ=255)
                                                                 tatCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCAAACTCAGTAATAAGCTTTACTAATATATTg      <  1:361875/71‑1 (MQ=255)
                                                                     gagaTATTATATAGACTGCTTTTTTTTATTATTCGCAAACTCAGTAATAAGCTTTACTAATATATTGCAGc  >  1:1660591/1‑71 (MQ=255)
                                                                      |                                                                     
AATGGAGAGATATTCGCTGTGGTGCATAAACTTTCATATGCAACGCATGAATATAAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCAAACTCAGTAATAAGCTTTACTAATATATTGCAGC  >  minE/1622177‑1622316

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: