Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 1,196,798 (C)6→7 coding (354/1911 nt) yoaA ← conserved hypothetical protein with nucleoside triphosphate hydrolase domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,196,7921.C97.6% 137.4 / ‑3.4 41coding (360/1911 nt)yoaAconserved hypothetical protein with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (14/26);  total (15/26)
Fisher's exact test for biased strand distribution p-value = 3.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

CGACTGTTTGATTAGACCAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTTTCAGCA  >  minE/1196728‑1196856
                                                                 |                                                                
cGACTGTTTGATTAGACCAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGccccc                                                             >  1:1133897/1‑70 (MQ=255)
        tGATTAGACCAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAg                                                       >  1:280069/1‑69 (MQ=255)
         gATTAGACCAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAgcg                                                     <  1:2894350/70‑1 (MQ=255)
          aTTAGACCAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAgcgc                                                    <  1:2497679/70‑1 (MQ=255)
          aTTAGACCAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAgcgc                                                    <  1:246532/70‑1 (MQ=255)
          aTTAGACCAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGcc                                                   >  1:860083/1‑71 (MQ=255)
            tAGACCAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCct                                                  >  1:691811/1‑70 (MQ=255)
                ccAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGt                                             <  1:635875/71‑1 (MQ=255)
                ccAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGt                                             <  1:2101789/71‑1 (MQ=255)
                 cAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGtt                                            >  1:995008/1‑71 (MQ=255)
                          agcagGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGtt                                   <  1:2228102/71‑1 (MQ=255)
                           gcagGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTc                                  <  1:1702667/71‑1 (MQ=255)
                            cagGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTc                                  >  1:742330/1‑70 (MQ=255)
                            cagGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCg                                 <  1:622037/71‑1 (MQ=255)
                            cagGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCg                                 <  1:510941/71‑1 (MQ=255)
                            cagGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCg                                 <  1:1408542/71‑1 (MQ=255)
                             agGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGa                                <  1:3018187/71‑1 (MQ=255)
                               gATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGtt                                   <  1:835455/66‑1 (MQ=255)
                               gATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAgg                              <  1:1818669/71‑1 (MQ=255)
                                 tCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGATGTTCGAGACGTTCGAgg                              >  1:2896129/1‑69 (MQ=255)
                                  cacaTCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGa                                >  1:1988026/1‑66 (MQ=255)
                                   acaTCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCaga                          <  1:1700633/71‑1 (MQ=255)
                                   acaTCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCag                           <  1:1085535/70‑1 (MQ=255)
                                   acaTCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCag                           <  1:654933/70‑1 (MQ=255)
                                    caTCGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCaga                          <  1:569806/70‑1 (MQ=255)
                                       cGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGtt                                   >  1:2203199/1‑58 (MQ=255)
                                       cGCTTAAGATTTGTACCGGCAGATCGCCCCCCCGCAAGCGCCTGCTGTTCGAGACGtt                                   >  1:936606/1‑58 (MQ=255)
                                        gCTTAAGATTTGTACCGGCAGATCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGtt                   >  1:843811/1‑71 (MQ=255)
                                         cTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCg                                          <  1:1085958/49‑1 (MQ=38)
                                         cTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCg                                          <  1:392662/49‑1 (MQ=38)
                                         cTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCg                                          <  1:1730588/49‑1 (MQ=38)
                                         cTTAAGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGt                    <  1:933780/71‑1 (MQ=255)
                                            aaGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCGGAGGTAGtt                   <  1:804436/69‑1 (MQ=255)
                                            aaGATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGttt                  >  1:2477214/1‑70 (MQ=255)
                                              gATTTGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAg               >  1:645321/1‑71 (MQ=255)
                                                  tGTACCGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCaga                          >  1:1695622/1‑56 (MQ=255)
                                                      ccGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCtt        <  1:1260691/70‑1 (MQ=255)
                                                      ccGGCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCtt        <  1:1641591/70‑1 (MQ=255)
                                                        ggCAGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCagag                         >  1:1729858/1‑51 (MQ=39)
                                                           aGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTTTcagca  <  1:341548/71‑1 (MQ=255)
                                                           aGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTTTcagca  <  1:17428/71‑1 (MQ=255)
                                                           aGATCGCCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTTTcag    >  1:1316456/1‑69 (MQ=255)
                                                                 |                                                                
CGACTGTTTGATTAGACCAGGAGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTTTCAGCA  >  minE/1196728‑1196856

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: