Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,912,429 C→T 100% A135T (GCG→ACG)  pyrI ← aspartate carbamoyltransferase, regulatory subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,912,4290CT98.7% 154.3 / ‑6.0 47A135T (GCG→ACG) pyrIaspartate carbamoyltransferase, regulatory subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (33/13);  minor base G (0/1);  total (33/14)
Fisher's exact test for biased strand distribution p-value = 2.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGA  >  minE/2912387‑2912488
                                          |                                                           
caTTTTGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:1596319/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:2838236/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:2335365/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:2362400/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:2449811/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:2494403/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:2496950/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:2569435/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:2681947/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:2824336/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:1061604/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:344503/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:636842/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:65063/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:67782/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:746798/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:972563/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:986842/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:1652029/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:1029059/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:1087109/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:1156472/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:1271533/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:1502255/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:1645408/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:1741755/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:1830666/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:1859065/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:2087890/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:2058615/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:2043265/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGTGata                                                         >  1:1982076/1‑47 (MQ=255)
caTTATGGGAAAACTCTTTTTCACAGTATTTACATTTGAGCGTGata                                                         >  1:771014/1‑47 (MQ=255)
                                cttttGAGCGGGATATCATTGGCGCGTTTTCGCACGGCTAAGCTGGATGAAACCGGTTCGGCATGTCTGa  <  1:866371/68‑1 (MQ=255)
                                cATTTGAGCGTGATATCATTGGCTCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGa  <  1:2889921/70‑1 (MQ=255)
                                cATTTGAGCGTGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGa  <  1:2189836/70‑1 (MQ=255)
                                cATTTGAGCGTGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGa  <  1:2770380/70‑1 (MQ=255)
                                cATTTGAGCGTGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGa  <  1:2464814/70‑1 (MQ=255)
                                cATTTGAGCGTGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGa  <  1:130678/70‑1 (MQ=255)
                                cATTTGAGCGTGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGa  <  1:2495495/70‑1 (MQ=255)
                                cATTTGAGCGTGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGa  <  1:793639/70‑1 (MQ=255)
                                cATTTGAGCGTGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGa  <  1:883231/70‑1 (MQ=255)
                                cATTTGAGCGTGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGa  <  1:957030/70‑1 (MQ=255)
                                cATTTGAGCGTGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGa  <  1:2358580/70‑1 (MQ=255)
                                cATTTGAGCGTGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGa  <  1:2556934/70‑1 (MQ=255)
                                cATCTGAGCGTGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGa  <  1:1534761/70‑1 (MQ=255)
                                 aTTTGAGCGTGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGa  <  1:2704729/69‑1 (MQ=255)
                                          |                                                           
CATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGA  >  minE/2912387‑2912488

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: