Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,622,299	G→A	T1036T (ACG→ACA)	rhsB →	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	3,622,299	0	G→A	100.0%	35.0	22	T1036T (ACG→ACA)	rhsB	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  new base (8/12):  ref base (0/1):  total (9/13)

CCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGG  >  NC_000913/3622176‑3622494                                                                                                                            |                                                                                                                                                                                                    cccGGTGGGCCGCAGGGTGGC‑AAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGagc                                                                                                                                                                                            <  3:466284/132‑3 (MQ=12)                                                  cTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACAATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCt                                                                        >  1:407718/1‑200 (MQ=14)                                                       gggCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACAATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTcagc                                                                   <  1:347148/200‑1 (MQ=14)                                                        ggCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCTACCGGCTGACCACCATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTcagca                                                                  >  1:101193/1‑200 (MQ=12)                                                            ggATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACAATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTcc                                                              <  1:682810/200‑1 (MQ=14)                                                              aTGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACAATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCcgg                                                            >  1:598467/1‑200 (MQ=14)                                                                        cACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACAATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGAc                                                  >  1:258084/1‑200 (MQ=14)                                                                                  gCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACAATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAgtg                                        <  1:191282/200‑1 (MQ=14)                                                                                           caggTACGGCTGGGAGCGCGAACGGCTGACCACAATGCAGAACGACAGAAACCGCATCCAGACGATTTAGCAGCCGGGGAGCTTCACGCCACTCAAAAGGGTCGAAGCTGCCACCGGTGAGCTGGCGAAAAc                                                                                                   <  3:542467/130‑1 (MQ=11)                                                                                             ggTATGGCTGGGACGGCGACCGGCTGACCACAATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCccg                             >  1:575035/1‑200 (MQ=9)                                                                                                  ggCTGGGACGGCGACCGGCTGACCACAATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGt                        >  1:83399/1‑200 (MQ=14)                                                                                                      gggACGGCGACCGGCTGACCACAATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTc                                                                      <  2:19703/150‑1 (MQ=9)                                                                                                         actgCGACCGGCTGACCACAATACAGAACTACATAAACCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGTGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGAGCCGCAGCCTGGCGGATGCCCTTcagc                                                                   <  2:198977/147‑1 (MQ=11)                                                                                                         aCGGCGACCGGCTGCCCACAATACAGAACGAAAGAACCCGCATCCAGAAGATTTAACAGCCGGAGAGCTTCACGCCAATCATAAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGACTGGCGGATGCCCTTcagc                                                                   <  2:85416/150‑1 (MQ=11)                                                                                                         aCGGCGACCGGCTGACCACAATACAGATCGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTcagc                                                                   <  2:533846/150‑1 (MQ=9)                                                                                                         aCGGCGACCGGCTGACCACAATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTcagc                                                                   <  2:523401/150‑1 (MQ=9)                                                                                                         aCGGCGACCGGCTGACCACAATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTcagc                                                                   <  2:603987/150‑1 (MQ=9)                                                                                                              gACCGGCTGACCACAATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGACGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGAt            <  1:254617/200‑1 (MQ=9)                                                                                                              gACCGGCTGACCACAATACAGAACGACAGAACCCGCATCCAGACGAGTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTcc                                                              <  2:661827/150‑1 (MQ=9)                                                                                                                      gACCACAATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGAcc    >  1:617665/1‑200 (MQ=14)                                                                                                                        ccACAATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCcg                                                             >  4:255807/1‑141 (MQ=9)                                                                                                                        ccACAATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCgg  >  1:685007/1‑200 (MQ=14)                                                                                                                            |                                                                                                                                                                                                    CCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGG  >  NC_000913/3622176‑3622494

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