Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,049,641 G→A 18.5% S238S (AGC→AGT glnE ← fused deadenylyltransferase and adenylyltransferase for glutamine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,049,6410GA18.5% 90.2 / 11.2 59S238S (AGC→AGTglnEfused deadenylyltransferase and adenylyltransferase for glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base G (22/26);  new base A (4/7);  total (26/33)
Fisher's exact test for biased strand distribution p-value = 7.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATCGCGTAACGCTCCCAGTCGCGCCCCTGCTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACTTTCGCCAAACGGACGCAGCCGCATATCCACGCGATAGACGAAGCCATC  >  minE/2049580‑2049709
                                                             |                                                                    
ccATCGCGTAACGCTCCCAGTCGCGCCCCTGCTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGc                                                                <  1:2610900/68‑1 (MQ=255)
ccATCGCGTAACGCTCCCAGTCGCGCCCCTGCTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGc                                                                <  1:273689/68‑1 (MQ=255)
ccATCGCGTAACGCTCCCAGTCGCGCCCCTGCTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGc                                                                <  1:276642/68‑1 (MQ=255)
ccATCGCGTAACGCTCCCAGTCGCGCCCCTGCTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGc                                                                <  1:2817222/68‑1 (MQ=255)
  aTCGCGTAACGCTCCCAGTCGCGCCCCTGCTCCTGGTAATAATCTTCCAGCGCCGCAAAACTCAGCAcc                                                             >  1:124946/1‑69 (MQ=255)
  aTCGCGTAACGCTCCCAGTCGCGCCCCTGCTCCTGGTAATAATCTTCCAGCGCCGCAAAACTCAGCAcc                                                             >  1:2715330/1‑69 (MQ=255)
  aTCGCGTAACGCTCCCAGTCGCGCCCCTGCTCCTGGTAATAATCTTCCAGCGCCGCAAAACTCAGCAcc                                                             >  1:617744/1‑69 (MQ=255)
  aTCGCGTAACGCTCCCAGTCGCGCCCCTGCTCCTGGTAATAATCTTCCAGCGCCGCAAAACTCAGCAc                                                              >  1:1706349/1‑68 (MQ=255)
    cgcgTAACGCTCCCAGTCGCGCCCCTGCTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAg                                                           >  1:2022140/1‑69 (MQ=255)
    cgcgTAACGCTCCCAGTCGCGCCCCTGCTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAg                                                           >  1:3314475/1‑69 (MQ=255)
    cgcgTAACGCTCCCAGTCGCGCCCCTGCTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAg                                                           >  1:2968977/1‑69 (MQ=255)
        tAACGCTCCCAGTCGCGCCCCTGCTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAGCg                                                         >  1:2894397/1‑67 (MQ=255)
        tAACGCTCCCAGTCGCGACCCTGCTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAGCgg                                                        >  1:2208460/1‑68 (MQ=255)
             cTCCCAGTCGCGCCCCTGCTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACt                                                  <  1:3311155/69‑1 (MQ=255)
             cTCCCAGTCGCGCCCCTGCTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACt                                                  <  1:1808892/69‑1 (MQ=255)
                      gcgcCCCTGCTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACTTTCGCCaaa                                         >  1:1744495/1‑69 (MQ=255)
                      gcgcCCCTGCTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACTTTCGCCaaa                                         >  1:2946800/1‑69 (MQ=255)
                      gcgcCCCTGCTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACTTTCGCCaaa                                         >  1:117676/1‑69 (MQ=255)
                      gcgcCCCTGCTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACTTTCGCCaaa                                         >  1:1825954/1‑69 (MQ=255)
                      gcgcCCCTGCTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACTTTCGCCaaa                                         >  1:1838952/1‑69 (MQ=255)
                        gcCCCTGCTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACTTTCGCCAAACg                                       <  1:284085/69‑1 (MQ=255)
                        gcCCCTGCTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACTTTCGCCAAACg                                       <  1:2861112/69‑1 (MQ=255)
                        gcCCCTGCTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACTTTCGCCAAACg                                       <  1:2864829/69‑1 (MQ=255)
                        gcCCCTGCTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACTTTCGCCAAACg                                       <  1:3124067/69‑1 (MQ=255)
                        gcCCCTGCTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACTTTCGCCAAACg                                       <  1:2372972/69‑1 (MQ=255)
                        gcCCCTGCTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACTTTCGCCAAACg                                       <  1:3236446/69‑1 (MQ=255)
                        gcCCCTGCTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACTTTCGCCAAACg                                       <  1:1051644/69‑1 (MQ=255)
                        gcCCCTGCTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACTTTCGCCAAACg                                       <  1:1559897/69‑1 (MQ=255)
                        gcCCCTGCTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACTTTCGCCAAACg                                       <  1:683957/69‑1 (MQ=255)
                        gcCCCTGCTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACTTTCGCCAAACg                                       <  1:1388499/69‑1 (MQ=255)
                        gcCCCTGCTCCTGGTAATAATCTTACAACGCGGCAAAGCTCAGCACCAGCGGGCCACTTTCGCCAAACg                                       <  1:434920/69‑1 (MQ=255)
                               cTCCTGGTAATAATCTTCCAACGGGGCAAAGCTCAGCAcc                                                             >  1:1845994/1‑40 (MQ=255)
                               cTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCAcc                                                             >  1:3402165/1‑40 (MQ=255)
                               cTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCAcc                                                             >  1:753493/1‑40 (MQ=255)
                               cTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCAcc                                                             >  1:358342/1‑40 (MQ=255)
                               cTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCAcc                                                             >  1:2747700/1‑40 (MQ=255)
                               cTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCAcc                                                             >  1:3354211/1‑40 (MQ=255)
                               cTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCAcc                                                             >  1:1318150/1‑40 (MQ=255)
                                 ccTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACTTTCGCCAAACGGACGCAGc                               <  1:1530856/68‑1 (MQ=255)
                                 ccTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACTTTCGCCAAACGGACGCAGc                               <  1:1123328/68‑1 (MQ=255)
                                 ccTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACTTTCGCCAAACGGACGCAGc                               <  1:1317451/68‑1 (MQ=255)
                                 ccTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACTTTCGCCAAACGGACGCAGc                               <  1:1367662/68‑1 (MQ=255)
                                 ccTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACTTTCGCCAAACGGACGCAGc                               <  1:1441562/68‑1 (MQ=255)
                                 ccTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACTTTCGCCAAACGGACGCAGc                               <  1:2088288/68‑1 (MQ=255)
                                 ccTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACTTTCGCCAAACGGACGCAGc                               <  1:1582458/68‑1 (MQ=255)
                                 ccTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACTTTCGCCAAACGGACGCAGc                               <  1:1745430/68‑1 (MQ=255)
                                      taataaTCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACTTTCGCCAAACGGACGCAGCCGCata                         >  1:3366045/1‑69 (MQ=255)
                                      taataaTCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACTTTCGCCAAACGGACGCAGCCGCata                         >  1:3405475/1‑69 (MQ=255)
                                      taataaTCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACTTTCGCCAAACGGACGCAGCCGCata                         >  1:866715/1‑69 (MQ=255)
                                      taataaTCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACTTTCGCCAAACGGACGCAGCCGCata                         >  1:1096827/1‑69 (MQ=255)
                                      taataaTCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACTTTCGCCAAACGGACGCAGCCGCata                         >  1:1818196/1‑69 (MQ=255)
                                           atCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGcc                                                     <  1:1090748/36‑1 (MQ=255)
                                                       cgCAAAACTCAGCACCAGCGGGCCACTTTCACCAAACGGACGCAGCCGCATATCCACGCGATAGACGaa        <  1:3062011/68‑1 (MQ=255)
                                                       cgCAAAACTCAGCACCAGCGGGCCACTTTCACCAAACGGACGCAGCCGCATATCCACGCGATAGACGaa        <  1:1515502/68‑1 (MQ=255)
                                                       cgCAAAACTCAGCACCAGCGGGCCACTTTCACCAAACGGACGCAGCCGCATATCCACGCGATAGACGaa        <  1:508876/68‑1 (MQ=255)
                                                       cgCAAAACTCAGCACCAGCGGGCCACTTTCACCAAACGGACGCAGCCGCATATCCACGCGATAGACGaa        <  1:174997/68‑1 (MQ=255)
                                                       cgCAAAACTCAGCACCAGCGGGCCACTTTCACCAAACGGACGCAGCCGCATATCCACGCGATAGACGaa        <  1:391174/68‑1 (MQ=255)
                                                       cgCAAAACTCAGCACCAGCGGGCCACTTTCACCAAACGGACGCAGCCGCATATCCACGCGATAGACGaa        <  1:2785825/68‑1 (MQ=255)
                                                        gCAAAACTCAGCACCAGCGGGCCACTTTCACCAAACGGACGCAGCCGCATATCCACGCGATAGACGaa        <  1:567579/68‑1 (MQ=255)
                                                             gCTCAGCACCAGCGGGCCACTTTCGCCAAACGGACGCAGCCGCATATCCACGCGATAGACGAAGCCATc  <  1:451583/69‑1 (MQ=255)
                                                             gCTCAGCACCAGCGGGCCACTTTCGCCAAACGGACGCAGCCGCATATCCACGCGATAGACGAAGCCATc  <  1:1897664/69‑1 (MQ=255)
                                                             gCTCAGCACCAGCGGGCCACTTTCGCCAAACGGACGCAGCCGCATATCCACGCGATAGACGAAGCCATc  <  1:2110360/69‑1 (MQ=255)
                                                             gCTCAGCACCAGCGGGCCACTTTCGCCAAACGGACGCAGCCGCATATCCACGCGATAGACGAAGCCATc  <  1:1484962/69‑1 (MQ=255)
                                                             gCTCAGCACCAGCGGGCCACTTTCGCCAAACGGACGCAGCCGCATATCCACGCGATAGACGAAGCCATc  <  1:734801/69‑1 (MQ=255)
                                                             gCTCAGCACCAGCGGGCCACTTTCGCCAAACGGACGCAGCCGCATATCCACGCGATAGACGAAGCCATc  <  1:778602/69‑1 (MQ=255)
                                                             gCTCAGCACCAGCGGGCCACTTTCGCCAAACGGACGCAGCCGCATATCCACGCGATAGACGAAGCCATc  <  1:3083736/69‑1 (MQ=255)
                                                             |                                                                    
CCATCGCGTAACGCTCCCAGTCGCGCCCCTGCTCCTGGTAATAATCTTCCAACGCGGCAAAGCTCAGCACCAGCGGGCCACTTTCGCCAAACGGACGCAGCCGCATATCCACGCGATAGACGAAGCCATC  >  minE/2049580‑2049709

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: