Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,188,658 G→A 76.0% A67A (GCC→GCT nanA ← N‑acetylneuraminate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,188,6580GA76.0% 106.2 / 29.2 87A67A (GCC→GCTnanAN‑acetylneuraminate lyase
Reads supporting (aligned to +/- strand):  ref base G (5/16);  new base A (16/50);  total (21/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GGCTTTCGGCGGTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAA  >  minE/2188594‑2188704
                                                                |                                              
ggCTTTCGGCGGTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGcg                                             >  1:3240806/1‑68 (MQ=255)
       ggcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGAtt                                       <  1:2958118/67‑1 (MQ=255)
       ggcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGAtt                                       <  1:1858030/67‑1 (MQ=255)
       ggcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGAtt                                       <  1:1874101/67‑1 (MQ=255)
       ggcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGAtt                                       <  1:3325512/67‑1 (MQ=255)
       ggcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGAtt                                       <  1:3164245/67‑1 (MQ=255)
       ggcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGAtt                                       <  1:3160121/67‑1 (MQ=255)
       ggcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGAtt                                       <  1:3128660/67‑1 (MQ=255)
       ggcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGAtt                                       <  1:2191818/67‑1 (MQ=255)
       ggcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGAtt                                       <  1:2201400/67‑1 (MQ=255)
       ggcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGAtt                                       <  1:2028501/67‑1 (MQ=255)
       ggcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGAtt                                       <  1:2938942/67‑1 (MQ=255)
       ggcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGAtt                                       <  1:2362181/67‑1 (MQ=255)
       ggcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGAtt                                       <  1:2420755/67‑1 (MQ=255)
       ggcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGAtt                                       <  1:2422627/67‑1 (MQ=255)
       ggcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGAtt                                       <  1:284619/67‑1 (MQ=255)
       ggcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGAtt                                       <  1:2474555/67‑1 (MQ=255)
       ggcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGAtt                                       <  1:2479301/67‑1 (MQ=255)
       ggcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGAtt                                       <  1:2572171/67‑1 (MQ=255)
       ggcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGAtt                                       <  1:1164859/67‑1 (MQ=255)
       ggcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGAtt                                       <  1:1186278/67‑1 (MQ=255)
       ggcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGAtt                                       <  1:893464/67‑1 (MQ=255)
       ggcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGAtt                                       <  1:841111/67‑1 (MQ=255)
       ggcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGAtt                                       <  1:1561711/67‑1 (MQ=255)
       ggcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGAtt                                       <  1:1566407/67‑1 (MQ=255)
       ggcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGAtt                                       <  1:1659998/67‑1 (MQ=255)
       ggcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGAtt                                       <  1:1768034/67‑1 (MQ=255)
       ggcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGAtt                                       <  1:1707048/67‑1 (MQ=255)
        gcggTGCTGACGCAACCGATGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTc                                     <  1:2978013/68‑1 (MQ=255)
        gcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTc                                     <  1:332560/68‑1 (MQ=255)
        gcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTc                                     <  1:1642384/68‑1 (MQ=255)
        gcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTc                                     <  1:2827792/68‑1 (MQ=255)
        gcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTc                                     <  1:348936/68‑1 (MQ=255)
        gcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTc                                     <  1:1281838/68‑1 (MQ=255)
        gcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTc                                     <  1:677120/68‑1 (MQ=255)
        gcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTc                                     <  1:2411595/68‑1 (MQ=255)
        gcggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTc                                     <  1:2474245/68‑1 (MQ=255)
         cggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTcc                                    <  1:1508592/68‑1 (MQ=255)
         cggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTc                                     <  1:3309827/67‑1 (MQ=255)
         cggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTc                                     <  1:781889/67‑1 (MQ=255)
         cggTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGAtt                                       <  1:1795718/65‑1 (MQ=255)
         cggTGCTGACCCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTc                                     <  1:3373227/67‑1 (MQ=255)
              cTGACGCAACCGACGTGGTCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACTAtt                                       <  1:783615/60‑1 (MQ=255)
                gACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTCCAGTACCTg                            >  1:2009726/1‑69 (MQ=255)
                  cGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTc                                     <  1:2717127/58‑1 (MQ=255)
                   gCAACCGCCGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTTCCTGtt                          >  1:2510672/1‑68 (MQ=255)
                   gCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTCCAGTACCTGtt                          <  1:93718/68‑1 (MQ=255)
                   gCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTCCAGTACCTGtt                          <  1:2557303/68‑1 (MQ=255)
                   gCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGTTc                         >  1:2679150/1‑69 (MQ=255)
                   gCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGTTc                         >  1:1378441/1‑69 (MQ=255)
                   gCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGTTc                         >  1:823378/1‑69 (MQ=255)
                   gCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGTTc                         >  1:526339/1‑69 (MQ=255)
                   gCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGTTc                         >  1:3396485/1‑69 (MQ=255)
                   gCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGTTc                         >  1:3327143/1‑69 (MQ=255)
                   gCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGTTc                         >  1:1863601/1‑69 (MQ=255)
                   gCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGTTc                         >  1:3313399/1‑69 (MQ=255)
                   gCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGTTc                         >  1:2169171/1‑69 (MQ=255)
                   gCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGTTc                         >  1:2623357/1‑69 (MQ=255)
                   gCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGTTc                         >  1:2991773/1‑69 (MQ=255)
                   gCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGTTc                         >  1:3061237/1‑69 (MQ=255)
                   gCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGTTc                         >  1:2282637/1‑69 (MQ=255)
                   gCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGAGTTCCAGTACCTGTTc                         >  1:1817876/1‑69 (MQ=255)
                   gCAACCGACGTGGGCGATGAGTATAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGTTc                         >  1:245053/1‑69 (MQ=255)
                     aaCCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGAtt                                       <  1:476505/53‑1 (MQ=255)
                        cGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGAtt                                       <  1:1798140/50‑1 (MQ=255)
                               ggCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAg             >  1:2140372/1‑69 (MQ=255)
                               ggCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAg             >  1:2165622/1‑69 (MQ=255)
                                  gATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGtt                          <  1:217779/53‑1 (MQ=255)
                                  gATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGtt                          <  1:3302097/53‑1 (MQ=255)
                                  gATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGtt                          <  1:2699268/53‑1 (MQ=255)
                                  gATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGtt                          <  1:885662/53‑1 (MQ=255)
                                  gATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGtt                          <  1:1411548/53‑1 (MQ=255)
                                  gATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGtt                          <  1:1006278/53‑1 (MQ=255)
                                  gATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGtt                          <  1:1623597/53‑1 (MQ=255)
                                  gATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGtt                          <  1:1643006/53‑1 (MQ=255)
                                  gATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGtt                          <  1:703397/53‑1 (MQ=255)
                                  gATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGtt                          <  1:563517/53‑1 (MQ=255)
                                  gATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGtt                          <  1:358011/53‑1 (MQ=255)
                                  gATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGtt                          <  1:2285036/53‑1 (MQ=255)
                                  gATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGtt                          <  1:1843936/53‑1 (MQ=255)
                                  gATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGtt                          <  1:3361673/53‑1 (MQ=255)
                                  gATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGtt                          <  1:2151340/53‑1 (MQ=255)
                                  gATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGtt                          <  1:2959642/53‑1 (MQ=255)
                                  gATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGtt                          <  1:3325390/53‑1 (MQ=255)
                                  gATGAGTTTAATCTTACCTTTCGCCTCTTCAGCGACCATTTCCAGTACCTGtt                          <  1:64295/53‑1 (MQ=39)
                                  gATGACTTTAATCTTACCTTTCGCCTCTTCAGCGACGATTTCCAGTACCTGtt                          <  1:2140391/53‑1 (MQ=39)
                                          tAATCTTACCTTTCGCCTCTTCGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACaa  >  1:2258280/1‑69 (MQ=255)
                                                                |                                              
GGCTTTCGGCGGTGCTGACGCAACCGACGTGGGCGATGAGTTTAATCTTACCTTTCGCCTCTTCGGCGACGATTTCCAGTACCTGTTCACGCTCGGAAAGGCTTTGTACAA  >  minE/2188594‑2188704

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: