Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,475,562	(G)6→7	coding (126/1815 nt)	yjgL →	SopA‑central‑domain‑like hexapeptide repeat protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,475,556	1	.	G	100.0%	178.5 / NA	51	A40A (GCC→GCG)	yjgL	SopA‑central‑domain‑like hexapeptide repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (27/24);  total (27/24)

CACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑G‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTCA  >  NC_000913/4475420‑4475687                                                                                                                                          |                                                                                                                                     cACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑g‑g                                                                                                                                   <  2:206239/139‑1 (MQ=255)      aaGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑g‑gggggg                                                                                                                              >  2:807406/1‑138 (MQ=255)      aaGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑g‑gggggg                                                                                                                              <  1:959788/139‑2 (MQ=255)      aaGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑g‑gggggg                                                                                                                              >  1:838604/1‑138 (MQ=255)       aGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑g‑ggggggc                                                                                                                             >  1:47794/1‑137 (MQ=255)         gAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑g‑ggggggcag                                                                                                                           <  2:995686/139‑5 (MQ=255)              cGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATaa                                                                                                                       >  2:840395/1‑139 (MQ=255)              cGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATaa                                                                                                                       >  2:256464/1‑139 (MQ=255)              cGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATaa                                                                                                                       >  2:229751/1‑139 (MQ=255)              cGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATaa                                                                                                                       >  1:641668/1‑139 (MQ=255)              cGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑g‑gggggg                                                                                                                              >  2:851876/1‑130 (MQ=255)              cGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑g‑gggggg                                                                                                                              <  1:851876/131‑2 (MQ=255)                    gAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAGAAg                                                                                                                 <  1:984959/139‑1 (MQ=255)                                    ttAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCt                                                                                                 >  1:103117/1‑139 (MQ=255)                                        attattCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAg                                                                                             <  1:229751/139‑1 (MQ=255)                                            ttCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGAATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCAtt                                                                                         >  2:588896/1‑139 (MQ=255)                                                      attacattaGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGttt                                                                               >  1:38188/1‑139 (MQ=255)                                                      attacattaGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑g‑gggggg                                                                                                                              <  2:749758/91‑2 (MQ=255)                                                      attacattaGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑g‑gggggg                                                                                                                              >  1:749758/1‑90 (MQ=255)                                                        tacattaGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGttttt                                                                             >  1:19792/1‑139 (MQ=255)                                                        tacattaGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGttttt                                                                             >  1:236541/1‑139 (MQ=255)                                                             taGCCTACAATTTTAAAAAAAAAAGTGAATTTTTAAATTCCGGGCGTGTTGGAATGAATGATTCTGCCCGTATTGCCGG‑GGGGGCAGGATAACAGAATGAATATTCTTTCTCCTGGAGCATTTTGAGAGTGTTTGGCaa                                                                        <  2:581730/139‑1 (MQ=255)                                                             taGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCaa                                                                        >  1:447737/1‑139 (MQ=255)                                                               gCCGACAATTTTAAACAAAAAAGTGTAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑g‑gggggg                                                                                                                              <  2:413366/82‑2 (MQ=255)                                                               gCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGCGTTGGAATGAATGATTTTGCCCGTATTGCC‑g‑gggggg                                                                                                                              >  1:413366/1‑81 (MQ=255)                                                                ccGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAAc                                                                                                                      <  2:606417/90‑1 (MQ=255)                                                                ccGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAAc                                                                                                                      >  1:606417/1‑90 (MQ=255)                                                                ccGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCaa                                                                        >  2:845757/1‑136 (MQ=255)                                                                ccGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCaa                                                                        <  1:845757/136‑1 (MQ=255)                                                                ccGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGa                                                                     >  2:786390/1‑139 (MQ=255)                                                                  gACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGAATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGata                                                                   <  1:588896/139‑1 (MQ=255)                                                                   aCAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGatat                                                                  <  1:786390/139‑1 (MQ=255)                                                                   aCAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGatat                                                                  <  2:641668/139‑1 (MQ=255)                                                                   aCAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGatat                                                                  <  1:368952/139‑1 (MQ=255)                                                                                 aaaaGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTAt                                                    >  2:166409/1‑139 (MQ=255)                                                                                 aaaaGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTAt                                                    <  2:838604/139‑1 (MQ=255)                                                                                 aaaaGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTAt                                                    >  1:531889/1‑139 (MQ=255)                                                                                        aaGTTTTAAATACATGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGAtt                                             >  2:712276/1‑139 (MQ=255)                                                                                                aaTACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCaaa                                     <  2:47794/139‑1 (MQ=255)                                                                                                       ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGa                                                                                                                          >  2:501287/1‑47 (MQ=255)                                                                                                       ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGa                                                                                                                          <  1:884137/47‑1 (MQ=255)                                                                                                       ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGa                                                                                                                          >  2:884137/1‑47 (MQ=255)                                                                                                       ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGa                                                                                                                          <  1:501287/47‑1 (MQ=255)                                                                                                       ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAg                                                                                                                    <  1:798800/53‑1 (MQ=255)                                                                                                       ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAg                                                                                                                    >  2:798800/1‑53 (MQ=255)                                                                                                       ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATAttt                                                                <  1:415915/105‑1 (MQ=255)                                                                                                       ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATAttt                                                                >  2:415915/1‑105 (MQ=255)                                                                                                       ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCa                              >  2:188740/1‑139 (MQ=255)                                                                                                            gtTGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGaa                         >  2:734967/1‑139 (MQ=255)                                                                                                             ttGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAg                        <  1:188740/139‑1 (MQ=255)                                                                                                              tGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCa                                       >  2:901518/1‑123 (MQ=255)                                                                                                              tGGAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCa                                       <  1:901518/123‑1 (MQ=255)                                                                                                               ggAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAg                                                                                                                    <  1:1052458/45‑1 (MQ=255)                                                                                                               ggAATGAATGATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAg                                                                                                                    >  2:1052458/1‑45 (MQ=255)                                                                                                                 aatgaatgATTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCggg                    <  1:166409/139‑1 (MQ=255)                                                                                                                  atgaatgaTTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGt                   <  2:531889/139‑1 (MQ=255)                                                                                                                   tgaatgaTTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTa                  <  1:807406/139‑1 (MQ=255)                                                                                                                     aatgaTTTTGCCCGTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGa                >  2:393653/1‑139 (MQ=255)                                                                                                                             tGCCCGTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATatgat        <  1:393653/139‑1 (MQ=255)                                                                                                                               cccGTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGc                       <  1:266368/122‑1 (MQ=255)                                                                                                                               cccGTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGc                       >  2:266368/1‑122 (MQ=255)                                                                                                                               cccGTATTGCC‑GTCGGGGCAAGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTGTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAatatcaggg      <  2:386249/139‑6 (MQ=255)                                                                                                                                  gTATTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAATATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTc   <  1:420108/139‑1 (MQ=255)                                                                                                                                   aaTTGCCGG‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTACTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTCa  >  1:649194/2‑139 (MQ=255)                                                                                                                                          |                                                                                                                                     CACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑G‑GGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTCA  >  NC_000913/4475420‑4475687

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