Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,302,610 G→T 55.0% A208A (GCC→GCA metB ← cystathionine gamma‑synthase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,302,6100GT55.0% 5.7 / 91.6 80A208A (GCC→GCAmetBcystathionine gamma‑synthase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (18/18);  new base T (29/15);  total (47/33)
Fisher's exact test for biased strand distribution p-value = 1.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

ACGCCAATATTGTTTGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCT  >  minE/2302542‑2302673
                                                                    |                                                               
aCGCCAATATTGTTTGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCg                                                                 >  1:2887750/1‑69 (MQ=255)
aCGCCAATATTGTTTGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCg                                                                 >  1:2806983/1‑69 (MQ=255)
aCGCCAATATTGTTTGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCg                                                                 >  1:1762133/1‑69 (MQ=255)
aCGCCAATATTGTTTGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCg                                                                 >  1:2298690/1‑69 (MQ=255)
  gCCAATATTGTTTGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGc                                                               <  1:1193303/69‑1 (MQ=255)
       tatTGTTTGCCCACCAGGCAAGTTCAGTGACAACGTCGGGGTCTTTAGCAATCACTACGCCTGCCAcca                                                          >  1:1975364/1‑67 (MQ=255)
       tatTGTTTGCCCACCAGGCAAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACTACGCCTGCCAcca                                                          >  1:2198634/1‑67 (MQ=255)
       tatTGTTTGCCCACCAGGCAAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACTACGCCTGCCAcca                                                          >  1:525298/1‑67 (MQ=255)
       tatTGTTTGCCCACCAGGCAAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACTACGCCTGCCAcca                                                          >  1:2122374/1‑67 (MQ=255)
       tatTGTTTGCCCACCAGGCAAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACTACGCCTGCCAcca                                                          >  1:514233/1‑67 (MQ=255)
       tatTGTTTGCCCACCAGGCAAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACTACGCCTGCCAcca                                                          >  1:2275590/1‑67 (MQ=255)
       tatTGTTTGCCCACCAGGCAAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACTACGCCTGCCAcca                                                          >  1:1751196/1‑67 (MQ=255)
       tatTGTTTGCCCACCAGGCAAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACTACGCCTGCCAcca                                                          >  1:668695/1‑67 (MQ=255)
       tatTGTTTGCCCACCAGGCAAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACTACGCCTGCCAcca                                                          >  1:3317897/1‑67 (MQ=255)
       tatTGTTTGCCCACCAGGCAAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACTACGCCTGCCAcca                                                          >  1:1158236/1‑67 (MQ=255)
       tatTGTTTGCCCACCAGGCAAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACTACGCCTGCCAcca                                                          >  1:835486/1‑67 (MQ=255)
       tatTGTTTGCCCACCAGGCAAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACTACGCCTGCCAcca                                                          >  1:3359634/1‑67 (MQ=255)
       tatTGTTTGCCCACCAGGCAAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACTACGCCTGCCAcc                                                           >  1:3190145/1‑67 (MQ=255)
       tatTGTTTGCCCACCAGGCAAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACTACGCCTGCCAcc                                                           >  1:1831937/1‑67 (MQ=255)
        atTGTTTGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTa                                                          >  1:608183/1‑68 (MQ=255)
        atTGTTTGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTa                                                          >  1:1279625/1‑68 (MQ=255)
        atTGTTTGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTAc                                                         >  1:3400173/1‑69 (MQ=255)
        atTGTTTGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTAc                                                         >  1:1877847/1‑69 (MQ=255)
        atTGTTTGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTAc                                                         >  1:1522647/1‑69 (MQ=255)
         ttGTTTGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACg                                                        >  1:959847/1‑69 (MQ=255)
              tGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTg                                                    <  1:1551276/68‑1 (MQ=255)
              tGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTg                                                    <  1:672303/68‑1 (MQ=255)
              tGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTg                                                    <  1:2003141/68‑1 (MQ=255)
              tGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTg                                                    <  1:1068124/68‑1 (MQ=255)
                           aGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGa                                      <  1:1896637/69‑1 (MQ=255)
                           aGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGa                                      <  1:2128774/69‑1 (MQ=255)
                                   gACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCg                               >  1:3386673/1‑68 (MQ=255)
                                    aCAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTg                             >  1:2361909/1‑69 (MQ=255)
                                    aCAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTg                             >  1:2505632/1‑69 (MQ=255)
                                    aCAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTg                             >  1:2214748/1‑69 (MQ=255)
                                    aCAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTg                             >  1:570209/1‑69 (MQ=255)
                                    aCAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTg                             >  1:85533/1‑69 (MQ=255)
                                    aCAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTg                             >  1:1294200/1‑69 (MQ=255)
                                    aCAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTg                             >  1:809708/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:1007093/69‑1 (MQ=255)
                                          tCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGaa                       <  1:2474539/69‑1 (MQ=255)
                                          tCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGaa                       <  1:1470227/69‑1 (MQ=255)
                                          tCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGaa                       <  1:365693/69‑1 (MQ=255)
                                          tCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGaa                       <  1:1985153/69‑1 (MQ=255)
                                              ggTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGc                    <  1:1187981/68‑1 (MQ=255)
                                              ggTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGc                    <  1:2755793/68‑1 (MQ=255)
                                              ggTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGc                    <  1:864537/68‑1 (MQ=255)
                                                           cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      <  1:340119/69‑1 (MQ=255)
                                                           cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      <  1:3322982/69‑1 (MQ=255)
                                                           cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      <  1:575625/69‑1 (MQ=255)
                                                           cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      <  1:1012483/69‑1 (MQ=255)
                                                           cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      <  1:538081/69‑1 (MQ=255)
                                                           cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      <  1:2856865/69‑1 (MQ=255)
                                                           cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      <  1:2107339/69‑1 (MQ=255)
                                                           cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      <  1:2217494/69‑1 (MQ=255)
                                                           cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      <  1:2305187/69‑1 (MQ=255)
                                                           cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      <  1:1154749/69‑1 (MQ=255)
                                                           cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      <  1:2428135/69‑1 (MQ=255)
                                                           cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      <  1:1126621/69‑1 (MQ=255)
                                                           cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      <  1:1061382/69‑1 (MQ=255)
                                                           cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      <  1:2707068/69‑1 (MQ=255)
                                                           cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      <  1:2754873/69‑1 (MQ=255)
                                                           caccacGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      >  1:857440/1‑69 (MQ=255)
                                                           caccacGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      >  1:2830562/1‑69 (MQ=255)
                                                           caccacGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      >  1:1955804/1‑69 (MQ=255)
                                                           caccacGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      >  1:751189/1‑69 (MQ=255)
                                                               acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt  >  1:136818/1‑69 (MQ=255)
                                                               acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt  >  1:1436260/1‑69 (MQ=255)
                                                               acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt  >  1:999427/1‑69 (MQ=255)
                                                               acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt  >  1:735527/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:1686191/1‑69 (MQ=255)
                                                               acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt  >  1:2161636/1‑69 (MQ=255)
                                                               acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt  >  1:416201/1‑69 (MQ=255)
                                                               acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt  >  1:2330078/1‑69 (MQ=255)
                                                               acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt  >  1:3388987/1‑69 (MQ=255)
                                                               acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt  >  1:2555854/1‑69 (MQ=255)
                                                               acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt  >  1:2671489/1‑69 (MQ=255)
                                                               acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt  >  1:3088166/1‑69 (MQ=255)
                                                                    |                                                               
ACGCCAATATTGTTTGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCT  >  minE/2302542‑2302673

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