Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,125,105 A→G 54.8% R177R (CGT→CGC pheT ← phenylalanine tRNA synthetase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,125,1050AG54.8% 17.4 / 66.8 84R177R (CGT→CGCpheTphenylalanine tRNA synthetase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (21/17);  new base G (28/18);  total (49/35)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

GGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATG  >  minE/1125042‑1125165
                                                               |                                                            
ggTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACgcgc                                                         >  1:2180302/1‑69 (MQ=255)
ggTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACgcgc                                                         >  1:2172309/1‑69 (MQ=255)
ggTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACgcgc                                                         >  1:799620/1‑69 (MQ=255)
ggTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACgcgc                                                         >  1:723087/1‑69 (MQ=255)
ggTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACgcgc                                                         >  1:2210105/1‑69 (MQ=255)
ggTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACgcgc                                                         >  1:2798365/1‑69 (MQ=255)
ggTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACgcgc                                                         >  1:3286779/1‑69 (MQ=255)
ggTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACgcgc                                                         >  1:2807485/1‑69 (MQ=255)
  tCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCaa                                                       >  1:2253249/1‑69 (MQ=255)
      aCCAACCGGAACGATTTCCGGTTCAACCAGCGGCAGCTGGTTGAACACGGCAACGTCGCGCGCAACAcc                                                   >  1:3241539/1‑69 (MQ=255)
      aCCAACCGGAACGATTTCCGGTTCAACCAGCGGCAGCTGGTTCAACACGGCAACGTCGCGCGCAacac                                                    >  1:2187984/1‑68 (MQ=255)
      aCCAACCGGAACGATTTCCGGTTCAACCAGCGGCAGCTGGTTCAACACGGCAACGTCGCGCGCAACAcc                                                   >  1:88472/1‑69 (MQ=255)
      aCCAACCGGAACGATTTCCGGTTCAACCAGCGGCAGCTGGTTCAACACGGCAACGTCGCGCGCAACAcc                                                   >  1:2167796/1‑69 (MQ=255)
      aCCAACCGGAACGATTTCCGGTTCAACCAGCGGCAGCTGGTTCAACACGGCAACGTCGCGCGCAACAcc                                                   >  1:2934512/1‑69 (MQ=255)
      aCCAACCGGAACGATTTCCGGTTCAACCAGCGGCAGCTGGTTCAACACGGCAACGTCGCGCGCAACAcc                                                   >  1:1812798/1‑69 (MQ=255)
      aCCAACCGGAACGATTTCCGGTTCAACCAGCGGCAGCTGGTTCAACACGGCAACGTCGCGCGCAACAcc                                                   >  1:3357169/1‑69 (MQ=255)
      aCCAACCGGAACGATTTCCGGTTCAACCAGCGGCAGCTGGTTCAACACGGCAACGTCGCGCGCAACAcc                                                   >  1:2798925/1‑69 (MQ=255)
      aCCAACCGGAACGATTTCCGGTTCAACCAGCGGCAGCTGGTTCAACACGGCAACGTCGCGCGCAACAcc                                                   >  1:1659409/1‑69 (MQ=255)
      aCCAACCGGAACGATTTCCGGTTCAACCAGCGGCAGCTGGTTCAACACGGCAACGTCGCGCGCAACAcc                                                   >  1:2327500/1‑69 (MQ=255)
      aCCAACCGGAACGATTTCCGGTTCAACCAGCGGCAGCTGGTTCAACACGGCAACGTCGCGCGCAACAcc                                                   >  1:1404104/1‑69 (MQ=255)
      aCCAACCGGAACGATTTCCGGTTCAACCAGCGGCAGCTGGTTCAACACGGCAACGTCGCGCGCAACAcc                                                   >  1:616876/1‑69 (MQ=255)
      aCCAACCGGAACGATTTCCGGTTCAACCAGCGGCAGCTGGTTCAACACGGCAACGTCGCGCGCAACAcc                                                   >  1:1197273/1‑69 (MQ=255)
         agccGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAat                                                <  1:1884758/67‑1 (MQ=255)
         aaCCGGAACGATTTCCTGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAat                                                <  1:2819441/69‑1 (MQ=255)
         aaCCGGAACGATTTCCGGTTTAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAat                                                <  1:2180433/69‑1 (MQ=255)
         aaCCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAat                                                <  1:2817411/69‑1 (MQ=255)
         aaCCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAat                                                <  1:2331151/69‑1 (MQ=255)
         aaCCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAat                                                <  1:2730926/69‑1 (MQ=255)
         aaCCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAat                                                <  1:300631/69‑1 (MQ=255)
         aaCCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAat                                                <  1:3016695/69‑1 (MQ=255)
         aaCCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAat                                                <  1:305769/69‑1 (MQ=255)
         aaCCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAat                                                <  1:196456/69‑1 (MQ=255)
         aaCCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAat                                                <  1:3298333/69‑1 (MQ=255)
         aaCCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAat                                                <  1:337432/69‑1 (MQ=255)
         aaCCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAat                                                <  1:393408/69‑1 (MQ=255)
         aaCCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAat                                                <  1:1483432/69‑1 (MQ=255)
         aaCCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAat                                                <  1:1156618/69‑1 (MQ=255)
         aaCCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGGCACGCGCAACACCAat                                                <  1:1127110/69‑1 (MQ=255)
                            tGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAat                                                <  1:1764894/50‑1 (MQ=255)
                            tGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCa                             >  1:22615/1‑69 (MQ=255)
                            tGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCa                             >  1:2019938/1‑69 (MQ=255)
                                     ggCAGCTGGTTCAACACGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACGGTTTGGc                    <  1:1110785/69‑1 (MQ=255)
                                     ggCAGCTGGTTCAACACGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACGGTTTGGc                    <  1:2530924/69‑1 (MQ=255)
                                     ggCAGCTGGTTCAACACGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACGGTTTGGc                    <  1:728402/69‑1 (MQ=255)
                                     ggCAGCTGGTTCAACACGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACGGTTTGGc                    <  1:1293628/69‑1 (MQ=255)
                                     ggCAGCTGGTTCAACACGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACGGTTTGGc                    <  1:3380884/69‑1 (MQ=255)
                                     ggCAGCTGGTTCAACACGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACGGTTTGGc                    <  1:3262722/69‑1 (MQ=255)
                                     ggCAGCTGGTTCAACACGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACGGTTTGGc                    <  1:2133777/69‑1 (MQ=255)
                                     ggCAGCTGGTTCAACACGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACGGTTTGGc                    <  1:3137979/69‑1 (MQ=255)
                                     ggCAGCTGGTTCAACACGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACGGTTTGGc                    <  1:3103916/69‑1 (MQ=255)
                                     ggCAGCTGGTTCAACACGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACGGTTTGGc                    <  1:2115104/69‑1 (MQ=255)
                                      gCAGCTGGTTCAACACGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACGGTTTGGc                    <  1:2575739/68‑1 (MQ=255)
                                                cAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGa          >  1:509062/1‑68 (MQ=255)
                                                cAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGAt         >  1:605953/1‑69 (MQ=255)
                                                cAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGAt         >  1:801611/1‑69 (MQ=255)
                                                cAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGAt         >  1:2827767/1‑69 (MQ=255)
                                                cAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGAt         >  1:475504/1‑69 (MQ=255)
                                                cAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGAt         >  1:1347890/1‑69 (MQ=255)
                                                cAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGAt         >  1:3354816/1‑69 (MQ=255)
                                                cAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGAt         >  1:2302479/1‑69 (MQ=255)
                                                cAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGAt         >  1:1758570/1‑69 (MQ=255)
                                                 aacACGGCAACGTCGCGCGCAACACCAATGATACCTGAGCAGTCGGCACTGTTTGGCGTCACGCTAAtt        >  1:3366780/3‑69 (MQ=255)
                                                 aacACGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACTGTTTGGCGTCACGCTAAtt        >  1:226630/3‑69 (MQ=255)
                                                 aacACGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACTGTTTGGCGTCACGCTAAtt        >  1:421248/3‑69 (MQ=255)
                                                 aacACGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACTGTTTGGCGTCACGCTAAtt        >  1:1555673/3‑69 (MQ=255)
                                                 aacACGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACTGTTTGGCGTCACGCTAAtt        >  1:613496/3‑69 (MQ=255)
                                                 aacACGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACTGTTTGGCGTCACGCTAAtt        >  1:3367547/3‑69 (MQ=255)
                                                 aacACGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACTGTTTGGCGTCACGCTAAtt        >  1:330599/3‑69 (MQ=255)
                                                 aacACGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACTGTTTGGCGTCACGCTAAtt        >  1:1846718/3‑69 (MQ=255)
                                                 aacACGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACTGTTTGGCGTCACGCTAAtt        >  1:1891516/3‑69 (MQ=255)
                                                 aacACGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACTGTTTGGCGTCACGCTAAtt        >  1:697509/3‑69 (MQ=255)
                                                 aacACGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACTGTTTGGCGTCACGCTAAtt        >  1:118739/3‑69 (MQ=255)
                                                 aacACGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACTGTTTGGCGTCACGCTAAtt        >  1:2909791/3‑69 (MQ=255)
                                                 aacACGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACTGTTTGGCGTCACGCTAAtt        >  1:19766/3‑69 (MQ=255)
                                                 aacACGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACTGTTTGGCGTCACGCTAAt         >  1:2554910/3‑68 (MQ=255)
                                                 aacACGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACTGTTTGGCGTCACGCTAAt         >  1:3420264/3‑68 (MQ=255)
                                                  acACGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACGGTTTGGCGTCACGCTAAtt        <  1:356624/67‑1 (MQ=255)
                                                  acACGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACGGTTTGGCGTCACGCTAAtt        <  1:2760630/67‑1 (MQ=255)
                                                  acACGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACGGTTTGGCGTCACGCTAAtt        <  1:1937142/67‑1 (MQ=255)
                                                  acACGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACGGTTTGGCGTCACGCTAAtt        <  1:3067396/67‑1 (MQ=255)
                                                  acACGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACGGTTTGGCGTCACGCTAAtt        <  1:806047/67‑1 (MQ=255)
                                                  acACGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACGGTTTGGCGTCACGCTAAtt        <  1:2073350/67‑1 (MQ=255)
                                                    aCGGCAACGTCGCGCGCAACACCAATGATACCTAAGCAGTCGGCACGGTTTGGc                    <  1:2200353/54‑1 (MQ=255)
                                                       gCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATg  >  1:2777923/1‑69 (MQ=255)
                                                               |                                                            
GGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATG  >  minE/1125042‑1125165

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