Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,521,023 A→G 58.6% S293S (TCT→TCC usg ← predicted semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,521,0230AG58.6% 37.7 / 66.7 94S293S (TCT→TCCusgpredicted semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (25/14);  new base G (38/17);  total (63/31)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.19e-01

ACCGACCAGAACTGGACTTGCTCCGGCATACCGTAGTCATTACGCACGCAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCACCTACCTGAGTTGGGAATTCGTTCTCTTCAGAG  >  minE/1520967‑1521083
                                                        |                                                            
aCCGACCAGAACTGAACTTGCTCCGGCATACCGTAGTCATTACGTATGCAGCCAACGGa                                                            >  1:2867053/1‑59 (MQ=255)
aCCGACCAGAACTGAACTTGCTCCGGCATACCGTAGTCATTACGTATGCAGCCAACGGa                                                            >  1:2853986/1‑59 (MQ=255)
aCCGACCAGAACTGAACTTGCTCCGGCATACCGTAGTCATTACGTATGCAGCCAACGGa                                                            >  1:2755707/1‑59 (MQ=255)
aCCGACCAGAACTGAACTTGCTCCGGCATACCGTAGTCATTACGTATGCAGCCAACGGa                                                            >  1:114799/1‑59 (MQ=255)
aCCGACCAGAACTGAACTTGCTCCGGCATACCGTAGTCATTACGTATGCAGCCAACGGa                                                            >  1:2750446/1‑59 (MQ=255)
aCCGACCAGAACTGAACTTGCTCCGGCATACCGTAGTCATTACGTATGCAGCCAACGGa                                                            >  1:2523842/1‑59 (MQ=255)
aCCGACCAGAACTGAACTTGCTCCGGCATACCGTAGTCATTACGTATGCAGCCAACGGa                                                            >  1:2458587/1‑59 (MQ=255)
aCCGACCAGAACTGAACTTGCTCCGGCATACCGTAGTCATTACGTATGCAGCCAACGGa                                                            >  1:1347959/1‑59 (MQ=255)
aCCGACCAGAACTGAACTTGCTCCGGCATACCGTAGTCATTACGTATGCAGCCAACGGa                                                            >  1:3214426/1‑59 (MQ=255)
aCCGACCAGAACTGAACTTGCTCCGGCATACCGTAGTCATTACGTATGCAGCCAACGGa                                                            >  1:241743/1‑59 (MQ=255)
aCCGACCAGAACTGAACTTGCTCCGGCATACCGTAGTCATTACGTATGCAGCCAACGGa                                                            >  1:2389550/1‑59 (MQ=255)
aCCGACCAGAACTGAACTTGCTCCGGCATACCGTAGTCATTACGTATGCAGCCAACGGa                                                            >  1:1624303/1‑59 (MQ=255)
aCCGACCAGAACTGAACTTGCTCCGGCATACCGTAGTCATTACGTATGCAGCCAACGGa                                                            >  1:345202/1‑59 (MQ=255)
aCCGACCAGAACTGAACTTGCTCCGGCATACCGTAGTCATTACGTATGCAGCCAACGGa                                                            >  1:2141247/1‑59 (MQ=255)
aCCGACCAGAACTGAACTTCCTCCGGCATACCGTAGTCATTACGTATGCAGCCAACGGa                                                            >  1:396003/1‑59 (MQ=255)
       aGAACTGGACTTGCTCCGGCATACCGTAGTCATTACGCACGCAGCCAACAGAAAGATGCGGCGTAc                                              >  1:3083295/1‑66 (MQ=255)
                cTTGCTCCGGCATACCGTAGTCATTACGCACGCAGCCAACAGAAAGATGCGGCGTACCCGAAGCATc                                    <  1:1019196/67‑1 (MQ=255)
                cTTGCTCCGGCATACCGTAGTCATTACGCACGCAGCCAACAGAAAGATGCGGCGTACCCGAAGCATc                                    <  1:188840/67‑1 (MQ=255)
                cTTGCTCCGGCATACCGTAGTCATTACGCACGCAGCCAACAGAAAGATGCGGCGTACCCGAAGCATc                                    <  1:2324935/67‑1 (MQ=255)
                cTTGCTCCGGCATACCGTAGTCATTACGCACGCAGCCAACAGAAAGATGCGGCGTACCCGAAGCATc                                    <  1:2732666/67‑1 (MQ=255)
                              ccGTAGTCATTACGCACGCAGCCAACAGAAAGATGCGGCGTACCCGACGCATCa                                   <  1:440425/54‑1 (MQ=255)
                              ccGTAGTCATTACGCACGCAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCa                                   <  1:3324898/54‑1 (MQ=255)
                              ccGTAGTCATTACGCACGCAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCa                                   <  1:2367269/54‑1 (MQ=255)
                              ccGTAGTCATTACGCACGCAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCa                                   <  1:2348780/54‑1 (MQ=255)
                              ccGTAGTCATTACGCACGCAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCa                                   <  1:2318018/54‑1 (MQ=255)
                              ccGTAGTCATTACGCACGCAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCa                                   <  1:812593/54‑1 (MQ=255)
                              ccGTAGTCATTACGCACGCAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCa                                   <  1:822343/54‑1 (MQ=255)
                              ccGTAGTCATTACGCACGCAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCa                                   <  1:1183177/54‑1 (MQ=255)
                              ccGTAGTCATTACGCACGCAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCa                                   <  1:2834605/54‑1 (MQ=255)
                              ccGTAGTCATTACGCACGCAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCACCTACCTGAGTTggg                    >  1:812344/1‑69 (MQ=255)
                              ccGTAGTCATTACGCACGCAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCACCTACCTGAGTTggg                    >  1:507469/1‑69 (MQ=255)
                              ccGTAGTCATTACGCACGCAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCACCTACCTGAGTTggg                    >  1:3306948/1‑69 (MQ=255)
                              ccGTAGTCATTACGCACGCAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCACCTACCTGAGTTggg                    >  1:2316729/1‑69 (MQ=255)
                              ccGTAGTCATTACGCACGCAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCACCTACCTGAGTTggg                    >  1:2968388/1‑69 (MQ=255)
                              ccGTAGTCATTACGCACGCAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCACCTACCTGAGTTggg                    >  1:2964536/1‑69 (MQ=255)
                              ccGTAGTCATTACGCACGCAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCACCTACCTGAGTTggg                    >  1:1449471/1‑69 (MQ=255)
                                      aTTACGCACGCAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCACCTACCTGAGTTGGGAATTCGtt            >  1:13738/1‑69 (MQ=255)
                                      aTTACGCACGCAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCACCTACCTGAGTTGGGAATTCGtt            >  1:8543/1‑69 (MQ=255)
                                      aTTACGCACGCAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCACCTACCTGAGTTGGGAATTCGtt            >  1:2341537/1‑69 (MQ=255)
                                      aTTACGCACGCAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCACCTACCTGAGTTGGGAATTCGtt            >  1:688983/1‑69 (MQ=255)
                                      aTTACGCACGCAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCACCTACCTGAGTTGGGAATTCGtt            >  1:231634/1‑69 (MQ=255)
                                      aTTACGCACGCAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCACCTACCTGAGTTGGGAATTCGtt            >  1:2880008/1‑69 (MQ=255)
                                      aTTACGCACGCAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCACCTACCTGAGTTGGGAATTCGtt            >  1:1608720/1‑69 (MQ=255)
                                      aTTACGCACGCAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCACCTACCTGAGTTGGGAATTCGtt            >  1:190093/1‑69 (MQ=255)
                                      aTTACGCACGCAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCACCTACCTGAGTTGGGAATTCGta            >  1:3248429/1‑68 (MQ=255)
                                         acgcacgcAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCACCTACCTGAGTTGGGAATTCGtt            >  1:2887839/1‑66 (MQ=255)
                                         acgcacgcAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCACCTACCTGAGTTGGGAATTCGtt            >  1:2501395/1‑66 (MQ=255)
                                         acgcacgcAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCACCTACCTGAGTTGGGAATTCGTtc           >  1:2127598/1‑67 (MQ=255)
                                         acgcacgcAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCACCTACCTGAGTTGGGAATTCGTtc           >  1:2056321/1‑67 (MQ=255)
                                         acgcacgcAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCACCTACCTGAGTTGGGAATTCGTtc           >  1:554177/1‑67 (MQ=255)
                                         acgcacgcAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCACCTACCTGAGTTGGGAATTCGTtc           >  1:1068948/1‑67 (MQ=255)
                                         acgcacgcAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCACCTACCTGAGTTGGGAATTCGTtc           >  1:1840170/1‑67 (MQ=255)
                                         acgcacgcAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCACCTACCTGAGTTGGGAATTCGTtc           >  1:645991/1‑67 (MQ=255)
                                             acgcaGCCAACAGAAAGATGCGGCGTACCCGAAGCATCACCTACCTGAGTTGGGAATTCGTTCTCTTCa     <  1:586602/69‑1 (MQ=255)
                                                caGCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGTCGGGAATTCGTTCTCTTCagag  <  1:104477/69‑1 (MQ=255)
                                                caGCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGTCGGGAATTCGTTCTCTTCagag  <  1:2942372/69‑1 (MQ=255)
                                                caGCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGTCGGGAATTCGTTCTCTTCagag  <  1:2946196/69‑1 (MQ=255)
                                                caGCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGTCGGGAATTCGTTCTCTTCagag  <  1:3345194/69‑1 (MQ=255)
                                                caGCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGTCGGGAATTCGTTCTCTTCagag  <  1:673860/69‑1 (MQ=255)
                                                caGCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGTCGGGAATTCGTTCTCTTCagag  <  1:83760/69‑1 (MQ=255)
                                                caGCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGTCGGGAATTCGTTCTCTTCagag  <  1:2937442/69‑1 (MQ=255)
                                                caGCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGTCGGGAATTCGTTCTCTTCagag  <  1:1780092/69‑1 (MQ=255)
                                                caGCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGTCGGGAATTCGTTCTCTTCagag  <  1:1723841/69‑1 (MQ=255)
                                                caGCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGTCGGGAATTCGTTCTCTTCagag  <  1:1720444/69‑1 (MQ=255)
                                                caGCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGTCGGGAATTCGTTCTCTTCagag  <  1:1637491/69‑1 (MQ=255)
                                                caGCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGTCGGGAATTCGTTCTCTTCagag  <  1:1286154/69‑1 (MQ=255)
                                                caGCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGTCGGGAATTCGTTCTCTTCagag  <  1:2279450/69‑1 (MQ=255)
                                                caGCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGTCGGGAATTCGTTCTCTTCagag  <  1:2368990/69‑1 (MQ=255)
                                                caGCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGTCGGGAATTCGTTCTCTTCagag  <  1:1051645/69‑1 (MQ=255)
                                                caGCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGTCGGGAATTCGTTCTCTTCagag  <  1:2636568/69‑1 (MQ=255)
                                                 aGCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGTCGGGAATTCGTTCTCTTCagag  <  1:2821913/68‑1 (MQ=255)
                                                  gCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGtc                       >  1:1241144/1‑45 (MQ=255)
                                                  gCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGtc                       >  1:881412/1‑45 (MQ=255)
                                                  gCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGtc                       >  1:83741/1‑45 (MQ=255)
                                                  gCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGtc                       >  1:894319/1‑45 (MQ=255)
                                                  gCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGtc                       >  1:959225/1‑45 (MQ=255)
                                                  gCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGtc                       >  1:748881/1‑45 (MQ=255)
                                                  gCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGtc                       >  1:287370/1‑45 (MQ=255)
                                                  gCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGtc                       >  1:1792673/1‑45 (MQ=255)
                                                  gCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGtc                       >  1:631463/1‑45 (MQ=255)
                                                  gCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGtc                       >  1:624308/1‑45 (MQ=255)
                                                  gCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGtc                       >  1:1963662/1‑45 (MQ=255)
                                                  gCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGtc                       >  1:2000350/1‑45 (MQ=255)
                                                  gCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGtc                       >  1:2122933/1‑45 (MQ=255)
                                                  gCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGtc                       >  1:400520/1‑45 (MQ=255)
                                                  gCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGtc                       >  1:2164771/1‑45 (MQ=255)
                                                  gCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGtc                       >  1:2222749/1‑45 (MQ=255)
                                                  gCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGtc                       >  1:2431650/1‑45 (MQ=255)
                                                  gCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGtc                       >  1:2431918/1‑45 (MQ=255)
                                                  gCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGtc                       >  1:2434335/1‑45 (MQ=255)
                                                  gCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGtc                       >  1:3013989/1‑45 (MQ=255)
                                                  gCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGtc                       >  1:276878/1‑45 (MQ=255)
                                                  gCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGt                        >  1:1893761/1‑45 (MQ=255)
                                                  gCCAACGGAAAGATGCGGCGTACCAGAAGCATCTCCTACCTGAGt                        >  1:914394/1‑45 (MQ=255)
                                                        |                                                            
ACCGACCAGAACTGGACTTGCTCCGGCATACCGTAGTCATTACGCACGCAGCCAACAGAAAGATGCGGCGTACCCGAAGCATCACCTACCTGAGTTGGGAATTCGTTCTCTTCAGAG  >  minE/1520967‑1521083

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: