Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,766,316	G→A	*1378* (TAG→TAA)	rhsA →	Rhs protein with putative toxin 55 domain, putative polysaccharide synthesis/export protein, putative neighboring cell growth inhibitor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,766,316	0	G	A	100.0%	230.8 / NA	71	*1378* (TAG→TAA)	rhsA	Rhs protein with putative toxin 55 domain, putative polysaccharide synthesis/export protein, putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (33/38);  total (33/38)

ATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCACA  >  NC_000913/3766180‑3766442                                                                                                                                         |                                                                                                                               aTGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAAcc                                                                                                                              >  2:1422859/1‑139 (MQ=255)  tGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAAccc                                                                                                                             <  1:2149/139‑1 (MQ=255)    ggAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCtt                                                                                                                           >  1:606723/1‑139 (MQ=255)         aaaCCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGa                                                                                                                      <  2:421333/139‑1 (MQ=255)                aCAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGa                                                                                                                      >  2:61036/1‑132 (MQ=255)                aCAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGa                                                                                                                      <  1:61036/132‑1 (MQ=255)                   gATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                            <  1:331647/139‑1 (MQ=255)                   gATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                            <  1:1422859/139‑1 (MQ=255)                    aTGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCaa                                                                                                           >  2:1396554/1‑139 (MQ=255)                        gACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTc                                                                                                       <  2:606723/139‑1 (MQ=255)                         aCGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTc                                                                                                       >  1:1218477/1‑138 (MQ=255)                         aCGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTc                                                                                                       <  2:1218477/138‑1 (MQ=255)                          cGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCaa                                                                                                     <  2:1375246/139‑1 (MQ=255)                               ttCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                            >  2:50789/1‑127 (MQ=255)                               ttCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                            <  1:50789/127‑1 (MQ=255)                                      aCTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCaga                                                                                         >  1:456759/1‑139 (MQ=255)                                         aaTTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACATAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCagaaaa                                                                                      >  1:760017/1‑139 (MQ=255)                                         aaTTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCagaaaa                                                                                      >  1:1429155/1‑139 (MQ=255)                                         aaTTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCagaaaa                                                                                      >  2:682965/1‑139 (MQ=255)                                         aaTTGATTGTGGTGAGATTAGTGCTAAAGATGCAAAACGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAACAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCagaaaa                                                                                      >  2:1008924/1‑139 (MQ=255)                                           ttgattgTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCagaaaaga                                                                                    <  2:456759/139‑1 (MQ=255)                                             gattgTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATGGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAATCCATGTGGAGGTTCCTGCAATGTCAAATACATACCtgaaaagaaa                                                                                  <  1:1396554/139‑1 (MQ=37)                                                tgtgGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGc                                                                               <  1:816635/139‑1 (MQ=255)                                                    gtgATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGt                                                                           >  2:808372/1‑139 (MQ=255)                                                       atatTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTaaa                                                                        <  1:682965/139‑1 (MQ=255)                                                           tAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                                                    <  2:1429155/139‑1 (MQ=255)                                                           tAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                                                    <  2:760017/139‑1 (MQ=255)                                                            aGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCa                                                                                           >  1:1245566/1‑115 (MQ=255)                                                            aGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCa                                                                                           <  2:1245566/115‑1 (MQ=255)                                                               gCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAg                                                                                                                     >  2:3535/1‑86 (MQ=255)                                                               gCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAg                                                                                                                     <  1:3535/86‑1 (MQ=255)                                                               gCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATgg                                                                >  2:491873/1‑139 (MQ=255)                                                                   aaGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAagag                                                                     >  1:67789/1‑130 (MQ=255)                                                                   aaGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAagag                                                                     <  2:67789/130‑1 (MQ=255)                                                                        gCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAata                                                                                                  <  2:1312734/96‑1 (MQ=255)                                                                        gCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAata                                                                                                  >  1:1312734/1‑96 (MQ=255)                                                                         cAAAAAGCACACAGACAGCCTGGAATTGTCGTCACTCCAGTCAGTCAAACGATAAAAAAAGATAACCCTTGCGGAGGTTCCTGCAATGTCAAATACAAACCCGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAAt                                                      <  1:1008924/139‑1 (MQ=18)                                                                         cAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAAt                                                      <  2:808088/139‑1 (MQ=255)                                                                                    caGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTc                                                                                                       <  1:1042694/79‑1 (MQ=255)                                                                                    caGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTc                                                                                                       >  2:1042694/1‑79 (MQ=255)                                                                                     aGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGaaa                                          <  1:491873/139‑1 (MQ=255)                                                                                      gAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGAAAc                                         >  1:1180367/1‑139 (MQ=255)                                                                                               gAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTtgtg                                                                                                                        <  2:600049/51‑1 (MQ=255)                                                                                               gAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTtgtg                                                                                                                        >  1:600049/1‑51 (MQ=255)                                                                                                     tcgtcACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGc                                                                               <  1:306449/86‑1 (MQ=255)                                                                                                     tcgtcACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGc                                                                               >  2:306449/1‑86 (MQ=255)                                                                                                     tcgtcACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGt                                                               <  1:330849/102‑1 (MQ=255)                                                                                                     tcgtcACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGt                                                               >  2:330849/1‑102 (MQ=255)                                                                                                     tcgtcACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGAAACTAAATGATGATGAat                          >  2:1464441/1‑139 (MQ=255)                                                                                                          aCTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGAAACTAAATGATGATGAATTAttt                     >  2:1219999/1‑139 (MQ=255)                                                                                                            tCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                            <  2:728307/50‑1 (MQ=255)                                                                                                            tCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                            >  1:728307/1‑50 (MQ=255)                                                                                                               aGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTta                                                            <  1:476662/95‑1 (MQ=255)                                                                                                               aGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTta                                                            >  2:476662/1‑95 (MQ=255)                                                                                                                     tCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                            >  1:371619/1‑41 (MQ=255)                                                                                                                     tCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                            <  2:371619/41‑1 (MQ=255)                                                                                                                     tCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCa                                                                                           >  2:1002851/1‑58 (MQ=255)                                                                                                                     tCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCa                                                                                           <  1:1002851/58‑1 (MQ=255)                                                                                                                     tCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGAAACTAAATGATGATGAATTATTTCg                   <  2:977437/130‑1 (MQ=255)                                                                                                                     tCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGAAACTAAATGATGATGAATTATTCCg                   >  1:977437/1‑130 (MQ=255)                                                                                                                      cAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTa                                                                          <  2:1270866/74‑1 (MQ=255)                                                                                                                      cAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTa                                                                          >  1:1270866/1‑74 (MQ=255)                                                                                                                      cAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAAt                                                      <  2:401501/94‑1 (MQ=255)                                                                                                                      cAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAAt                                                      >  1:401501/1‑94 (MQ=255)                                                                                                                      cAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTtgat         <  2:1180367/139‑1 (MQ=255)                                                                                                                      cAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTtgat         <  1:1219999/139‑1 (MQ=255)                                                                                                                             tAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCagaaaagaaa                                                                                  <  2:1022190/59‑1 (MQ=255)                                                                                                                             tAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCagaaaagaaa                                                                                  >  1:1022190/1‑59 (MQ=255)                                                                                                                             tAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTtg           >  1:768559/1‑130 (MQ=255)                                                                                                                             tAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTtg           <  2:768559/130‑1 (MQ=255)                                                                                                                             tAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATcaca  <  1:808372/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                               ATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCACA  >  NC_000913/3766180‑3766442

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑