Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,569,793 A→G 11.0% G21G (GGA→GGG yibI → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,569,7930AG11.0% 220.5 / 15.4 82G21G (GGA→GGGyibIpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (37/36);  new base G (6/3);  total (43/39)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.84e-01

AGGCAGATATGTTTCTAAACTATTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGAC  >  minE/2569723‑2569856
                                                                      |                                                               
aGGCAGATATGTTTCTAAACTATTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGa                                                                 <  1:3706631/71‑1 (MQ=255)
     aatatGTTCCTGGACTATTTTGCACTGGGAGTGCTTATTTTTGTATTTCTGGTGATCTTCTATGGGatcat                                                            >  1:1099372/2‑71 (MQ=255)
       tatGTTTCTAAACTATTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATcg                                                          >  1:239493/1‑71 (MQ=255)
          gTTTCTAAACTATTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGa                                                       >  1:1036736/1‑71 (MQ=255)
           tttCTAAACTATTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATcgc                                                         <  1:3007544/68‑1 (MQ=255)
           tttCTAAACTATTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATcgc                                                         <  1:960544/68‑1 (MQ=255)
             tCTAAACTATTTCGCGCTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATcgc                                                         <  1:1639939/66‑1 (MQ=255)
               tAAACTATTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACAt                                                  >  1:476931/1‑71 (MQ=255)
                aaaCTATTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACa                                                   >  1:1031908/1‑69 (MQ=255)
                 aaCTATTTCGCGTTGGTAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACAt                                                  <  1:2338748/69‑1 (MQ=255)
                 aaCTATTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACa                                                   <  1:1647543/68‑1 (MQ=255)
                 aaCTATTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACAt                                                  <  1:2121342/69‑1 (MQ=255)
                  aCTATTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGAc                                               <  1:3658050/71‑1 (MQ=255)
                      tttCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAatcatc                                                           <  1:2773618/55‑1 (MQ=255)
                      tttCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATcc                                           >  1:3021482/1‑71 (MQ=255)
                        tCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATccc                                          <  1:1708663/70‑1 (MQ=255)
                           cgTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTAt                                      >  1:960427/1‑71 (MQ=255)
                             ttGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCtt                                        <  1:101214/67‑1 (MQ=255)
                              tGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCt                                    >  1:3114019/1‑70 (MQ=255)
                              tGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCt                                    >  1:587957/1‑70 (MQ=255)
                              tGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCt                                    >  1:894829/1‑70 (MQ=255)
                              tGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCt                                    >  1:1007652/1‑70 (MQ=255)
                               gggAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACAt                                             <  1:2919653/60‑1 (MQ=255)
                               gggAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGAAATCATCGCGATACATGAc                                               <  1:86651/58‑1 (MQ=255)
                                ggAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACa                                                   >  1:1890939/1‑53 (MQ=255)
                                ggAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACa                                                   >  1:1870730/1‑53 (MQ=255)
                                 gAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGAt                                 >  1:467429/1‑70 (MQ=255)
                                  aGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTg                               <  1:3123798/71‑1 (MQ=255)
                                   gTGCTTATTTTTGTATTTCTGGTGATCTTCTATGGGATCATTATTTTACATGATATTCCCTACCTGATTg                               >  1:1295997/1‑70 (MQ=255)
                                   gTGCTTATTTTTGTATTTCTGGTGATCTTCTATGGGATCATTATTTTACATGATATTCCCTACCTGATTg                               >  1:1642966/1‑70 (MQ=255)
                                   gTGCTTATTTTTGTATTTCTGGTGATCTTCTATGGGATCATTATTTTACATGATATTCCCTACCTGATTGc                              >  1:3934834/1‑71 (MQ=255)
                                    tGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATc                                            <  1:1674716/56‑1 (MQ=255)
                                     gCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATAc                                                    >  1:2033589/1‑47 (MQ=255)
                                     gCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGa                                  <  1:1467725/65‑1 (MQ=255)
                                      cTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGa                                                       >  1:1790345/1‑43 (MQ=255)
                                      cTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCa                            >  1:462838/1‑70 (MQ=255)
                                      cTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCa                            >  1:1497526/1‑70 (MQ=255)
                                       tGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGAtt                                >  1:3436943/1‑65 (MQ=255)
                                       tGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGAt                                 <  1:3059743/64‑1 (MQ=255)
                                         aTCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaa                          <  1:3897787/69‑1 (MQ=255)
                                         aTCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaa                          <  1:2674449/69‑1 (MQ=255)
                                           cTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaaaa                        >  1:2114060/1‑69 (MQ=255)
                                           cTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAgc                      >  1:2045046/1‑71 (MQ=255)
                                           cTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAg                       >  1:3757642/1‑70 (MQ=255)
                                             tCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGAt                                 <  1:3266974/58‑1 (MQ=255)
                                              cGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCt                                    >  1:1502056/1‑54 (MQ=255)
                                              cGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaaaa                        <  1:3618881/66‑1 (MQ=255)
                                              cGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaaaa                        <  1:2397959/66‑1 (MQ=255)
                                              cGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaaaa                        <  1:995146/66‑1 (MQ=255)
                                              cGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaaaa                        <  1:3864142/66‑1 (MQ=255)
                                              cGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaaaa                        <  1:3868828/66‑1 (MQ=255)
                                              cGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAgcgc                    <  1:3941678/70‑1 (MQ=255)
                                              cGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAg                       >  1:3363620/1‑67 (MQ=255)
                                               gTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCaa                  <  1:478405/71‑1 (MQ=255)
                                                tCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCt                                    <  1:2798126/52‑1 (MQ=255)
                                                tCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCa                   >  1:2025948/1‑69 (MQ=255)
                                                  tttCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCa               >  1:798974/1‑71 (MQ=255)
                                                   ttCTGGTGATCTTCTATGGGATCATTATTTTACATGATATTCCCTACCTGATTGCCAAAAAACGTaatc                >  1:2529465/1‑67 (MQ=255)
                                                   ttCTGGTGATCTTCTATGGGATCATTATTTTACATGATATTCCCTACCTGATTGCCAAAAAACGTaat                 >  1:2307023/1‑67 (MQ=255)
                                                    tCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAAcc                <  1:1955310/68‑1 (MQ=255)
                                                    tCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCAt              <  1:995434/70‑1 (MQ=255)
                                                    tCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATc             >  1:1074034/1‑71 (MQ=255)
                                                    tCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATc             >  1:2266703/1‑71 (MQ=255)
                                                    tCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGACAAAAAgcg                     >  1:3170979/1‑63 (MQ=255)
                                                     cTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATcc            >  1:2047506/1‑71 (MQ=255)
                                                      tGGTGATCTTCTATGGGATCATTATTTTACATGATATTCCCTACCTGATTGCCaaaaa                        <  1:696601/58‑1 (MQ=255)
                                                       ggTGATCTTCTATGGGATCATTATTTTACATGATATTCCCTACCTGATTGCCAAAAAacgta                   <  1:1893992/62‑3 (MQ=255)
                                                        gTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaaa                         <  1:1866585/55‑1 (MQ=255)
                                                        gTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCaa                  <  1:748539/62‑1 (MQ=255)
                                                        gTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATccc           <  1:1767787/69‑1 (MQ=255)
                                                        gTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCa         <  1:2789212/71‑1 (MQ=255)
                                                        gTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCa         >  1:2533292/1‑71 (MQ=255)
                                                        gTGATCTTCTATGGGATCATTATTTTACATGATATTCCCTACCTGATTGCCAAAAAACGTAATCATcct           <  1:3379106/69‑2 (MQ=255)
                                                         tGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaa                           <  1:3122978/52‑1 (MQ=255)
                                                           aTTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGAt                                 >  1:1876429/1‑44 (MQ=255)
                                                           aTTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGAt                                 >  1:3765969/1‑44 (MQ=255)
                                                            ttttttATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGcc     >  1:2100605/1‑71 (MQ=255)
                                                             tttttATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGcc     >  1:2744494/1‑70 (MQ=255)
                                                             tttttATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGcc     >  1:3350040/1‑70 (MQ=255)
                                                              ttttATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCg    <  1:249353/70‑1 (MQ=255)
                                                               tttATGGAATCATCGCGATACATGACATCCCTTATCTg                                   >  1:572273/1‑38 (MQ=255)
                                                               tttATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGAc  >  1:571351/1‑71 (MQ=255)
                                                                    ggAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaaa                         <  1:3220887/43‑1 (MQ=255)
                                                                      |                                                               
AGGCAGATATGTTTCTAAACTATTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGAC  >  minE/2569723‑2569856

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: