Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,775,8420T→A100.0% 0.1 3intergenic (‑12/+237)insJ‑3/cysHIS150 hypothetical protein/phosphoadenosine phosphosulfate reductase
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (1/0):  total (3/0)

GGGTATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGC  >  REL606/2775811‑2775871
                                |                             
gggTATTTTGGCTTTGACATATTAGGGGCTAATct                            >  1:19531613/1‑34 (MQ=0)
  gTATTTTGGCTTTGACATATTAGGGGCTAATTCAt                          >  1:8769005/1‑35 (MQ=1)
  gTATTTTGGCTTTGACATATTAGGGGCTAATTCAt                          >  1:3803619/1‑35 (MQ=1)
  gTATTTTGGCTTTGACATATTAGGGGCTAATTCAt                          >  1:3066492/1‑35 (MQ=1)
  gTATTTTGGCTTTGACATATTAGGGGCTAATTCAt                          >  1:20462309/1‑35 (MQ=1)
  gTATTTTGGCTTTGACATATTAGGGGCTAATTCAt                          >  1:19539415/1‑35 (MQ=1)
  gTATTTTGGCTTTGACATATTAGGGGCTAATTCAt                          >  1:10196218/1‑35 (MQ=1)
  gTATTTTGGCTTTGACATATTAGGGGCTAATTCAt                          >  1:11226991/1‑35 (MQ=1)
  gTATTTTGGCTTTGACATATTAGGGGCTAATTCAt                          >  1:12898767/1‑35 (MQ=1)
  gTATTTTGGCTTTGACATATTAGGGGCTAATTCAt                          >  1:12919419/1‑35 (MQ=1)
  gTATTTTGGCTTTGACATATTAGGGGCTAATACAt                          >  1:3610114/1‑35 (MQ=1)
  gTATTTTGGCTTTGACATATTAGGGGCTAAACCAt                          >  1:6634321/1‑35 (MQ=1)
     ttttGGCTTTGACATATTAGGGGCTAATTCCTTTc                       >  1:6658228/1‑35 (MQ=1)
     ttttGGCTTTGACATATTAGGGGCTAATTCCTTTc                       >  1:15329371/1‑35 (MQ=1)
     ttttGGCTTTGACATATTAGGGGCTAATTCCTTTc                       >  1:14319577/1‑35 (MQ=1)
      tttGGCTTTGACATATTAGGGGCTAATCCActtta                      >  1:19625523/1‑33 (MQ=1)
      tttGGCTTTGACATATTAGGGGCTAATACAATTct                      >  1:10990167/1‑34 (MQ=1)
       ttGGCTTTGACATATTAGGGGCTAAATCATTTctt                     >  1:13044533/1‑33 (MQ=1)
       ttGGCTTCGACATATTAGGGGCTAATCCACTTTAt                     >  1:20486401/1‑35 (MQ=0)
        tGGCTTTGACATATTAGGGGCTAATTCATTTTAtt                    >  1:4306105/1‑34 (MQ=1)
        tGGCTTTGACATATTAGGGGCTAATTCATTTTAtt                    >  1:8726934/1‑34 (MQ=1)
        tGGCTTTGACATATTAGGGGCTAATGCACTTTATc                    >  1:15182614/1‑35 (MQ=0)
         ggCTTTGACATATTAGGGGCTAATTCATTTTAttg                   >  1:15927886/1‑33 (MQ=1)
         ggCTTTGACATATTAGGGGCTAATCCCTTTTAttg                   >  1:5750152/1‑33 (MQ=1)
          gCTTTGACATATTAGGGGCTAATCCCATTTCTTcc                   >  1:12383780/1‑34 (MQ=0)
              tGACATATTAGGGGCTAATCCATTTCATCGCccca              >  1:11441091/1‑33 (MQ=1)
                          ggCTATTCCATTTCATCGTCCAACAAAATCGGGGc  >  1:13894126/1‑35 (MQ=1)
                                |                             
GGGTATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGC  >  REL606/2775811‑2775871

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ 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.