Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,291,520 C→T 31.0% D744D (GAC→GAT ppc → phosphoenolpyruvate carboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,291,5200CT31.0% 69.0 / 34.0 74D744D (GAC→GATppcphosphoenolpyruvate carboxylase
Reads supporting (aligned to +/- strand):  ref base C (23/28);  new base T (12/11);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 6.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

TGATGCTCCCCGCCTGGCTGGGTGCAGGTACGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCCGCGATTGGCCATTCTTC  >  minE/2291465‑2291568
                                                       |                                                
tGATGCTCCCCGCCTGGCTGGGTGCAGGTACGGCGCTGCAAAAAGTGGTCGAAGATGGCAAACAgagcg                                     >  1:3408791/1‑69 (MQ=255)
tGATGCTCCCCGCCTGGCTGGGTGCAGGTACGGCGCTGCAAAAAGTGGTCGAAGATGGCAAACAgagcg                                     >  1:2628736/1‑69 (MQ=255)
tGATGCTCCCCGCCTGGCTGGGTGCAGGTACGGCGCTGCAAAAAGTGGTCGAAGATGGCAAACAgagcg                                     >  1:2603959/1‑69 (MQ=255)
tGATGCTCCCCGCCTGGCTGGGTGCAGGTACGGCGCTGCAAAAAGTGGTCGAAGATGGCAAACAgagcg                                     >  1:25276/1‑69 (MQ=255)
tGATGCTCCCCGCCTGGCTGGGTGCAGGTACGGCGCTGCAAAAAGTGGTCGAAGATGGCAAACAgagcg                                     >  1:1057080/1‑69 (MQ=255)
tGATGCTCCCCGCCTGGCTGGGTGCAGGTACGGCGCTGCAAAAAGTGGTCGAAGATGGCAAACAgagcg                                     >  1:2862858/1‑69 (MQ=255)
tGATGCTCCCCGCCTGGCTGGGTGCAGGTACGGCGCTGCAAAAAGTGGTCGAAGATGGCAAACAgagcg                                     >  1:1715618/1‑69 (MQ=255)
tGATGCTCCCCGCCTGGCTGGGTGCAGGTACGGCGCTGCAAAAAGTGGTCGAAGATGGCAAACAgagcg                                     >  1:49278/1‑69 (MQ=255)
tGATGCTCCCCGCCTGGCTGGGTGCAGGTACGGCGCTGCAAAAAGTGGTCGAAGATGGCAAACAgagcg                                     >  1:1098745/1‑69 (MQ=255)
tGATGCTCCCCGCCTGGCTGGGTGCAGGTACGGCGCTGCAAAAAGTGGTCGAAGATGGCAAACAgagcg                                     >  1:784504/1‑69 (MQ=255)
tGATGCTCCCCGCCTGGCTGGGTGCAGGTACGGCGCTGCAAAAAGTGGTCGAAGATGGCAAACAgagcg                                     >  1:1348284/1‑69 (MQ=255)
tGATGCTCCCCGCCTGGCTGGGTGCAGGTACGGCGCTGCAAAAAGTGGTCGAAGATGGCAAACAgagc                                      >  1:806642/1‑68 (MQ=255)
            ccTGGCTGGGTGCAGGTACGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTa                         <  1:2596737/69‑1 (MQ=255)
                      tGCAGGTACGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCcgcg                <  1:1074251/68‑1 (MQ=255)
                      tGCAGGTACGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCcgcg                <  1:770551/68‑1 (MQ=255)
                      tGCAGGTACGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCcgcg                <  1:231008/68‑1 (MQ=255)
                      tGCAGGTACGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCcgcg                <  1:206225/68‑1 (MQ=255)
                      tGCAGGTACGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCcgcg                <  1:2858849/68‑1 (MQ=255)
                      tGCAGGTACGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCcgcg                <  1:2892597/68‑1 (MQ=255)
                      tGCAGGTACGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCcgcg                <  1:1867232/68‑1 (MQ=255)
                      tGCAGGTACGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCcgcg                <  1:569604/68‑1 (MQ=255)
                      tGCAGGTACGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCcgcg                <  1:1523895/68‑1 (MQ=255)
                           gTACGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCCGCGATTGGc          <  1:1811363/69‑1 (MQ=255)
                           gTACGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCCGCGATTGGc          <  1:2075945/69‑1 (MQ=255)
                           gTACGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCCGCGATTGGc          <  1:1163153/69‑1 (MQ=255)
                           gTACGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCCGCGATTGGc          <  1:2476481/69‑1 (MQ=255)
                           gTACGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCCGCGATTGGc          <  1:265337/69‑1 (MQ=255)
                             accgcgcTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCCGCGATTGGCCa        >  1:2708880/4‑69 (MQ=255)
                             aCGGCGCTTCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTa                         <  1:434854/52‑1 (MQ=255)
                             aCGGCGCTGCAAAAAGTGGTCGAAGACGGCTAACAGAGCGAGCTGGAGGCTATGTGCCGCGAt              >  1:735283/1‑63 (MQ=255)
                             aCGGCGCTGCAAAAAGTGGTCGAAGACGGCAATCAGAGCGAGCTGGAGGCTATGTGCCGCGATTGGcc         >  1:2908785/1‑68 (MQ=255)
                             aCGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTg                                >  1:2088038/1‑45 (MQ=255)
                             aCGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTa                         <  1:2672571/52‑1 (MQ=255)
                             aCGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTa                         >  1:415632/1‑52 (MQ=255)
                             aCGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTa                         <  1:844097/52‑1 (MQ=255)
                             aCGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTa                         <  1:1241865/52‑1 (MQ=255)
                             aCGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTa                         <  1:2057705/52‑1 (MQ=255)
                             aCGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTa                         >  1:1622799/1‑52 (MQ=255)
                             aCGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTa                         <  1:180450/52‑1 (MQ=255)
                             aCGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTa                         <  1:1833797/52‑1 (MQ=255)
                             aCGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTa                         <  1:1863598/52‑1 (MQ=255)
                             aCGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCCGCGATTGGcc         >  1:943448/1‑68 (MQ=255)
                             aCGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCCGCGATTGGcc         >  1:586173/1‑68 (MQ=255)
                             aCGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCCGCGATTGGcc         >  1:54044/1‑68 (MQ=255)
                             aCGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCCGCGATTGGcc         >  1:1799668/1‑68 (MQ=255)
                             aCGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCCGCGATTGGc          <  1:261250/67‑1 (MQ=255)
                             aCGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCCGCGATTGGc          <  1:3063415/67‑1 (MQ=255)
                             aCGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCCGCGATTGGCCa        >  1:2803485/1‑69 (MQ=255)
                             aCGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCCGCGATTGGCCa        >  1:1126382/1‑69 (MQ=255)
                             aCGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCCGCGATTGGCCa        >  1:1383582/1‑69 (MQ=255)
                             aCGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCCGCGATTGGCCa        >  1:2416953/1‑69 (MQ=255)
                             aCGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCCGCGATTGGCCa        >  1:2745789/1‑69 (MQ=255)
                             aCGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCCGCGATTGGCCa        >  1:548278/1‑69 (MQ=255)
                             aCGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCCGCGATTGGCCa        >  1:2014625/1‑69 (MQ=255)
                             aCGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCCGCGATTGGCCa        >  1:297581/1‑69 (MQ=255)
                             aCGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCCGCGATTGGCCa        >  1:3053410/1‑69 (MQ=255)
                             aCGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCCGCGATTGGCCa        >  1:2025793/1‑69 (MQ=255)
                             aCGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCCGCGATTGGCCa        >  1:327579/1‑69 (MQ=255)
                             aCGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCCGCGATTGGCCa        >  1:2835052/1‑69 (MQ=255)
                              cGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCcgcg                <  1:843124/60‑1 (MQ=255)
                                gcgcTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCCGCGATTGGc          <  1:369294/64‑1 (MQ=255)
                                          aaGTGGTCGAAGATGGCAAACAGAGCGAGCTGGAAGCCATGTGCcgc                 <  1:3095189/47‑1 (MQ=38)
                                          aaGTGGTCGAAGATGGCAAACAGAGCGAGCTGGAAGCCATGTGCcgc                 <  1:2018474/47‑1 (MQ=38)
                                          aaGTGGTCGAAGATGGCAAACAGAGCGAGCTGGAAGCCATGTGCcgc                 <  1:2562183/47‑1 (MQ=38)
                                          aaGTGGTCGAAGATGGCAAACAGAGCGAGCTGGAAGCCATGTGCcgc                 <  1:1501202/47‑1 (MQ=38)
                                          aaGTGGTCGAAGATGGCAAACAGAGCGAGCTGGAAGCCATGTGCcgc                 <  1:1539627/47‑1 (MQ=38)
                                          aaGTGGTCGAAGATGGCAAACAGAGCGAGCTGGAAGCCATGTGCcgc                 <  1:674986/47‑1 (MQ=38)
                                          aaGTGGTCGAAGATGGCAAACAGAGCGAGCTGGAAGCCATGTGCcgc                 <  1:556001/47‑1 (MQ=38)
                                          aaGTGGTCGAAGATGGCAAACAGAGCGAGCTGGAAGCCATGTGCcgc                 <  1:540665/47‑1 (MQ=38)
                                          aaGTGGTCGAAGATGGCAAACAGAGCGAGCTGGAAGCCATGTGCcgc                 <  1:528573/47‑1 (MQ=38)
                                          aaGTGGTCGAAGATGGCAAACAGAGCGAGCTGGAAGCCATGTGCcgc                 <  1:1629367/47‑1 (MQ=38)
                                          aaGTGGTCGAAGATGGCAAACAGAGCGAGCTGGAAGCCATGTGCcgc                 <  1:2810648/47‑1 (MQ=38)
                                          aaGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCcgcg                <  1:818118/48‑1 (MQ=255)
                                                tCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCCGCGATTGGCCAttcttc  >  1:501306/1‑56 (MQ=255)
                                                       |                                                
TGATGCTCCCCGCCTGGCTGGGTGCAGGTACGGCGCTGCAAAAAGTGGTCGAAGACGGCAAACAGAGCGAGCTGGAGGCTATGTGCCGCGATTGGCCATTCTTC  >  minE/2291465‑2291568

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