Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,144,080 G→T 39.0% A117A (GCG→GCT argG → argininosuccinate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,144,0800GT39.0% 51.6 / 47.5 64A117A (GCG→GCTargGargininosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base G (19/20);  new base T (4/21);  total (23/41)
Fisher's exact test for biased strand distribution p-value = 8.77e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTATTTCAACACGACGCCGCTGGGCCGCGCCGTGACTGGTACCATGCTGGTTGCTGCGATGAAAGAAGATGGCGTGAATATCTGGGGTGACGGTAGCACCTACAAAGGAAACGATATCG  >  minE/2144023‑2144141
                                                         |                                                             
cTATTTCAACACGACGCCGCTGGGCCGCGCCGTGACTGGTACCATGCTGGTTGTTGCGATGAaagaa                                                      <  1:1244348/67‑1 (MQ=255)
cTATTTCAACACGACGCCGCTGGGCCGCGCCGTGACTGGTACCATGCTGGTTGCTGCGTTGAaagaa                                                      <  1:1974269/67‑1 (MQ=255)
cTATTTCAACACGACGCCGCTGGGCCGCGCCGTGACTGGTACCATGCTGGTTGCTGCGATGAaagaa                                                      <  1:773254/67‑1 (MQ=255)
cTATTTCAACACGACGCCGCTGGGCCGCGCCGTGACTGGTACCATGCTGGTTGCTGCGATGAaagaa                                                      <  1:3009025/67‑1 (MQ=255)
cTATTTCAACACGACGCCGCTGGGCCGCGCCGTGACTGGTACCATGCTGGTTGCTGCGATGAaagaa                                                      <  1:2652208/67‑1 (MQ=255)
cTATTTCAACACGACGCCGCTGGGCCGCGCCGTGACTGGTACCATGCTGGTTGCTGCGATGAaagaa                                                      <  1:1551113/67‑1 (MQ=255)
     tCAACACGACGCCGCTGGGCCGCGCCGTGACTGGTACCATGCTGGTTGCTGCGATGAaagaa                                                      <  1:2230357/62‑1 (MQ=255)
     tCAACACGACGCCGCTGGGCCGCGCCGTGACTGGTACCATGCTGGTTGCTGCGATGAAAGAAGATGGc                                                >  1:3190904/1‑68 (MQ=255)
     tCAACACGACGCCGCTGGGCCGCGCCGTGACCGGCACCATGCTGGTTGCTGCTATGAaag                                                        <  1:1624884/60‑1 (MQ=255)
     tCAACACGACGCCGCTGGGCCGCGCCGTGACCGGCACCATGCTGGTTGCTGCTATGAaag                                                        <  1:2011778/60‑1 (MQ=255)
     tCAACACGACGCCGCTGGGCCGCGCCGTGACCGGCACCATGCTGGTTGCTGCTATGAaag                                                        <  1:3076906/60‑1 (MQ=255)
     tCAACACGACGCCGCTGGGCCGCGCCGTGACCGGCACCATGCTGGTTGCTGCTATGAaag                                                        <  1:1797970/60‑1 (MQ=255)
     tCAACACGACGCCGCTGGGCCGCGCCGTGACCGGCACCATGCTGGTTGCTGCTATGAaag                                                        <  1:877111/60‑1 (MQ=255)
     tCAACACGACGCCGCTGGGCCGCGCCGTGACCGGCACCATGCTGGTTGCTGCTATGAaag                                                        <  1:2380063/60‑1 (MQ=255)
        acacGACGCCGCTGGGCCGCGCCGTGACTGGTACCATGCTGGTTGCTGCGATGAAAGAAGATGGCGTg                                             <  1:2844686/68‑1 (MQ=255)
        acacGACGCCGCTGGGCCGCGCCGTGACTGGTACCATGCTGGTTGCTGCGATGAAAGAAGATGGCGTg                                             <  1:283399/68‑1 (MQ=255)
        acacGACGCCGCTGGGCCGCGCCGTGACTGGTACCATGCTGGTTGCTGCGATGAAAGAAGATGGCGTg                                             <  1:294338/68‑1 (MQ=255)
         cacGACGCCGCTGGGCCGCGCCGTGACTGGTACCATGCTGGTTGCTGCGATGAAAGAAGATGGCGTg                                             <  1:244073/67‑1 (MQ=255)
                 cgcTGGGCCGCGCCGTGACCGGCACCATGCTGGTTGCTGCTATGAAAGAAGATGGCGTGAATATCTggg                                   <  1:1964643/69‑1 (MQ=255)
                 cgcTGGGCCGCGCCGTGACCGGCACCATGCTGGTTGCTGCTATGAAAGAAGATGGCGTGAATATCTggg                                   <  1:2199669/69‑1 (MQ=255)
                 cgcTGGGCCGCGCCGTGACCGGCACCATGCTGGTTGCTGCTATGAAAGAAGATGGCGTGAATATCTggg                                   <  1:2545385/69‑1 (MQ=255)
                 cgcTGGGCCGCGCCGTGACCGGCACCATGCTGGTTGCTGCTATGAAAGAAGATGGCGTGAATATCTggg                                   <  1:2591391/69‑1 (MQ=255)
                 cgcTGGGCCGCGCCGTGACCGGCACCATGCTGGTTGCTGCTATGAAAGAAGATGGCGTGAATATCTggg                                   <  1:325636/69‑1 (MQ=255)
                 cgcTGGGCCGCGCCGTGACCGGCACCATGCTGGTTGCTGCTATGAAAGAAGATGGCGTGAATATCTggg                                   <  1:855328/69‑1 (MQ=255)
                 cgcTGGGCCGCGCCGTGACCGGCACCATGCTGGTTGCTGCTATGAAAGAAGATGGCGTGAATATCTggg                                   <  1:3432133/69‑1 (MQ=255)
                 cgcTGGGCCGCGCCGTGACCGGCACCATGCTGGTTGCTGCTATGAAAGAAGATGGCGTGAATATCTggg                                   <  1:141164/69‑1 (MQ=255)
                 cgcTGGGCCGCGCCGTGACCGGCACCATGCTGGTTGCTGCTATGAAAGAAGATGGCGTGAATATCTggg                                   <  1:998510/69‑1 (MQ=255)
                 cgcTGGGCCGCGCCGTGACCGGCACCATGCTGGTTGCTGCTATGAAAGAAGATGGCGTGAATATCTggg                                   <  1:1917956/69‑1 (MQ=255)
                 cgcTGGGCCGCGCCGTGACCGGCACCATGCTGGTTGCTGCTATGAAAGAAGATGGCGTAAATATCTggg                                   <  1:1961665/69‑1 (MQ=255)
                  gcTTGGCCGCGCCGTGACCGGCACCATGCTGGTTGCTGCTATGAAAGAAGATGGCGTGAATATCTggg                                   <  1:2511519/68‑1 (MQ=255)
                     gggCCGCGCCGTGACTGGTACCATGCTGGTTGCTGCGATGAAAGAAGATGGCGTGAATATCTGGGGTGa                               >  1:2068318/1‑69 (MQ=255)
                      ggCCGCGCCGTGACCGGCACCATGCTGGTTGCTGCTATGAtag                                                        <  1:1716358/43‑1 (MQ=25)
                           cgcCGTGACCGGCACCATGCTGGTTGCTGCTATGAaaga                                                       >  1:1206981/1‑39 (MQ=37)
                           cgcCGTGACCGGCACCATGCTGGTTGCTGCTATGAAAGAAGATGGCGTg                                             >  1:972206/1‑49 (MQ=255)
                           cgcCGTGACCGGCACCATGCTGGTTGCTGCTATGAAAGAAGATGGCGTGAATATCt                                      >  1:224418/1‑56 (MQ=255)
                           cgcCGTGACCGGCACCATGCTGGTTGCTGCTATGAAAGAAGATGGCGTGAATATCTGGGGTgctg                             >  1:1689584/1‑62 (MQ=255)
                            gcCGTGACTGGTACCATGCTGGTTGCTGCGATGAAAGAAGATGGCGTGAATATCTGGGGTGACGGTAGc                        <  1:2850778/69‑1 (MQ=255)
                            gcCGTGACTGGTACCATGCTGGTTGCTGCGATGAAAGAAGATGGCGTGAATATCTGGGGTGACGGTAGc                        <  1:584123/69‑1 (MQ=255)
                            gcCGTGACTGGTACCATGCTGGTTGCTGCGATGAAAGAAGATGGCGTGAATATCTGGGGTGACGGTAGc                        <  1:2477991/69‑1 (MQ=255)
                            gcCGTGACTGGTACCATGCTGGTTGCTGCGATGAAAGAAGATGGCGTGAATATCTGGGGTGACGGTAGc                        <  1:2489416/69‑1 (MQ=255)
                             ccGTGACTGGTACCATGCTGGTTGCTGCGATGAAAGAAGATGGCGTGAATATCTGGGGTGACGGTAGc                        <  1:226024/68‑1 (MQ=255)
                             ccGTGACTGGTACCATGCTGGTTGCTGCGATGAAAGAAGATGGCGTGAATATCTGGGGTGACGGTAGc                        <  1:1873490/68‑1 (MQ=255)
                                  aCTGGTACCATGCTGGTTGCTGCGATGAAAGAAGATGGCGTGAATATCTGGGGTGACgg                            >  1:127178/1‑59 (MQ=255)
                                  aCTGGTACCATGCTGGTTGCTGCGATGAAAGAAGATGGCGTGAATATCTGGGGTGACGGTAGCACCTa                   >  1:406478/1‑68 (MQ=255)
                                  aCTGGTACCATGCTGGTTGCTGCGATGAAAGAAGATGGCGTGAATATCTGGGGTGACGGTAGCACCTa                   >  1:2177696/1‑68 (MQ=255)
                                  aCTGGTACCATGCTGGTTGCTGCGATGAAAGAAGATGGCGTGAATATCTGGGGTGACGGTAGCACCTa                   >  1:3249592/1‑68 (MQ=255)
                                           aTGCTGGTTGCTGCGATGAAAGAAGATGGCGTGAATATCTGGGGTGACGGTAGCACCTACAAAGGaaa          >  1:1182502/1‑68 (MQ=255)
                                           aTGCTGGTTGCTGCGATGAAAGAAGATGGCGTGAATATCTGGGGTGACGGTAGCACCTACAAAGGAAAc         >  1:1833979/1‑69 (MQ=255)
                                           aTGCTGGTTGCTGCGATGAAAGAAGATGGCGTGAATATCTGGGGTGACGGTAGCACCTACAAAGGAAAc         >  1:1082356/1‑69 (MQ=255)
                                           aTGCTGGTTGCTGCGATGAAAGAAGATGGCGTGAATATCTGGGGTGACGGTAGCACCTACAAAGGAAAc         >  1:1301194/1‑69 (MQ=255)
                                           aTGCTGGTTGCTGCGATGAAAGAAGATGGCGTGAATATCTGGGGTGACGGTAGCACCTACAAAGGAAAc         >  1:1829835/1‑69 (MQ=255)
                                           aTGCTGGTTGCTGCGATGAAAGAAGATGGCGTGAATATCTGGGGTGACGGTAGCACCTACAAAGGAAAc         >  1:2360237/1‑69 (MQ=255)
                                              cTGGTTGCTGCTATGAAAGAAGATGGCGTGAATATCTggt                                   <  1:2677028/40‑2 (MQ=39)
                                              cTGGTTGCTGCGATGAAAGAAGATGGCGTGAATATCTGGGGTGACGGTAGCACCTACAAAGGAAACGat      >  1:2107671/1‑69 (MQ=255)
                                              cTGGTTGCTGCGATGAAAGAAGATGGCGTGAATATCTGGGGTGACGGTAGCACCTACAAAGGAAACGat      >  1:153060/1‑69 (MQ=255)
                                              cTGGTTGCTGCGATGAAAGAAGATGGCGTGAATATCTGGGGTGACGGTAGCACCTACAAAGGAAACGat      >  1:3261653/1‑69 (MQ=255)
                                              cTGGTTGCTGCGATGAAAGAAGATGGCGTGAATATCTGGGGTGACGGTAGCACCTACAAAGGAAACGat      >  1:2276985/1‑69 (MQ=255)
                                              cTGGTTGCTGCGATGAAAGAAGATGGCGTGAATATCTGGGGTGACGGTAGCACCTACAAAGGAAACGat      >  1:92747/1‑69 (MQ=255)
                                              cTGGTTGCTGCGATGAAAGAAGATGGCGTGAATATCTGGGGTGACGGTAGCACCTACAAAGGAAACGat      >  1:2676271/1‑69 (MQ=255)
                                              cTGGTTGCTGCGATGAAAGAAGATGGCGTGAATATCTGGGGTGACGGTAGCACCTACAAAGGAAACGa       >  1:1181627/1‑68 (MQ=255)
                                                  ttGCTGCTATGAAAGAAGATGGCGTGAATATCTggt                                   <  1:2232346/36‑2 (MQ=38)
                                                   tgctgcGATGAAAGAAGATGGCGTGAATATCTGGGGTGACGGTAGCACCTACAAAGGAAACGATATCg  <  1:2795621/68‑1 (MQ=255)
                                                   tgctgcGATGAAAGAAGATGGCGTGAATATCTGGGGTGACGGTAGCACCTACAAAGGAAACGATATCg  <  1:2302541/68‑1 (MQ=255)
                                                   tgctgcGATGAAAGAAGATGGCGTGAATATCTGGGGTGACGGTAGCACCTACAAAGGAAACGATATCg  <  1:1453128/68‑1 (MQ=255)
                                                         |                                                             
CTATTTCAACACGACGCCGCTGGGCCGCGCCGTGACTGGTACCATGCTGGTTGCTGCGATGAAAGAAGATGGCGTGAATATCTGGGGTGACGGTAGCACCTACAAAGGAAACGATATCG  >  minE/2144023‑2144141

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