Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0009131,396,3110A→C100.0% 1.5 4Q79P (CAG→CCG) insH1IS5 transposase and trans‑activator
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (0/1):  total (3/1)

GGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCCGCTGGAAACCATGCTACGCATTCACTGCATGCAGCATTGGTACAACCTGAGCGATGGCG  >  NC_000913/1396018‑1396337
                                                                                                                                                                                                                                                                                                     |                          
gggAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCCGCTGGAAACCCTGCTACGCCTTCCCTGCATGCCCCATTgg                     >  2:107116/1‑301 (MQ=255)
   aaggaagTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGATATCATGTTTTTCATCGGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACGTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCCGCTGGAAACCATGCTACGCATTCACTGCATGCAGCATTGGTa                   <  2:152033/300‑1 (MQ=25)
               aCATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCACGGCTGGTAATGGCCGGCGACCTTATCCCATGGAAACCATGCTACCCATTCACTGCAGCCGGCATTTGTCCACccgaagagag      >  1:230676/1‑293 (MQ=11)
                    aGCTAAGGAAGGTGCGAACAAGTCACTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCCACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGCAAAGAGAGTTTCTTGTCGCGCCTGGAGCGGATTCTTCCATGGAAAAACATGGTGGGCGTCATCGAGCCGTGTTACCCCAAGGCTGGTAATGGCGGCCGACCTTATCCGCTGGCACCCCTGTCACGCCTTCACTCCATGCCGCGTTGGTCAAGACAGAGCGATGccg  >  2:360760/1‑300 (MQ=11)
                                                                                                                                                                                                                                                                                                     |                          
GGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCCGCTGGAAACCATGCTACGCATTCACTGCATGCAGCATTGGTACAACCTGAGCGATGGCG  >  NC_000913/1396018‑1396337

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