Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 217,941 A→T 9.4% intergenic (‑33/‑73) 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,9410AT 9.4% 238.1 / 15.0 85intergenic (‑33/‑73)yafU/yafFpredicted inner membrane protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (36/41);  new base T (4/4);  total (40/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

GCCAAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTG  >  minE/217873‑218010
                                                                    |                                                                     
gCCAAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCag                                                                      <  1:2198570/70‑1 (MQ=255)
    aGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGaa                                                                  <  1:884661/70‑1 (MQ=255)
    aGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAg                                                                 <  1:2941092/71‑1 (MQ=255)
    aGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAg                                                                 <  1:4112173/71‑1 (MQ=255)
    aGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAg                                                                 <  1:3752923/71‑1 (MQ=255)
    aGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAg                                                                 >  1:273513/1‑71 (MQ=255)
      aaaaTTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGtt                                                               >  1:421166/1‑71 (MQ=255)
        aaTTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCagag                                                                    <  1:2983655/64‑1 (MQ=255)
        aaTTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCg                                                             >  1:4090969/1‑71 (MQ=255)
         aTTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGc                                                            <  1:812661/71‑1 (MQ=255)
              ccAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAg                                                       >  1:1662292/1‑71 (MQ=255)
                 gATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCagag                                                                    <  1:2701962/55‑1 (MQ=255)
                 gATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCaa                                                    <  1:3897863/71‑1 (MQ=255)
                   tctcGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGaa                                                                  >  1:1153410/1‑55 (MQ=255)
                   tctcGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAg                                                                 >  1:3820880/1‑56 (MQ=255)
                    ctcGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGcc                                                           >  1:2854829/1‑61 (MQ=255)
                    ctcGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATcc                                                 >  1:773493/1‑71 (MQ=255)
                       gCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGa                                               <  1:1239906/70‑1 (MQ=255)
                       gCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGa                                               >  1:988750/1‑70 (MQ=255)
                        ctctGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGaaa                                             >  1:1249923/1‑71 (MQ=255)
                           tGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGaa                                                                  >  1:1988392/1‑47 (MQ=255)
                           tGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAAc                                            <  1:3178897/69‑1 (MQ=255)
                           tGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACAGCagag                                                                    <  1:8359/45‑1 (MQ=255)
                            gAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAAcc                                           <  1:3476126/69‑1 (MQ=255)
                            gAACTCATAATAAATTCTATGGTTCTGCTCATTTAACCGCAgaga                                                                   <  1:2796176/45‑1 (MQ=255)
                             aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCca                                                          >  1:573681/1‑53 (MQ=255)
                             aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCa                                          <  1:3991144/69‑1 (MQ=255)
                               cTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTgg                                      >  1:2931494/1‑71 (MQ=255)
                                tCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAAc                                            <  1:394632/64‑1 (MQ=255)
                                 cATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGaaa                                             >  1:2664722/1‑62 (MQ=255)
                                 cATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCa                                    >  1:988695/1‑71 (MQ=255)
                                 cATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCa                                    >  1:388499/1‑71 (MQ=255)
                                  ataataAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGtt                                                               >  1:933284/1‑43 (MQ=255)
                                  ataataAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGc                                     <  1:1423329/69‑1 (MQ=255)
                                  ataataAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCa                                    <  1:2323429/70‑1 (MQ=255)
                                   taataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCg                                                             <  1:1708743/44‑1 (MQ=255)
                                   taataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACt                                        <  1:2302957/65‑1 (MQ=255)
                                   taataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACg                                  <  1:2386851/71‑1 (MQ=255)
                                     ataaATTCTATGGTTCTGCTGATTTAACCGCAgaga                                                                   <  1:1076021/36‑1 (MQ=255)
                                     ataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAg                                                                 <  1:166028/38‑1 (MQ=255)
                                     ataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCa                                          >  1:3191486/1‑61 (MQ=255)
                                     ataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTg                                <  1:2953796/71‑1 (MQ=255)
                                       aaaTTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTg                                >  1:3918264/1‑69 (MQ=255)
                                       aaaTTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTg                                >  1:10688/1‑69 (MQ=255)
                                         aTTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGgaga                            <  1:2564520/71‑1 (MQ=255)
                                          ttCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCAc                                   <  1:1307870/63‑1 (MQ=255)
                                           tCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGa                                               <  1:558680/50‑1 (MQ=255)
                                           tCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGaa                           >  1:3966259/1‑70 (MQ=255)
                                               tGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATaa                        <  1:975312/69‑1 (MQ=255)
                                               tGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATa                         <  1:245991/68‑1 (MQ=255)
                                               cagTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGaaa                                             >  1:849788/3‑48 (MQ=255)
                                                agTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGc                      <  1:134594/69‑1 (MQ=255)
                                                  ttCTGCTGATTTTACCGCAGAGAAGTTCGCCACAGCAAt                                                   <  1:1140689/39‑1 (MQ=255)
                                                  ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAAt                                                   <  1:4005408/39‑1 (MQ=255)
                                                  ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAAt                                                   <  1:2432146/39‑1 (MQ=255)
                                                  ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                         <  1:1022923/49‑1 (MQ=255)
                                                  ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                         <  1:3856120/49‑1 (MQ=255)
                                                  ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                         <  1:3905891/49‑1 (MQ=255)
                                                  ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTg                    <  1:2774097/70‑1 (MQ=255)
                                                  ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTg                    <  1:1339770/70‑1 (MQ=255)
                                                    ctgctgATTTAACCGCTGAGAAATTCGCCACAGCGATc                                                  >  1:118034/1‑38 (MQ=25)
                                                    ctgctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTg                                       >  1:784773/1‑49 (MQ=255)
                                                    ctgctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAg                       >  1:18459/1‑65 (MQ=255)
                                                     tgctgATTTAACCGCTGAGAAATTCGCCACAGCGATCCGaaat                                            >  1:3563026/1‑42 (MQ=25)
                                                      gctgATTTAAGCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTg                                >  1:2471184/1‑54 (MQ=255)
                                                       ctgATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGcact                <  1:4125316/69‑2 (MQ=255)
                                                       ctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGa                                               <  1:2732487/38‑1 (MQ=255)
                                                       ctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAg               <  1:3776470/70‑1 (MQ=255)
                                                         gATTTAACCGCAGAGAAGTTCGCCAGAGCAATCCGAAACCACTGGCACGt                                 <  1:1817019/50‑1 (MQ=255)
                                                         gATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCg            >  1:326242/1‑71 (MQ=255)
                                                          aTTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGt           >  1:512236/1‑71 (MQ=255)
                                                          aTTTAACCGCTGAGAAATCCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGc                   <  1:3992780/63‑1 (MQ=255)
                                                          aTTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAAcc                                           <  1:1639405/39‑1 (MQ=255)
                                                          aTTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                         >  1:618077/1‑41 (MQ=255)
                                                          aTTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACt                                        >  1:2921357/1‑42 (MQ=255)
                                                          aTTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACt                                        >  1:3464803/1‑42 (MQ=255)
                                                           tttAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTc          >  1:2025755/1‑71 (MQ=255)
                                                           tttAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTc          <  1:1901136/71‑1 (MQ=255)
                                                           tttAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTg                                <  1:1959724/49‑1 (MQ=255)
                                                            ttAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                         >  1:2563916/1‑39 (MQ=255)
                                                            ttAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGca                  >  1:1873956/1‑62 (MQ=255)
                                                            ttAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCt         >  1:496185/1‑71 (MQ=255)
                                                            ttAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCt         >  1:215941/1‑71 (MQ=255)
                                                            ttAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCt         >  1:2440319/1‑71 (MQ=255)
                                                               aCCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGga                              >  1:3381260/1‑47 (MQ=255)
                                                               aCCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTg        <  1:1041706/69‑1 (MQ=255)
                                                               aCCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGa      <  1:1979510/71‑1 (MQ=255)
                                                               aCCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTCCACAGGCGTCTgg       <  1:1929880/70‑1 (MQ=255)
                                                                   cAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTg  >  1:3055515/1‑71 (MQ=255)
                                                                    |                                                                     
GCCAAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTG  >  minE/217873‑218010

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: