Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 157,729 A→G 34.2% L84L (TTG→CTG)  pfs ← 5'‑methylthioadenosine/S‑adenosylhomocysteine nucleosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE157,7290AG34.2% 50.3 / 40.4 73L84L (TTG→CTG) pfs5'‑methylthioadenosine/S‑adenosylhomocysteine nucleosidase
Reads supporting (aligned to +/- strand):  ref base A (22/26);  new base G (5/20);  total (27/46)
Fisher's exact test for biased strand distribution p-value = 4.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.01e-01

ACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAAC  >  minE/157661‑157793
                                                                    |                                                                
aCCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCaa                                                                  <  1:2974186/69‑1 (MQ=255)
aCCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCaa                                                                  <  1:20468/69‑1 (MQ=255)
   aaaTGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGt                                                               <  1:2863678/69‑1 (MQ=255)
   aaaTGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGt                                                               <  1:1195611/69‑1 (MQ=255)
   aaaTGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGt                                                               <  1:2707928/69‑1 (MQ=255)
   aaaTGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGt                                                               <  1:1958100/69‑1 (MQ=255)
    aaTGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGtt                                                              <  1:1390358/69‑1 (MQ=255)
                  cgcgTCGTTATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGt                                                               <  1:2079404/54‑1 (MQ=255)
                  cgcgTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGt                                                               <  1:1001551/54‑1 (MQ=255)
                  cgcgTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGt                                                               <  1:3235663/54‑1 (MQ=255)
                  cgcgTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGt                                                               <  1:3112329/54‑1 (MQ=255)
                  cgcgTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGt                                                               <  1:3319510/54‑1 (MQ=255)
                  cgcgTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGt                                                               <  1:2196006/54‑1 (MQ=255)
                  cgcgTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGt                                                               <  1:1820346/54‑1 (MQ=255)
                  cgcgTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGt                                                               <  1:1679402/54‑1 (MQ=255)
                  cgcgTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGt                                                               <  1:1612760/54‑1 (MQ=255)
                  cgcgTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGt                                                               <  1:621037/54‑1 (MQ=255)
                  cgcgTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGt                                                               <  1:689066/54‑1 (MQ=255)
                    cgtcgtGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGt                                                               <  1:2902999/52‑1 (MQ=255)
                        gtGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCg                                               >  1:2114857/1‑64 (MQ=255)
                        gtGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGa                                          >  1:1844255/1‑69 (MQ=255)
                        gtGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGa                                          >  1:3087728/1‑69 (MQ=255)
                        gtGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGa                                          >  1:2925442/1‑69 (MQ=255)
                        gtGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGa                                          >  1:734573/1‑69 (MQ=255)
                           aTAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCaa                                                                  <  1:3122804/42‑1 (MQ=255)
                                gcgCTTCGTCCGAGACATCGATATCGCCCACTTTCAGCGTTCGTGCCAGGCCACCGGCAGAACCgg                                     >  1:1260654/3‑66 (MQ=255)
                                gcgCTTCGTCCGAGACAACGATATCGCCCACTTTCAGCGTTGGTGCCAGGCCACCGGCAGAACCGgtg                                   >  1:321194/3‑68 (MQ=255)
                                gcgCTTCGTCCGAGACAACGATATCGCCCACTTTCAGCGTTGGTGCCAGGCCACCGGCAGAACCGgtg                                   >  1:2507543/3‑68 (MQ=255)
                                gcgCTTCGTCCGAGACAACGATATCGCCCACTTTCAGCGTTGGTGCCAGGCCACCGGCAGAACCGgt                                    >  1:1753819/3‑67 (MQ=255)
                                gcgCTTCGTCCGAGACAACGATATCGCCCACTTTCAGCGTTGGTGCCAGGCCACCGGCAGAACCGgt                                    >  1:3001061/3‑67 (MQ=255)
                                     tCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTtaat                              <  1:1583557/68‑1 (MQ=255)
                                     tCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTtaat                              <  1:317702/68‑1 (MQ=255)
                                     tCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTtaat                              <  1:1181939/68‑1 (MQ=255)
                                            agaCAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATcaca                       <  1:744838/68‑1 (MQ=255)
                                            agaCAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATcaca                       <  1:2147291/68‑1 (MQ=255)
                                            agaCAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATcaca                       <  1:2680132/68‑1 (MQ=255)
                                            agaCAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATcaca                       <  1:2628875/68‑1 (MQ=255)
                                                            ctcTTTCAGCGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGtt       <  1:322275/66‑1 (MQ=255)
                                                             aCTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGTCTTGCAGTGTTcc     <  1:3316314/69‑1 (MQ=255)
                                                             aCTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGTCTTGCAGTGTTcc     <  1:3312242/69‑1 (MQ=255)
                                                                ttCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAg           >  1:2342714/1‑60 (MQ=255)
                                                                ttCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAAc  >  1:2994537/1‑69 (MQ=255)
                                                                ttCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAAc  >  1:697848/1‑69 (MQ=255)
                                                                ttCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAAc  >  1:2282311/1‑69 (MQ=255)
                                                                ttCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAAc  >  1:1903546/1‑69 (MQ=255)
                                                                ttCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAAc  >  1:1081938/1‑69 (MQ=255)
                                                                ttCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAAc  >  1:1967525/1‑69 (MQ=255)
                                                                ttCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAAc  >  1:3273717/1‑69 (MQ=255)
                                                                ttCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAAc  >  1:2594424/1‑69 (MQ=255)
                                                                ttCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAAc  >  1:456551/1‑69 (MQ=255)
                                                                ttCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAAc  >  1:3290046/1‑69 (MQ=255)
                                                                ttCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAAc  >  1:2911661/1‑69 (MQ=255)
                                                                ttCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAAc  >  1:1412672/1‑69 (MQ=255)
                                                                ttCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAAc  >  1:2558182/1‑69 (MQ=255)
                                                                ttCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAAc  >  1:2833261/1‑69 (MQ=255)
                                                                ttCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAAc  >  1:2573891/1‑69 (MQ=255)
                                                                ttCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAAAAATCACATCTGGCTTGCAGTGTTCCAAc  >  1:3070270/1‑69 (MQ=255)
                                                                 tCAGCGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGtt       <  1:597712/63‑1 (MQ=255)
                                                                 tCAGCGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGtt       <  1:808327/63‑1 (MQ=255)
                                                                 tCAGCGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGtt       <  1:350210/63‑1 (MQ=255)
                                                                 tCAGCGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGtt       <  1:3208997/63‑1 (MQ=255)
                                                                 tCAGCGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGtt       <  1:3198449/63‑1 (MQ=255)
                                                                 tCAGCGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGtt       <  1:2912571/63‑1 (MQ=255)
                                                                 tCAGCGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGtt       <  1:225440/63‑1 (MQ=255)
                                                                 tCAGCGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGtt       <  1:2240952/63‑1 (MQ=255)
                                                                 tCAGCGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGtt       <  1:2239927/63‑1 (MQ=255)
                                                                 tCAGCGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGtt       <  1:2085322/63‑1 (MQ=255)
                                                                 tCAGCGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGtt       <  1:2026004/63‑1 (MQ=255)
                                                                 tCAGCGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGtt       <  1:2023584/63‑1 (MQ=255)
                                                                 tCAGCGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGtt       <  1:145933/63‑1 (MQ=255)
                                                                 tCAGCGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGtt       <  1:1332970/63‑1 (MQ=255)
                                                                 tCAGCGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGtt       <  1:1321194/63‑1 (MQ=255)
                                                                 tCAGCGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGtt       <  1:129603/63‑1 (MQ=255)
                                                                 tCAGCGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGtt       <  1:1103183/63‑1 (MQ=255)
                                                                 tCAGCGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGtt       <  1:1057948/63‑1 (MQ=255)
                                                                 tCAGCGTTGGTGCCAGGCCACCGGCAGAACCGGGGTTAATAATCACATCTGGCTTGCAGTGtt       <  1:3075452/63‑1 (MQ=255)
                                                                    |                                                                
ACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAAC  >  minE/157661‑157793

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: