Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 267,465 A→G 100% T394T (ACT→ACC dxs ← 1‑deoxyxylulose‑5‑phosphate synthase, thiamine‑requiring, FAD‑requiring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE267,4650AG84.4% 96.4 / 6.9 57T394T (ACT→ACCdxs1‑deoxyxylulose‑5‑phosphate synthase, thiamine‑requiring, FAD‑requiring
Reads supporting (aligned to +/- strand):  ref base A (0/9);  new base G (40/8);  total (40/17)
Fisher's exact test for biased strand distribution p-value = 2.70e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCGAACAGGACCGGAAGCTTTTGAATCGCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAAGTGGAGTAAATCGCGACAATGGGTTTGTACCCACCAATCGCCAGACCCGCAGCA  >  minE/267397‑267519
                                                                    |                                                      
gCGAACAGGACCGGAAGCTTTTGAATCGCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGaaa                                                        <  1:437717/69‑1 (MQ=255)
gCGAACAGGACCGGAAGCTTTTGAATCGCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGaaa                                                        <  1:3174573/69‑1 (MQ=255)
gCGAACAGGACCGGAAGCTTTTGAATCGCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGaaa                                                        <  1:2877522/69‑1 (MQ=255)
gCGAACAGGACCGGAAGCTTTTGAATCGCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGaaa                                                        <  1:1756273/69‑1 (MQ=255)
 cGAACAGGACCGGAAGCTTTTGAATCGCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGaaa                                                        <  1:1896382/68‑1 (MQ=255)
        gACCGGAAGCTTTTGAATCGCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAAGTGGAGTa                                                <  1:1091924/69‑1 (MQ=255)
        gACCGGAAGCTTTTGAATCGCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAAGTGGAGTa                                                <  1:1496290/69‑1 (MQ=255)
        gACCGGAAGCTTTTGAATCGCCACGTCATGCAGCACCTGATCATAGGCGCCTTGCAGGAAAGTGGAGTa                                                <  1:1984572/69‑1 (MQ=255)
                       aaTCGCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAAGTGGAGTAAATCGCGACAATgg                                  <  1:1099099/68‑1 (MQ=255)
                       aaTCGCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAAGTGGAGTAAATCGCGACAATgg                                  <  1:470047/68‑1 (MQ=255)
                       aaTCGCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAAGTGGAGTAAATCGCGACAATgg                                  <  1:46182/68‑1 (MQ=255)
                       aaTCGCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAAGTGGAGTAAATCGCGACAATgg                                  <  1:1969040/68‑1 (MQ=255)
                         tCGCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAGGTGGAGt                                                 >  1:995895/1‑51 (MQ=255)
                         tCGCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAAt                                    >  1:1607685/1‑64 (MQ=255)
                         tCGCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATg                                   >  1:1689191/1‑65 (MQ=255)
                         tCGCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATg                                   >  1:3306147/1‑65 (MQ=255)
                         tCGCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATg                                   >  1:3282900/1‑65 (MQ=255)
                         tCGCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATg                                   >  1:328034/1‑65 (MQ=255)
                         tCGCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATg                                   >  1:1461589/1‑65 (MQ=255)
                         tCGCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATg                                   >  1:313625/1‑65 (MQ=255)
                         tCGCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATg                                   >  1:2832444/1‑65 (MQ=255)
                         tCGCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATg                                   >  1:280863/1‑65 (MQ=255)
                         tCGCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATg                                   >  1:263637/1‑65 (MQ=255)
                         tCGCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATg                                   >  1:2041280/1‑65 (MQ=255)
                         tCGCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATg                                   >  1:1684524/1‑65 (MQ=255)
                           gCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCTCAACa                                      >  1:1565067/1‑60 (MQ=255)
                           gCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCTACAATGGGtt                               >  1:2306560/1‑67 (MQ=255)
                           gCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATGGGtt                               >  1:784107/1‑67 (MQ=255)
                           gCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATGGGtt                               >  1:2293556/1‑67 (MQ=255)
                           gCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATGGGtt                               >  1:92928/1‑67 (MQ=255)
                           gCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATGGGtt                               >  1:506181/1‑67 (MQ=255)
                           gCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATGGGtt                               >  1:1650724/1‑67 (MQ=255)
                           gCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATGGGtt                               >  1:2669036/1‑67 (MQ=255)
                           gCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATGGGtt                               >  1:1427720/1‑67 (MQ=255)
                           gCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATGGGtt                               >  1:2825784/1‑67 (MQ=255)
                           gCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATGGGtt                               >  1:1874212/1‑67 (MQ=255)
                           gCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATGGGtt                               >  1:3190060/1‑67 (MQ=255)
                           gCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATGGGtt                               >  1:3419236/1‑67 (MQ=255)
                           gCCACGTCATGCAGCACCTGATAATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATGGGtt                               >  1:869655/1‑67 (MQ=255)
                                         caCCTGATCATAGGCGCGTTGCAGGAAAGTGGAGTAAATCGCGACAATGGGTTTGTACCCACCAATCGc               <  1:823779/69‑1 (MQ=255)
                                         caCCTGATCATAGGCGCGTTGCAGGAAAGTGGAGTAAATCGCGACAATGGGTTTGTACCCACCAATCGc               <  1:20259/69‑1 (MQ=255)
                                          aCCTGATCATAGGCGCGTTTCAGGAAGGTGGAGTAGATCGCAACAATGGGTTTGTAGCCTCCAATc                 >  1:1432310/1‑66 (MQ=255)
                                          aCCTGATCATAGGCGCGTTGCAGGTAGGTGGAGTa                                                >  1:337641/1‑35 (MQ=37)
                                          aCCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATGGGTTTGTAGcctcc                     >  1:396958/1‑62 (MQ=255)
                                          aCCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATGGGTTTGTAGCCTCCAATc                 >  1:432012/1‑66 (MQ=255)
                                          aCCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATGGGTTTGTAGCCTCCAATc                 >  1:1897281/1‑66 (MQ=255)
                                          aCCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATGGGTTTGTAGCCTCCAATc                 >  1:1270586/1‑66 (MQ=255)
                                          aCCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATGGGTTTGTAGCCTCCAATc                 >  1:717394/1‑66 (MQ=255)
                                          aCCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATGGGTTTGTAGCCTCCAATc                 >  1:1875378/1‑66 (MQ=255)
                                          aCCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATGGGTTTGTAGCCTCCAATc                 >  1:1475585/1‑66 (MQ=255)
                                          aCCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATGGGTTTGTAGCCTCCAATc                 >  1:193886/1‑66 (MQ=255)
                                          aCCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATGGGTTTGTAGCCTCCAATc                 >  1:2337692/1‑66 (MQ=255)
                                          aCCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCGCAACAATGGGTTTGTAGCCTCCAATc                 >  1:2208084/1‑66 (MQ=255)
                                          aCCTGATCATAGGCGCGTTGCAGGAAGGTGGAGTAGATCCCAACAATGGGTTTGTAGCCTCCAATc                 >  1:2351589/1‑66 (MQ=255)
                                                       cgcgTTGCAGGAAGGTGGAGTAGATCGCAACAATGGGTTTGTAGCCTCCAATCGCCAGACCCgcagca  <  1:3206777/68‑1 (MQ=255)
                                                       cgcgTTGCAGGAAGGTGGAGTAGATCGCAACAATGGGTTTGTAGCCTCCAATCGCCAGACCCgcagca  <  1:90811/68‑1 (MQ=255)
                                                       cgcgTTGCAGGAAGGTGGAGTAGATCGCAACAATGGGTTTGTAGCCTCCAATCGCCAGACCCgcagca  <  1:1226057/68‑1 (MQ=255)
                                                       cgcgTTGCAGGAAGGTGGAGTAGATCGCAACAATGGGTTTGTAGCCTCCAATCGCCAGACCCgcagca  <  1:3110022/68‑1 (MQ=255)
                                                       cgcgTTGCAGGAAGGTGGAGTAGATCGCAACAATGGGTTTGTAGCCTCCAATCGCCAGACCCgcagca  <  1:157190/68‑1 (MQ=255)
                                                       cgcgTTGCAGGAAGGTGGAGTAGATCGCAACAATGGGTTTGTAGCCTCCAATCGCCAGACCCgcagca  <  1:2321842/68‑1 (MQ=255)
                                                       cgcgTTGCAGGAAGGTGGAGTAGATCGCAACAATGGGTTTGTAGCCTCCAATCGCCAGACCCgcagca  <  1:200264/68‑1 (MQ=255)
                                                        gcgTTGCAGGAAGGTGGAGTAGATCGCAACAATCGGTTTGTAGCCTCCAATCGCCAGACCCgcagca  <  1:2828582/67‑1 (MQ=255)
                                                                    |                                                      
GCGAACAGGACCGGAAGCTTTTGAATCGCCACGTCATGCAGCACCTGATCATAGGCGCGTTGCAGGAAAGTGGAGTAAATCGCGACAATGGGTTTGTACCCACCAATCGCCAGACCCGCAGCA  >  minE/267397‑267519

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: