Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,557,6290T→G100.0% 21.9 27T53P (ACC→CCC) yjiWhypothetical protein
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (17/0):  ref base (1/9):  total (18/9)
Fisher's exact test strand distribution p-value = 2.13e-06
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

ATCACTTTGACATCTACCGCCGTGCCGGTGGTAAAACCGGCGGCTTCCAGCCACT  >  REL606/4557598‑4557652
                               |                       
aTCACTTTGACATCTACCGCCGTGCCGGTGGGaaa                      >  1:18659954/1‑35 (MQ=255)
aTCACTTTGACATCTACCGCCGTGCCGGTGGGaaa                      >  1:6093991/1‑35 (MQ=255)
aTCACTTTGACATCTACCGCCGTGCCGGTGGGaaa                      >  1:19163292/1‑35 (MQ=255)
 tCACTTTGACATCTACCGCCGTGCCGGTGGGaaaa                     >  1:5136784/1‑35 (MQ=255)
 tCACTTTGACATCTACCGCCGTGCCGGTGGGaaaa                     >  1:12839008/1‑35 (MQ=255)
 tCACTTTGACATCTACCGCCGTGCCGGTGGGaaaa                     >  1:13827665/1‑35 (MQ=255)
  cACTTTGACATCTACCGCCGTGCCGGTGGTAAAAc                    >  1:12214194/1‑35 (MQ=255)
  cACTTTGACATCTACCGCCGTGCCGGTGGGaaaaa                    >  1:2015973/1‑34 (MQ=255)
  cACTTTGACATCTACCGCCGTGCCGGTGGGAAAAc                    >  1:11711235/1‑35 (MQ=255)
  cACTTTGACATCTACCGCCGTGCCGGTGGGAAAAc                    >  1:20237183/1‑35 (MQ=255)
   aCTTTGACATCTACCGCCGTGCCGGTGGGaaaaac                   >  1:17265330/1‑33 (MQ=37)
   aCTTTGACATCTACCGCCGTGCCGGTGGGaaaaac                   >  1:15264169/1‑33 (MQ=37)
   aCTTTGACATCTACCGCCGTGCCGGTGGGaaaaac                   >  1:19591200/1‑33 (MQ=37)
   aCTTTGACATCTACCGCCGTGCCGGTGGGAAAAcc                   >  1:13347954/1‑35 (MQ=255)
     tttGACATCTACCGCCGTGCCGGTGGGAAAACcgg                 >  1:16077790/1‑35 (MQ=255)
      ttGACATCTACCGCCGTGCCGGTGGGAAAAcccgg                >  1:10879095/1‑32 (MQ=25)
      ttGACATCTACCGCCGTGCCGGTGGGAAAACcggc                >  1:1402512/1‑35 (MQ=255)
      ttGACATCTACCGCCGTGCCGGTGGGAAAACcggc                >  1:18694735/1‑35 (MQ=255)
       tGACATCTACCGCCGTGCCGGTGGGAAAACcgggg               >  1:18670817/1‑33 (MQ=37)
        gACATCTACCGCCGTGCCGGTGGGGAAAAcggggg              >  1:4657888/1‑35 (MQ=16)
        gACATCTACCGCCGTGCCGGTGGGAAAACccgcgg              >  1:17393501/1‑35 (MQ=25)
          cATCTACCGCCGTGCCGGTGGGAAAACcggcgggg            >  1:8332214/1‑33 (MQ=37)
           aTCTACCGCCGTGCCGGTGGGGAAAAccgcggggt           >  1:19234587/1‑32 (MQ=12)
             cccccgccgTGCCGGTGGTAAAACCGGCGGCTTcc         <  1:7134782/32‑1 (MQ=255)
              ccccgccgTGCCGGTGGTAAAACCGGCGGCTTCCa        <  1:15484695/33‑1 (MQ=255)
              ccccgccgTGCCGGTGGTAAAACCGGCGGCTTCCa        <  1:15108456/33‑1 (MQ=255)
              ccccgccgTGCCGGTGGTAAAACCGGCGGCTTCCa        <  1:769560/33‑1 (MQ=255)
              ccccgccgTGCCGGTGATAAAACCGGCGGCTTCCa        <  1:14430849/33‑1 (MQ=37)
               gccgccgTGCCGGGGGTAAAACCGGCGGCTTCCAg       <  1:9889178/34‑1 (MQ=255)
               cccgccgTGCCGGTGGTAAAACCGGCGGCTTCCAg       <  1:5457324/34‑1 (MQ=255)
                ccgccgTGCCGGTGGTAAAACCGGCGGCTTCCAGc      <  1:11292922/35‑1 (MQ=255)
                 cgccgTGCCGGTGGTAAACCCGGCGGCTTCCAGcc     <  1:12170104/35‑1 (MQ=255)
                    cgTGCCGGTGTTAAAACCGGCGGGCTCCCGCCCCt  >  1:15309437/1‑35 (MQ=12)
                               |                       
ATCACTTTGACATCTACCGCCGTGCCGGTGGTAAAACCGGCGGCTTCCAGCCACT  >  REL606/4557598‑4557652

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.