Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,622,314	G→A	R1041R (AGG→AGA)	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,314	0	G→A	100.0%	74.9	39	R1041R (AGG→AGA)	rhsB	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  new base (23/16):  ref base (0/0):  total (23/16)

CTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATC  >  NC_000913/3622225‑3622509                                                                                          |                                                                                                                                                                                                    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:603987/150‑1 (MQ=9)                                                        aCGGCGACCGGCTGACCACAATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTcagc                                                                                  <  2:523401/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)                                                                          caaTACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGATGGCGAAAACGCAGCGCCGCGGCCTGGCGGATGCCCTTCAGCAGTCcggcg                                                                         >  4:236851/3‑141 (MQ=255)                                                                          caaTACAGAACGACAGAACCCGCAGCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTcagc                                                                                  <  3:562947/130‑1 (MQ=12)                                                                            aTACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCTCGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCggcggc                                                                       >  4:503410/1‑141 (MQ=9)                                                                            aTACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGt                                                              >  2:113048/1‑150 (MQ=14)                                                                            aTACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTggc            >  1:730051/1‑199 (MQ=14)                                                                            aTACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCACCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGa            >  1:500599/1‑200 (MQ=14)                                                                               cAGAACGACAGAACCCGCATCCAGACGATTTATGAGCCGGGTAGCTTCAGTCCCCTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCATCCTGGCGGATGCCCTTCAGCCGTCCGGCGGCGAAGACGGTGGc                                                           >  2:706088/1‑150 (MQ=11)                                                                               cAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGc                                                           <  2:317873/150‑1 (MQ=14)                                                                               cAGAACGACAGAACCCGCATACAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGc                                                           <  2:513075/150‑1 (MQ=9)                                                                                 gAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCCCCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAg                                                         >  2:380496/1‑150 (MQ=9)                                                                                 gAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAg                                                         <  2:79589/150‑1 (MQ=14)                                                                                 gAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAg                                                         >  2:800384/1‑150 (MQ=14)                                                                                    cGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCcggcgg                                                                        >  3:34899/1‑132 (MQ=17)                                                                                    cGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACgg                                                               >  4:80730/1‑141 (MQ=14)                                                                                    cGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAgtgt                                                      >  2:572975/1‑150 (MQ=14)                                                                                      aCAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGgcggcg                                                                      >  3:521907/1‑132 (MQ=17)                                                                                      aCAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATc  >  1:483667/1‑200 (MQ=14)                                                                                      aCAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATc  >  1:248671/1‑200 (MQ=14)                                                                                          |                                                                                                                                                                                                    CTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATC  >  NC_000913/3622225‑3622509

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.