Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,572,711 T→C 25.9% intergenic (+1479/‑300) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,572,7110TC25.9% 270.7 / 105.3 166intergenic (+1479/‑300)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (70/53);  new base C (27/16);  total (97/69)
Fisher's exact test for biased strand distribution p-value = 5.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 188 total aligned reads displayed.

GGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAG  >  minE/2572641‑2572780
                                                                      |                                                                     
ggTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCacctacct                                                                       >  1:3298856/1‑71 (MQ=255)
  tGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTc                                                                     <  1:1357980/71‑1 (MQ=255)
    gCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTc                                                                     >  1:1142815/1‑69 (MQ=255)
    gCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTc                                                                     >  1:4964614/1‑69 (MQ=255)
    gCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCa                                                                   <  1:1422235/71‑1 (MQ=255)
      tAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGc                                                                 >  1:2094602/1‑71 (MQ=255)
        aTCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGccc                                                               <  1:5210802/71‑1 (MQ=255)
         tCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCa                                                              <  1:2969901/71‑1 (MQ=255)
           ccAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAAt                                                            >  1:1848093/1‑71 (MQ=255)
            cAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTc                                                                     <  1:1032664/61‑1 (MQ=255)
            cAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTc                                                                     <  1:3296740/61‑1 (MQ=255)
            cAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATc                                                           <  1:2416703/71‑1 (MQ=255)
             aaGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTc                                                                     <  1:917203/60‑1 (MQ=255)
             aaGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCg                                                          >  1:1465919/1‑71 (MQ=255)
             aaGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCg                                                          >  1:4834487/1‑71 (MQ=255)
               ggggATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCgaa                                                        <  1:724763/71‑3 (MQ=255)
               ggggATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCg                                                          >  1:3756246/1‑69 (MQ=255)
                  gATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGcc                                                     >  1:1086747/1‑71 (MQ=255)
                          ccAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCa                                                                   >  1:1602066/1‑49 (MQ=255)
                          ccAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCa                                                                   >  1:4658891/1‑49 (MQ=255)
                           cAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                            <  1:2017498/71‑1 (MQ=255)
                             gagaTCGGATTCAACGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGatat                                           <  1:3519381/67‑1 (MQ=255)
                             gAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCa                                                              <  1:2400665/51‑1 (MQ=255)
                             gAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtg                                               <  1:3999812/66‑1 (MQ=255)
                             gAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTtata                                          <  1:1731037/71‑1 (MQ=255)
                              aaaTCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTg                                                   >  1:2153154/1‑61 (MQ=255)
                                aTCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgt                                                <  1:4897259/62‑1 (MQ=255)
                                aTCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCg                                       <  1:1734648/71‑1 (MQ=255)
                                 ccGGATTCAGCGGATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAgtga                                              <  1:4494340/62‑2 (MQ=255)
                                    ggTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGcccc                                   <  1:2673676/71‑1 (MQ=255)
                                     attCAGCGGATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCgg                                                         <  1:5664344/47‑1 (MQ=255)
                                     attCAGCGGATACGCATACAGACTCCATCCCCCCTCCAGCCCTATCGGGTCCTGAGTGATATACCGCCCct                                  <  1:2205840/70‑1 (MQ=255)
                                     attCAACGGATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGat                                                       >  1:320179/2‑49 (MQ=255)
                                      ttCAGCGGATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATACCTCCCctg                                 >  1:4336450/1‑71 (MQ=255)
                                      ttCAACGTATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                            >  1:2577749/1‑60 (MQ=255)
                                      ttCAACGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGccccag                                 >  1:1136220/1‑69 (MQ=255)
                                      ttCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCaa                                                             <  1:1049083/43‑1 (MQ=255)
                                      ttCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCaa                                                             <  1:2808058/43‑1 (MQ=255)
                                      ttCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                            >  1:2413257/1‑60 (MQ=255)
                                       tCAACGGATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATACCGCCCctg                                 >  1:27694/1‑70 (MQ=255)
                                       tCAACGGATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATACCGCCCct                                  >  1:4019555/1‑69 (MQ=255)
                                        cAACGGATACTGATACATATTCCATCCACCTCTCAACCCAATCGGGTCCGGGGTGATATACCGCCCctgc                                <  1:5249456/70‑1 (MQ=255)
                                        cAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAAt                                                            >  1:4909787/1‑42 (MQ=255)
                                        cAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgct                               >  1:3994771/1‑71 (MQ=255)
                                         aaCGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTg                                                   >  1:3002682/1‑50 (MQ=255)
                                         aaCGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAAt                                                            >  1:4998908/1‑41 (MQ=255)
                                         aaCGGATACTGATAAAAATTCCACCTACCTTCCAg                                                                  >  1:1215970/1‑35 (MQ=255)
                                         aaCGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCg                                                          >  1:671114/1‑43 (MQ=255)
                                          aCGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCcagcag                              >  1:3815513/1‑68 (MQ=255)
                                          aCGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCcagcag                              >  1:554642/1‑68 (MQ=255)
                                          aCGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATAc                                         >  1:610003/1‑59 (MQ=255)
                                           cGGATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATACCGCCCCTGCAGCg                            <  1:3995756/71‑1 (MQ=255)
                                            ggATACTGATAAAAATTCCATCCCCCTTCAGCCCAATCGGATCCTGAGTGATATACCTCCCCTGCAGCggg                          >  1:2878568/1‑70 (MQ=255)
                                            ggATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATACCGCCCCTGCAGCgg                           >  1:5171066/1‑71 (MQ=255)
                                            ggATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                              <  1:1911403/52‑1 (MQ=255)
                                            ggATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATAc                                         >  1:3432235/1‑57 (MQ=255)
                                            ggATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATAc                                         >  1:5633958/1‑57 (MQ=255)
                                             gATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCCAGCAGCgg                           >  1:5212955/1‑70 (MQ=255)
                                             gATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAgtgat                                             >  1:1654182/1‑50 (MQ=255)
                                             gATACTGATAAAAATTCCACCTACCTTCCAGCCCAATc                                                           >  1:428937/1‑38 (MQ=255)
                                              aTACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGa                                                        >  1:3232991/1‑40 (MQ=255)
                                              aTACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGt                          >  1:2729850/1‑70 (MQ=255)
                                               tACTGATAAAAATTCCATCCCCACTTCAGCCCAATCGGATCCTGAg                                                 >  1:4576623/1‑46 (MQ=255)
                                               tACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGa                                                        <  1:1103184/39‑1 (MQ=255)
                                                 cTGATAAAAATTCCATCCCCCCTTCAGCCCAATCg                                                          <  1:2042505/35‑1 (MQ=16)
                                                 cTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                            >  1:76943/1‑49 (MQ=255)
                                                 cTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGccc                                    >  1:4294827/1‑57 (MQ=255)
                                                  tGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCt                                                    <  1:3630668/40‑1 (MQ=255)
                                                  tGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCt                                                    <  1:2519766/40‑1 (MQ=255)
                                                   gATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCCAGCAGCgg                           >  1:4979797/1‑64 (MQ=255)
                                                   gATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCt                                                    <  1:1263718/39‑1 (MQ=16)
                                                   gATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCt                                                    <  1:5332975/39‑1 (MQ=16)
                                                    aTAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGAtatat                                         >  1:5594015/1‑48 (MQ=255)
                                                    aTAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGAtata                                          >  1:3729263/1‑48 (MQ=255)
                                                    aTAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGa                                                  >  1:4320881/1‑40 (MQ=255)
                                                    aTAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGAta                                            >  1:2750914/1‑46 (MQ=255)
                                                    aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATAc                                         >  1:4383217/1‑49 (MQ=255)
                                                     tAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCCagcagc                             <  1:1807372/60‑1 (MQ=255)
                                                     tAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCtgctgc                             <  1:2339800/60‑1 (MQ=255)
                                                      aaaaaTTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTAc                 >  1:1167399/1‑71 (MQ=255)
                                                       aaaaTTCCATCACCCCTTCAGCCCAATCGGATCCTGAGTGAtata                                          >  1:5468939/1‑45 (MQ=255)
                                                       aaaaTTCCACCTACCTTCCAGCCCAATCGGAGCATGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTAcc                >  1:4425874/1‑71 (MQ=255)
                                                         aaTTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGAta                                            >  1:1699164/1‑41 (MQ=255)
                                                         aaTTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATAcctc                                      <  1:4132582/47‑3 (MQ=255)
                                                           ttCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGat             >  1:3171790/1‑70 (MQ=255)
                                                             ccACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataa           <  1:5065269/70‑1 (MQ=255)
                                                             ccACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataa           >  1:4211728/1‑70 (MQ=255)
                                                              cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataat          <  1:967637/70‑1 (MQ=255)
                                                                cctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGt       <  1:2116581/71‑1 (MQ=255)
                                                                 ctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCgg                           >  1:2143794/1‑50 (MQ=255)
                                                                 ctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGt       >  1:2631853/1‑70 (MQ=255)
                                                                  tacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgc                                <  1:158107/44‑1 (MQ=255)
                                                                  tacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgc                                <  1:5145287/44‑1 (MQ=255)
                                                                  tacctTCCAGCCCAATCGGAGCCTGAGAGTTATACCGCCCCAGCTGCGGTTAAAGTTACCGAAAATAGAta     >  1:2432791/1‑71 (MQ=255)
                                                                   acctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta     >  1:2439551/1‑70 (MQ=255)
                                                                    cctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTTCGGTTAATGTTACCGATAATAGTta     <  1:3063471/69‑1 (MQ=255)
                                                                    cctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta     <  1:4763221/69‑1 (MQ=255)
                                                                    cctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta     <  1:5310379/69‑1 (MQ=255)
                                                                     ctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAg  >  1:1943673/1‑71 (MQ=255)
                                                                      ttCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTa                  >  1:815782/1‑54 (MQ=255)
                                                                      |                                                                     
GGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAG  >  minE/2572641‑2572780

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: