Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 870,646 G→A 61.7% V215V (GTG→GTA topA → DNA topoisomerase I, omega subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE870,6460GA61.7% 34.2 / 47.5 73V215V (GTG→GTAtopADNA topoisomerase I, omega subunit
Reads supporting (aligned to +/- strand):  ref base G (13/15);  new base A (35/10);  total (48/25)
Fisher's exact test for biased strand distribution p-value = 1.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

GCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGC  >  minE/870599‑870699
                                               |                                                     
gCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTc                                  <  1:1384528/69‑1 (MQ=255)
gCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTc                                  <  1:279077/69‑1 (MQ=255)
gCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTc                                  <  1:1720075/69‑1 (MQ=255)
gCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTc                                  <  1:1644660/69‑1 (MQ=255)
gCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTc                                  <  1:1946407/69‑1 (MQ=255)
gCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTc                                  <  1:1473652/69‑1 (MQ=255)
gCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTc                                  <  1:1383579/69‑1 (MQ=255)
gCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTc                                  <  1:3362775/69‑1 (MQ=255)
gCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTc                                  <  1:818749/69‑1 (MQ=255)
gCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTc                                  <  1:2248125/69‑1 (MQ=255)
gCGGTCCTCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTc                                  >  1:2973917/1‑69 (MQ=255)
gCGGTCCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTc                                  >  1:2061901/1‑69 (MQ=255)
gCGGTCCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTc                                  >  1:2721921/1‑69 (MQ=255)
gCGGTCCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTc                                  >  1:2141379/1‑69 (MQ=255)
gCGGTCCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTc                                  >  1:2900636/1‑69 (MQ=255)
gCGGTCCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTc                                  >  1:2746982/1‑69 (MQ=255)
gCGGTCCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTc                                  >  1:2013208/1‑69 (MQ=255)
gCGGTCCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTc                                  >  1:2212614/1‑69 (MQ=255)
gCGGTCCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTc                                  >  1:1941321/1‑69 (MQ=255)
gCGGTCCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTc                                  >  1:1928502/1‑69 (MQ=255)
gCGGTCCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTc                                  >  1:2439265/1‑69 (MQ=255)
gCGGTCCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTc                                  >  1:2273898/1‑69 (MQ=255)
gCGGTCCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTc                                  >  1:1674564/1‑69 (MQ=255)
gCGGTCCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTc                                  >  1:293947/1‑69 (MQ=255)
gCGGTCCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTc                                  >  1:3217464/1‑69 (MQ=255)
gCGGTCCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTc                                  >  1:419613/1‑69 (MQ=255)
gCGGTCCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTc                                  >  1:810345/1‑69 (MQ=255)
gCGGTCCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTc                                  >  1:1065947/1‑69 (MQ=255)
gCGGTCCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTc                                  >  1:868485/1‑69 (MQ=255)
 cGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTc                                  <  1:1082910/68‑1 (MQ=255)
     tCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTg                                         >  1:73851/1‑57 (MQ=255)
     tCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGa                                >  1:3178806/1‑66 (MQ=255)
     tCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGAt                               >  1:2768253/1‑67 (MQ=255)
     tCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGAt                               >  1:1022194/1‑67 (MQ=255)
     tCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGAt                               >  1:2943577/1‑67 (MQ=255)
     tCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGAt                               >  1:822859/1‑67 (MQ=255)
     tCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGAt                               >  1:257216/1‑67 (MQ=255)
     tCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGAt                               >  1:2282286/1‑67 (MQ=255)
     tCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGAt                               >  1:1531057/1‑67 (MQ=255)
     tCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGAt                               >  1:1015844/1‑67 (MQ=255)
     tCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGAt                               >  1:1363776/1‑67 (MQ=255)
     tCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGAt                               >  1:1604688/1‑67 (MQ=255)
     tCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGAt                               >  1:1667/1‑67 (MQ=255)
      cGCCTTGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTACCGGAAGAGTTCTggg                                       <  1:1658122/58‑1 (MQ=255)
      cGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTACCGGAAGAGTTCTggg                                       <  1:16983/58‑1 (MQ=255)
      cGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTACCGGAAGAGTTCTggg                                       <  1:3200746/58‑1 (MQ=255)
      cGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTACCGGAAGAGTTCTggg                                       <  1:1313768/58‑1 (MQ=255)
      cGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTACCGGAAGAGTTCTggg                                       <  1:1689065/58‑1 (MQ=255)
      cGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTACCGGAAGAGTTCTggg                                       <  1:2052639/58‑1 (MQ=255)
      cGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTACCGGAAGAGTTCTggg                                       <  1:562392/58‑1 (MQ=255)
      cGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTACCGGAAGAGTTCTggg                                       <  1:571883/58‑1 (MQ=255)
      cGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTACCGGAAGAGTTCTggg                                       <  1:837272/58‑1 (MQ=255)
      cGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTACCGGAAGAGTTCTggg                                       <  1:2264359/58‑1 (MQ=255)
                          gcgtgaAATTAAAGCGTTCGTACCGGAAGAGTTCTg                                         >  1:1098617/1‑36 (MQ=255)
                          gcgtgaAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTg         >  1:2248828/1‑68 (MQ=255)
                          gcgtgaAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTg         >  1:530133/1‑68 (MQ=255)
                          gcgtgaAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTg         >  1:1218787/1‑68 (MQ=255)
                          gcgtgaAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTg         >  1:1191876/1‑68 (MQ=255)
                          gcgtgaAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTg         >  1:401490/1‑68 (MQ=255)
                          gcgtgaAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTg         >  1:3203432/1‑68 (MQ=255)
                          gcgtgaAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTg         >  1:1589317/1‑68 (MQ=255)
                          gcgtgaAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTg         >  1:2966005/1‑68 (MQ=255)
                          gcgtgaAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTg         >  1:1639667/1‑68 (MQ=255)
                          gcgtgaAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTg         >  1:2101693/1‑68 (MQ=255)
                          gcgtgaAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTg         >  1:2639905/1‑68 (MQ=255)
                          gcgtgaAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTg         >  1:2572265/1‑68 (MQ=255)
                          gcgtgaAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTg         >  1:223423/1‑68 (MQ=255)
                          gcgtgaAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTg         >  1:2544859/1‑68 (MQ=255)
                          gcgtgaAATTAAAGCGTTCGTACCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTg         >  1:2525558/1‑68 (MQ=255)
                                 aTTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGc  <  1:3119992/68‑1 (MQ=255)
                                 aTTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGc  <  1:590260/68‑1 (MQ=255)
                                 aTTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGc  <  1:2601163/68‑1 (MQ=255)
                                 aTTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGc  <  1:1103407/68‑1 (MQ=255)
                                               |                                                     
GCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGC  >  minE/870599‑870699

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: