Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,325,522 T→G 21.2% K71T (AAA→ACA)  wcaL ← predicted glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,325,5220TG21.2% 34.9 / 21.7 33K71T (AAA→ACA) wcaLpredicted glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base T (15/11);  new base G (2/5);  total (17/16)
Fisher's exact test for biased strand distribution p-value = 2.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

AGGTTGAGCGCCTGCCAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCAT  >  minE/1325453‑1325590
                                                                     |                                                                    
aGGTTGAGCGCCTGCCAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGtt                                                                      <  1:1899144/70‑1 (MQ=255)
      aGCGCCTGCCATGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCa                                                                <  1:178423/70‑1 (MQ=255)
           cTGCCAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTGTCGCCa                                                                <  1:1440669/65‑1 (MQ=255)
            tGCGAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTGTCGCCa                                                                <  1:313614/64‑1 (MQ=255)
                aGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTGTCGCCACTTTGCCCGTa                                                     <  1:914505/71‑1 (MQ=255)
                aGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTGTCGCCACTTTGCCCGTa                                                     <  1:1737922/71‑1 (MQ=255)
                   tATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGt                                                  >  1:307799/1‑71 (MQ=255)
                   tATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTGTCGCCACTTTGCCCGTAGGt                                                  >  1:488094/1‑71 (MQ=255)
                     tttttACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGccc                                                        >  1:1524933/1‑63 (MQ=255)
                     tttttACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGtt                                                 <  1:348701/70‑1 (MQ=255)
                     tttttACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTc                                                <  1:2119053/71‑1 (MQ=255)
                     tttttACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTc                                                <  1:1587248/71‑1 (MQ=255)
                     tttttACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTc                                                <  1:1522916/71‑1 (MQ=255)
                       tttACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAg                                                    >  1:586379/1‑65 (MQ=255)
                       tttACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCg                                               <  1:11566/70‑1 (MQ=255)
                       tttACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGt                                              >  1:1809817/1‑71 (MQ=255)
                       tttACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTGTCGCCACTTTGCCCGTAGGTTCGt                                              >  1:1199108/1‑71 (MQ=255)
                               gAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAAc                                      >  1:90659/1‑71 (MQ=255)
                                   gCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGt                                              <  1:390323/59‑1 (MQ=255)
                                     cgcgCAAGGTCTGGCTGGCTCGGTGGCTCAGTTTCGCCACTTTGc                                                          >  1:1008760/1‑45 (MQ=255)
                                     cgcgCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTcc                                            >  1:153726/1‑59 (MQ=255)
                                     cgcgCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTcc                                            >  1:698547/1‑59 (MQ=255)
                                     cgcgCAAGGTCTGGCTGGCTCGGTCGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCgg                                >  1:332743/1‑71 (MQ=255)
                                           aGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTg                                                           >  1:1652531/1‑38 (MQ=255)
                                           aGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGtt                             >  1:227277/1‑68 (MQ=255)
                                             gTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCg                                 <  1:1163386/62‑1 (MQ=255)
                                               ctggctggCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGgcag                      >  1:750795/1‑71 (MQ=255)
                                                    tggctCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTg                          <  1:1890501/62‑1 (MQ=255)
                                                          ggTGGCGCAGTGTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTa            <  1:266286/70‑1 (MQ=255)
                                                              gcgcAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCg       >  1:1887256/1‑71 (MQ=255)
                                                              gcgcAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCg       >  1:407504/1‑71 (MQ=255)
                                                               cgcAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGt      >  1:851000/1‑71 (MQ=255)
                                                                 cAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTcc    <  1:582727/71‑1 (MQ=255)
                                                                   gTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCAt  <  1:1256637/71‑1 (MQ=255)
                                                                     |                                                                    
AGGTTGAGCGCCTGCCAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCAT  >  minE/1325453‑1325590

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: