Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,302,616 T→C 46.2% V206V (GTA→GTG metB ← cystathionine gamma‑synthase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,302,6160TC46.2% 4.2 / 54.3 65V206V (GTA→GTGmetBcystathionine gamma‑synthase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (18/17);  new base C (15/15);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGTTTGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCT  >  minE/2302550‑2302673
                                                                  |                                                         
atTGTTTGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTa                                                          >  1:1279625/1‑68 (MQ=255)
atTGTTTGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTa                                                          >  1:608183/1‑68 (MQ=255)
atTGTTTGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTAc                                                         >  1:1877847/1‑69 (MQ=255)
atTGTTTGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTAc                                                         >  1:3400173/1‑69 (MQ=255)
atTGTTTGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTAc                                                         >  1:1522647/1‑69 (MQ=255)
 ttGTTTGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACg                                                        >  1:959847/1‑69 (MQ=255)
      tGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTg                                                    <  1:2003141/68‑1 (MQ=255)
      tGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTg                                                    <  1:1068124/68‑1 (MQ=255)
      tGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTg                                                    <  1:1551276/68‑1 (MQ=255)
      tGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTg                                                    <  1:672303/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:809708/1‑69 (MQ=255)
                            aCAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTg                             >  1:1294200/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:2505632/1‑69 (MQ=255)
                            aCAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTg                             >  1:2361909/1‑69 (MQ=255)
                                cGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATg                         <  1:1598652/69‑1 (MQ=255)
                                cGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATg                         <  1:11469/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:1985153/69‑1 (MQ=255)
                                  tCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGaa                       <  1:365693/69‑1 (MQ=255)
                                  tCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGaa                       <  1:1470227/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:1012483/69‑1 (MQ=255)
                                                   cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      <  1:2754873/69‑1 (MQ=255)
                                                   cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      <  1:2856865/69‑1 (MQ=255)
                                                   cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      <  1:3322982/69‑1 (MQ=255)
                                                   cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      <  1:340119/69‑1 (MQ=255)
                                                   cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      <  1:538081/69‑1 (MQ=255)
                                                   cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      <  1:575625/69‑1 (MQ=255)
                                                   cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      <  1:1061382/69‑1 (MQ=255)
                                                   cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      <  1:2428135/69‑1 (MQ=255)
                                                   cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      <  1:1126621/69‑1 (MQ=255)
                                                   cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      <  1:2305187/69‑1 (MQ=255)
                                                   cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      <  1:2217494/69‑1 (MQ=255)
                                                   cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      <  1:1154749/69‑1 (MQ=255)
                                                   cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      <  1:2107339/69‑1 (MQ=255)
                                                   cactacGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      <  1:2707068/69‑1 (MQ=255)
                                                   caccacGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      >  1:2830562/1‑69 (MQ=255)
                                                   caccacGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      >  1:857440/1‑69 (MQ=255)
                                                   caccacGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      >  1:1955804/1‑69 (MQ=255)
                                                   caccacGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGc      >  1:751189/1‑69 (MQ=255)
                                                       acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt  >  1:1436260/1‑69 (MQ=255)
                                                       acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt  >  1:416201/1‑69 (MQ=255)
                                                       acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt  >  1:735527/1‑69 (MQ=255)
                                                       acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt  >  1:136818/1‑69 (MQ=255)
                                                       acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt  >  1:999427/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:3388987/1‑69 (MQ=255)
                                                       acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt  >  1:3088166/1‑69 (MQ=255)
                                                       acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt  >  1:2330078/1‑69 (MQ=255)
                                                       acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt  >  1:2555854/1‑69 (MQ=255)
                                                       acGCCTGCCACCACGTCAGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCt  >  1:2671489/1‑69 (MQ=255)
                                                                  |                                                         
ATTGTTTGCCCACCAGGCCAGTTCAGTGACAACGTCCGGGTCTTTAGCAATCACCACGCCGGCCACTACGTCTGAGTGACCGTTCAGATATTTCGTGCATGAATGCAACACCAGATCGGCACCT  >  minE/2302550‑2302673

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: