Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	3,622,043	0	T→G	100.0%	1.0	7	L951R (CTC→CGC)	rhsB	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (0/4):  ref base (1/1):  total (1/6)

CGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGGACTCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATC  >  NC_000913/3622003‑3622336                                         |                                                                                                                                                                                                                                                                                                      cgggaggaccgtgcaggaagacTGAAAGAGAAAACCGGCCGGCTGCAGGAAGGGGGTACCGGCACGGATGAGGAGCGGACGCACAGGGACCGTTACGCCGGGCAGGACCGGTGGGTGCCCTACACGCGGACACAATAGGAAGAGGCGCTGGTGGGAAGTCGCTGTACTGACGACCCGCGGGGACGGAGGGTGGAAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATAc                                   <  2:237935/280‑1 (MQ=2) cggcaaaggtgtataaaagacaggACAGAGAAACCCGACCGCATCCCGCAAGGGGTTATCCGCACGGATGATGAGCGGACTCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATAc                                   <  2:207544/278‑1 (MQ=11)    aaTGACCGGCACGGCAGGCTGACAGAGAAAACCGACCGCATCCCGGAAGGGGGTATCCGCACGGATGATGAGCGGACTCACCGGTACCATTACGACAGGCAGCACCGGCTGGGGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGa                                <  1:495198/300‑1 (MQ=2)            ccACGGGAGGCTGACAGGGAAAACCGCCCCCCGCCCGGACGGGGTTAGCCGCACGGATGATGAGCGGCATCACCGGTACCACTACGACAGTCAGCAGCGGCGGGTGCACTACACGCGGACACAATATGAAGGGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAgg                        <  2:670697/300‑1 (MQ=11)                        gACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGGACTCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGCCCCGGATAcgtggctcttataccaatctcggg            >  1:207551/1‑278 (MQ=11)                           aaagaAACCCCACCGCAGCCCGGCAGGGGTTATCCCGAAGGCGGTCGGTCGGACTCAGTGGGACCATTACGACATTCATCCGCGGCGGGTGCACTACACGCGGACACAATAGGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCAGCAGGGTGGCAAAACGGGGGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACg         <  2:664344/299‑1 (MQ=11)                                  aCCGACCTCATCCCGGAAGGGGTGATCCGAACGGATGATGAGCGGACTCACCGGTACCATTAAGACAGTCAGCACCGGCGGGTGCACGACACGCGGACACAAAATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGAGGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATc  <  1:155570/301‑1 (MQ=255)                                         |                                                                                                                                                                                                                                                                                                      CGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGGACTCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATC  >  NC_000913/3622003‑3622336

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