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%	189.2 / NA	59	*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 (27/32);  total (27/32)

CCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCACAATTCCTTG  >  NC_000913/3766191‑3766450                                                                                                                              |                                                                                                                                       cccGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGt                                                                                                                           <  1:39516/139‑1 (MQ=255)        gATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                                    <  2:1174607/139‑1 (MQ=255)        gATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                                    <  2:533464/139‑1 (MQ=255)        gATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                                    <  1:476630/139‑1 (MQ=255)         aTGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGAT‑AAAAAAGATAACCCTTGTGGAGGTTCCTGCAAt                                                                                                                  <  2:463448/139‑1 (MQ=255)                   cTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAatac                                                                                                         >  2:656435/1‑139 (MQ=255)                    ttCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTtgtg                                                                                                                                >  2:595392/1‑115 (MQ=255)                    ttCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTtgtg                                                                                                                                <  1:595392/115‑1 (MQ=255)                    ttCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAataca                                                                                                        >  1:470029/1‑139 (MQ=255)                          aaCTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCag                                                                                                  <  2:1135397/139‑1 (MQ=255)                              aaTTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCaaaaaa                                                                                              >  1:235882/1‑139 (MQ=255)                                    ttgTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAgg                                                                                        <  2:235882/139‑1 (MQ=255)                                      gtggtgATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCa                                                                                      <  2:839766/139‑1 (MQ=255)                                      gtggtgATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCa                                                                                      <  2:470029/139‑1 (MQ=255)                                       tggtgATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCaa                                                                                     >  1:604785/1‑139 (MQ=255)                                       tggtgATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCaa                                                                                     <  1:656435/139‑1 (MQ=255)                                          tgATATTAGTACTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTa                                                                                  <  2:631945/139‑1 (MQ=255)                                            atatTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTaaa                                                                                <  2:1029565/139‑1 (MQ=255)                                             tatTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAag                                                                               >  2:242688/1‑139 (MQ=255)                                                tAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAGTACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                                                            <  1:779785/139‑1 (MQ=255)                                                 aGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAAGCAAGTAAAGAGTa                                                                           >  2:527745/1‑139 (MQ=255)                                                  gTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTAt                                                                          >  2:999550/1‑139 (MQ=255)                                                   tGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCagaaaagaaa                                                                                          >  2:573931/1‑122 (MQ=255)                                                   tGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCagaaaagaaa                                                                                          <  1:573931/122‑1 (MQ=255)                                                   tGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATg                                                                         >  2:934221/1‑139 (MQ=255)                                                    gCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATgg                                                                        <  1:999550/139‑1 (MQ=255)                                                      tAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTc                                                                                                               <  1:10275/98‑1 (MQ=255)                                                      tAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTc                                                                                                               >  2:10275/1‑98 (MQ=255)                                                        aaGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                                    <  1:914444/91‑1 (MQ=255)                                                        aaGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                                    >  2:914444/1‑91 (MQ=255)                                                         aGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTtat                                                                   >  2:592467/1‑139 (MQ=255)                                                             gCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATaa                                                               >  2:743803/1‑139 (MQ=255)                                                              cAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAAt                                                              <  1:934221/139‑1 (MQ=255)                                                                      acacaGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                                    <  2:274102/77‑1 (MQ=255)                                                                      acacaGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                                    >  1:274102/1‑77 (MQ=255)                                                                             aaGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCagaaaaga                                                                                            <  2:1072256/94‑1 (MQ=255)                                                                             aaGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCagaaaaga                                                                                            >  1:1072256/1‑94 (MQ=255)                                                                             aaGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTaaa                                                                                <  2:282692/106‑1 (MQ=255)                                                                             aaGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTaaa                                                                                >  1:282692/1‑106 (MQ=255)                                                                             aaGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTta                                                                    >  2:264503/1‑118 (MQ=255)                                                                             aaGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTta                                                                    <  1:264503/118‑1 (MQ=255)                                                                               gCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                                                            <  1:144626/108‑1 (MQ=255)                                                                               gCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                                                            >  2:144626/1‑108 (MQ=255)                                                                                         gtcgtcACTCCAGACAGTCAAACGATAAAAAAAGATAACTCTTGTGGAGGTTCCTGCa                                                                                                                    <  1:501653/58‑1 (MQ=18)                                                                                         gtcgtcACTCCAGACAGTCAAACGATAAAAAAAGATAACTCTTGTGGAGGTTCCTGCa                                                                                                                    >  2:501653/1‑58 (MQ=21)                                                                                          tcgtcACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTc                                                                                                               >  2:568746/1‑62 (MQ=255)                                                                                          tcgtcACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTc                                                                                                               <  1:568746/62‑1 (MQ=255)                                                                                           cgtcACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGAAACTAAATGATGATGAatt                                 <  1:743803/139‑1 (MQ=255)                                                                                              cACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGc                                                                                       >  2:689219/1‑82 (MQ=255)                                                                                              cACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGc                                                                                       <  1:689219/82‑1 (MQ=255)                                                                                                    aGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                                                            >  1:749350/1‑87 (MQ=255)                                                                                                    aGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGAAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAagagc                                                                            <  2:749350/87‑2 (MQ=255)                                                                                                           cacaCGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                                                            >  2:971676/4‑80 (MQ=255)                                                                                                           cAAACGATAAAAAAAGATCACCCTTGTGGAGGTTCCTGCa                                                                                                                    <  1:733043/40‑1 (MQ=255)                                                                                                           cAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                                    >  2:733043/1‑40 (MQ=255)                                                                                                           cAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                                    >  2:531569/1‑40 (MQ=255)                                                                                                           cAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                                                                    <  1:531569/40‑1 (MQ=255)                                                                                                           cAAAAGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                                                            <  1:971676/80‑1 (MQ=255)                                                                                                                          gATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCACAATTCCTTg  >  2:156435/1‑139 (MQ=255)                                                                                                                              |                                                                                                                                       CCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGAAACTAAATGATGATGAATTATTTCGCTTACTTGATGATCACAATTCCTTG  >  NC_000913/3766191‑3766450

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