Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 67,361 G→A 27.4% G303G (GGC→GGT leuA ← 2‑isopropylmalate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE67,3610GA27.4% 76.6 / 25.7 69G303G (GGC→GGTleuA2‑isopropylmalate synthase
Reads supporting (aligned to +/- strand):  ref base G (19/31);  new base A (8/11);  total (27/42)
Fisher's exact test for biased strand distribution p-value = 7.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TTTGGTTCAGACCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCTTTGT  >  minE/67300‑67419
                                                             |                                                          
tttGGTTCAGCCCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATcc                                                      >  1:1316520/1‑68 (MQ=255)
     ttCAGACCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGgtg                                                 <  1:2300113/68‑1 (MQ=255)
     ttCAGACCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGgtg                                                 <  1:385648/68‑1 (MQ=255)
     ttCAGACCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGgtg                                                 <  1:1903406/68‑1 (MQ=255)
          aCCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATAc                                            <  1:2104952/68‑1 (MQ=255)
               tAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATATCggagg                                      <  1:441469/69‑1 (MQ=255)
               tAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCggagg                                      <  1:1769173/69‑1 (MQ=255)
               tAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCggagg                                      <  1:2360249/69‑1 (MQ=255)
               tAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCggagg                                      <  1:2475742/69‑1 (MQ=255)
               tAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCggagg                                      <  1:2255160/69‑1 (MQ=255)
               tAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCggagg                                      <  1:2077608/69‑1 (MQ=255)
               tAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCggagg                                      <  1:2572708/69‑1 (MQ=255)
               tAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCggagg                                      <  1:2845526/69‑1 (MQ=255)
               tAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCggagg                                      <  1:1938951/69‑1 (MQ=255)
               tAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCggagg                                      <  1:2957173/69‑1 (MQ=255)
               tAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCggagg                                      <  1:3244746/69‑1 (MQ=255)
               tAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCggagg                                      <  1:1494236/69‑1 (MQ=255)
               tAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCggagg                                      <  1:390579/69‑1 (MQ=255)
               tAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCggagg                                      <  1:444822/69‑1 (MQ=255)
               aataTTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCggagg                                      <  1:981394/66‑1 (MQ=255)
                 gATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGgtg                                                 >  1:2671777/1‑56 (MQ=255)
                 gATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGgtg                                                 >  1:3031831/1‑56 (MQ=255)
                 gATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGgtg                                                 >  1:281081/1‑56 (MQ=255)
                 gATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGgtg                                                 >  1:267223/1‑56 (MQ=255)
                 gATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGgtg                                                 >  1:261652/1‑56 (MQ=255)
                 gATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGgtg                                                 >  1:3292122/1‑56 (MQ=255)
                 gATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGgtg                                                 >  1:907419/1‑56 (MQ=255)
                 gATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGgtg                                                 >  1:1126634/1‑56 (MQ=255)
                 gATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGgtg                                                 >  1:1135237/1‑56 (MQ=255)
                 gATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGgtg                                                 >  1:2057547/1‑56 (MQ=255)
                 gATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGgtg                                                 >  1:1256654/1‑56 (MQ=255)
                 gATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGgtg                                                 >  1:2030217/1‑56 (MQ=255)
                                aTTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCt                                                     <  1:3394464/37‑1 (MQ=255)
                                aTTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCt                                                     <  1:2989737/37‑1 (MQ=255)
                                      tAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGcc                <  1:2407626/68‑1 (MQ=255)
                                      tAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGcc                <  1:1333375/68‑1 (MQ=255)
                                      tAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGcc                <  1:3208649/68‑1 (MQ=255)
                                      tAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGcc                <  1:3177830/68‑1 (MQ=255)
                                      tAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGcc                <  1:1727194/68‑1 (MQ=255)
                                      tAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGcc                <  1:1988453/68‑1 (MQ=255)
                                      tAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGcc                <  1:2595626/68‑1 (MQ=255)
                                      tAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGcc                <  1:2002410/68‑1 (MQ=255)
                                        gTTTTCGCGGTTTTTCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCcaac             >  1:2282342/1‑69 (MQ=255)
                                        gTTTTCGCGGTTTTTCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCcaac             >  1:1151966/1‑69 (MQ=255)
                                        gTTTTCGCGGTTTTTCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCcaac             >  1:2920317/1‑69 (MQ=255)
                                        gTTTTCGCGGTTTTTCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCcaac             >  1:1714183/1‑69 (MQ=255)
                                        gTTTTCGCGGTTTTTCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCcaac             >  1:113701/1‑69 (MQ=255)
                                        gTTTTCGCGGTTTTTCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCcaac             >  1:1146017/1‑69 (MQ=255)
                                        gTTTTCGCGGTTTTTCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCcaac             >  1:40713/1‑69 (MQ=255)
                                        gTTTTCGCGGTTTTTCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCcaac             >  1:3189514/1‑69 (MQ=255)
                                        gTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGCcaac             <  1:1084642/69‑1 (MQ=255)
                                        gTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGCcaac             <  1:3207759/69‑1 (MQ=255)
                                           ttCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGCcaacaa           >  1:1658292/1‑68 (MQ=255)
                                           ttCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAAt          >  1:2942100/1‑69 (MQ=255)
                                           ttCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAAt          >  1:44899/1‑69 (MQ=255)
                                           ttCGCGGTTTTTCAGCACGCAATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGCcaacaa           >  1:3133332/1‑68 (MQ=255)
                                                   ttttCAGTACACTATCCTGGTGGATACCGGAGGAGTGTGCGAAt                           <  1:2807979/44‑1 (MQ=25)
                                                   ttttCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAAt                           <  1:2340772/44‑1 (MQ=37)
                                                   ttttCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAAt                           <  1:791427/44‑1 (MQ=37)
                                                   ttttCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAAt                           <  1:1993936/44‑1 (MQ=37)
                                                   ttttCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAAt                           <  1:3428923/44‑1 (MQ=37)
                                                   ttttCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAAt                           <  1:2666574/44‑1 (MQ=37)
                                                   ttttCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAAt                           <  1:3147708/44‑1 (MQ=37)
                                                   ttttCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAAt                           <  1:3120243/44‑1 (MQ=37)
                                                   ttttCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAAt                           <  1:2795113/44‑1 (MQ=37)
                                                   ttttCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAAt                           <  1:2852562/44‑1 (MQ=37)
                                                   ttttCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCtttgt  >  1:2769140/1‑69 (MQ=255)
                                                   ttttCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCtttgt  >  1:2848113/1‑69 (MQ=255)
                                                      tCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAAt                           <  1:2901649/41‑1 (MQ=37)
                                                             |                                                          
TTTGGTTCAGACCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCTTTGT  >  minE/67300‑67419

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: