Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 218,166 C→T 36.1% Y51Y (TAC→TAT yafF → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE218,1660CT36.1% 93.6 / 91.3 97Y51Y (TAC→TATyafFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/32);  new base T (14/21);  total (44/53)
Fisher's exact test for biased strand distribution p-value = 5.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 101 total aligned reads displayed.

GACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAAC  >  minE/218097‑218236
                                                                     |                                                                      
gACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTAc                                                                        <  1:1429119/70‑1 (MQ=255)
 aCGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCt                                                                      >  1:642518/1‑71 (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)
                ttCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGc                                                       <  1:1626290/71‑1 (MQ=255)
                ttCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGc                                                       <  1:2383706/71‑1 (MQ=255)
                 tCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTcgcg                                                                  >  1:5689/1‑59 (MQ=255)
                 tCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACg                                                      >  1:2276626/1‑71 (MQ=255)
                  cAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCtt                                                         >  1:1861789/1‑67 (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)
                   aaGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACggg                                                    >  1:2057439/1‑71 (MQ=255)
                    aGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTAcc                                                                       <  1:2795499/51‑1 (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)
                          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:2674328/71‑1 (MQ=255)
                            ttAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCt                                           <  1:1114153/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)
                               aGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCa                                                               >  1:2054281/1‑48 (MQ=255)
                               aGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAGCGGGCtt                                          <  1:2100372/69‑1 (MQ=255)
                               aGACGTAAGATGCGAAAAGCAGCAATGGACAGAAACTACCTGGCGTCAGTcc                                                           >  1:1627727/1‑52 (MQ=255)
                                gACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCg                                       >  1:1700406/1‑71 (MQ=255)
                                gACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCg                                       >  1:114188/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)
                                  cGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgt                                      >  1:3072249/1‑70 (MQ=255)
                                  cGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTgg                                                                    >  1:2171587/1‑40 (MQ=255)
                                  cGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCttac                                                       <  1:3455665/53‑3 (MQ=255)
                                   gTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgtaa                                    >  1:1832434/1‑70 (MQ=255)
                                   gTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGc                                                                   <  1:1968357/40‑1 (MQ=255)
                                   gTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCtt                                                         >  1:3102107/1‑50 (MQ=255)
                                   gTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCtt                                                         >  1:2046771/1‑50 (MQ=255)
                                   gTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCtt                                                         >  1:3712929/1‑50 (MQ=255)
                                    tAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGCGGGGAGc                                                >  1:560724/1‑58 (MQ=255)
                                     aaGATGCGAAAAGCAGCCATGGACAGAAACTATCt                                                                      <  1:1557778/35‑1 (MQ=255)
                                     aaGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGt                                                                 >  1:1415825/1‑40 (MQ=255)
                                     aaGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGTGGAGCGGGCTTTCgtag                                    >  1:3456885/1‑69 (MQ=255)
                                     aaGATCCCAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCATACgtagt                                   >  1:2704951/1‑70 (MQ=255)
                                      aGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTg                                                        <  1:2649544/48‑1 (MQ=255)
                                      aGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTaa                                 >  1:2145273/1‑71 (MQ=255)
                                      aGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTaa                                 <  1:1763557/71‑1 (MQ=255)
                                       gATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGCGGGGAGCGGGCTTTc                                        >  1:3837798/1‑63 (MQ=255)
                                       gATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGc                                                       <  1:219190/48‑1 (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)
                                             aaaaGCAGCCATGGACAGAAACTACCTGGCGTCAGt                                                             <  1:1122428/36‑1 (MQ=255)
                                              aaaGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgta                                     >  1:2859295/1‑59 (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:167262/39‑1 (MQ=255)
                                                 gcagcCATGGACAGAAACTACCTGGCGTCAGTCCTTGCg                                                      <  1:3156852/39‑1 (MQ=255)
                                                  cagcCATGAACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGActct                     <  1:338180/71‑1 (MQ=255)
                                                   agcCATGGATAGAAACTATCTCGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgtaat                                   <  1:1388687/56‑3 (MQ=255)
                                                   agcCATGGATAGAAACTATCTCGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgtaat                                   <  1:2745131/56‑3 (MQ=255)
                                                       aTGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGCGCGGGCTTTCGTAGTAAt                                >  1:2516777/1‑55 (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)
                                                               aaaCTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTc              >  1:2100415/1‑65 (MQ=255)
                                                               aaaCTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTaa         <  1:2332396/70‑1 (MQ=255)
                                                                aaCTACCTGGCGTCAGTCCTTACGGGGAGCGGGCttt                                         <  1:1687960/37‑1 (MQ=255)
                                                                 aCTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGGCCTTAAaca      >  1:856117/1‑71 (MQ=255)
                                                                    aCCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATAcccc                           <  1:2816994/47‑1 (MQ=255)
                                                                    aCCTGGCGTCAGTCCTTGCGGGGAGCGGGATTTCGTAGTaa                                 >  1:1626310/1‑41 (MQ=255)
                                                                     ccTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAAc  <  1:2506768/71‑1 (MQ=255)
                                                                     |                                                                      
GACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAAC  >  minE/218097‑218236

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: