Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,654,281 A→G 45.5% S133S (TCA→TCG purM → phosphoribosylaminoimidazole synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,654,2810AG45.5% 4.5 / 43.9 55S133S (TCA→TCGpurMphosphoribosylaminoimidazole synthetase
Reads supporting (aligned to +/- strand):  ref base A (11/19);  new base G (13/12);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGATTAC  >  minE/1654213‑1654344
                                                                    |                                                               
aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCa                                                                 <  1:963121/69‑1 (MQ=255)
aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCa                                                                 <  1:956039/69‑1 (MQ=255)
aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCa                                                                 <  1:758955/69‑1 (MQ=255)
aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCa                                                                 <  1:533837/69‑1 (MQ=255)
aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCa                                                                 <  1:3065052/69‑1 (MQ=255)
aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCa                                                                 <  1:2720461/69‑1 (MQ=255)
aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCa                                                                 <  1:1009093/69‑1 (MQ=255)
aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCa                                                                 <  1:1862861/69‑1 (MQ=255)
 ccGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCa                                                                 <  1:2551581/68‑1 (MQ=255)
       aaCTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGtt                                                          <  1:2995661/69‑1 (MQ=255)
       aaCTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGtt                                                          <  1:120791/69‑1 (MQ=255)
       aaCTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGtt                                                          <  1:2719768/69‑1 (MQ=255)
                     aCCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGgtgggtgga                                            <  1:1254876/69‑2 (MQ=255)
                        gCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGa                                          <  1:1448157/68‑1 (MQ=255)
                        gCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGa                                          <  1:116505/68‑1 (MQ=255)
                        gCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGa                                          <  1:376433/68‑1 (MQ=255)
                        gCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGa                                          <  1:475694/68‑1 (MQ=255)
                            cAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGGAATCGGGCTGTTCACTGGTGGGTGGCGAAAc                                       >  1:2927943/1‑67 (MQ=255)
                            cAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAAc                                       >  1:409887/1‑67 (MQ=255)
                            cAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAAc                                       >  1:2632510/1‑67 (MQ=255)
                            cAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAAc                                       >  1:1246116/1‑67 (MQ=255)
                            cAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAAc                                       >  1:3324030/1‑67 (MQ=255)
                            cAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAAc                                       >  1:3392501/1‑67 (MQ=255)
                            cAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAAc                                       >  1:2332358/1‑67 (MQ=255)
                            cAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAAc                                       >  1:1301203/1‑67 (MQ=255)
                            cAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAAc                                       >  1:2232599/1‑67 (MQ=255)
                            cAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAAc                                       >  1:1647953/1‑67 (MQ=255)
                            cAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAAc                                       >  1:782039/1‑67 (MQ=255)
                            cAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAAc                                       >  1:2291997/1‑67 (MQ=255)
                            cAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAAc                                       >  1:1824735/1‑67 (MQ=255)
                                ggTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGa                                 <  1:2420300/69‑1 (MQ=255)
                                ggTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGa                                 <  1:3013107/69‑1 (MQ=255)
                                      cAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGcc                           >  1:837375/1‑69 (MQ=255)
                                                gCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGa                                 <  1:828471/53‑1 (MQ=255)
                                                gCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGa                                 <  1:1708266/53‑1 (MQ=255)
                                                gCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGa                                 <  1:673840/53‑1 (MQ=255)
                                                gCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGa                                 <  1:2084155/53‑1 (MQ=255)
                                                gCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGcc                           >  1:393137/1‑59 (MQ=255)
                                                gCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTa                  >  1:3111596/1‑68 (MQ=255)
                                                gCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTAt                 >  1:3107818/1‑69 (MQ=255)
                                                gCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTAt                 >  1:107914/1‑69 (MQ=255)
                                                gCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTAt                 >  1:121884/1‑69 (MQ=255)
                                                gCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTAt                 >  1:2294011/1‑69 (MQ=255)
                                                        ttGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACTGTGa         <  1:958428/69‑1 (MQ=255)
                                                        ttGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGa         <  1:2523487/69‑1 (MQ=255)
                                                        ttGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGa         <  1:2497895/69‑1 (MQ=255)
                                                        ttGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGa         <  1:3409329/69‑1 (MQ=255)
                                                        ttGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGa         <  1:2444346/69‑1 (MQ=255)
                                                        ttGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGa         <  1:233567/69‑1 (MQ=255)
                                                        ttGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGa         <  1:1118381/69‑1 (MQ=255)
                                                        ttGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGa         <  1:1751732/69‑1 (MQ=255)
                                                        ttGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGa         <  1:1885840/69‑1 (MQ=255)
                                                        ttGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGa         <  1:826946/69‑1 (MQ=255)
                                                        ttGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGa         <  1:1949255/69‑1 (MQ=255)
                                                        ttGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGa         <  1:1150875/69‑1 (MQ=255)
                                                          gTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAgg       >  1:283121/1‑69 (MQ=255)
                                                          gTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAgg       >  1:1266144/1‑69 (MQ=255)
                                                          gTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAgg       >  1:452220/1‑69 (MQ=255)
                                                          gTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAgg       >  1:2780492/1‑69 (MQ=255)
                                                                 aTCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGATTAc  <  1:1898464/67‑1 (MQ=255)
                                                                 aTCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGATTAc  <  1:61218/67‑1 (MQ=255)
                                                                 aTCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGATTAc  <  1:1960256/67‑1 (MQ=255)
                                                                 aTCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGATTAc  <  1:2204569/67‑1 (MQ=255)
                                                                 aTCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGATTAc  <  1:2883684/67‑1 (MQ=255)
                                                                    |                                                               
ACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGATTAC  >  minE/1654213‑1654344

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: