Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,148,960 T→A 100% A68A (GCA→GCT glmM ← phosphoglucosamine mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,148,9600TA99.1% 366.9 / ‑4.5 105A68A (GCA→GCTglmMphosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (0/104);  total (0/105)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
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.

Only 100 of 105 total aligned reads displayed.

CCACGGCCGGTGTTGGCATCGGGCCAGTGAAGAGTGCGGAAAGGCCCGCTGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATATAGCCAGAAATACGCGTG  >  minE/2148911‑2149014
                                                 |                                                      
ccACGGCCGGTGTTGGCATCGGGCCAGTGAAGAGGGCGGAAAGGCCCGCTGCCGCCAGACCCGAtt                                        <  1:39459/66‑1 (MQ=255)
      ccGGTGTTGTCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:156642/67‑1 (MQ=255)
      ccGGTGTTGGCATCGTACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:2175867/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGCCCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:2067071/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCTCTTCCAGTGc                                 <  1:2257758/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCTCTTCCAGTGc                                 <  1:2255144/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:1785953/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:1012006/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:1816748/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:1828994/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:1907493/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:1920758/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:1936157/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:2033398/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:2068817/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:2074603/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:221003/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:443772/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:558740/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:594793/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:605525/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:637516/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:658424/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:737137/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:775455/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:838649/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:846933/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:902597/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:1586476/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:1069417/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:1076534/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:1573840/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:1564192/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:1133684/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:1159973/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:1200534/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:1524970/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:1306052/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:1389530/67‑1 (MQ=255)
      ccGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:1379299/67‑1 (MQ=255)
      ccGGTCTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:30010/67‑1 (MQ=255)
       cGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:1608425/66‑1 (MQ=255)
       cGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:917980/66‑1 (MQ=255)
       cGGTGTTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:2262165/66‑1 (MQ=255)
        ggTGTTGGCATCGGACCAGTGAAGAGTGCTGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:1325877/65‑1 (MQ=255)
         gtgtTGGCATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:1014941/64‑1 (MQ=255)
                cATCGGACCAGTGAAGAGTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:467354/57‑1 (MQ=255)
                                 gTGCGGAAAGGCCCGCAGCCGCCAGACCCGCTTCCAGTGc                                 <  1:591952/40‑1 (MQ=255)
                                      atcaGGCCCGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:119468/63‑1 (MQ=255)
                                           gCCCGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:788942/61‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCGAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:363484/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATTTAGCCAGAAATACGCgtg  <  1:2029341/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:517719/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:505126/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:638301/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:538837/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:1852235/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:489389/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:1557549/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:110548/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:610035/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:612571/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:1078055/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:4238/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:1766115/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:1041689/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:812764/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:812998/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:837637/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:1033456/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:857176/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:1001189/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:908826/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:1628159/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:2210422/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:1989019/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:2020723/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:2022204/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:193701/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:1323350/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:1318588/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:2088351/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:2145434/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:215708/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:1923652/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:1279154/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:1189253/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:2220544/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:1265966/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:1537483/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:1905550/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:2314800/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:236389/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:300074/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:1852330/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:389581/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:39928/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:1975086/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACCCGCTTCAAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:1693711/59‑1 (MQ=255)
                                             ccGCAGCCGCCAGACACGCTTCCAGTGCTGACTCCAGCATGTAGCCAGAAATACGCgtg  <  1:421779/59‑1 (MQ=255)
                                                 |                                                      
CCACGGCCGGTGTTGGCATCGGGCCAGTGAAGAGTGCGGAAAGGCCCGCTGCCGCCAGACCCGCTTCCAGTGCTGACTCCAGCATATAGCCAGAAATACGCGTG  >  minE/2148911‑2149014

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: