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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,569,7840TC11.7% 203.2 / 15.3 77I18I (ATT→ATCyibIpredicted 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

GTGAATCCGAGGCAGATATGTTTCTAAACTATTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCC  >  minE/2569714‑2569853
                                                                      |                                                                     
gTGAATCCGAGGCAGATATGTTTCTAAACTATTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGAtt                                                                       <  1:3331806/71‑1 (MQ=255)
gTGAATCCGAGGCAGATATGTTTCTAAACTATTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGAtt                                                                       <  1:1211800/71‑1 (MQ=255)
  gAATCCGAGGCAGATATGTTTCTAAACTATTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGAttt                                                                      >  1:389360/1‑70 (MQ=255)
    aTCCGAGGCAGATATGTTTCTAAACTATTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGAtttt                                                                     <  1:2171483/69‑1 (MQ=255)
     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: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)
                              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: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: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:995146/66‑1 (MQ=255)
                                                       cGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaaaa                     <  1:3868828/66‑1 (MQ=255)
                                                       cGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaaaa                     <  1:3864142/66‑1 (MQ=255)
                                                       cGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaaaa                     <  1:2397959/66‑1 (MQ=255)
                                                       cGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaaaa                     <  1:3618881/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:3350040/1‑70 (MQ=255)
                                                                      tttttATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGcc  >  1:2744494/1‑70 (MQ=255)
                                                                      |                                                                     
GTGAATCCGAGGCAGATATGTTTCTAAACTATTTCGCGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCC  >  minE/2569714‑2569853

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: