Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,565,635 C→T 100% T137T (ACG→ACA glk ← glucokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,565,6350CT79.3% 28.9 / 2.5 24T137T (ACG→ACAglkglucokinase
Reads supporting (aligned to +/- strand):  ref base C (5/0);  new base T (19/0);  total (24/0)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGCCTGGCAAGCTTACCCAACGCTTATCGACATGGACCAGATGCGCAACCCCAAGCCCCGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATG  >  minE/1565577‑1565703
                                                          |                                                                    
cGCCTGGCAAGCTTACCCAACGCTTATCGACATGGACCAGATGCGCAACCCCAAGCCCCGTTCCGGCAc                                                            >  1:2381690/1‑69 (MQ=255)
                ccAACGCTTATCGACATGGACCAGATGCGCAACCCCAAGCCCCGTTCCGGCACCGTAAACCGCAATAgg                                            >  1:3097141/1‑69 (MQ=255)
                         aTCGACATGGACCAGATGCGCAACCCCGAGCCCTGTTCCGGCACCGTAAACCGCAATAGGCtt                                         >  1:687964/1‑63 (MQ=255)
                         aTCGACATGGACCAGATGCGCAACCCCAAGCCCTGTTCCGGCACCGTAAACCGCAATAGGCtt                                         >  1:2125297/1‑63 (MQ=255)
                         aTCGACATGGACCAGATGCGCAACCCCAAGCCCTGTTCCGGCACCGTAAACCGCAATAGGCt                                          >  1:960137/1‑62 (MQ=255)
                         aTCGACATGGACCAGATGCGCAACCCCAAGCCCTGTTCCGGCACCGTAAACCGCAATAGGCt                                          >  1:1846943/1‑62 (MQ=255)
                                  gACCAGATGCGCAACCCCAAGCCCCGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCTGTTc                          >  1:518060/1‑69 (MQ=255)
                                  gACCAGATGCGCAACCCCAAGCCCCGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGtt                           >  1:2896575/1‑68 (MQ=255)
                                  gACCAGATGCGCAACCCCAAGCCCCGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGTTc                          >  1:328179/1‑69 (MQ=255)
                                          gcgcAACCCCAAGCCCTGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCTAcc                               >  1:2407332/1‑56 (MQ=255)
                                          gcgcAACCCCAAGCCCTGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCTAcc                               >  1:619540/1‑56 (MQ=255)
                                          gcgcAACCCCAAGCCCTGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCTAcc                               >  1:3271471/1‑56 (MQ=255)
                                          gcgcAACCCCAAGCCCTGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCTAcc                               >  1:3195015/1‑56 (MQ=255)
                                          gcgcAACCCCAAGCCCTGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCTAcc                               >  1:292950/1‑56 (MQ=255)
                                          gcgcAACCCCAAGCCCTGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCTAcc                               >  1:2463883/1‑56 (MQ=255)
                                          gcgcAACCCCAAGCCCTGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCTAcc                               >  1:2352263/1‑56 (MQ=255)
                                          gcgcAACCCCAAGCCCTGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCTAcc                               >  1:21665/1‑56 (MQ=255)
                                          gcgcAACCCCAAGCCCTGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCTAcc                               >  1:196746/1‑56 (MQ=255)
                                          gcgcAACCCCAAGCCCTGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCTAcc                               >  1:1796469/1‑56 (MQ=255)
                                          gcgcAACCCCAAGCCCTGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCTAcc                               >  1:1382287/1‑56 (MQ=255)
                                          gcgcAACCCCAAGCCCTGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCTAcc                               >  1:1271615/1‑56 (MQ=255)
                                          gcgcAACCCCAAGCCCTGTTCCGGCACCGTAAACCGCAATAGGCTTAACTTCTAcc                               >  1:2905941/1‑56 (MQ=255)
                                          gcgcAACCCCAAGCCCTGTTCCGGCACCGTAAACAGCAATAGGCTTACCTTcta                                 >  1:769283/1‑52 (MQ=255)
                                          gcgcAACCCAAGCCCTGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTc                                   >  1:2706048/1‑51 (MQ=255)
                                                          cGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGCCACCAAACTGAATCATATg  <  1:3059752/69‑1 (MQ=255)
                                                          cGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATg  <  1:2345259/69‑1 (MQ=255)
                                                          cGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATg  <  1:23268/69‑1 (MQ=255)
                                                          cGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATg  <  1:1866578/69‑1 (MQ=255)
                                                          cGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATg  <  1:1621031/69‑1 (MQ=255)
                                                          cGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATg  <  1:1504294/69‑1 (MQ=255)
                                                          cGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGTGCTGCGCCACCAAACTGAATCAGATg  <  1:1263423/69‑1 (MQ=255)
                                                          |                                                                    
CGCCTGGCAAGCTTACCCAACGCTTATCGACATGGACCAGATGCGCAACCCCAAGCCCCGTTCCGGCACCGTAAACCGCAATAGGCTTACCTTCGACCGGTTCTGCGCCACCAAACTGAATCAGATG  >  minE/1565577‑1565703

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: