Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,017,320 T→C 100% L352L (TTG→CTG)  manA → mannose‑6‑phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,017,3200TC80.5% 106.7 / 15.5 71L352L (TTG→CTG) manAmannose‑6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base T (14/0);  new base C (41/16);  total (55/16)
Fisher's exact test for biased strand distribution p-value = 2.95e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AAGAAACCACCATTAGCCAGCAGAGTGCCGCCATTTTGTTCTGCGTCGAAGGCGATGCAACGTTGTGGAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGAATCAGCGTTTATTGCCGCCAA  >  minE/1017258‑1017381
                                                              |                                                             
aaGAAACCACCATTAGCCAGCAGAGTGCCGCCATTTTGTTCTGGGTCGAAGGCGATGCAACGTtgtg                                                           >  1:792856/1‑67 (MQ=255)
aaGAAACCACCATTAGCCAGCAGAGTGCCGCCATTTTGTTCTGCGTCGAAGGCGATGCAACGTtgtg                                                           >  1:3027753/1‑67 (MQ=255)
aaGAAACCACCATTAGCCAGCAGAGTGCCGCCATTTTGTTCTGCGTCGAAGGCGATGCAACGTtgtg                                                           >  1:2709647/1‑67 (MQ=255)
aaGAAACCACCATTAGCCAGCAGAGTGCCGCCATTTTGTTCTGCGTCGAAGGCGATGCAACGTtgtg                                                           >  1:1895492/1‑67 (MQ=255)
        accaTCAGCCAGCAGAGTGCCTCCATTTTGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGtt                                                   >  1:3185909/1‑67 (MQ=255)
        accaTCAGCCAGCAGAGTGCCGCCATTTTGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGtt                                                   >  1:1114126/1‑67 (MQ=255)
        accaTCAGCCAGCAGAGTGCCGCCATTTTGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGtt                                                   >  1:671625/1‑67 (MQ=255)
        accaTCAGCCAGCAGAGTGCCGCCATTTTGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGtt                                                   >  1:3147641/1‑67 (MQ=255)
        accaTCAGCCAGCAGAGTGCCGCCATTTTGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGtt                                                   >  1:3010770/1‑67 (MQ=255)
        accaTCAGCCAGCAGAGTGCCGCCATTTTGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGtt                                                   >  1:2821634/1‑67 (MQ=255)
        accaTCAGCCAGCAGAGTGCCGCCATTTTGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGtt                                                   >  1:274915/1‑67 (MQ=255)
        accaTCAGCCAGCAGAGTGCCGCCATTTTGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGtt                                                   >  1:2628497/1‑67 (MQ=255)
        accaTCAGCCAGCAGAGTGCCGCCATTTTGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGtt                                                   >  1:2290361/1‑67 (MQ=255)
        accaTCAGCCAGCAGAGTGCCGCCATTTTGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGtt                                                   >  1:2198716/1‑67 (MQ=255)
        accaTCAGCCAGCAGAGTGCCGCCATTTTGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGtt                                                   >  1:2141680/1‑67 (MQ=255)
        accaTCAGCCAGCAGAGTGCCGCCATTTTGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGtt                                                   >  1:2065535/1‑67 (MQ=255)
        accaTCAGCCAGCAGAGTGCCGCCATTTTGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGtt                                                   >  1:1086144/1‑67 (MQ=255)
        accaTCAGCCAGCAGAGTGCCGCCATTTTGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGtt                                                   >  1:1462394/1‑67 (MQ=255)
        accaTCAGCCAGCAGAGTGCCGCCATTTTGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGt                                                    >  1:2568952/1‑66 (MQ=255)
            ttAGCCAGCAGAGTGCCGCCATTTTGTTCTGCGTCGAAGGCGATGCAACGTTGTGGAAAGGTTCTcagc                                             >  1:301578/1‑69 (MQ=255)
            ttAGCCAGCAGAGTGCCGCCATTTTGTTCTGCGTCGAAGGCGATGCAACGTTGTGGAAAGGTTCTcagc                                             >  1:3039201/1‑69 (MQ=255)
            ttAGCCAGCAGAGTGCCGCCATTTTGTTCTGCGTCGAAGGCGATGCAACGTTGTGGAAAGGTTCTcagc                                             >  1:2312722/1‑69 (MQ=255)
            ttAGCCAGCAGAGTGCCGCCATTTTGTTCTGCGTCGAAGGCGATGCAACGTTGTGGAAAGGTTCTcagc                                             >  1:378741/1‑69 (MQ=255)
            ttAGCCAGCAGAGTGCCGCCATTTTGTTCTGCGTCGAAGGCGATGCAACGTTGTGGAAAGGTTCTcagc                                             >  1:2045819/1‑69 (MQ=255)
            ttAGCCAGCAGAGTGCCGCCATTTTGTTCTGCGTCGAAGGCGATGCAACGTTGTGGAAAGGTTCTcagc                                             >  1:1659852/1‑69 (MQ=255)
            ttAGCCAGCAGAGTGCCGCCATTTTGTTCTGCGTCGAAGGCGATGCAACGTTGTGGAAAGGTTCTcag                                              >  1:1488524/1‑68 (MQ=255)
            ttAGCCAGCAGAGTGCCGCCATTTCGTTCTGCGTCGAAGGCGATGCAACGTTGTGGAAAGGTTCTcagc                                             >  1:1319895/1‑69 (MQ=255)
                             gccATTTTGTTCTGCGTCGAAGGCGATGCAACGTTGTGGAAAGGTTCTCAGCAGTTACATCTTAAACCg                            >  1:3413459/1‑69 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTa                                 >  1:2311150/1‑59 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTg                        >  1:285073/1‑68 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTg                        >  1:1936795/1‑68 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTg                        >  1:2252072/1‑68 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTg                        >  1:1115681/1‑68 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTg                        >  1:234609/1‑68 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       <  1:3153135/69‑1 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       >  1:1371857/1‑69 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       >  1:1232774/1‑69 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       >  1:3122922/1‑69 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       >  1:3086596/1‑69 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       >  1:689961/1‑69 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       <  1:553551/69‑1 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       <  1:1169424/69‑1 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       <  1:3329352/69‑1 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       >  1:3342622/1‑69 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       >  1:3431737/1‑69 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       <  1:591498/69‑1 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       <  1:405142/69‑1 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       <  1:471526/69‑1 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       <  1:587199/69‑1 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       <  1:2955101/69‑1 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       >  1:1572901/1‑69 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       <  1:1643585/69‑1 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       >  1:1728747/1‑69 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       >  1:1524096/1‑69 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       <  1:1941737/69‑1 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       >  1:2028725/1‑69 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       <  1:1503532/69‑1 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       <  1:2086501/69‑1 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       <  1:2124365/69‑1 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       >  1:2221328/1‑69 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       <  1:2354427/69‑1 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       >  1:2611612/1‑69 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       >  1:1391412/1‑69 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       >  1:2634463/1‑69 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       >  1:2654644/1‑69 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       <  1:1545209/69‑1 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       >  1:1346987/1‑69 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       >  1:1341095/1‑69 (MQ=255)
                                  tttGTTCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGCTGa                       >  1:1412216/1‑69 (MQ=255)
                                  tttGTCCTGCGTCGAAGGCGATGCAACGCTGTGTAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGa                       >  1:1362445/1‑69 (MQ=255)
                                                        gCAACGTTGTGGAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGAATCAGCGTTTATTGCCGCCaa  >  1:2796289/1‑68 (MQ=255)
                                                              |                                                             
AAGAAACCACCATTAGCCAGCAGAGTGCCGCCATTTTGTTCTGCGTCGAAGGCGATGCAACGTTGTGGAAAGGTTCTCAGCAGTTACAGCTTAAACCGGGTGAATCAGCGTTTATTGCCGCCAA  >  minE/1017258‑1017381

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: