Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 195,616 C→T 36.4% L71L (CTG→CTA yaeF ← predicted lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE195,6160CT36.4% 30.5 / 30.8 33L71L (CTG→CTAyaeFpredicted lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (6/15);  new base T (5/7);  total (11/22)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

ACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAGAGCTTTGA  >  minE/195548‑195685
                                                                    |                                                                     
aCGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAg                                                                     <  1:1671530/71‑1 (MQ=255)
 tgTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAg                                                                     >  1:2146545/2‑70 (MQ=255)
 cGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGc                                                                    <  1:2213689/71‑1 (MQ=255)
   tGACTCAGAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAATAGCaa                                                                  <  1:1802621/71‑1 (MQ=255)
       tCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACa                                                                      <  1:737580/63‑1 (MQ=255)
         acaGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGcc                                                            >  1:549062/1‑71 (MQ=255)
               gAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAATAGCAAATCGCCGGGttt                                                      <  1:1534558/71‑1 (MQ=255)
                  gTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTgag                                                   <  1:356799/71‑1 (MQ=255)
                    gCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAATAGCAAATCGCCGGGttt                                                      >  1:262305/1‑66 (MQ=255)
                     cTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAATAGCAAATCGCCGGGTTTGAGAt                                                 >  1:1704077/1‑70 (MQ=255)
                      tGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTga                                                    <  1:117046/66‑1 (MQ=255)
                       gAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAg                                              >  1:1003289/1‑71 (MQ=255)
                        aaGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAATAGCAAATCGCCGGGTTTGAGATCAgg                                             <  1:1527036/71‑1 (MQ=255)
                        aaGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCa                                               <  1:1238006/69‑1 (MQ=255)
                                          gAGGTTACCCCAAGGCTTGAGGAGAATAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCttt                              >  1:2032927/1‑68 (MQ=255)
                                           aGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCgg                                                          >  1:1565297/1‑39 (MQ=255)
                                           aGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTatt                          <  1:474688/71‑1 (MQ=255)
                                            ggTTACCCCAAGGCTTGAGGAGAATAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTatt                          <  1:2050221/70‑1 (MQ=255)
                                            ggTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGttt                                                      <  1:2041951/42‑1 (MQ=255)
                                            ggTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGttt                                                      <  1:2065603/42‑1 (MQ=255)
                                             gTTACCCCAAGGCTTGAGGAGAATAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTattga                        <  1:564712/71‑1 (MQ=255)
                                              ttACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGt                                                        <  1:2077104/38‑1 (MQ=255)
                                              ttACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGt                                                        <  1:1257988/38‑1 (MQ=255)
                                              ttACCCCAAGGCTTCAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTAtt                                   <  1:1132015/59‑1 (MQ=255)
                                               tACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTAtt                                   <  1:1311392/58‑1 (MQ=255)
                                                    cAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTga                                                    <  1:1794605/36‑1 (MQ=255)
                                                      aGGCTTGAGGAGAATAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCtt                               >  1:913822/1‑55 (MQ=255)
                                                       ggCTTGAGGAGAATAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCtttt                             <  1:1258993/56‑1 (MQ=255)
                                                       ggCTTGAGGAGAATAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCtttt                             <  1:1006853/56‑1 (MQ=255)
                                                        gCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTg              >  1:1021382/1‑70 (MQ=255)
                                                         cTTGAGGAGAATAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTat                           >  1:1527137/1‑56 (MQ=255)
                                                              ggagAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAGAGc       <  1:1644319/71‑1 (MQ=255)
                                                                   aCAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAGAGCTTTGa  >  1:883854/1‑71 (MQ=255)
                                                                    cAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCgg                >  1:336735/1‑56 (MQ=255)
                                                                    |                                                                     
ACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTATTGATTGTTCGGTGAAAGAGCTTTGA  >  minE/195548‑195685

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: