Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,279,191 C→T 7.8% L324L (CTA→TTA)  yeeJ → adhesin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,279,1910CT 7.8% 158.8 / 3.5 51L324L (CTA→TTA) yeeJadhesin
Reads supporting (aligned to +/- strand):  ref base C (22/25);  new base T (2/2);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

TGAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGT  >  minE/1279127‑1279260
                                                                |                                                                     
tGAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCg                                                                 >  1:1757504/1‑71 (MQ=255)
 gAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGc                                                                >  1:2566779/1‑71 (MQ=255)
  aaCTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCg                                                                 <  1:4109024/69‑1 (MQ=255)
    cTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTg                                                             <  1:1985323/71‑1 (MQ=255)
        aCAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCg                                                         <  1:2761269/71‑1 (MQ=255)
        aCAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCg                                                         <  1:2104505/71‑1 (MQ=255)
          aaCGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTg                                                             <  1:956121/65‑1 (MQ=255)
          aaCGATTATGAAGCACGACCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCt                                                              >  1:3196251/1‑64 (MQ=255)
           aCGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCAc                                                      >  1:1625532/1‑71 (MQ=255)
            ggATTATCAGGAACGCCCGGCGAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCagc                                                     >  1:3126947/2‑69 (MQ=255)
             gATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGcc                                                               >  1:581064/1‑60 (MQ=255)
             gATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCt                                                    >  1:910732/1‑71 (MQ=255)
             gATTATCAGGAACGCCCGGCGAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCt                                                    <  1:2662854/71‑1 (MQ=255)
              aTTATCAGGAACGCCCGGCGAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGc                                                        <  1:1419477/66‑1 (MQ=255)
                     aGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTa                                            <  1:485619/71‑1 (MQ=255)
                     aGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTa                                            <  1:2608193/71‑1 (MQ=255)
                     aGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTa                                            <  1:2513920/71‑1 (MQ=255)
                      gCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTa                                            <  1:2702910/70‑1 (MQ=255)
                               gCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTaaa                                          >  1:2985413/1‑63 (MQ=255)
                                ccAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTa                                            <  1:2165892/60‑1 (MQ=255)
                                 cAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTAt                                 >  1:2391338/1‑70 (MQ=255)
                                 cAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATg                                <  1:3805810/71‑1 (MQ=255)
                                   aTGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTa                                            <  1:639173/57‑1 (MQ=255)
                                    tGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCagc                                                     >  1:3804651/1‑45 (MQ=255)
                                     ggCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTaaa                                          >  1:2665002/1‑57 (MQ=255)
                                     ggCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTa                                  >  1:3931117/1‑65 (MQ=255)
                                      gCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAg                           <  1:2589805/71‑1 (MQ=255)
                                      gCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAg                           <  1:2032599/71‑1 (MQ=255)
                                       cTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGt                          <  1:792457/71‑1 (MQ=255)
                                           gATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGt                                     >  1:4045751/1‑56 (MQ=255)
                                           gATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGtatta                      >  1:3454334/1‑71 (MQ=255)
                                            aTGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAAc                             >  1:539588/1‑63 (MQ=255)
                                             tGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGta                         <  1:1465955/66‑1 (MQ=255)
                                                aCGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGc                                                >  1:2928907/1‑40 (MQ=255)
                                                aCGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGt                                     >  1:1348917/1‑51 (MQ=255)
                                                 cgcgCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGa                >  1:438868/1‑71 (MQ=255)
                                                  gcgcAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGAt               <  1:3917097/71‑1 (MQ=255)
                                                  gcgcAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGAt               >  1:2203185/1‑71 (MQ=255)
                                                   cgcAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTc                                    <  1:1461261/49‑1 (MQ=255)
                                                   cgcAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATg              <  1:3566021/71‑1 (MQ=255)
                                                    gcAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATg              <  1:1852795/70‑1 (MQ=255)
                                                    gcAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATg              <  1:1401545/70‑1 (MQ=255)
                                                       gAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAAc                             <  1:1865654/52‑1 (MQ=255)
                                                         aaGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGc                  >  1:1910476/1‑61 (MQ=255)
                                                         aaGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTgg        <  1:1878218/71‑1 (MQ=255)
                                                         aaGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTgg        >  1:2354644/1‑71 (MQ=255)
                                                            cTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGccc     >  1:921553/1‑71 (MQ=255)
                                                             tGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGAt               <  1:247926/60‑1 (MQ=255)
                                                               gCTACCCGCCTGGCCGCATCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGt  >  1:2034695/1‑71 (MQ=255)
                                                               gCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGt  <  1:3745945/71‑1 (MQ=255)
                                                               gCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGt  >  1:2211100/1‑71 (MQ=255)
                                                                cTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGa             >  1:2062383/1‑59 (MQ=255)
                                                                |                                                                     
TGAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGT  >  minE/1279127‑1279260

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: ‑