Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 218,182 G→A 12.7% A57T (GCG→ACG)  yafF → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE218,1820GA12.7% 166.8 / 12.2 63A57T (GCG→ACG) yafFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/28);  new base A (4/4);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

TTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACA  >  minE/218113‑218250
                                                                     |                                                                    
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)
 tCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACg                                                                    >  1:2276626/1‑71 (MQ=255)
  cAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTTGCGTCAGTCCTTGCgg                                                                   >  1:2889748/1‑71 (MQ=255)
  cAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCgg                                                                   >  1:3418360/1‑71 (MQ=255)
  cAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACgg                                                                   >  1:2222874/1‑71 (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)
            ttAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAGCGGGCt                                                         >  1:3535584/1‑71 (MQ=255)
            ttAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCt                                                         <  1:2674328/71‑1 (MQ=255)
            ttAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCt                                                         <  1:1114153/71‑1 (MQ=255)
             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)
                gACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCg                                                     >  1:1700406/1‑71 (MQ=255)
                gACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCg                                                     >  1:114188/1‑71 (MQ=255)
                gACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGa                                                                <  1:2466753/60‑1 (MQ=255)
                  cGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgt                                                    >  1:3072249/1‑70 (MQ=255)
                   gTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgtaa                                                  >  1:1832434/1‑70 (MQ=255)
                    tAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGCGGGGAGc                                                              >  1:560724/1‑58 (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)
                          gCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTc                                                      >  1:2038019/1‑60 (MQ=255)
                           cGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACGGGGAGCGGGCTTTCgtaatc                                                <  1:3786105/65‑4 (MQ=255)
                            gAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATAccc                                          <  1:2576615/70‑1 (MQ=255)
                              aaaGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgta                                                   >  1:2859295/1‑59 (MQ=255)
                                ggcagcTATGGTTAGAAACTACCTGACGGCAGTCCTTGCGGGGAGCAGGCTTTCgtaa                                                  <  1:2911956/57‑2 (MQ=255)
                                 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)
                                                             cAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGAct                                     >  1:3807805/1‑42 (MQ=255)
                                                              aGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGActct                                   <  1:1072645/43‑1 (MQ=255)
                                                                 ccTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGAct                                     >  1:2784837/1‑38 (MQ=255)
                                                                  cTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCgc                                       >  1:3336290/1‑34 (MQ=255)
                                                                  cTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTcaccac   <  1:3050594/71‑1 (MQ=255)
                                                                  cTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTcacca    <  1:676791/70‑1 (MQ=255)
                                                                    tGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCg                              >  1:3123493/1‑42 (MQ=255)
                                                                     gCGGGGAGCGGGCTTTCGTAGTAATACCCCGActc                                    >  1:287213/1‑35 (MQ=255)
                                                                     gCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTcccc                               >  1:3231989/1‑40 (MQ=255)
                                                                     gCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTcc                           >  1:3090294/1‑44 (MQ=255)
                                                                     gCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACcaca  <  1:1624581/69‑1 (MQ=255)
                                                                     |                                                                    
TTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACA  >  minE/218113‑218250

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: