Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,569,799 C→T 10.5% I23I (ATC→ATT yibI → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,569,7990CT10.5% 213.1 / 12.3 76I23I (ATC→ATTyibIpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (36/32);  new base T (5/3);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

TATGTTTCTAAACTATTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCA  >  minE/2569730‑2569864
                                                                     |                                                                 
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:1007652/1‑70 (MQ=255)
                       tGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCt                                            >  1:587957/1‑70 (MQ=255)
                       tGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCt                                            >  1:894829/1‑70 (MQ=255)
                       tGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCt                                            >  1:3114019/1‑70 (MQ=255)
                        gggAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACAt                                                     <  1:2919653/60‑1 (MQ=255)
                        gggAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGAAATCATCGCGATACATGAc                                                       <  1:86651/58‑1 (MQ=255)
                         ggAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACa                                                           >  1:1870730/1‑53 (MQ=255)
                         ggAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACa                                                           >  1:1890939/1‑53 (MQ=255)
                          gAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGAt                                         >  1:467429/1‑70 (MQ=255)
                           aGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTg                                       <  1:3123798/71‑1 (MQ=255)
                            gTGCTTATTTTTGTATTTCTGGTGATCTTCTATGGGATCATTATTTTACATGATATTCCCTACCTGATTg                                       >  1:1642966/1‑70 (MQ=255)
                            gTGCTTATTTTTGTATTTCTGGTGATCTTCTATGGGATCATTATTTTACATGATATTCCCTACCTGATTg                                       >  1:1295997/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:1497526/1‑70 (MQ=255)
                               cTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCa                                    >  1:462838/1‑70 (MQ=255)
                                tGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGAtt                                        >  1:3436943/1‑65 (MQ=255)
                                tGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGAt                                         <  1:3059743/64‑1 (MQ=255)
                                  aTCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaa                                  <  1:2674449/69‑1 (MQ=255)
                                  aTCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaa                                  <  1:3897787/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:2397959/66‑1 (MQ=255)
                                       cGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaaaa                                <  1:3864142/66‑1 (MQ=255)
                                       cGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaaaa                                <  1:995146/66‑1 (MQ=255)
                                       cGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaaaa                                <  1:3618881/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:2533292/1‑71 (MQ=255)
                                                 gTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCa                 <  1:2789212/71‑1 (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:3350040/1‑70 (MQ=255)
                                                      tttttATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGcc             >  1:2744494/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)
                                                                atcatcGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTc   <  1:1982544/70‑1 (MQ=255)
                                                                 tcatcGCGATACATGACATCCCTTATCTGATTGCCaaaaa                                >  1:2465508/1‑40 (MQ=255)
                                                                 tcatcGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCa  >  1:57758/1‑70 (MQ=255)
                                                                 tcatcGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCa  >  1:75986/1‑70 (MQ=255)
                                                                 tcatcGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCa  >  1:1920899/1‑70 (MQ=255)
                                                                     |                                                                 
TATGTTTCTAAACTATTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCA  >  minE/2569730‑2569864

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: