Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,239,208 C→T 40.6% L341L (CTG→CTA zwf ← glucose‑6‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,239,2080CT40.6% 26.1 / 42.8 64L341L (CTG→CTAzwfglucose‑6‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (4/34);  new base T (12/14);  total (16/48)
Fisher's exact test for biased strand distribution p-value = 2.52e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TTCTTTAAACAGATTCAGTTCAGGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGTTTACCAGTACGCAGGTAGAATGGCACACCGGCCCAGCGCCAGTTATCAATGTCG  >  minE/1239145‑1239268
                                                               |                                                            
ttCTTTAAACAGATTCAGTTCAGGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACg                                                         <  1:1725034/69‑1 (MQ=255)
 tCTTTAAACAGATTCAGTTCAGGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGt                                                        <  1:1168361/69‑1 (MQ=255)
 tCTTTAAACAGATTCAGTTCAGGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGt                                                        <  1:2788838/69‑1 (MQ=255)
 tCTTTAAACAGATTCAGTTCAGGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGt                                                        <  1:3407270/69‑1 (MQ=255)
 tCTTTAAACAGATTCAGTTCAGGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGt                                                        <  1:2096884/69‑1 (MQ=255)
  cTTTAAACAGATTCAGTTCAGGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGt                                                        <  1:2426714/68‑1 (MQ=255)
  cTTTAAACAGATTCAGTTCAGGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGt                                                        <  1:201772/68‑1 (MQ=255)
         cAGATTCAGTTCAGGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGTTTACCAGt                                                <  1:876510/69‑1 (MQ=255)
         cAGATTCAGTTCAGGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGTTTACCAGt                                                <  1:2424741/69‑1 (MQ=255)
         cAGATTCAGTTCAGGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGTTTACCAGt                                                <  1:2330842/69‑1 (MQ=255)
         cAGATTCAGTTCAGGTGGTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGTTTACCAGt                                                <  1:2265074/69‑1 (MQ=255)
                   tcagGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGTTTACCAGTACGCAGGTAg                                      <  1:1027921/69‑1 (MQ=255)
                   tcagGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGTTTACCAGTACGCAGGTAg                                      <  1:2663200/69‑1 (MQ=255)
                   tcagGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGTTTACCAGTACGCAGGTAg                                      <  1:2667282/69‑1 (MQ=255)
                   tcagGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGTTTACCAGTACGCAGGTAg                                      <  1:2196972/69‑1 (MQ=255)
                   tcagGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGTTTACCAGTACGCAGGTAg                                      >  1:2196092/1‑69 (MQ=255)
                   tcagGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGTTTACCAGTACGCAGGTAg                                      <  1:851632/69‑1 (MQ=255)
                   tcagGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGTTTACCAGTACGCAGGTAg                                      <  1:1324071/69‑1 (MQ=255)
                   tcagGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGTTTACCAGTACGCAGGTAg                                      <  1:1716368/69‑1 (MQ=255)
                   tcagGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGTTTACCAGTACGCAGGTAg                                      <  1:1588770/69‑1 (MQ=255)
                    cagGTGTTTTGAAATATACCACGACTTCAGAACATTTGGTCGGTAGACGTTTACCAg                                                 <  1:3213126/57‑1 (MQ=255)
                    cagGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGTAGACGTTTACCAg                                                 <  1:1923860/57‑1 (MQ=255)
                    cagGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGTAGACGTTTACCAg                                                 <  1:2501644/57‑1 (MQ=255)
                    cagGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGTAGACGTTTACCAg                                                 <  1:121240/57‑1 (MQ=255)
                    cagGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGTAGACGTTTACCAg                                                 <  1:3424868/57‑1 (MQ=255)
                    cagGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGTAGACGTTTACCAg                                                 <  1:2208782/57‑1 (MQ=255)
                    cagGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGTAGACGTTTACCAg                                                 <  1:1597026/57‑1 (MQ=255)
                    cagGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGTAGACGTTTACCAg                                                 <  1:3074962/57‑1 (MQ=255)
                    cagGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGTAGACGTTTACCAg                                                 <  1:314079/57‑1 (MQ=255)
                    cagGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGTAGACGTTTACCAg                                                 <  1:2061033/57‑1 (MQ=255)
                    cagGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGTAGACGTTTACCAg                                                 <  1:3161002/57‑1 (MQ=255)
                    cagGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGTAGACGTTTACCAg                                                 <  1:2032274/57‑1 (MQ=255)
                     agGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCAGTAGACGTTTACCAg                                                 <  1:3271782/56‑1 (MQ=255)
                            ttGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGTTTACCAGTACGCATGTAGAATgg                                 <  1:2364134/65‑1 (MQ=255)
                            ttGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGTTTACCAGTACGCAGGTAGAATgg                                 <  1:3403947/65‑1 (MQ=255)
                            ttGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGTTTACCAGTACGCAGGTAGAATgg                                 <  1:291193/65‑1 (MQ=255)
                            ttGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGTTTACCAGTACGCAGGTAGAATgg                                 <  1:2903351/65‑1 (MQ=255)
                            ttGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGTTTACCAGTACGCAGGTAGAATgg                                 <  1:922956/65‑1 (MQ=255)
                            ttGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGTTTACCAGTACGCAGGTAGAATgg                                 <  1:1115582/65‑1 (MQ=255)
                            ttGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGTTTACCAGTACGCAGGTAGAATgg                                 <  1:2446488/65‑1 (MQ=255)
                            ttGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGTTTACCAGTACGCAGGTAGAATgg                                 <  1:2320257/65‑1 (MQ=255)
                            ttGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGTTTACCAGTACGCAGGTAGAATgg                                 <  1:2049145/65‑1 (MQ=255)
                            ttGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGTTTACCAGTACGCAGGTAGAATGGcaca                             >  1:943904/1‑69 (MQ=255)
                             tGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGTTTACCAGTACGCAGGTAGAATgg                                 <  1:2613479/64‑1 (MQ=255)
                                aaTAGACCACGACTTCAGAACATTTGGTCGGTAGACGTTTACCAg                                                 <  1:4197/45‑1 (MQ=255)
                                  tAGACCACGACTTCAGAACATTTGGTCGGCAGACGttt                                                      <  1:52686/38‑1 (MQ=255)
                                  tAGACCACGACTTCAGAACATTTGGTCGGCAGACGttt                                                      <  1:2156823/38‑1 (MQ=255)
                                  tAGACCACGACTTCAGAACATTTGGTCGGCAGACGttt                                                      <  1:508232/38‑1 (MQ=255)
                                  tAGACCACGACTTCAGAACATTTGGTAGGCAGACGTTTACCAGTACGCAGGTAGAATGGCACACCgg                         >  1:1276224/1‑67 (MQ=255)
                                  tAGACCACGACTTCAGAACATTAGGTCGGCAGACGTTTACCAGTACGCAGGTAGAATGGCACACCgg                         >  1:1616264/1‑67 (MQ=255)
                                             ttCAGAACATTTGGTCGGCAGACGTTTACCAGTACGCAGGTAGAATGGCACACCGGCCCAGCGCCAGtt            <  1:955326/69‑1 (MQ=255)
                                             ttCAGAACATTTGGTCGGCAGACGTTTACCAGTACGCAGGTAGAATGGCACACCGGCCCAGCGCCAGtt            <  1:212855/69‑1 (MQ=255)
                                                       ttGGTCGGTAGACGTTTACCAGTACGCAGGTAGAATGGCACACCGGCCCAGGGCCAGTTATCAATGTCg  >  1:1921209/1‑69 (MQ=255)
                                                       ttGGTCGGTAGACGTTTACCAGTACGCAGGTAGAATGGCACACCGGCCCAGCGCCAGTTATCAATGTc   >  1:365602/1‑68 (MQ=255)
                                                       ttGGTCGGTAGACGTTTACCAGTACGCAGGTAGAATGGCACACCGGCCCAGCGCCAGTTATCAATGTCg  >  1:1239862/1‑69 (MQ=255)
                                                       ttGGTCGGTAGACGTTTACCAGTACGCAGGTAGAATGGCACACCGGCCCAGCGCCAGTTATCAATGTCg  >  1:1151115/1‑69 (MQ=255)
                                                       ttGGTCGGTAGACGTTTACCAGTACGCAGGTAGAATGGCACACCGGCCCAGCGCCAGTTATCAATGTCg  >  1:1457844/1‑69 (MQ=255)
                                                       ttGGTCGGTAGACGTTTACCAGTACGCAGGTAGAATGGCACACCGGCCCAGCGCCAGTTATCAATGTCg  >  1:583292/1‑69 (MQ=255)
                                                       ttGGTCGGTAGACGTTTACCAGTACGCAGGTAGAATGGCACACCGGCCCAGCGCCAGTTATCAATGTCg  >  1:1455803/1‑69 (MQ=255)
                                                       ttGGTCGGTAGACGTTTACCAGTACGCAGGTAGAATGGCACACCGGCCCAGCGCCAGTTATCAATGTCg  >  1:1856511/1‑69 (MQ=255)
                                                       ttGGTCGGTAGACGTTTACCAGTACGCAGGTAGAATGGCACACCGGCCCAGCGCCAGTTATCAATGTCg  >  1:2941740/1‑69 (MQ=255)
                                                       ttGGTCGGTAGACGTTTACCAGTACGCAGGTAGAATGGCACACCGGCCCAGCGCCAGTTATCAATGTCg  >  1:2569572/1‑69 (MQ=255)
                                                       ttGGTCGGTAGACGTTTACCAGTACGCAGGTAGAATGGCACACCGGCCCAGCGCCAGTTATCAAAGTc   >  1:839689/1‑68 (MQ=255)
                                                       ttGGTCGGTAGACGTTTACCAGTACGCAGGTAGAATGGCACACAGGCCCAGCGCCAGTTATCAATGTc   >  1:21141/1‑68 (MQ=255)
                                                               |                                                            
TTCTTTAAACAGATTCAGTTCAGGTGTTTTGAAATAGACCACGACTTCAGAACATTTGGTCGGCAGACGTTTACCAGTACGCAGGTAGAATGGCACACCGGCCCAGCGCCAGTTATCAATGTCG  >  minE/1239145‑1239268

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: