Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,569,787 T→C 11.7% F19F (TTT→TTC yibI → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,569,7870TC11.7% 203.4 / 15.9 77F19F (TTT→TTCyibIpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (33/35);  new base C (6/3);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCGAGGCAGATATGTTTCTAAACTATTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGAC  >  minE/2569719‑2569856
                                                                    |                                                                     
tCCGAGGCAGATATGTTTCTAAACTATTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGAtttttt                                                                      <  1:2747569/70‑1 (MQ=255)
    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)
               tttCTAAACTATTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTa                                                                     >  1:1519263/1‑56 (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)
                         aTTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGAttttt                                                                       <  1:3561090/44‑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)
                                 ttGGGAGTGCTGATCTTCGTCTTTCTGGTGAttttt                                                                       <  1:227984/36‑1 (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:3114019/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: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: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: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:2397959/66‑1 (MQ=255)
                                                  cGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaaaa                        <  1:3868828/66‑1 (MQ=255)
                                                  cGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaaaa                        <  1:3864142/66‑1 (MQ=255)
                                                  cGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaaaa                        <  1:3618881/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: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)
                                                                    |                                                                     
TCCGAGGCAGATATGTTTCTAAACTATTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGAC  >  minE/2569719‑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: