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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE218,1820GA 5.5% 243.9 / 2.8 73A57T (GCG→ACG) yafFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (34/35);  new base A (2/2);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

ATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACA  >  minE/218112‑218234
                                                                      |                                                    
aTTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTg                                                      >  1:3813829/1‑71 (MQ=255)
aTTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTg                                                      >  1:4119509/1‑71 (MQ=255)
aTTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTg                                                      >  1:807176/1‑71 (MQ=255)
aTTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTg                                                      >  1:1203054/1‑71 (MQ=255)
   cAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCg                                                    >  1:187658/1‑70 (MQ=255)
   cAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCTTTGCgg                                                   >  1:675992/1‑71 (MQ=255)
      ggCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCggg                                                  >  1:1901107/1‑69 (MQ=255)
        cAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGc                                              <  1:833421/71‑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)
             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)
                aGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCggg                                                  <  1:1726861/59‑1 (MQ=255)
                aGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCgg                                                   >  1:377940/1‑58 (MQ=255)
                aGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCgg                                                   >  1:1368511/1‑58 (MQ=255)
                 gACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCg                                     <  1:3654419/71‑1 (MQ=255)
                   cGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGTAGCGGGCTTTCgta                                   <  1:3869391/71‑1 (MQ=255)
                   cGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgt                                    <  1:2551760/70‑1 (MQ=255)
                   cGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgt                                    <  1:1015421/70‑1 (MQ=255)
                   cGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCg                                     >  1:1410887/1‑69 (MQ=255)
                    gTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCg                                                    <  1:3732799/53‑1 (MQ=255)
                     tAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGa                                                >  1:753863/1‑56 (MQ=255)
                     tAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgtagt                                 <  1:1559773/71‑1 (MQ=255)
                     tAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgtagt                                 <  1:2904104/71‑1 (MQ=255)
                     tAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgtagt                                 <  1:308399/71‑1 (MQ=255)
                     tAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgtag                                  >  1:3508145/1‑70 (MQ=255)
                      aaGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgtagt                                 >  1:851487/1‑70 (MQ=255)
                      aaGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCg                                     >  1:665057/1‑66 (MQ=255)
                       aGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAGc                                              <  1:849033/56‑1 (MQ=255)
                       aGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGGGGGAGCGGGCTTTCGTAGTaa                               >  1:457266/1‑70 (MQ=255)
                       aGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACgg                                                   >  1:3920696/1‑51 (MQ=255)
                        gATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCggg                                                  >  1:434918/1‑51 (MQ=255)
                        gATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCg                                     <  1:3550700/64‑1 (MQ=255)
                           gCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgta                                   >  1:3530399/1‑63 (MQ=255)
                             gAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACGGGGAGCGGGCTTTCgtaatc                                <  1:3379625/64‑4 (MQ=255)
                              aaaaGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGc                                                     <  1:2527603/42‑1 (MQ=255)
                              aaaaGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGc                                                     <  1:2598712/42‑1 (MQ=255)
                              aaaaGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAg                                               <  1:772944/48‑1 (MQ=255)
                               aaaGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgt                                    >  1:2063674/1‑58 (MQ=255)
                                aaGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGCg                                                    >  1:2681255/1‑41 (MQ=255)
                                aaGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGa                       >  1:1611346/1‑70 (MQ=255)
                                   cagcCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgt                                    <  1:2399621/54‑1 (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:264221/1‑71 (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)
                                              aGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgtag                                  >  1:1670170/1‑45 (MQ=255)
                                              aGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTccc                >  1:4109688/1‑63 (MQ=255)
                                              aGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCtt         >  1:2025505/1‑70 (MQ=255)
                                              aGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTa        <  1:1098499/71‑1 (MQ=255)
                                               gAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATAc                            <  1:3293331/50‑1 (MQ=255)
                                               gAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTa        <  1:1166993/70‑1 (MQ=255)
                                                 aaCTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATAccc                          >  1:3564837/1‑50 (MQ=255)
                                                  aCTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTaa                               >  1:3507135/1‑44 (MQ=255)
                                                  aCTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTcccc               <  1:24954/60‑1 (MQ=255)
                                                  aCTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTc            >  1:2463527/1‑63 (MQ=255)
                                                  aCTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAaca    <  1:1809961/71‑1 (MQ=255)
                                                   cTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAacac   <  1:42574/71‑1 (MQ=255)
                                                      ccTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAAcaca  <  1:473885/69‑1 (MQ=255)
                                                           cGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGActc                    >  1:3313205/1‑46 (MQ=255)
                                                                gTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTaaa      >  1:2267528/1‑55 (MQ=255)
                                                                  ccTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTcc                 >  1:1547503/1‑42 (MQ=255)
                                                                   cTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGCCTCTCCCCGTCCTTAAaca    >  1:987514/1‑54 (MQ=255)
                                                                     tGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTccc                >  1:299251/1‑40 (MQ=255)
                                                                      |                                                    
ATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACA  >  minE/218112‑218234

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: