Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 218,169 G→C 15.2% L52L (CTG→CTC yafF → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE218,1690GC15.2% 206.3 / 23.9 86L52L (CTG→CTCyafFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (34/38);  new base C (4/9);  total (39/47)
Fisher's exact test for biased strand distribution p-value = 3.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

GAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAAC  >  minE/218100‑218236
                                                                     |                                                                   
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)
             ttCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTg                                                        <  1:3432835/70‑1 (MQ=255)
             ttCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGc                                                       <  1:2383706/71‑1 (MQ=255)
             ttCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGc                                                       <  1:1626290/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)
                   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: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)
                            aGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCa                                                               >  1:2054281/1‑48 (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)
                             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:2046771/1‑50 (MQ=255)
                                gTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCtt                                                         >  1:3102107/1‑50 (MQ=255)
                                gTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCtt                                                         >  1:3712929/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)
                                   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)
                                                  ccATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgtagta                                  <  1:917290/55‑1 (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)
                                                                   cTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACaa   <  1:3752986/69‑1 (MQ=255)
                                                                     |                                                                   
GAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAAC  >  minE/218100‑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: