Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,525,940 G→C 54.2% A78A (GCG→GCC pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,525,9400GC54.2% 15.4 / 94.8 83A78A (GCG→GCCpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (21/17);  new base C (4/41);  total (25/58)
Fisher's exact test for biased strand distribution p-value = 7.08e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.47e-01

GCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGC  >  minE/2525872‑2525994
                                                                    |                                                      
gCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCg                                                        >  1:1197598/1‑69 (MQ=255)
  aTTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGAc                                                      <  1:1010818/69‑1 (MQ=255)
  aTTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGAc                                                      <  1:626076/69‑1 (MQ=255)
  aTTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGAc                                                      <  1:2227028/69‑1 (MQ=255)
  aTTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGAc                                                      <  1:2083943/69‑1 (MQ=255)
  aTTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGAc                                                      <  1:2044697/69‑1 (MQ=255)
  aTTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGAc                                                      <  1:7851/69‑1 (MQ=255)
  aTTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGAc                                                      <  1:914738/69‑1 (MQ=255)
  aTTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGAc                                                      <  1:1411461/69‑1 (MQ=255)
   ttGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGAcc                                                     >  1:300633/1‑69 (MQ=255)
   ttGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGAcc                                                     >  1:982977/1‑69 (MQ=255)
   ttGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGAcc                                                     >  1:888389/1‑69 (MQ=255)
   ttGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGAc                                                      <  1:663811/68‑1 (MQ=255)
     ggtggATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCAc                                                   >  1:342358/1‑69 (MQ=255)
                     gAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCt                                    >  1:238265/1‑68 (MQ=255)
                     gAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCt                                    >  1:74856/1‑68 (MQ=255)
                     gAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTg                                   >  1:1367985/1‑69 (MQ=255)
                     gAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTg                                   >  1:3145048/1‑69 (MQ=255)
                     gAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTg                                   >  1:992265/1‑69 (MQ=255)
                     gAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTg                                   >  1:1092736/1‑69 (MQ=255)
                     gAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTg                                   >  1:1096708/1‑69 (MQ=255)
                     gAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTg                                   >  1:294271/1‑69 (MQ=255)
                     gAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTg                                   >  1:714956/1‑69 (MQ=255)
                     gAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTg                                   >  1:955128/1‑69 (MQ=255)
                     gAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTg                                   >  1:781302/1‑69 (MQ=255)
                     gAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTg                                   >  1:1665647/1‑69 (MQ=255)
                     gAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTg                                   >  1:1738275/1‑69 (MQ=255)
                     gAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTg                                   >  1:887605/1‑69 (MQ=255)
                     gAGTTCGATCTGCTGTTTGGACCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTg                                   >  1:1193016/1‑69 (MQ=255)
                                      tGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCAtgac                    <  1:2288078/67‑1 (MQ=255)
                                      tGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCAtgac                    <  1:341493/67‑1 (MQ=255)
                                      tGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCAtgac                    <  1:2096562/67‑1 (MQ=255)
                                         cccTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCAtgaccttg                <  1:838904/68‑1 (MQ=255)
                                         cccTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCAtgaccttg                <  1:3167350/68‑1 (MQ=255)
                                         cccTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCAtgaccttg                <  1:2072089/68‑1 (MQ=255)
                                          ccTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCAtgaccttga               <  1:1966702/68‑1 (MQ=255)
                                          ccTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCAtgaccttg                <  1:929681/67‑1 (MQ=255)
                                              cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCAtgacctt                 >  1:1573139/1‑62 (MQ=255)
                                                   aaaGTGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa      <  1:2839297/68‑1 (MQ=255)
                                                   aaaGTGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa      <  1:3378464/68‑1 (MQ=255)
                                                   aaaGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCTTa      <  1:663915/68‑1 (MQ=255)
                                                   aaaGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa      <  1:977098/68‑1 (MQ=255)
                                                   aaaGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa      <  1:657526/68‑1 (MQ=255)
                                                   aaaGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa      <  1:917290/68‑1 (MQ=255)
                                                   aaaGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa      <  1:491979/68‑1 (MQ=255)
                                                   aaaGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa      <  1:457429/68‑1 (MQ=255)
                                                   aaaGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa      <  1:400728/68‑1 (MQ=255)
                                                   aaaGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa      <  1:364686/68‑1 (MQ=255)
                                                   aaaGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa      <  1:895972/68‑1 (MQ=255)
                                                   aaaGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa      <  1:693568/68‑1 (MQ=255)
                                                   aaaGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa      <  1:783485/68‑1 (MQ=255)
                                                   aaaGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa      <  1:3025609/68‑1 (MQ=255)
                                                   aaaGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa      <  1:363380/68‑1 (MQ=255)
                                                   aaaGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa      <  1:2145780/68‑1 (MQ=255)
                                                   aaaGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa      <  1:1776260/68‑1 (MQ=255)
                                                   aaaGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa      <  1:1767784/68‑1 (MQ=255)
                                                   aaaGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa      <  1:2780098/68‑1 (MQ=255)
                                                   aaaGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa      <  1:1193702/68‑1 (MQ=255)
                                                   aaaGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa      <  1:1171016/68‑1 (MQ=255)
                                                   aaaGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa      <  1:3089867/68‑1 (MQ=255)
                                                   aaaGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa      <  1:3105705/68‑1 (MQ=255)
                                                   aaaGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa      <  1:1096302/68‑1 (MQ=255)
                                                   aaaGGGATCCCGATTGCCACAACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa      <  1:48935/68‑1 (MQ=255)
                                                     aGGGATCCCGATTTCCACCACTACCGCCGTGGCGCTGGCGGAGca                           <  1:972539/45‑1 (MQ=38)
                                                     aGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGca                           <  1:1184095/45‑1 (MQ=255)
                                                     aGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGca                           <  1:14481/45‑1 (MQ=255)
                                                     aGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGca                           <  1:1461425/45‑1 (MQ=255)
                                                     aGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGca                           <  1:1933895/45‑1 (MQ=255)
                                                     aGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGca                           <  1:3113796/45‑1 (MQ=255)
                                                     aGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGca                           <  1:1967055/45‑1 (MQ=255)
                                                     aGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGca                           <  1:2341897/45‑1 (MQ=255)
                                                     aGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGca                           <  1:708732/45‑1 (MQ=255)
                                                     aGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGca                           <  1:2805844/45‑1 (MQ=255)
                                                     aGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGca                           <  1:2957617/45‑1 (MQ=255)
                                                     aGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGca                           <  1:3072900/45‑1 (MQ=255)
                                                     aGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGca                           <  1:383238/45‑1 (MQ=255)
                                                     aGGGATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa      <  1:1044931/66‑1 (MQ=255)
                                                     aGGGATCCCGATTGCCACCACTACAGCCGTGGCGCTGGCGGAGca                           <  1:1199108/45‑1 (MQ=38)
                                                       ggATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGc  >  1:2693033/1‑68 (MQ=255)
                                                       ggATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGt       >  1:808437/1‑63 (MQ=255)
                                                       ggATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGc  >  1:775990/1‑68 (MQ=255)
                                                       ggATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGc  >  1:1338174/1‑68 (MQ=255)
                                                       ggATCCCGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGc  >  1:2996615/1‑68 (MQ=255)
                                                            ccGATTGCCACCACTACCGCCGTGGCGCTGGCGGAGca                           <  1:2854676/38‑1 (MQ=255)
                                                                    |                                                      
GCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGC  >  minE/2525872‑2525994

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: