Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,569,811 C→T 9.7% D27D (GAC→GAT yibI → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,569,8110CT 9.7% 236.5 / 12.0 82D27D (GAC→GATyibIpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (42/32);  new base T (5/3);  total (47/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

ACTATTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGA  >  minE/2569741‑2569881
                                                                      |                                                                      
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:1007652/1‑70 (MQ=255)
            tGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCt                                                             >  1:894829/1‑70 (MQ=255)
            tGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCt                                                             >  1:587957/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)
               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)
                   gCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGa                                                           <  1:1467725/65‑1 (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: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:3868828/66‑1 (MQ=255)
                            cGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaaaa                                                 <  1:3618881/66‑1 (MQ=255)
                            cGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaaaa                                                 <  1:995146/66‑1 (MQ=255)
                            cGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaaaa                                                 <  1:3864142/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:3765969/1‑44 (MQ=255)
                                         aTTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGAt                                                          >  1:1876429/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:1920899/1‑70 (MQ=255)
                                                      tcatcGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCa                   >  1:75986/1‑70 (MQ=255)
                                                      tcatcGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCa                   >  1:57758/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)
                                                              aTACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTAtt                     <  1:3783623/60‑1 (MQ=255)
                                                               tACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATc                                      <  1:1830793/42‑1 (MQ=255)
                                                                 cATGACATCCCTTATCTGATTGCCAAAAAGCGCAAc                                          >  1:3943071/1‑36 (MQ=255)
                                                                 cATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATa                 >  1:1921031/1‑61 (MQ=255)
                                                                 cATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATa                 >  1:990938/1‑61 (MQ=255)
                                                                 cATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCt        <  1:124014/70‑1 (MQ=255)
                                                                   tGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACg                          >  1:1094361/1‑50 (MQ=255)
                                                                    gACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGt    >  1:2138577/1‑71 (MQ=255)
                                                                    gACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGt    >  1:1269066/1‑71 (MQ=255)
                                                                      cATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGcc                              <  1:276661/43‑1 (MQ=255)
                                                                      cATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGt    >  1:3620947/1‑69 (MQ=255)
                                                                      cATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGa  <  1:1503116/71‑1 (MQ=255)
                                                                      cATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGa  <  1:3176992/71‑1 (MQ=255)
                                                                      |                                                                      
ACTATTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGA  >  minE/2569741‑2569881

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: