Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,526,291 G→A 100% P195P (CCG→CCA pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,526,2910GA88.0% 151.0 / 10.8 58P195P (CCG→CCApyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (7/0);  new base A (0/51);  total (7/51)
Fisher's exact test for biased strand distribution p-value = 3.33e-09
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.

GATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGG  >  minE/2526240‑2526316
                                                   |                         
gATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACAt            >  1:1027648/1‑67 (MQ=255)
          atcaCCCTGAAAGAGCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:2038212/53‑1 (MQ=37)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGTCGGa                <  1:297786/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:591171/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:210939/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:2206461/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:22771/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:282018/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:292593/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:329176/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:361393/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:392499/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:433689/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:452836/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:2070249/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:616912/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:63315/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:682887/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:726999/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:734449/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:747705/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:757932/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:814511/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:849184/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:889038/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:906916/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:922041/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:990229/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:1908466/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:1056588/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:1057059/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:1275945/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:1329866/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:1342225/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:1479055/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:1487980/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:1509377/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:1558143/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:1568480/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:1721194/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:1903155/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:1863778/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:1852383/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:1845918/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:1840401/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:1834703/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:1753400/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:1645195/53‑1 (MQ=255)
          atcaCCCTGAAAGACCTGATTGCCTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:1072568/53‑1 (MQ=38)
          atcaCCCCGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:300237/53‑1 (MQ=255)
           tcaCCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:65452/52‑1 (MQ=255)
            caCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGcggcgg  >  1:1483198/1‑65 (MQ=255)
            caCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGcggcgg  >  1:1402430/1‑65 (MQ=255)
            caCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGcggcgg  >  1:59769/1‑65 (MQ=255)
            caCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGcggcgg  >  1:155638/1‑65 (MQ=255)
            caCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGcggcgg  >  1:1106882/1‑65 (MQ=255)
            caCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGcggcgg  >  1:1922219/1‑65 (MQ=255)
                     aGACCTGATTGCTTACCTGGAAGAGAAGCCAGAAATGGCGGa                <  1:1264489/42‑1 (MQ=25)
                                                   |                         
GATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGG  >  minE/2526240‑2526316

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