Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 218,030 A→G 73.7% Y6C (TAC→TGC)  yafF → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE218,0300AG73.7% 326.1 / 126.4 201Y6C (TAC→TGC) yafFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/25);  new base G (98/50);  total (126/75)
Fisher's exact test for biased strand distribution p-value = 9.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 201 total aligned reads displayed.

CCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTT  >  minE/217962‑218096
                                                                    |                                                                  
ccGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTAc                                                                   >  1:3945071/1‑70 (MQ=255)
 cGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCa                                                                  >  1:3648158/1‑70 (MQ=255)
   aaaCCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACaaaa                                                               >  1:1000416/1‑71 (MQ=255)
     atcACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAAt                                                              >  1:4553976/3‑70 (MQ=255)
     aCCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCaaaa                                                               >  1:716916/1‑69 (MQ=39)
     aCCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAAt                                                              >  1:4678397/1‑70 (MQ=255)
                gTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGaaa                                                  >  1:5445949/1‑71 (MQ=255)
                 tGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGaaat                                                 >  1:3074972/1‑70 (MQ=255)
                 tGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGaaat                                                 >  1:4662215/1‑70 (MQ=255)
                 tGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGa                                                    <  1:4260102/68‑1 (MQ=255)
                   gagaATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGaaatg                                                >  1:3893140/1‑68 (MQ=255)
                   gagaATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGaaatg                                                >  1:5649601/1‑68 (MQ=255)
                   gagaATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACg                                                >  1:4017662/1‑70 (MQ=255)
                      aaTAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACgcag                                             >  1:311949/1‑70 (MQ=255)
                        tAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATaa                                                            >  1:1531927/1‑53 (MQ=255)
                        tAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATaa                                                            >  1:3849266/1‑53 (MQ=255)
                          aGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAAt                                                              >  1:3687921/1‑49 (MQ=38)
                          aGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATaag                                                           >  1:5621552/1‑52 (MQ=255)
                          aGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATaa                                                            >  1:2125525/1‑51 (MQ=255)
                          aGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGaaatg                                                >  1:103484/1‑61 (MQ=255)
                          aGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGaaatg                                                >  1:1722551/1‑61 (MQ=255)
                          aGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGaa                                        >  1:183671/1‑71 (MQ=255)
                          aGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGaa                                        >  1:3431688/1‑71 (MQ=255)
                             tGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAatta                                     >  1:702751/1‑71 (MQ=255)
                             agCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAatt                                      >  1:2921597/2‑70 (MQ=255)
                                actggCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTAttt                                  >  1:2476369/4‑71 (MQ=255)
                                 ctggCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGTAAATGCAGCAGAATTAtttt                                 >  1:379680/3‑71 (MQ=255)
                                 ctggCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTAtttt                                 >  1:2305285/3‑71 (MQ=255)
                                 caGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt                                 <  1:261775/71‑1 (MQ=255)
                                 caGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt                                 <  1:4452388/71‑1 (MQ=255)
                                 caGGCGGCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAgg                                                     <  1:381832/51‑1 (MQ=255)
                                  tggCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGa                                                    <  1:1621040/50‑1 (MQ=255)
                                    gCGTCTGGACGTGGTAATGAATGAAGACGACTGCaaa                                                                >  1:4178891/1‑37 (MQ=255)
                                    gCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAg                              >  1:4531804/1‑71 (MQ=255)
                                    gCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAg                              >  1:3619865/1‑71 (MQ=255)
                                     cGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGaaa                                                  >  1:3002787/1‑50 (MQ=255)
                                     cGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAgg                             <  1:3748484/71‑1 (MQ=255)
                                       tCTGGACGTGGTAATGAATGAAGACGACTGCAAAATaag                                                           <  1:872150/39‑1 (MQ=255)
                                       tCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAatta                                     >  1:2613107/1‑61 (MQ=255)
                                        cTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGa                           >  1:3031516/1‑70 (MQ=255)
                                         tGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAagag                                                      >  1:4237334/1‑42 (MQ=255)
                                          ggACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGaa                                        <  1:3307548/55‑1 (MQ=255)
                                           gACGTGGTATTGAATGAAGACGACTGCAAAATAAGAAGAGGa                                                    <  1:1100447/42‑1 (MQ=38)
                                           gACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGa                           >  1:1587740/1‑67 (MQ=255)
                                           gACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGa                           >  1:4419242/1‑67 (MQ=255)
                                               tGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGca                    >  1:1214694/1‑70 (MQ=255)
                                                ggTAATGAATGAAGACGACTGCAAAATAAGAAGAGGaaa                                                  >  1:2072664/1‑39 (MQ=255)
                                                ggTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGcaca                  >  1:3224214/1‑71 (MQ=255)
                                                 gTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACAt                 >  1:5314869/1‑71 (MQ=255)
                                                 gTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCa                               >  1:5633730/1‑57 (MQ=255)
                                                 gTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACAt                 >  1:5524303/1‑71 (MQ=255)
                                                 gTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACAt                 >  1:2329927/1‑71 (MQ=255)
                                                  tattgattgaaGACGACTGCAAAATAAGTAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGcaca                  <  1:3340743/66‑1 (MQ=255)
                                                  tAATGAATGAAGACGACTGCAAAATAAGAAGAGGa                                                    >  1:1373102/1‑35 (MQ=255)
                                                   aatgaatgAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGcaca                  >  1:90337/1‑68 (MQ=255)
                                                   aatgaatgAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt                                 <  1:1145339/53‑1 (MQ=255)
                                                    atgaatgaAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCa                               >  1:1223351/1‑54 (MQ=39)
                                                    atgaatgaAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACgg                      >  1:2680293/1‑63 (MQ=255)
                                                    atgaatgaAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACAtt                >  1:4311200/1‑69 (MQ=255)
                                                    atgaatgaAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGc              >  1:4188913/1‑71 (MQ=255)
                                                     tgagtgaaGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGc              >  1:4580661/1‑70 (MQ=255)
                                                      gaatgaaGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTAtttt                                 >  1:2553238/1‑50 (MQ=255)
                                                       aatgaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCa                               >  1:580294/1‑51 (MQ=39)
                                                       aatgaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTAt           >  1:5424152/1‑71 (MQ=255)
                                                       aatgaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAattat                                    >  1:4847072/1‑46 (MQ=255)
                                                        ttgatgaCGACTGCAAAATAAGAAGAGGAAATgca                                              <  1:412162/34‑1 (MQ=21)
                                                        atgaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTAtttt                                 <  1:1299706/48‑1 (MQ=38)
                                                        atgaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTAttt                                  <  1:4750907/47‑1 (MQ=38)
                                                        atgaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAg                              >  1:479685/1‑51 (MQ=39)
                                                        atgaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCa                               >  1:5344283/1‑50 (MQ=255)
                                                         tgaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTAt           >  1:4403200/1‑69 (MQ=255)
                                                         tgaaGACGACCGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTAtt          >  1:5295270/1‑70 (MQ=255)
                                                          caaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt                                 <  1:3349400/45‑1 (MQ=255)
                                                           aaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTaa        >  1:4026744/1‑70 (MQ=255)
                                                           aaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTaa        >  1:4932436/1‑70 (MQ=255)
                                                            cgacgacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTaa        <  1:2274405/68‑1 (MQ=255)
                                                            acacgacTGCAAAATAAGAAGAGGAAACGCCGCAGAATTAtttt                                 <  1:5632639/42‑1 (MQ=38)
                                                            aGACGACTGCAAAATAAGAAGAGGAAATGcagcag                                          <  1:2019559/35‑1 (MQ=25)
                                                            aGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAatt                                      <  1:2441357/39‑1 (MQ=37)
                                                            aGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTc                                <  1:2969280/45‑1 (MQ=38)
                                                            aGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTAtt          <  1:5099815/67‑1 (MQ=255)
                                                            aGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAatt                                      <  1:2166000/39‑1 (MQ=37)
                                                            aGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTAtttt                                 <  1:2795275/44‑1 (MQ=38)
                                                            aGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAatt                                      <  1:3176870/39‑1 (MQ=255)
                                                            aGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAata      <  1:2389810/71‑1 (MQ=255)
                                                             gacgacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGa                           <  1:831704/49‑1 (MQ=38)
                                                              acgacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGa                           <  1:3765346/48‑1 (MQ=38)
                                                              acgacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGa                           <  1:2519896/48‑1 (MQ=38)
                                                              acgacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGca                    <  1:3524767/55‑1 (MQ=255)
                                                              acgacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATAtt    >  1:450137/1‑71 (MQ=255)
                                                              acgaaTGCAAAATAAGAAGAGGAAATGCAGCAGAATTAttt                                  <  1:5676934/41‑1 (MQ=25)
                                                               cgacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTAtttt                                 <  1:408413/41‑1 (MQ=38)
                                                               cgacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGa                           <  1:5400852/47‑1 (MQ=38)
                                                               cgacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGca                    <  1:2262752/54‑1 (MQ=39)
                                                                gacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTAttt                                  <  1:1357884/39‑1 (MQ=37)
                                                                gacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCa                               >  1:3481111/1‑42 (MQ=38)
                                                                gacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATAtttt  <  1:2239943/71‑1 (MQ=255)
                                                                gacTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt                                 <  1:3239038/40‑1 (MQ=255)
                                                                gacTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt                                 <  1:149770/40‑1 (MQ=255)
                                                                gacTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATAtttt  <  1:2872822/71‑1 (MQ=255)
                                                                    |                                                                  
CCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTT  >  minE/217962‑218096

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: