Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,572,941 G→A 49.0% intergenic (+1709/‑70) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,572,9410GA49.0% ‑0.4 / 69.0 51intergenic (+1709/‑70)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (10/16);  new base A (8/17);  total (18/33)
Fisher's exact test for biased strand distribution p-value = 7.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAAATCATTTGGTTGGCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAACCAGACGGTAATTTAA  >  minE/2572873‑2573011
                                                                    |                                                                      
aTAAATCATTTGGTTGGCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGAt                                                                      <  1:3295118/71‑1 (MQ=255)
aTAAATCATTTGGTTGGCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGAt                                                                      <  1:3188118/71‑1 (MQ=255)
aTAAATCATTTGGTTGGCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGAt                                                                      >  1:187977/1‑71 (MQ=255)
aTAAATCATTTGGTTGGCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAAAt                                                                      <  1:619320/71‑1 (MQ=255)
aTAAATCATTTGGTTGGCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAAAt                                                                      >  1:844976/1‑71 (MQ=255)
    aTCATTTGGTTGGCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGt                                                                  <  1:1356765/71‑1 (MQ=255)
     tCATTTGGTTGGCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGt                                                                  >  1:1835705/1‑70 (MQ=255)
     tCATTTGGTTGGCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTa                                                                 <  1:3258565/71‑1 (MQ=255)
       aTTTGGTTGGCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAg                                                               <  1:3378349/71‑1 (MQ=255)
           ggttggCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGat                                                             <  1:2177753/69‑1 (MQ=255)
           ggttggCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGat                                                             <  1:1533305/69‑1 (MQ=255)
           ggttggCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTAGTAAGata                                                            <  1:3149796/70‑1 (MQ=255)
            gttggCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAAATTCATAAGATATc                                                          <  1:2446983/71‑1 (MQ=255)
             ttggCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAAATTc                                                                    >  1:259625/1‑60 (MQ=255)
              tggCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAg                                                        >  1:967142/1‑71 (MQ=255)
                gcgcACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAAATTCGTAAg                                                               >  1:3002976/1‑62 (MQ=255)
                      tGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAAATTCGTAAGATATCAGCCACtat                                                 >  1:1929913/1‑70 (MQ=255)
                        gAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATAcc                                              >  1:1452939/1‑71 (MQ=255)
                        gAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAAATTCGTAAGatat                                                           <  1:1585914/58‑1 (MQ=255)
                        gAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAAATTCGTAAGatat                                                           <  1:1979414/58‑1 (MQ=255)
                           taatCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGgaga                                                                       <  1:1288950/43‑1 (MQ=255)
                           taatCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTaa                                                                >  1:77168/1‑50 (MQ=255)
                           taatCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGa                                            >  1:1230610/1‑70 (MQ=255)
                           taatCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAAATTCGt                                                                  <  1:703616/48‑1 (MQ=255)
                            aatCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACta                                                  >  1:1390285/1‑63 (MQ=255)
                             atCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGatat                                         <  1:977178/71‑1 (MQ=255)
                                  aGTTTTATCACTACAGCTTATAGAGGCTTAAGGAAATTCGTAAGATATCAGCCACTATACCGAtata                                        <  1:1743912/67‑1 (MQ=255)
                                    ttttATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAaataat                                  <  1:357918/71‑1 (MQ=255)
                                      ttATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAaataat                                  >  1:2394901/1‑69 (MQ=255)
                                      ttATCACTACAGCTTATAGAGGCTTAAGGAAATTCGTAAGATATCAGCCACta                                                  >  1:3761166/1‑53 (MQ=255)
                                      ttATCACTACAGCTTATAGAGGCTTAAGGAAATTCGTAAGATATCAGCCACTATACCGATATAAataata                                 <  1:3121627/70‑1 (MQ=255)
                                       tATCACTACAGCTTATAGAGGCTTAAGGAAATTCGTAAGATATCAGCCACTATACCGATATAAAtaataa                                <  1:1447203/70‑1 (MQ=255)
                                        aTCACTACAGCTTATAGAGGCTTAAGGAAATTCGTAAGat                                                             <  1:3673999/40‑1 (MQ=255)
                                        aTCACTACAGCTTATAGAGGCTTAAGGAAATTCGTAAGat                                                             <  1:2996673/40‑1 (MQ=255)
                                        aTCACTACAGCTTATAGAGGCTTAAGGAAATTCGTAAGat                                                             <  1:1323217/40‑1 (MQ=255)
                                          cACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACt                            >  1:1377222/1‑71 (MQ=255)
                                          cACTACAGCTTATAGAGGCTTAAGGAAATTCGTAAGatat                                                           >  1:1882068/1‑40 (MQ=255)
                                           aCTACAGCTTATAGAGGCTTAAGGAAATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTc                           <  1:491364/71‑1 (MQ=255)
                                           aCTACAGCTTATAGAGGCTTAAGGAAATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTc                           <  1:2833291/71‑1 (MQ=255)
                                             tACAGCTTATAGAGGCTTAAGGAAATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCAc                         >  1:3801722/1‑71 (MQ=255)
                                               cAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACt                        <  1:3208275/70‑1 (MQ=255)
                                                 gCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACt                        <  1:2518800/68‑1 (MQ=255)
                                                     ataGGGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAAc                 >  1:2736402/1‑71 (MQ=255)
                                                      taGAGGCTTAAGGAAATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACCTGCAAAcc                <  1:925975/71‑1 (MQ=255)
                                                        gagGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATataaataa                                   <  1:3305894/50‑1 (MQ=255)
                                                        gagGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATataaataa                                   <  1:1871633/50‑1 (MQ=255)
                                                        gagGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATataaataa                                   <  1:1579516/50‑1 (MQ=255)
                                                        gagGCTTAAGGAAATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACCTGCAAACCAg              <  1:2379440/71‑1 (MQ=255)
                                                          ggCTTAAGGAAATTCGTAAGATATCAGCCACTATACCGat                                           <  1:599604/40‑1 (MQ=255)
                                                            cTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAACCAGACgg          <  1:1317361/71‑1 (MQ=255)
                                                             ttAAGGAAATTCGTAAGATATCAGCCACTATAGCGATATAAATATTAAGACt                            <  1:599544/52‑1 (MQ=255)
                                                                 ggAAATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACCTGCAAACCAGACGGTAAtt     >  1:1195105/1‑71 (MQ=255)
                                                                   agaTTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAACCAGACGGTAATTTa   <  1:1640416/71‑1 (MQ=255)
                                                                    gaTTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAACCAGACGGTAATTTaa  <  1:985088/71‑1 (MQ=255)
                                                                    |                                                                      
ATAAATCATTTGGTTGGCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAACCAGACGGTAATTTAA  >  minE/2572873‑2573011

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: