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%	207.7 / NA	64	*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 (28/36);  total (28/36)

AAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGAAACTAAATGA  >  NC_000913/3766187‑3766410                                                                                                                                  |                                                                                               aaaaCCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTgg                                                                                       <  2:832246/139‑1 (MQ=255)       cGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTc                                                                                 >  2:1173183/1‑139 (MQ=255)          cAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCGAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTg                                                                              <  2:1270996/139‑1 (MQ=255)            gttgtgACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                            <  2:1049740/137‑1 (MQ=255)            gATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                            <  1:899174/139‑1 (MQ=255)            gATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                            <  1:1188457/139‑1 (MQ=255)                   cGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCaa                                                                     <  1:735542/139‑1 (MQ=255)                    gTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCaaa                                                                    >  2:274838/1‑139 (MQ=255)                             gAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTc                                                                       >  1:1224194/1‑127 (MQ=255)                             gAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTc                                                                       <  2:1224194/127‑1 (MQ=255)                                cTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACAagaa                                                        >  1:368275/1‑139 (MQ=255)                                  aaTTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                            <  2:1231645/117‑1 (MQ=255)                                  aaTTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                            >  1:1231645/1‑117 (MQ=255)                                  aaTTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCagaaaa                                                      >  1:1208881/1‑139 (MQ=255)                                  aaTTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCagaaaa                                                      >  1:91675/1‑139 (MQ=255)                                    ttgattgTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCagaaaaga                                                    <  1:1173183/139‑1 (MQ=255)                                        ttgTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAgg                                                <  2:10609/139‑1 (MQ=255)                                               gATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTAAAATACATACCAGAAAAGAAAGGCAAGTaa                                         >  1:1321652/1‑139 (MQ=255)                                                atatTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTaaa                                        <  2:1208881/139‑1 (MQ=255)                                                    taatGCTAAAGATTCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAATCCGTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                    <  1:904326/136‑1 (MQ=21)                                                    tAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                    <  2:1101219/139‑1 (MQ=255)                                                     aGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGa                                                                                      >  1:995666/1‑88 (MQ=255)                                                     aGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGa                                                                                      <  2:131398/88‑1 (MQ=255)                                                     aGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGa                                                                                      <  2:995666/88‑1 (MQ=255)                                                     aGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGa                                                                                      >  1:131398/1‑88 (MQ=255)                                                      gTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                    <  1:130816/137‑1 (MQ=255)                                                      gTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                    >  2:130816/1‑137 (MQ=255)                                                       tGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAataca                                                                <  1:897134/108‑1 (MQ=255)                                                       tGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAataca                                                                >  2:897134/1‑108 (MQ=255)                                                       tGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTaaa                                        >  2:1129737/1‑132 (MQ=255)                                                       tGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTaaa                                        <  1:1129737/132‑1 (MQ=255)                                                        gCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGa                                                                                      <  1:316937/85‑1 (MQ=255)                                                        gCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGa                                                                                      >  2:316937/1‑85 (MQ=255)                                                         cTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGt                               <  2:1130327/139‑1 (MQ=255)                                                            aaGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTtgtg                                                                                        <  1:472160/79‑1 (MQ=255)                                                            aaGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTtgtg                                                                                        >  2:472160/1‑79 (MQ=255)                                                                 gCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTa                                          >  2:34652/1‑120 (MQ=255)                                                                 gCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTa                                          <  1:34652/120‑1 (MQ=255)                                                                  cAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAAt                      <  1:205639/139‑1 (MQ=255)                                                                  cAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAAt                      <  1:274838/139‑1 (MQ=255)                                                                         cacacaGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                            <  2:589244/78‑1 (MQ=255)                                                                         cacacaGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                            >  1:589244/1‑78 (MQ=255)                                                                          acacaGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                    >  1:396730/1‑117 (MQ=255)                                                                          acacaGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGt                                    <  2:396730/117‑1 (MQ=255)                                                                           cacaGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAg             <  2:1321652/139‑1 (MQ=255)                                                                            acaGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                            >  2:1005189/1‑75 (MQ=255)                                                                            acaGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                            <  1:1005189/75‑1 (MQ=255)                                                                            acaGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCagaaaaga                                                    <  1:949927/99‑1 (MQ=255)                                                                            acaGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCagaaaaga                                                    >  2:949927/1‑99 (MQ=255)                                                                                aaaGCCTGGAATTGTCGTCACTCCAGACAGTCAAACTATAAAAAAAGTTAACCCTTGTGGa                                                                                      >  2:869904/1‑61 (MQ=37)                                                                                aaaGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGa                                                                                      <  1:869904/61‑1 (MQ=255)                                                                                    ccTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAAt                      >  1:413122/1‑121 (MQ=255)                                                                                    ccTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAAt                      <  2:413122/121‑1 (MQ=255)                                                                                      tGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGAAACTAAatga  >  1:694734/1‑139 (MQ=255)                                                                                             gtcgtcACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGt                                           <  2:1369674/91‑1 (MQ=255)                                                                                             gtcgtcACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGt                                           >  1:1369674/1‑91 (MQ=255)                                                                                                  cACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTaaa                                        >  2:785947/1‑89 (MQ=255)                                                                                                  cACTCCAGACAGTCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTaaa                                        <  1:785947/89‑1 (MQ=255)                                                                                                              tCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAAt                                                                          <  1:1165186/43‑1 (MQ=255)                                                                                                              tCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAAt                                                                          >  2:1165186/1‑43 (MQ=255)                                                                                                              tCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTa                                   <  1:1235149/82‑1 (MQ=255)                                                                                                              tCAAACGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTa                                   >  2:1235149/1‑82 (MQ=255)                                                                                                                   cGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                            <  2:399680/36‑1 (MQ=255)                                                                                                                   cGATAAAAAAAGATAACCCTTGTGGAGGTTCCTGCa                                                                            >  1:399680/1‑36 (MQ=255)                                                                                                                                  |                                                                                               AAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATCAATGTAAGAAACTAAATGA  >  NC_000913/3766187‑3766410

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: ‑