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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE218,1570CT15.9% 292.4 / 39.1 119D48D (GAC→GATyafFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (51/49);  new base T (5/14);  total (56/63)
Fisher's exact test for biased strand distribution p-value = 7.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 119 total aligned reads displayed.

ATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCG  >  minE/218090‑218222
                                                                   |                                                                 
atatTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAg                                                                 >  1:1397117/1‑70 (MQ=255)
atatTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAg                                                                 >  1:3891349/1‑70 (MQ=255)
atatTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAg                                                                 <  1:1122330/70‑1 (MQ=255)
 tatTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACa                                                                  <  1:3529269/68‑1 (MQ=255)
 tatTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGa                                                                <  1:3083070/70‑1 (MQ=255)
 tatTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGa                                                                <  1:3439000/70‑1 (MQ=255)
 tatTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGATAGaa                                                               <  1:1368838/71‑1 (MQ=255)
  atTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATAGACAGaaa                                                              >  1:2105796/1‑71 (MQ=255)
   ttttGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAAc                                                             >  1:1118617/1‑71 (MQ=255)
   ttttGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAAc                                                             <  1:1755981/71‑1 (MQ=255)
     ttGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACt                                                            <  1:2634741/70‑1 (MQ=255)
     ttGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTa                                                           >  1:3356191/1‑71 (MQ=255)
     ttGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTa                                                           >  1:2226934/1‑71 (MQ=255)
       gACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTAc                                                          <  1:1429119/70‑1 (MQ=255)
        aCGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCt                                                        <  1:486881/71‑1 (MQ=255)
          gAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTg                                                       >  1:1157819/1‑70 (MQ=255)
               agaAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCa                                                 <  1:2567343/71‑1 (MQ=255)
                taaGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAg                                                <  1:2192812/70‑1 (MQ=255)
                  aGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGc                                                     <  1:3423784/64‑1 (MQ=255)
                  aGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTc                                              <  1:2304032/71‑1 (MQ=255)
                  aGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTc                                              <  1:2308547/71‑1 (MQ=255)
                   ggTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACATAAACTACCTGGCGTCAGTcc                                             <  1:3012776/71‑1 (MQ=255)
                    gTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCt                                            >  1:2977425/1‑71 (MQ=255)
                     tATTCAAGGCAGGGTTAAGACGTAATATGCGAAAAGCAGCCATGGACAGAAACTATCTg                                                       <  1:1249339/59‑1 (MQ=255)
                     tATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTg                                                       <  1:915348/59‑1 (MQ=255)
                     tATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCtt                                           <  1:1822875/71‑1 (MQ=255)
                       ttCAAGGCAGTGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTg                                          <  1:1133043/70‑1 (MQ=255)
                       ttCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCt                                                        <  1:2825384/56‑1 (MQ=255)
                       ttCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCt                                                        <  1:2829873/56‑1 (MQ=255)
                       ttCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTg                                          <  1:3432835/70‑1 (MQ=255)
                       ttCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGa                                                                <  1:1511386/48‑1 (MQ=255)
                       ttCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGa                                                                <  1:3141568/48‑1 (MQ=255)
                       ttCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGc                                         <  1:2383706/71‑1 (MQ=255)
                       ttCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGc                                         <  1:1626290/71‑1 (MQ=255)
                        tCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTcgcg                                                    >  1:5689/1‑59 (MQ=255)
                         cAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCtt                                           >  1:1861789/1‑67 (MQ=255)
                         cAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAg                                                                 >  1:1485514/1‑45 (MQ=255)
                         cAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTTGCGTCAGTCCTTGCgg                                       >  1:2889748/1‑71 (MQ=255)
                         cAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCgg                                       >  1:3418360/1‑71 (MQ=255)
                         cAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACgg                                       >  1:2222874/1‑71 (MQ=255)
                         cAAGGCAGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTg                                                       >  1:225888/1‑54 (MQ=255)
                          aaGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGaaa                                                              >  1:555751/1‑47 (MQ=255)
                          aaGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGaaa                                                              >  1:1200662/1‑47 (MQ=255)
                          aaGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACggg                                      >  1:2057439/1‑71 (MQ=255)
                           aGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTAcc                                                         <  1:2795499/51‑1 (MQ=255)
                            ggCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTat                                                          >  1:488574/1‑48 (MQ=255)
                            ggCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTAc                                                          >  1:2202170/1‑49 (MQ=255)
                             gCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAg                                   >  1:1061903/1‑71 (MQ=255)
                              cAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGCGGGGAGc                                  <  1:1069311/71‑1 (MQ=255)
                                gggTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACGGGGAGCgg                                <  1:1027957/71‑1 (MQ=255)
                                 ggTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTa                                                           >  1:713361/1‑43 (MQ=255)
                                 ggTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGc                                         <  1:1948180/61‑1 (MQ=255)
                                 ggTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGt                                               <  1:2095755/55‑1 (MQ=255)
                                   ttAAGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTcgcg                                                    <  1:1745777/48‑1 (MQ=255)
                                   ttAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAGCGGGCt                             >  1:3535584/1‑71 (MQ=255)
                                   ttAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCt                             <  1:1114153/71‑1 (MQ=255)
                                   ttAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCt                             <  1:2674328/71‑1 (MQ=255)
                                   ttAAGACGAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGt                                                   <  1:3814885/48‑1 (MQ=39)
                                    tAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACGGGTAGCGGGCtt                            >  1:1790473/1‑71 (MQ=255)
                                     aaGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAGCGGGCttt                           >  1:1243038/1‑71 (MQ=255)
                                     aaGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACGGGGAGCGGGCtt                            <  1:3607905/70‑1 (MQ=255)
                                      aGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAGCGGGCtt                            <  1:2100372/69‑1 (MQ=255)
                                      aGACGTAAGATGCGAAAAGCAGCAATGGACAGAAACTACCTGGCGTCAGTcc                                             >  1:1627727/1‑52 (MQ=255)
                                       gACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCg                         >  1:114188/1‑71 (MQ=255)
                                       gACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCg                         >  1:1700406/1‑71 (MQ=255)
                                       gACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCt                                                        <  1:887272/40‑1 (MQ=255)
                                       gACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTc                                              <  1:371921/50‑1 (MQ=255)
                                       gACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGa                                    <  1:2466753/60‑1 (MQ=255)
                                         cGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTgg                                                      >  1:2171587/1‑40 (MQ=255)
                                         cGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCttac                                         <  1:3455665/53‑3 (MQ=255)
                                          gTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGc                                                     <  1:1968357/40‑1 (MQ=255)
                                          gTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCtt                                           >  1:3102107/1‑50 (MQ=255)
                                          gTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCtt                                           >  1:2046771/1‑50 (MQ=255)
                                           tAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGCGGGGAGc                                  >  1:560724/1‑58 (MQ=255)
                                            aaGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGt                                                   >  1:1415825/1‑40 (MQ=255)
                                            aaGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGTGGAGCGGGCTTTCgtag                      >  1:3456885/1‑69 (MQ=255)
                                            aaGATCCCAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCATACgtagt                     >  1:2704951/1‑70 (MQ=255)
                                             aGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTaa                   >  1:2145273/1‑71 (MQ=255)
                                             aGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTaa                   <  1:1763557/71‑1 (MQ=255)
                                              gATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGCGGGGAGCGGGCTTTc                          >  1:3837798/1‑63 (MQ=255)
                                               aTGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGt                                               >  1:689513/1‑41 (MQ=255)
                                                tGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGt                                                   >  1:2997168/1‑36 (MQ=255)
                                                tGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTcc                                             <  1:1156780/42‑1 (MQ=255)
                                                 gCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTc                          >  1:2038019/1‑60 (MQ=255)
                                                  cGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCtt                                           >  1:264636/1‑42 (MQ=255)
                                                  cGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACGGGGAGCGGGCTTTCgtaatc                    <  1:3786105/65‑4 (MQ=255)
                                                   gAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATAccc              <  1:2576615/70‑1 (MQ=255)
                                                      aaGCAGCCATGGACAGAAACTACCTGGCGTCAGTcc                                             >  1:1308844/1‑36 (MQ=255)
                                                       ggcagcTATGGTTAGAAACTACCTGACGGCAGTCCTTGCGGGGAGCAGGCTTTCgtaa                      <  1:2911956/57‑2 (MQ=255)
                                                        gcagcCATGGATAGAAACTATCTCGCGTCAGTCCtt                                           <  1:2477675/36‑1 (MQ=37)
                                                        gcagcCATGGATAGAAACTATCTCGCGTCAGTCCTTGCggg                                      <  1:3587545/41‑1 (MQ=255)
                                                        gcagcCATGGACAGAAACTACCTGGCGTCAGTCCTTGCg                                        <  1:3156852/39‑1 (MQ=255)
                                                        gcagcCATGGACAGAAACTACCTGGCGTCAGTCCTTGCg                                        <  1:167262/39‑1 (MQ=255)
                                                         cagcCATGAACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGActct       <  1:338180/71‑1 (MQ=255)
                                                          agcCATGGATAGAAACTATCTCGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgtaat                     <  1:2745131/56‑3 (MQ=255)
                                                              aTGGACAGAAACTACCTGGCGTAAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCg  >  1:2602766/1‑71 (MQ=255)
                                                               tGGATAGAAACTATCTCGCGTCAGTCCTTGCGGGGAg                                   <  1:1398882/37‑1 (MQ=37)
                                                               tGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCg                         <  1:1676202/47‑1 (MQ=255)
                                                               tGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTcccc   >  1:904705/1‑69 (MQ=255)
                                                                 gACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCgg                                >  1:1173149/1‑38 (MQ=255)
                                                                   |                                                                 
ATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCG  >  minE/218090‑218222

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: