Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 693,202 G→A 27.8% E50E (GAG→GAA hyaE → protein involved in processing of HyaA and HyaB proteins

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE693,2020GA27.8% 72.9 / 25.5 68E50E (GAG→GAAhyaEprotein involved in processing of HyaA and HyaB proteins
Reads supporting (aligned to +/- strand):  ref base G (38/11);  new base A (13/6);  total (51/17)
Fisher's exact test for biased strand distribution p-value = 5.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATA  >  minE/693134‑693257
                                                                    |                                                       
gacgacTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCagag                                                         >  1:1509976/1‑69 (MQ=255)
   gacTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGCGATCCGAAACGCACGCCAGAAGt                                                       <  1:912131/68‑1 (MQ=255)
   gacTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGCGATCCGAAACGCACGCCAGAAGt                                                       <  1:770815/68‑1 (MQ=255)
   gacTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGCGATCCGAAACGCACGCCAGAAGt                                                       <  1:1242810/68‑1 (MQ=255)
   gacTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGCGATCCGAAACGCACGCCAGAAGt                                                       <  1:2552563/68‑1 (MQ=255)
   gacTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGCGATCCGAAACGCACGCCAGAAGt                                                       <  1:2597378/68‑1 (MQ=255)
   gacGGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGCGATCCGAAACGCACGCCAGAAGt                                                       <  1:2155542/68‑1 (MQ=255)
      tGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGc                                                   >  1:745150/1‑69 (MQ=255)
      tGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGc                                                   >  1:1383839/1‑69 (MQ=255)
       ggCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAgg                                                   >  1:2139725/1‑67 (MQ=255)
       ggCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGc                                                   >  1:1872029/1‑68 (MQ=255)
       ggCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGc                                                   >  1:2836956/1‑68 (MQ=255)
          ttACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGAt                                                <  1:2939475/68‑1 (MQ=255)
          ttACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGAt                                                <  1:336062/68‑1 (MQ=255)
          ttACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGAt                                                <  1:354177/68‑1 (MQ=255)
          ttACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGAt                                                <  1:820389/68‑1 (MQ=255)
          ttACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGAt                                                <  1:1152632/68‑1 (MQ=255)
          ttACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGAt                                                <  1:1637519/68‑1 (MQ=255)
          ttACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGAt                                                <  1:1702856/68‑1 (MQ=255)
          ttACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGAt                                                <  1:1506198/68‑1 (MQ=255)
                 aGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGt                                        >  1:52776/1‑69 (MQ=255)
                 aGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGt                                        >  1:1162462/1‑69 (MQ=255)
                 aGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGt                                        >  1:327556/1‑69 (MQ=255)
                 aGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGt                                        >  1:1816993/1‑69 (MQ=255)
                 aGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGt                                        >  1:2118277/1‑69 (MQ=255)
                 aGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGt                                        >  1:2768517/1‑69 (MQ=255)
                                  tgtTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTAc                       >  1:2642634/1‑69 (MQ=255)
                                  tgtTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTAc                       >  1:274360/1‑69 (MQ=255)
                                  tgtTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTAc                       >  1:585327/1‑69 (MQ=255)
                                  tgtTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTAc                       >  1:645129/1‑69 (MQ=255)
                                  tgtTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTAc                       >  1:1832995/1‑69 (MQ=255)
                                      attaAGCAGTGACGCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTgc                    >  1:1193961/1‑68 (MQ=255)
                                      attaAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTgc                    >  1:378489/1‑68 (MQ=255)
                                      attaAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTgc                    >  1:1461796/1‑68 (MQ=255)
                                      attaAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTgc                    >  1:3191938/1‑68 (MQ=255)
                                      attaAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTgc                    >  1:1826243/1‑68 (MQ=255)
                                      attaAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTgc                    >  1:1994538/1‑68 (MQ=255)
                                      attaAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTg                     >  1:321258/1‑67 (MQ=255)
                                      attaAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTg                     >  1:3224281/1‑67 (MQ=255)
                                      attaAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCAGTAATGATTGGCGAATTACTgc                    >  1:2163175/1‑68 (MQ=255)
                                         aaGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGcgcg                 >  1:468371/1‑68 (MQ=255)
                                         aaGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGa                >  1:2549116/1‑69 (MQ=255)
                                         aaGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGa                >  1:364589/1‑69 (MQ=255)
                                              gTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTc           >  1:511857/1‑69 (MQ=255)
                                              gTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTc           >  1:597995/1‑69 (MQ=255)
                                                 aCCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTccc         >  1:694577/1‑68 (MQ=255)
                                                 aCCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCg        >  1:124616/1‑69 (MQ=255)
                                                 aCCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCg        >  1:810849/1‑69 (MQ=255)
                                                 aCCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCg        >  1:3113472/1‑69 (MQ=255)
                                                 aCCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCg        >  1:1888816/1‑69 (MQ=255)
                                                 aCCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCg        >  1:2073245/1‑69 (MQ=255)
                                                 aCCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCg        >  1:2795330/1‑69 (MQ=255)
                                                 aCCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCg        >  1:1530409/1‑69 (MQ=255)
                                                   ccGAAACGCACGCCAGAAGTCAGCGATAACCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCg        >  1:1117548/1‑67 (MQ=255)
                                                   ccGAAACGCACGCCAGAAGTCAGCGATAACCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGAc      >  1:2463730/1‑69 (MQ=255)
                                                   ccGAAACGCACGCCAGAAGTCAGCGATAACCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGAc      >  1:986961/1‑69 (MQ=255)
                                                   ccGAAACGCACGCCAGAAGTCAGCGATAACCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGAc      >  1:1680821/1‑69 (MQ=255)
                                                   ccGAAACGCACGCCAGAAGTCAGCGATAACCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGAc      >  1:2246998/1‑69 (MQ=255)
                                                   ccGAAACGCACGCCAGAAGTCAGCGATAACCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGAc      >  1:2318411/1‑69 (MQ=255)
                                                   ccGAAACGCACGCCAGAAGTCAGCGATAACCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGAc      >  1:3190323/1‑69 (MQ=255)
                                                   ccGAAACGCACGCCAGAAGTCAGCGATAACCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGAc      >  1:2690789/1‑69 (MQ=255)
                                                   ccGAAACGCACGCCAGAAGTCAGCGATAACCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGAc      >  1:1505381/1‑69 (MQ=255)
                                                   ccGAAACGCACGCCAGAAGTCAGCGATAACCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGAc      >  1:364041/1‑69 (MQ=255)
                                                   ccGAAACGCACGCCAGAAGTCAGCGATAACCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGAc      >  1:3148236/1‑69 (MQ=255)
                                                   ccGAAACGCACGCCAGAAGTCAGCGATAACCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGAc      >  1:1500049/1‑69 (MQ=255)
                                                   ccGAAACGCACGCCAGAAGTCAGCGATAACCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGAc      >  1:3348315/1‑69 (MQ=255)
                                                       aaCGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACtata  <  1:3068033/69‑1 (MQ=255)
                                                       aaCGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACtata  <  1:1927516/69‑1 (MQ=255)
                                                       aaCGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACtata  <  1:1788751/69‑1 (MQ=255)
                                                                    |                                                       
GACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATA  >  minE/693134‑693257

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: