Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,844,848 A→T 38.0% A79A (GCT→GCA cysN ← sulfate adenylyltransferase, subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,844,8480AT38.0% 35.6 / 51.3 61A79A (GCT→GCAcysNsulfate adenylyltransferase, subunit 1
Reads supporting (aligned to +/- strand):  ref base A (21/16);  new base T (9/14);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 2.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.44e-01

TAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTG  >  minE/1844783‑1844913
                                                                 |                                                                 
tAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGcc                                                                <  1:1221712/69‑1 (MQ=255)
 aGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCaga                                                                   >  1:950700/1‑65 (MQ=255)
 aGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCaga                                                                   >  1:1178773/1‑65 (MQ=255)
 aGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCaga                                                                   >  1:1369663/1‑65 (MQ=255)
 aGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCaga                                                                   >  1:2482173/1‑65 (MQ=255)
 aGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCaga                                                                   >  1:2250083/1‑65 (MQ=255)
 aGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCaga                                                                   >  1:2194153/1‑65 (MQ=255)
 aGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCaga                                                                   >  1:1874136/1‑65 (MQ=255)
 aGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCaga                                                                   >  1:1874058/1‑65 (MQ=255)
             catCAATAGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGTGCCAGATCCAGCtttt                                                   >  1:1174150/3‑69 (MQ=255)
             catCAATAGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGTGCCAGATCCAGCtttt                                                   >  1:906691/3‑69 (MQ=255)
             catCAATAGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGTGCCAGATCCAGCtttt                                                   >  1:5442/3‑69 (MQ=255)
             catCAATAGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGTGCCAGATCCAGCtttt                                                   >  1:2955002/3‑69 (MQ=255)
             catCAATAGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGTGCCAGATCCAGCtttt                                                   >  1:2761820/3‑69 (MQ=255)
             catCAATAGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGTGCCAGATCCAGCtttt                                                   >  1:2063146/3‑69 (MQ=255)
             catCAATAGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGTGCCAGATCCAGCtttt                                                   >  1:1835966/3‑69 (MQ=255)
             catCAATAGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGTGCCAGATCCAGCttt                                                    >  1:3222408/3‑68 (MQ=255)
             catCAATAGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGTGCCAGATCCAGCttt                                                    >  1:1212396/3‑68 (MQ=255)
                         tGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGc                                       >  1:247796/1‑69 (MQ=255)
                         tGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGc                                       >  1:79926/1‑69 (MQ=255)
                              tGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGa                                  >  1:803481/1‑69 (MQ=255)
                              tGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGa                                  >  1:2349684/1‑69 (MQ=255)
                              tGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGa                                  >  1:2709132/1‑69 (MQ=255)
                              tGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGa                                  >  1:791345/1‑69 (MQ=255)
                              tGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGa                                  >  1:3336765/1‑69 (MQ=255)
                              tGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTAGCCCTGGGTGCCGTGa                                  >  1:2459123/1‑69 (MQ=255)
                              tGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCAAGCTTTTCGCCCTGGGTGCCGTGa                                  >  1:637024/1‑69 (MQ=255)
                                  cgcgCTCAGCTTGCAGGCTGTCCACCAGCAGTGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCt                              <  1:1282454/69‑1 (MQ=255)
                                  cgcgCTCAGCTTGCAGGCCGTCCACCAGCAGTGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCt                              <  1:1830639/69‑1 (MQ=255)
                                  cgcgCTCAGCTTGCAGGCCGTCCACCAGCAGTGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCt                              <  1:2724333/69‑1 (MQ=255)
                                  cgcgCTCAGCTTGCAGGCCGTCCACCAGCAGTGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCt                              <  1:3244068/69‑1 (MQ=255)
                                  cgcgCTCAGCTTGCAGGCCGTCCACCAGCAGTGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCt                              <  1:3232152/69‑1 (MQ=255)
                                  cgcgCTCAGCTTGCAGGCCGTCCACCAGCAGTGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCt                              <  1:3050109/69‑1 (MQ=255)
                                  cgcgCTCAGCTTGCAGGCCGTCCACCAGCAGTGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCt                              <  1:2397155/69‑1 (MQ=255)
                                  cgcgCTCAGCTTGCAGGCCGTCCACCAGCAGTGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCt                              <  1:1569383/69‑1 (MQ=255)
                                  cgcgCTCAGCTTGCAGGCCGTCCACCAGCAGTGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCt                              <  1:1616205/69‑1 (MQ=255)
                                  cgcgCTCAGCTTGCAGGCCGTCCACCAGCAGTGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCt                              <  1:748010/69‑1 (MQ=255)
                                  cgcgCTCAGCTTGCAGGCCGTCCACCAGCAGTGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCt                              <  1:1111117/69‑1 (MQ=255)
                                  cgcgCTCAGCTTGCAGGCCGTCCACCAGCAGTGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCt                              <  1:2011180/69‑1 (MQ=255)
                                  cgcgCTCAGCTTGCAGGCCGTCCACCAGCAGTGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCt                              <  1:1724171/69‑1 (MQ=255)
                                  cgcgCTCAGCTTGCAGGCCGTCCACCAGCAGCGCCAGATCCAGCTTTTCGCGCTGGGTGCCGTGACGCt                              >  1:3211313/1‑69 (MQ=255)
                                  cgcgCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCt                              >  1:1393266/1‑69 (MQ=255)
                                  cgcgCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCt                              >  1:1940645/1‑69 (MQ=255)
                                      cTCAGCTTGCAGGCCGTCCACCAGCAGTGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCt                              <  1:46085/65‑1 (MQ=255)
                                        cAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCg                                                 >  1:2623163/1‑44 (MQ=255)
                                        cAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGc                                                >  1:862251/1‑45 (MQ=255)
                                        cAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGc                                                >  1:1175534/1‑45 (MQ=255)
                                        cAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGc                                                >  1:975061/1‑45 (MQ=255)
                                        cAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGc                                                >  1:3065283/1‑45 (MQ=255)
                                        cAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGc                                                >  1:296034/1‑45 (MQ=255)
                                        cAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGc                                                >  1:2767874/1‑45 (MQ=255)
                                        cAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGc                                                >  1:1336835/1‑45 (MQ=255)
                                        cAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGc                                                >  1:2158239/1‑45 (MQ=255)
                                        cAGCTTGCAGGCCGTCCACCAGAAGAGCCAGATCCAGCtttt                                                   >  1:2733041/1‑42 (MQ=255)
                                                             gcagAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCt   <  1:3322669/69‑1 (MQ=255)
                                                             gcagAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCt   <  1:3169178/69‑1 (MQ=255)
                                                             gcagAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCt   <  1:2617374/69‑1 (MQ=255)
                                                             gcagAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCt   <  1:2035362/69‑1 (MQ=255)
                                                             gcagAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCt   <  1:2174954/69‑1 (MQ=255)
                                                             gcagAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCt   <  1:892319/69‑1 (MQ=255)
                                                             gcagAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCt   <  1:80657/69‑1 (MQ=255)
                                                              cagAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCt   <  1:2978651/68‑1 (MQ=255)
                                                              cagAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTg  <  1:179311/69‑1 (MQ=255)
                                                              cagAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTg  <  1:20583/69‑1 (MQ=255)
                                                              cagAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTg  <  1:1630359/69‑1 (MQ=255)
                                                              cagAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTg  <  1:2141874/69‑1 (MQ=255)
                                                              cagAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTg  <  1:2172865/69‑1 (MQ=255)
                                                              cagAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTg  <  1:667541/69‑1 (MQ=255)
                                                              cagAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTg  <  1:520810/69‑1 (MQ=255)
                                                                 |                                                                 
TAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTG  >  minE/1844783‑1844913

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: