Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,751,234 G→T 14.1% intergenic (‑20/‑187) yfiF ← / → trxC predicted methyltransferase/thioredoxin 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,751,2340GT14.1% 127.3 / 5.1 42intergenic (‑20/‑187)yfiF/trxCpredicted methyltransferase/thioredoxin 2
Reads supporting (aligned to +/- strand):  ref base G (15/21);  new base T (2/4);  total (17/25)
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.92e-01

GACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAGCGCGCAGATTACCCGAAGTCCGCGCGGTTCGCCATGATTTCGTACCAAAGCCTGCGAC  >  minE/1751164‑1751299
                                                                      |                                                                 
gACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTg                                                                   <  1:1723098/71‑1 (MQ=255)
gACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTg                                                                   <  1:2946305/71‑1 (MQ=255)
    tacatCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTGcgcg                                                               <  1:1093848/71‑1 (MQ=255)
      catCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTGCGCGaa                                                             <  1:1592095/71‑1 (MQ=255)
      catCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTGCGCGaa                                                             <  1:3953822/71‑1 (MQ=255)
      catCACTTTGACCTTGCCACTTTTACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTGCGCGa                                                              <  1:2649637/70‑1 (MQ=255)
       atCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAg                                                            >  1:846073/1‑71 (MQ=255)
       atCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAg                                                            <  1:2525412/71‑1 (MQ=255)
          aCTTTGCCCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTGCGCGaa                                                             >  1:2392036/1‑67 (MQ=255)
              tGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAgcgcg                                                        >  1:412823/1‑68 (MQ=255)
              tGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAgcgcg                                                        >  1:3696680/1‑68 (MQ=255)
               gACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAGCGCGCAg                                                     <  1:3742288/70‑1 (MQ=255)
                   ttGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAGCGCGCa                                                      <  1:3742242/65‑1 (MQ=255)
                   ttGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAGCGCGCa                                                      <  1:1776249/65‑1 (MQ=255)
                   ttGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAGCGCGCAGATTAc                                                >  1:3001297/1‑71 (MQ=255)
                    tGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAGCGCGCAGATTAcc                                               >  1:991126/1‑71 (MQ=255)
                     gCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAGCGCGCAGATTAccc                                              >  1:404084/1‑71 (MQ=255)
                       cGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAGCGCGCAGATTACCCg                                             <  1:1853147/70‑1 (MQ=255)
                           tttACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAGCGCGCAGATTACCCGAAGTc                                        >  1:1715031/1‑71 (MQ=255)
                            ttACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTTcgcg                                                               <  1:945922/47‑1 (MQ=255)
                            ttACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTTcgcg                                                               <  1:3701084/47‑1 (MQ=255)
                            ttACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTTcgcg                                                               <  1:1746722/47‑1 (MQ=255)
                            ttACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTTcgcg                                                               <  1:380911/47‑1 (MQ=255)
                            ttACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTGcgcg                                                               <  1:191080/47‑1 (MQ=255)
                            ttACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAGCGCGCAGATTACCCGAAGTc                                        <  1:3156055/70‑1 (MQ=255)
                              aCCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTTcgcgca                                                             >  1:2056772/1‑45 (MQ=255)
                              aCCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTTcgcgca                                                             >  1:880840/1‑45 (MQ=255)
                                cTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAGCGCGCAGATTACCCGAAGTCcgcgc                                   <  1:1452699/71‑1 (MQ=255)
                                  ttcatttcatCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAGCGCGCAGATTACCCGAAGTCCGCGCgg                                 >  1:2225610/1‑71 (MQ=255)
                                   tcatttcatCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAGCGCGCAGATTACCCGAAGTCCGCGCGGt                                >  1:222325/1‑71 (MQ=255)
                                     atttcatCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAGCGCGCAGATTACCCGaa                                           >  1:1671908/1‑58 (MQ=255)
                                     atttcatCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAGCGCGCAGATTACCCGAAGTCCGCGCGGtt                               <  1:4070977/70‑1 (MQ=255)
                                     atttcatCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAGCGCGCAGATTACCCGAAGTCCGCGCGGtt                               <  1:420168/70‑1 (MQ=255)
                                      tttcatCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAGCGCGCAGATTACCCGAAGTCcgcgc                                   <  1:717391/65‑1 (MQ=255)
                                       ttcatCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAGCGCGCAGATTACCCGAAGTCCGCGCGGTTCg                             <  1:2979130/70‑1 (MQ=255)
                                        tcatCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAGCGCGCAGATTAc                                                >  1:1991304/1‑50 (MQ=255)
                                         catCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAGCGCGCAGATTACCCGAAGTCCGCGCGGTTCGcc                           <  1:1017718/70‑1 (MQ=255)
                                         catCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAGCGCGCAGATTACCCGAAGTCCGCGCGGTTCGCCa                          >  1:2818948/1‑71 (MQ=255)
                                          atCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAGCGCGCAGATTACCCGAAGTCcg                                      <  1:811607/58‑1 (MQ=255)
                                             gTTCATGCTTTTCTCCACCAGCGCTGCGCGAAGCGCGCAGATTACCCGAAGTCCGCGCGGTTCGCCATGa                       <  1:613259/70‑1 (MQ=255)
                                             gTTCATGCTTTTCTCCACCAGCGCTGAGCGAAGCGCGCAGATTACCCGAAGTCCGCGCGGTTCGCCATGa                       <  1:3769050/70‑1 (MQ=255)
                                                 aTGCTTTTCTCCACCAGCGCTGCGCGAAGCGCGCAGATTACCCGAAGTCCGCGCGGTTCGCCATGATTTCg                  <  1:3861313/71‑1 (MQ=255)
                                                              ccaGCGCTGCGCGAAGCGCGCAGATTACCCGAAGTCCGCGCGGTTCGCCATGATTTCGTACCAAAGCCTGc     >  1:1541140/1‑71 (MQ=255)
                                                                  cgcTGCGCGAAGCGCGCAGATTACCCGAAGTCCGCGCGGTTCGCCATGATTTCGTACCAAAGCCTGCGAc  >  1:2575885/1‑70 (MQ=255)
                                                                      |                                                                 
GACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAGCGCGCAGATTACCCGAAGTCCGCGCGGTTCGCCATGATTTCGTACCAAAGCCTGCGAC  >  minE/1751164‑1751299

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: