Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,302,622 T→A 43.1% S204S (TCA→TCT metB ← cystathionine gamma‑synthase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,302,6220TA43.1% 23.8 / 74.9 72S204S (TCA→TCTmetBcystathionine gamma‑synthase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (24/17);  new base A (15/16);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.57e-01

TGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCTAATGCCAGCGGA  >  minE/2302556‑2302685
                                                                  |                                                               
tGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTg                                                                <  1:1551276/68‑1 (MQ=255)
tGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTg                                                                <  1:1068124/68‑1 (MQ=255)
tGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTg                                                                <  1:672303/68‑1 (MQ=255)
tGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTg                                                                <  1:2003141/68‑1 (MQ=255)
             aGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGa                                                  <  1:2128774/69‑1 (MQ=255)
             aGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGa                                                  <  1:1896637/69‑1 (MQ=255)
                     gACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCg                                           >  1:3386673/1‑68 (MQ=255)
                      aCAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTg                                         >  1:2214748/1‑69 (MQ=255)
                      aCAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTg                                         >  1:2505632/1‑69 (MQ=255)
                      aCAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTg                                         >  1:570209/1‑69 (MQ=255)
                      aCAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTg                                         >  1:1294200/1‑69 (MQ=255)
                      aCAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTg                                         >  1:2361909/1‑69 (MQ=255)
                      aCAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTg                                         >  1:809708/1‑69 (MQ=255)
                      aCAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTg                                         >  1:85533/1‑69 (MQ=255)
                          cGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATg                                     <  1:11469/69‑1 (MQ=255)
                          cGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATg                                     <  1:1598652/69‑1 (MQ=255)
                          cGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTAAGATATTTCGTGCATg                                     <  1:2831867/69‑1 (MQ=255)
                            tCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGaa                                   <  1:1985153/69‑1 (MQ=255)
                            tCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGaa                                   <  1:2474539/69‑1 (MQ=255)
                            tCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGaa                                   <  1:1007093/69‑1 (MQ=255)
                            tCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGaa                                   <  1:1470227/69‑1 (MQ=255)
                            tCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGaa                                   <  1:365693/69‑1 (MQ=255)
                                ggTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGc                                <  1:2755793/68‑1 (MQ=255)
                                ggTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGc                                <  1:1187981/68‑1 (MQ=255)
                                ggTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGc                                <  1:864537/68‑1 (MQ=255)
                                             cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc                  <  1:2856865/69‑1 (MQ=255)
                                             cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc                  <  1:2754873/69‑1 (MQ=255)
                                             cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc                  <  1:3322982/69‑1 (MQ=255)
                                             cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc                  <  1:2707068/69‑1 (MQ=255)
                                             cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc                  <  1:340119/69‑1 (MQ=255)
                                             cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc                  <  1:2428135/69‑1 (MQ=255)
                                             cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc                  <  1:538081/69‑1 (MQ=255)
                                             cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc                  <  1:575625/69‑1 (MQ=255)
                                             cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc                  <  1:2305187/69‑1 (MQ=255)
                                             cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc                  <  1:2217494/69‑1 (MQ=255)
                                             cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc                  <  1:1012483/69‑1 (MQ=255)
                                             cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc                  <  1:1061382/69‑1 (MQ=255)
                                             cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc                  <  1:2107339/69‑1 (MQ=255)
                                             cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc                  <  1:1126621/69‑1 (MQ=255)
                                             cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc                  <  1:1154749/69‑1 (MQ=255)
                                             caccacGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc                  >  1:857440/1‑69 (MQ=255)
                                             caccacGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc                  >  1:1955804/1‑69 (MQ=255)
                                             caccacGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc                  >  1:751189/1‑69 (MQ=255)
                                             caccacGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc                  >  1:2830562/1‑69 (MQ=255)
                                                 acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt              >  1:3088166/1‑69 (MQ=255)
                                                 acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt              >  1:999427/1‑69 (MQ=255)
                                                 acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt              >  1:735527/1‑69 (MQ=255)
                                                 acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt              >  1:416201/1‑69 (MQ=255)
                                                 acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt              >  1:136818/1‑69 (MQ=255)
                                                 acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt              >  1:3388987/1‑69 (MQ=255)
                                                 acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt              >  1:1436260/1‑69 (MQ=255)
                                                 acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt              >  1:2161636/1‑69 (MQ=255)
                                                 acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt              >  1:1502142/1‑69 (MQ=255)
                                                 acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt              >  1:1551920/1‑69 (MQ=255)
                                                 acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt              >  1:1577457/1‑69 (MQ=255)
                                                 acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt              >  1:2330078/1‑69 (MQ=255)
                                                 acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt              >  1:1686191/1‑69 (MQ=255)
                                                 acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt              >  1:2555854/1‑69 (MQ=255)
                                                 acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt              >  1:2671489/1‑69 (MQ=255)
                                                             aCGTCTGAGTGACCGTTCAGATGTTTCGTGCATGAATGCAACACCAGATCGGCACCTAATGCCAGCGGa  >  1:3031398/1‑69 (MQ=255)
                                                             aCGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCTAATGCCAGCGGa  >  1:2117792/1‑69 (MQ=255)
                                                             aCGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCTAATGCCAGCGGa  >  1:2315317/1‑69 (MQ=255)
                                                             aCGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCTAATGCCAGCGGa  >  1:185319/1‑69 (MQ=255)
                                                             aCGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCTAATGCCAGCGGa  >  1:115063/1‑69 (MQ=255)
                                                             aCGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCTAATGCCAGCGGa  >  1:1312527/1‑69 (MQ=255)
                                                             aCGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCTAATGCCAGCGGa  >  1:132811/1‑69 (MQ=255)
                                                             aCGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCTAATGCCAGCGGa  >  1:1627313/1‑69 (MQ=255)
                                                             aCGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCTAATGCCAGCGGa  >  1:2730513/1‑69 (MQ=255)
                                                             aCGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCTAATGCCAGCGGa  >  1:1430032/1‑69 (MQ=255)
                                                             aCGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCTAATGCCAGCGGa  >  1:1502899/1‑69 (MQ=255)
                                                             aCGTCTGAGTGACCGTTCAGATATTTAGTGCATGAATGCAACACCAGATCGGCACCTAATGCCAGCGGa  >  1:1499971/1‑69 (MQ=255)
                                                             aCGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc                  <  1:120472/53‑1 (MQ=255)
                                                                  |                                                               
TGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCTAATGCCAGCGGA  >  minE/2302556‑2302685

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: