Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,682,335 T→C 11.6% K459K (AAA→AAG yhgF ← predicted transcriptional accessory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,682,3350TC11.6% 123.5 / 7.4 43K459K (AAA→AAGyhgFpredicted transcriptional accessory protein
Reads supporting (aligned to +/- strand):  ref base T (23/15);  new base C (2/3);  total (25/18)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GTTTGCGGGCCAGTTGCGTCTGGCTGACGTCATGCTGATACTGACCTACGCCGATAGATTTCGGATCGATTTTCACCAGCTCCGCCAGCGGATCCTGCAAACGGCGGGCGATAGACACCGCGCCACGCAGCGAAACGTC  >  minE/2682265‑2682403
                                                                      |                                                                    
gTTTGCGGGCCAGTTGCGTCTGGCTGACGTCATGCTGATACTGACCTACGCCGATAGATTTCGGATCGAtt                                                                      <  1:2606632/71‑1 (MQ=255)
gTTTGCGGGCCAGTTGCGTCTGGCTGACGTCATGCTGATACTGACCTACGCCGATAGATTTCGGATCGAtt                                                                      <  1:2728658/71‑1 (MQ=255)
gTTTGCGGGCCAGTTGCGTCTGGCTGACGTCATGCTGATACTGACCTACGCCGATAGATTTCGGATCGAtt                                                                      <  1:3213142/71‑1 (MQ=255)
   tGCGGGCCAGTTGCGTCTGGCTGACGTCATGCTGATACTGACCTACGCCGATAGATTTCGGATCGAtttt                                                                    >  1:741153/1‑70 (MQ=255)
    gCGGGCCAGTTGCGTCTGGCTGACGTCATGCTGATACTGACCTACGCCGATAGATTTCGGATCGATTTTc                                                                   <  1:1801006/70‑1 (MQ=255)
      gggCCAGTTGCGTCTGGCTGACGTCATGCTGATACTGACCTACGCCGATAGATTTCGGATCGATTTTCAcc                                                                <  1:2251506/71‑1 (MQ=255)
        gCCAGTTGCGTCTGGCTGACGTCATGCTGATACTGACCTACGCCGATAGATTTCGGATCGATTTTCACCAg                                                              >  1:1933332/1‑71 (MQ=255)
               gCGTCTGGCTGACGTCATGCTGATACTGACCTACGCCGATAGATTTCGGATCGAtt                                                                      >  1:558308/1‑56 (MQ=255)
               gCGTCTGGCTGACGTCATGCTGATACTGACCTACGCCGATAGATTTCGGATCGAtt                                                                      >  1:64617/1‑56 (MQ=255)
               gCGTCTGGCTGACGTCATGCTGATACTGACCTACGCCGATAGATTTCGGATCGATTTTCACCAGCTCCGc                                                        >  1:2962185/1‑70 (MQ=255)
               gCGTCTGGCTGACGTCATGCTGATACTGACCTACGCCGATAGATTTCGGATCGATTTTCACCAGCTCCGc                                                        >  1:1841534/1‑70 (MQ=255)
               gCGTCTGGCTGACGTCATGCTGATACTGACCTACGCCGATAGATTTCGGATCGATTTTCACCAGCTCCGc                                                        >  1:188681/1‑70 (MQ=255)
               gCGTCTGGCTGACGTCATGCTGATACTGACCTACGCCGATAGATTTCGGATCGATTTTCACCAGCTCCGc                                                        <  1:2697400/70‑1 (MQ=255)
                   cTGGCTGACGTCATGCTGATACTGACCTACGCCGATAGATTTCGGATCGATTTTCACCAGCTc                                                           >  1:3309793/1‑63 (MQ=255)
                   cTGGCTGACGTCATGCTGATACTGACCTACGCCGATAGATTTCGGATCGATTTTCACCAGCTCCGCCAGCg                                                   >  1:2433256/1‑71 (MQ=255)
                    tGGCTGACGTCATGCTGATACTGACCTACGCCGATAGATTTCGGATCGATTTTCACCAGCTCCGCCAGCg                                                   >  1:1707292/1‑70 (MQ=255)
                    tGGCTGACGTCATGCTGATACTGACCTACGCCGATAGATTTCGGATCGATTTTCACCAGCTCCGCCAGCg                                                   >  1:2579544/1‑70 (MQ=255)
                      gCTGACGTCATGCTGATACTGACCTACGCCGATAGATTTCGGATCGATTTTCACCAGCTCCGCCAGCGGa                                                 <  1:3081906/70‑1 (MQ=255)
                       cTGACGTCATGCTGATACTGACCTACGCCGATAGATTTCGGATCGATTTTCACCAGCTCCGCCAGCGGATc                                               >  1:1913189/1‑71 (MQ=255)
                          aCGTCATGCTGATACTGACCTACGCCGATAGATTTCGGATCGATTTTCACCAGCTCCGCCAGCGGATCCTg                                            >  1:3245242/1‑71 (MQ=255)
                             tCATGCTGATACTGACCTACGCCGATAGATTTCGGATCGATTTTCACCAGCTCCGCCAGCGGATCCTGCa                                          >  1:1298807/1‑70 (MQ=255)
                               aTGCTGATACTGACCTACGCCGATAGATTTCGGATCGATTTTCACCAGCTCCGCCAGCGGATCCTGCaaa                                        <  1:2282698/70‑1 (MQ=255)
                               aTGCTGATACTGACCTACGCCGATAGATTTCGGATCGATTTTCACCAGCTCCGCCAGCGGATCCTGCaaa                                        <  1:1198282/70‑1 (MQ=255)
                               aTGCTGATACTGACCTACGCCGATAGATTTCGGATCGATTTTCACCAGCTCCGCCAGCGGATCCTGCAAAc                                       >  1:676111/1‑71 (MQ=255)
                                 gCTGATACTGACCTACGCCGATAGATTTCGGATCGATTTTCTCCAg                                                              >  1:2840550/1‑46 (MQ=255)
                                 gCTGATACTGACCTACGCCGATAGATTTCGGATCGATTTTCACCAg                                                              >  1:228192/1‑46 (MQ=255)
                                 gCTGATACTGACCTACGCCGATAGATTTCGGATCGATTTTCACCAg                                                              >  1:478935/1‑46 (MQ=255)
                                 gCTGATACTGACCTACGCCGATAGATTTCGGATCGATTTTCACCAg                                                              >  1:1761845/1‑46 (MQ=255)
                                  cTGATACTGACCTACGCCGATAGATTTCGGATCGATTTTCACCAGCTCCGCCAGCGGa                                                 <  1:1400733/58‑1 (MQ=255)
                                  cTGATACTGACCTACGCCGATAGATTTCGGATCGATTTTCACCAGCTCCGCCAGCGGATCCTGCAAAcggc                                    >  1:2338258/1‑71 (MQ=255)
                                  cTGATACTGACCTACGCCGATAGATTTCGGATCGATTTTCACCAGCTCCGCCAGCGGATCCTGCAAAcggc                                    >  1:1821034/1‑71 (MQ=255)
                                     aTACTGACCTACGCCGATAGATTTCGGATCGATTTTCACCAGCTCCGCCAGCGGATCCTGCAAAcggcg                                   <  1:1752186/69‑1 (MQ=255)
                                      tACTGACCTACGCCGATAGATTTCGGATCGATTTTCACCAGCTCCGCCAGCGGATCCTGCAAACggcgggc                                <  1:817519/71‑1 (MQ=255)
                                      tACTGACCTACGCCGATAGATTTCGGATCGATTTTCACCAGCTCCGCCAGCGGATCCTGCAAACggcgggc                                >  1:2055915/1‑71 (MQ=255)
                                        cTGACCTACGCCGATAGATTTCGGATCGATTTTCACCAGCTCCGCCAGCGGATCCTGCAAACGGCGGGCGa                              <  1:3139245/71‑1 (MQ=255)
                                        cTGACCTACGCCGATAGATTTCGGATCGATTTTCACCAGCTCCGCCAGCGGATCCTGCAAACGGCGGGCGa                              <  1:1481402/71‑1 (MQ=255)
                                         tGACCTACGCCGATAGATTTCGGATCGATTTTCACCAGCTCCGCCAGCGGATCCTGCAAACGGCGGGCGa                              >  1:1837115/1‑70 (MQ=255)
                                               aCGCCGATCGACTTCGGATCGATCTTGACCAGTTCGGCCAGCGGATcc                                              <  1:3331199/48‑1 (MQ=19)
                                               aCGCCGATCGACTTCGGATCGATCTTGACCAGTTCGGCCAGCGGATcc                                              <  1:2108225/48‑1 (MQ=19)
                                               aCGCCGATCGACTTCGGATCGATCTTGACCAGTTCGGCCAGCGGATcc                                              <  1:2169859/48‑1 (MQ=19)
                                                cGCCGATAGATTTCGGATCGATTTTCACCAGCTCCGCCAGCGGATCCTGCAAACGGCGGGCGATAGacac                       <  1:54263/70‑1 (MQ=255)
                                                      tAGATTTCGGATCGATTTTCACCAGCTCCGCCAGCGGATCCTg                                            >  1:1959605/1‑43 (MQ=255)
                                                      tAGATTTCGGATCGATTTTCACCAGCTCCGCCAGCGGATCCTGCAAACGGCGGGCGATAGACACCGCGcc                 <  1:679028/70‑1 (MQ=255)
                                                             cGGATCGATTTTCACCAGCTCCGCCAGCGGATCCTGCAAACGGCGGGCGATAGACACCGCGCCACGCAGc          >  1:3388262/1‑70 (MQ=255)
                                                             cGGATCGATTTTCACCAGCTCCGCCAGCGGATCCTGCAAACAGCGGGCGATAGACACCGCGCCACGCAGc          >  1:1670262/1‑70 (MQ=255)
                                                              ggATCGATTTTCACCAGCTCCGCCAGCGGATCCTGCAAACGGCGGGCGATAGACACCGCGCCACGCAGCGa        <  1:626333/71‑1 (MQ=255)
                                                                   gATCTTCACCAGCTCCGCCAGCGGGTCCTGCAGGCGACGGGCTATGGATACGGCGCCGCGCAGGGACACGt   >  1:1357987/1‑71 (MQ=255)
                                                                   gATCTTCACCAGCTCCGCCAGCGGGTCCTGCAGGCGACGGGCGATGGATACGGCGCCGCGCAGGGACACGt   >  1:55007/1‑71 (MQ=255)
                                                                    aTTTTCACCAGCTCCGCCAGCGGATCCTGCAAACGGCGGGCGATAGACACCGCGCCACGCAGCGAAACGTc  <  1:637323/71‑1 (MQ=255)
                                                                      |                                                                    
GTTTGCGGGCCAGTTGCGTCTGGCTGACGTCATGCTGATACTGACCTACGCCGATAGATTTCGGATCGATTTTCACCAGCTCCGCCAGCGGATCCTGCAAACGGCGGGCGATAGACACCGCGCCACGCAGCGAAACGTC  >  minE/2682265‑2682403

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