Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,028,416	A→G	intergenic (‑46/+33)	yodC ← / ← yedI	uncharacterized protein/DUF808 family inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,028,416	0	A	G	100.0%	210.9 / NA	56	intergenic (‑46/+33)	yodC/yedI	uncharacterized protein/DUF808 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (36/20);  total (36/20)

CGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACT  >  NC_000913/2028268‑2028563                                                                                                                                                     |                                                                                                                                                    cGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTg                                                                                                                                                   >  1:856988/1‑151 (MQ=255)      ttccAAAGCCACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGCCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGc                                                                                                                                                 >  2:43069/3‑148 (MQ=255)       aCCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTACTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAAt                                                                                                                                             >  2:398382/1‑151 (MQ=255)          aaCGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAAtttt                                                                                                                                          <  1:1048363/151‑1 (MQ=255)             gaaacTCAACAATACCGCTGGAGTATCAAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTAACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGc                                                                                                                                        >  1:372666/4‑150 (MQ=255)             gACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGc                                                                                                                                        >  1:372678/1‑150 (MQ=255)               cacTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCa                                                                                                                                       <  1:129965/149‑1 (MQ=255)                    aaCCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGc                                                                                                                                <  2:16195/151‑1 (MQ=255)                          aCCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTg                                                                                                                           <  2:288825/150‑1 (MQ=255)                          aCCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGc                                                                                                                          >  1:934330/1‑151 (MQ=255)                           ccGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGc                                                                                                                          <  1:665923/150‑1 (MQ=255)                            cGCTGGAGTATCCAGTGACAATCATCCGCGGGCAGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCg                                                                                                                         >  1:703854/1‑150 (MQ=255)                             gCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCg                                                                                                                         >  1:137366/1‑149 (MQ=255)                             gCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCg                                                                                                                         >  2:260901/1‑149 (MQ=255)                                ggAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGt                                                                                                                      <  1:740439/149‑1 (MQ=255)                                  aGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTa                                                                                                                  <  2:894023/151‑1 (MQ=255)                                  aGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTa                                                                                                                  <  1:644272/151‑1 (MQ=255)                                   gTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTaa                                                                                                                 <  2:985113/151‑1 (MQ=255)                                     aTCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGTGAGTATTAAt                                                                                                                <  2:243504/150‑1 (MQ=255)                                         aGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGcc                                                                                                           >  2:208541/1‑151 (MQ=255)                                         aGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGcc                                                                                                           >  1:876132/1‑151 (MQ=255)                                             aCAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGAc                                                                                                       >  2:672081/1‑151 (MQ=255)                                                     ccGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCAtt                                                                                               <  2:876873/151‑1 (MQ=255)                                                       gcgGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACCTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCAtttt                                                                                             >  2:859288/1‑151 (MQ=255)                                                          tgccgccCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTtctc                                                                                          >  2:348586/2‑151 (MQ=255)                                                             cgccCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTAc                                                                                       >  1:58581/1‑151 (MQ=255)                                                                ccTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGc                                                                                    >  2:276130/1‑151 (MQ=255)                                                                     ttAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCtt                                                                                  <  2:860864/148‑1 (MQ=255)                                                                         ccGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTcc                                                                           >  1:452319/1‑151 (MQ=255)                                                                                ttCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGcacc                                                                    >  2:580379/1‑151 (MQ=255)                                                                                   cTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGcaccac                                                                  <  2:892685/150‑1 (MQ=255)                                                                                     cACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGAt                                                               >  1:269293/1‑151 (MQ=255)                                                                                               aaaGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCgat                                                      <  1:672472/150‑1 (MQ=255)                                                                                               aaaGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGGGGCGAATTTTGCATAATAGCCAGTGCGATTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCgatg                                                     >  2:379954/1‑151 (MQ=255)                                                                                                aaGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCgatga                                                    <  1:759807/151‑1 (MQ=255)                                                                                                  gCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCgatgat                                                   >  1:835750/1‑150 (MQ=255)                                                                                                     cATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAAt                                               <  2:931506/151‑1 (MQ=255)                                                                                                     cATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAAt                                               >  2:854180/1‑151 (MQ=255)                                                                                                      aTACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATc                                              <  1:228758/151‑1 (MQ=255)                                                                                                        aCTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATcc                                             >  1:732209/1‑150 (MQ=255)                                                                                                               tcctTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATc                                      >  2:659908/1‑150 (MQ=255)                                                                                                                cctTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATc                                      >  2:279226/1‑149 (MQ=255)                                                                                                                 ctTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATc                                      >  1:1060379/1‑148 (MQ=255)                                                                                                                  ttttttCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATc                                      >  1:643817/1‑147 (MQ=255)                                                                                                                   tttttCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGa                                   <  1:349150/149‑1 (MQ=255)                                                                                                                      ttCGAGTGAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTa                               >  2:615313/1‑150 (MQ=255)                                                                                                                             gAAACTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGt                        <  2:446229/150‑1 (MQ=255)                                                                                                                               aaCTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCgg                     >  2:464316/1‑151 (MQ=255)                                                                                                                                aCTTGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGt                    >  1:692707/1‑151 (MQ=255)                                                                                                                                  ttGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTa                   >  1:621589/1‑150 (MQ=255)                                                                                                                                   tGTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAat                 >  2:473026/1‑151 (MQ=255)                                                                                                                                    gTTCACCTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCTCGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAata                <  2:934796/151‑1 (MQ=255)                                                                                                                                         ccTTAGTTGAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATatcatc           >  2:314247/1‑151 (MQ=255)                                                                                                                                                 gAAGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCactac   >  1:185170/1‑151 (MQ=255)                                                                                                                                                  aaGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCactact  >  1:917430/1‑151 (MQ=255)                                                                                                                                                  aaGGTGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCactact  >  1:754486/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCATAATAGCCAGTGCGAGTATTAATGTGCCTGACCGCGCATTTTCTCTACGGCTTTCACTCCCAGCACCACGATGCCGCCGATGATAAATCCAAGAATCAGATTTAAAACAGTCGGTAATATCATCGCCACTACT  >  NC_000913/2028268‑2028563

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