Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 218,157 C→T 9.7% D48D (GAC→GAT yafF → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE218,1570CT 9.7% 345.4 / 18.7 113D48D (GAC→GATyafFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (48/54);  new base T (4/7);  total (52/61)
Fisher's exact test for biased strand distribution p-value = 5.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 114 total aligned reads displayed.

ATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTC  >  minE/218090‑218224
                                                                   |                                                                   
atatTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAg                                                                   >  1:3181590/1‑70 (MQ=255)
atatTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAg                                                                   <  1:2291658/70‑1 (MQ=255)
atatTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGa                                                                  <  1:3355322/71‑1 (MQ=255)
atatTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGa                                                                  <  1:403665/71‑1 (MQ=255)
 tatTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAATATGCGAAAAGCAGCCATGGACAGaa                                                                 <  1:647353/71‑1 (MQ=255)
 tatTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGa                                                                  <  1:145670/70‑1 (MQ=255)
  atTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGAc                                                                     >  1:3642772/1‑66 (MQ=255)
  atTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGaaa                                                                >  1:2332996/1‑71 (MQ=255)
   ttttGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGa                                                                  <  1:3254819/68‑1 (MQ=255)
   ttttGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAAc                                                               <  1:3823550/71‑1 (MQ=255)
   ttttGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAAc                                                               >  1:789687/1‑71 (MQ=255)
    tttGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGa                                                                  <  1:1430416/67‑1 (MQ=255)
     ttaaCGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACt                                                              <  1:2144779/67‑1 (MQ=255)
      tGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTa                                                             <  1:291078/70‑1 (MQ=255)
       aaCGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTatc                                                           >  1:2302951/2‑69 (MQ=255)
       aaCGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTatc                                                           >  1:1913283/2‑69 (MQ=255)
        aCGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAAc                                                               <  1:3040248/66‑1 (MQ=255)
        aCGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAAc                                                               <  1:3789528/66‑1 (MQ=255)
           aaTGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGc                                                       <  1:1917859/71‑1 (MQ=255)
            aTGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCg                                                      <  1:1117020/71‑1 (MQ=255)
            aTGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCg                                                      >  1:3207457/1‑71 (MQ=255)
             tGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGt                                                     >  1:2223288/1‑71 (MQ=255)
             tGAGAAGGTATTAAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACt                                                              >  1:2022421/1‑62 (MQ=255)
                  aGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCt                                                          >  1:465618/1‑61 (MQ=255)
                   ggTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGaaa                                                                >  1:3102749/1‑54 (MQ=255)
                    gTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCt                                              >  1:2293620/1‑71 (MQ=255)
                    gTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGAACAGAAACTACCTGGCGTCAGTCCt                                              >  1:3230561/1‑71 (MQ=255)
                     tATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCtt                                             >  1:2291435/1‑71 (MQ=255)
                      aTTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTg                                            >  1:807176/1‑71 (MQ=255)
                      aTTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTg                                            >  1:1203054/1‑71 (MQ=255)
                      aTTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTg                                            >  1:4119509/1‑71 (MQ=255)
                       ttCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGa                                                                  <  1:1340411/48‑1 (MQ=255)
                       ttCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGa                                                                  <  1:4091948/48‑1 (MQ=255)
                       ttCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGa                                                                  <  1:3678727/48‑1 (MQ=255)
                         cAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGaaa                                                                >  1:2014657/1‑48 (MQ=255)
                         cAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGaaa                                                                >  1:3890734/1‑48 (MQ=255)
                         cAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCt                                                          <  1:3415655/54‑1 (MQ=255)
                         cAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCg                                          >  1:187658/1‑70 (MQ=255)
                          aaGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGaaa                                                                >  1:2562296/1‑47 (MQ=255)
                            ggCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCggg                                        >  1:1901107/1‑69 (MQ=255)
                              cAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGc                                    <  1:833421/71‑1 (MQ=255)
                               aGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGt                                                     <  1:2840579/53‑1 (MQ=255)
                               aGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCg                                   <  1:4057331/71‑1 (MQ=255)
                               aGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACGGGGAGCg                                   <  1:3230083/71‑1 (MQ=255)
                                gggTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCgggg                                       <  1:4075973/66‑1 (MQ=255)
                                 ggTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGc                                    <  1:3726520/68‑1 (MQ=255)
                                  gTTAAGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGCGGGGa                                      <  1:600666/65‑1 (MQ=255)
                                  gTTAAGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGCGGGGa                                      <  1:1968890/65‑1 (MQ=255)
                                  gTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGaaa                                                                >  1:2666520/1‑39 (MQ=255)
                                  gTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTAc                                                            >  1:4132895/1‑43 (MQ=255)
                                   ttAAGACGTAAGATGCGAAAAGCAGCCATGGACAGa                                                                  <  1:1003594/36‑1 (MQ=255)
                                   ttAAGACGTAAGATGCGAAAAGCAGCCATGGACAGa                                                                  <  1:122705/36‑1 (MQ=255)
                                   ttAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCt                               <  1:196569/71‑1 (MQ=255)
                                   ttAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACgg                                         >  1:1851928/1‑61 (MQ=255)
                                    tAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCtt                              >  1:1841244/1‑71 (MQ=255)
                                     aaGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCttt                             >  1:2480533/1‑71 (MQ=255)
                                      aGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACt                                                              <  1:2932131/37‑1 (MQ=255)
                                      aGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCggg                                        <  1:1726861/59‑1 (MQ=255)
                                      aGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCgg                                         >  1:1368511/1‑58 (MQ=255)
                                      aGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCgg                                         >  1:377940/1‑58 (MQ=255)
                                       gACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCg                           <  1:3654419/71‑1 (MQ=255)
                                         cGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCt                                                          >  1:2426795/1‑38 (MQ=38)
                                         cGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGTAGCGGGCTTTCgta                         <  1:3869391/71‑1 (MQ=255)
                                         cGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgt                          <  1:1015421/70‑1 (MQ=255)
                                         cGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgt                          <  1:2551760/70‑1 (MQ=255)
                                         cGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCg                           >  1:1410887/1‑69 (MQ=255)
                                          gTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCg                                          <  1:3732799/53‑1 (MQ=255)
                                          gTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAg                                                  <  1:3817410/45‑1 (MQ=255)
                                          gTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAg                                                  <  1:1195059/45‑1 (MQ=255)
                                          gTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAg                                                  <  1:1199974/45‑1 (MQ=255)
                                           tAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGa                                      >  1:753863/1‑56 (MQ=255)
                                           tAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgtagt                       <  1:1559773/71‑1 (MQ=255)
                                           tAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgtagt                       <  1:308399/71‑1 (MQ=255)
                                           tAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgtagt                       <  1:2904104/71‑1 (MQ=255)
                                           tAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgtag                        >  1:3508145/1‑70 (MQ=255)
                                            aaGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCg                           >  1:665057/1‑66 (MQ=255)
                                             aGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAGc                                    <  1:849033/56‑1 (MQ=255)
                                             aGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTc                                                    >  1:675423/1‑40 (MQ=255)
                                             aGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCtta                                            <  1:3404136/48‑2 (MQ=255)
                                             aGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGGGGGAGCGGGCTTTCGTAGTaa                     >  1:457266/1‑70 (MQ=255)
                                             aGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACgg                                         >  1:3920696/1‑51 (MQ=255)
                                              gATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCg                           <  1:3550700/64‑1 (MQ=255)
                                                 gCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgta                         >  1:3530399/1‑63 (MQ=255)
                                                   gAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCtta                                            <  1:2863869/42‑2 (MQ=255)
                                                   gAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACGGGGAGCGGGCTTTCgtaatc                      <  1:3379625/64‑4 (MQ=255)
                                                    aaaaGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGc                                           <  1:2598712/42‑1 (MQ=255)
                                                    aaaaGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGc                                           <  1:2527603/42‑1 (MQ=255)
                                                    aaaaGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAg                                     <  1:772944/48‑1 (MQ=255)
                                                      aaGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGCg                                          >  1:2681255/1‑41 (MQ=255)
                                                      aaGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGa             >  1:1611346/1‑70 (MQ=255)
                                                         cagcCATGGACAGAAACTACCTGGCGTCAGTACTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGActct         >  1:1312792/1‑71 (MQ=255)
                                                           gcCATGGACAGAAACTACCTGGCGTCAGTGCTTGCGGGGAGCGGGCTTTCgtagt                       <  1:3970781/55‑1 (MQ=255)
                                                           gcCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTcc       >  1:1153224/1‑71 (MQ=255)
                                                            ccATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATAcccc               >  1:1669971/1‑62 (MQ=255)
                                                            ccATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTccc      <  1:4069344/71‑1 (MQ=255)
                                                              aTGGATATAAACTATCTCGCGTCAGTCCTTGCGGGGa                                      <  1:2497878/37‑1 (MQ=25)
                                                                ggACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAg                                     <  1:3851405/36‑1 (MQ=255)
                                                                ggACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCg                           <  1:811416/46‑1 (MQ=255)
                                                                 gACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTc  <  1:1583840/70‑1 (MQ=255)
                                                                 gACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGGGGCTTTCGTAGTAATACCCCGa             >  1:567215/1‑58 (MQ=255)
                                                                   |                                                                   
ATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTC  >  minE/218090‑218224

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: