Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,642,808 T→C 15.6% Q66Q (CAA→CAG ypfJ ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,642,8080TC15.6% 77.0 / 7.6 32Q66Q (CAA→CAGypfJconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (15/12);  new base C (3/2);  total (18/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCA  >  minE/1642741‑1642875
                                                                   |                                                                   
ttGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAttgttgt                                                                  <  1:2376230/71‑1 (MQ=255)
 tGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAttgttgt                                                                  >  1:1070948/1‑70 (MQ=255)
  gCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAttgttgttg                                                                >  1:2448164/1‑71 (MQ=255)
  gCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAttgttgtt                                                                 <  1:1704458/70‑1 (MQ=255)
   ccAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAttgttgtt                                                                 <  1:855793/69‑1 (MQ=255)
    cAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTAGGGCTAATTGACCGCGTTGATTGCTGTTgg                                                               <  1:824030/70‑1 (MQ=255)
    cAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTAGGGCTAATTGACCGCGTTGATTGCTGTTgg                                                               <  1:3752530/70‑1 (MQ=255)
      gAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTAGGGCTAATTGACCGCGTTGATTGCTGTTGGGa                                                             >  1:3646630/1‑70 (MQ=255)
      gAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTAGGGCTAATTGACCGCGTTGATTGCTGTTGGGa                                                             >  1:2135943/1‑70 (MQ=255)
      gAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTAGGGCTAATTGACCGCGTTGATTGCTGTTGGGa                                                             >  1:351067/1‑70 (MQ=255)
           aCCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCgg                                                       >  1:853570/1‑71 (MQ=255)
            ccGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGCCCGCGTTGAttgttgt                                                                  >  1:2566813/1‑59 (MQ=255)
                ggTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGcc                                                  <  1:1482181/71‑1 (MQ=255)
                 gTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGccc                                                 <  1:2079174/71‑1 (MQ=255)
                   gAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCg                                                >  1:224524/1‑70 (MQ=255)
                        ttGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGtcat                                           >  1:3421468/1‑70 (MQ=255)
                        ttGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGtcat                                           >  1:2473393/1‑70 (MQ=255)
                         tGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGt                                              >  1:847699/1‑66 (MQ=255)
                                tCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCg                                    <  1:1826230/69‑1 (MQ=255)
                                tCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGt                                  >  1:3695849/1‑71 (MQ=255)
                                     ttCATTTGGGCTAATTGACCGCGTTGAttgttgttg                                                                >  1:3326141/1‑36 (MQ=255)
                                     ttCATTTGGGCTAATTGACCGCGTTGAttgttgttg                                                                >  1:3652310/1‑36 (MQ=255)
                                     ttCATTTGGGCTAATTGACCGCGTTGAttgttgttg                                                                >  1:1251970/1‑36 (MQ=255)
                                     ttCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCg                                                        <  1:3158454/44‑1 (MQ=255)
                                      tCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAAt                             >  1:1432580/1‑70 (MQ=255)
                                       cATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGc                                                      >  1:779590/1‑44 (MQ=255)
                                       cATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATc                            <  1:2681094/70‑1 (MQ=255)
                                       cATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATc                            <  1:3492781/70‑1 (MQ=255)
                                       cATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATc                            <  1:1098881/70‑1 (MQ=255)
                                        aTTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCaa                          <  1:922170/71‑1 (MQ=255)
                                         tttGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTg                                                    <  1:2821461/44‑1 (MQ=255)
                                                 aaTTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCAtag                  >  1:3537444/1‑70 (MQ=255)
                                                  aTTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCAtag                  >  1:1136578/1‑69 (MQ=255)
                                                   ttGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCAtagtag               >  1:3510671/1‑71 (MQ=255)
                                                     gACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACAcc                      >  1:3164573/1‑62 (MQ=255)
                                                     gACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGcc             >  1:312683/1‑71 (MQ=255)
                                                      aCCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCg                                    <  1:2890293/47‑1 (MQ=255)
                                                            ttGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAAcc      <  1:3789351/71‑1 (MQ=255)
                                                            ttGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAAcc      <  1:2322778/71‑1 (MQ=255)
                                                              gATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTa                                >  1:346820/1‑43 (MQ=255)
                                                              gATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTa                                >  1:1594696/1‑43 (MQ=255)
                                                                ttgttgttgGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCa  >  1:3628796/1‑71 (MQ=255)
                                                                ttgttgttgGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCa  >  1:2532230/1‑71 (MQ=255)
                                                                   |                                                                   
TTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCA  >  minE/1642741‑1642875

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: