Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,153,806 C→T 70.1% S270N (AGT→AAT)  astE ← succinylglutamate desuccinylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,153,8060CT70.1% 42.8 / 20.9 50S270N (AGT→AAT) astEsuccinylglutamate desuccinylase
Reads supporting (aligned to +/- strand):  ref base C (2/13);  new base T (27/8);  total (29/21)
Fisher's exact test for biased strand distribution p-value = 7.27e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGAT  >  minE/1153744‑1153863
                                                              |                                                         
gTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCa                                                     <  1:644422/69‑1 (MQ=255)
gTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCa                                                     <  1:2591555/69‑1 (MQ=255)
gTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCa                                                     <  1:1431788/69‑1 (MQ=255)
gTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCa                                                     <  1:3251575/69‑1 (MQ=255)
gTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGACATAAAATTCAGCGTGTCACTTGCCa                                                     <  1:2635380/69‑1 (MQ=255)
               gcgcCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAg                                       <  1:43965/68‑1 (MQ=255)
                             ccTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACggcgagtggcg                         >  1:1293532/1‑68 (MQ=255)
                                  ctcAAACGGCATAAAATTCAGCGTGTCTTTTGCCATATGCATTTc                                           >  1:1370526/1‑45 (MQ=38)
                                  ctcAAACGGCATAAAATTCAGCGTGTCATTTGCCat                                                    >  1:3296033/1‑36 (MQ=255)
                                  ctcAAACGGCATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACggcga                               >  1:2269079/1‑57 (MQ=255)
                                  ctcAAACGGCATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACggcga                               >  1:1009051/1‑57 (MQ=255)
                                  ctcAAACGGCATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACggcga                               >  1:2449735/1‑57 (MQ=255)
                                  ctcAAACGGCATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACggcga                               >  1:2624268/1‑57 (MQ=255)
                                  ctcAAACGGCATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACggcga                               >  1:2832692/1‑57 (MQ=255)
                                  ctcAAACGGCATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACggcga                               >  1:2874219/1‑57 (MQ=255)
                                  ctcAAACGGCATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACggcga                               >  1:2943769/1‑57 (MQ=255)
                                  ctcAAACGGCATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACggcga                               >  1:325557/1‑57 (MQ=255)
                                  ctcAAACGGCATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACggcga                               >  1:3351662/1‑57 (MQ=255)
                                  ctcAAACGGCATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACggcga                               >  1:3378165/1‑57 (MQ=255)
                                  ctcAAACGGCATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACggcga                               >  1:483061/1‑57 (MQ=255)
                                  ctcAAACGGCATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACggcga                               >  1:604009/1‑57 (MQ=255)
                                  ctcAAACGGCATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACggcga                               >  1:892976/1‑57 (MQ=255)
                                  ctcAAACGGCATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACggcga                               >  1:225950/1‑57 (MQ=255)
                                  ctcAAACGGCATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACggcga                               >  1:2258312/1‑57 (MQ=255)
                                  ctcAAACGGCATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACggcga                               >  1:2249430/1‑57 (MQ=255)
                                  ctcAAACGGCATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACggcga                               >  1:2113416/1‑57 (MQ=255)
                                  ctcAAACGGCATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACggcga                               >  1:1898986/1‑57 (MQ=255)
                                  ctcAAACGGCATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACggcga                               >  1:1000109/1‑57 (MQ=255)
                                  ctcAAACGGCATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACggcga                               >  1:1877620/1‑57 (MQ=255)
                                  ctcAAACGGCATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACggcga                               >  1:1817186/1‑57 (MQ=255)
                                  ctcAAACGGCATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACggcga                               >  1:1714872/1‑57 (MQ=255)
                                  ctcAAACGGCATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACggcga                               >  1:1364187/1‑57 (MQ=255)
                                  ctcAAACGGCATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACggcga                               >  1:137917/1‑57 (MQ=255)
                                  ctcAAACGGCATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACggcg                                >  1:1491413/1‑56 (MQ=255)
                                           cATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACGGCTAGTGGCGAGTAATTTGCGAAac          <  1:143110/69‑1 (MQ=255)
                                           cATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAac          <  1:941126/69‑1 (MQ=255)
                                           cATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAac          <  1:80140/69‑1 (MQ=255)
                                           cATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAac          <  1:489173/69‑1 (MQ=255)
                                           cATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAac          <  1:2827045/69‑1 (MQ=255)
                                           cATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAac          <  1:2755554/69‑1 (MQ=255)
                                           cATAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAac          <  1:2299461/69‑1 (MQ=255)
                                            aTAAAATTCAGCGTGTCATTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAac          <  1:1577451/68‑1 (MQ=255)
                                            aTAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAacc         >  1:1124219/1‑69 (MQ=255)
                                                   tCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGAt  <  1:2958784/69‑1 (MQ=255)
                                                   tCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGAt  <  1:3379070/69‑1 (MQ=255)
                                                   tCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGAt  <  1:411532/69‑1 (MQ=255)
                                                   tCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGAt  <  1:1267461/69‑1 (MQ=255)
                                                   tCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGAt  <  1:1216820/69‑1 (MQ=255)
                                                   tCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGAt  <  1:2903802/69‑1 (MQ=255)
                                                   tCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGAt  <  1:1068130/69‑1 (MQ=255)
                                                              |                                                         
GTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGAT  >  minE/1153744‑1153863

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: