Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,569,801 C→T 9.9% A24V (GCG→GTG) ‡ yibI → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,569,8010CT 9.9% 230.5 / 12.1 81A24V (GCG→GTG) ‡yibIpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (40/33);  new base T (5/3);  total (45/36)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

GTTTCTAAACTATTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCG  >  minE/2569733‑2569871
                                                                    |                                                                      
gTTTCTAAACTATTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGa                                                                      >  1:1036736/1‑71 (MQ=255)
 tttCTAAACTATTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATcgc                                                                        <  1:960544/68‑1 (MQ=255)
 tttCTAAACTATTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATcgc                                                                        <  1:3007544/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)
            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:1007652/1‑70 (MQ=255)
                    tGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCt                                                   >  1:894829/1‑70 (MQ=255)
                    tGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCt                                                   >  1:587957/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: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:3864142/66‑1 (MQ=255)
                                    cGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaaaa                                       <  1:3868828/66‑1 (MQ=255)
                                    cGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaaaa                                       <  1:995146/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:3765969/1‑44 (MQ=255)
                                                 aTTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGAt                                                >  1:1876429/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)
                                                             atcatcGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTc          <  1:1982544/70‑1 (MQ=255)
                                                              tcatcGCGATACATGACATCCCTTATCTGATTGCCaaaaa                                       >  1:2465508/1‑40 (MQ=255)
                                                              tcatcGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCa         >  1:1920899/1‑70 (MQ=255)
                                                              tcatcGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCa         >  1:57758/1‑70 (MQ=255)
                                                              tcatcGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCa         >  1:75986/1‑70 (MQ=255)
                                                                   gcgATACATGACATCCCTTATCTGATTGCCAAAAAg                                      >  1:1934114/1‑36 (MQ=255)
                                                                   gcgATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATAcggc   >  1:3639745/1‑71 (MQ=255)
                                                                    cgATACATGACATCCCTTATCTGATTGCCAAAAAGCGCACCCATCCCCATGCCGACGCTATTCATAcggag  >  1:916638/1‑69 (MQ=255)
                                                                    cgATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATAcggcg  <  1:1953051/71‑1 (MQ=255)
                                                                    |                                                                      
GTTTCTAAACTATTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCG  >  minE/2569733‑2569871

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: