Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	580,240	A→G	intergenic (+154/+11)	ybcW → / ← ylcI	DLP12 prophage, uncharacterized protein/DUF3950 family protein, DLP12 prophage

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	580,240	0	A	G	100.0%	230.6 / NA	65	intergenic (+154/+11)	ybcW/ylcI	DLP12 prophage, uncharacterized protein/DUF3950 family protein, DLP12 prophage
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (32/33);  total (32/33)

CATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTT  >  NC_000913/580104‑580370                                                                                                                                         |                                                                                                                                   cATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGcg                                                                                                                                  <  1:329397/139‑1 (MQ=255)    tCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCgg                                                                                                                               <  1:231069/139‑1 (MQ=255)     cAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCggg                                                                                                                              <  1:569810/139‑1 (MQ=255)                     tgACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTtcatc                                                                                                              <  1:1234383/139‑1 (MQ=255)                      gACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCgg                                                                                                                               <  1:863776/121‑1 (MQ=255)                      gACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCgg                                                                                                                               >  2:863776/1‑121 (MQ=255)                        cAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTTGTAAAACTGCAACTCAATTACTTCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACt                                                                                                           <  2:647805/139‑1 (MQ=255)                        cAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACt                                                                                                           <  2:496916/139‑1 (MQ=255)                           agCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGTGAAGAGCACGGGATGCTCATTCTTCATCACtttt                                                                                                        >  2:1311002/1‑139 (MQ=255)                           agCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTtca                                                                                                                >  2:872622/1‑131 (MQ=255)                           agCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTtca                                                                                                                <  1:872622/131‑1 (MQ=255)                           agCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACtttt                                                                                                        >  1:1008066/1‑139 (MQ=255)                           agCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACtttt                                                                                                        >  2:1182657/1‑139 (MQ=255)                               aaCACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAAtt                                                                                                    <  1:83850/139‑1 (MQ=255)                                   cGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATg                                                                                                <  1:1311002/139‑1 (MQ=255)                                      aGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCggg                                                                                                                              <  2:1067232/106‑1 (MQ=255)                                      aGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCggg                                                                                                                              >  1:1067232/1‑106 (MQ=255)                                      aGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTaa                                                                                                      >  1:294656/1‑130 (MQ=255)                                      aGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTaa                                                                                                      <  2:294656/130‑1 (MQ=255)                                                  tcaGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAAtt                                                                                                    <  2:827769/120‑1 (MQ=255)                                                  tcaGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAAtt                                                                                                    >  1:827769/1‑120 (MQ=255)                                                  tcaGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAAt                                                                                                     >  2:165658/1‑119 (MQ=255)                                                  tcaGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAAt                                                                                                     <  1:165658/119‑1 (MQ=255)                                                   caGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTc                                                                                <  2:28821/139‑1 (MQ=255)                                                           tGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTa                                                                        <  1:1182657/139‑1 (MQ=255)                                                                               tACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTCTCTGACGTTAGCCTCCGACGGCAGGCTTCa                                                    <  2:1008066/139‑1 (MQ=255)                                                                                    cAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTatat                                                                                           <  1:773498/95‑1 (MQ=255)                                                                                    cAATGCCCTCGTAATTAAGTAAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTatat                                                                                           >  2:773498/1‑95 (MQ=255)                                                                                      aTGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACtttt                                                                                                        >  2:329866/1‑80 (MQ=255)                                                                                      aTGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACtttt                                                                                                        <  1:329866/80‑1 (MQ=255)                                                                                      aTGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGAccc                                             >  1:875082/1‑139 (MQ=255)                                                                                      aTGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGAccc                                             >  1:990974/1‑139 (MQ=255)                                                                                        gCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAg                                           >  2:462843/1‑139 (MQ=255)                                                                                         cccTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTg                                                                                                   <  1:955954/82‑1 (MQ=255)                                                                                         cccTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTg                                                                                                   >  2:955954/1‑82 (MQ=255)                                                                                            tCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAgg                                                         >  1:221/1‑121 (MQ=255)                                                                                            tCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAgg                                                         <  2:221/121‑1 (MQ=255)                                                                                                     aGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATg                                                                                                                           >  2:411369/1‑46 (MQ=255)                                                                                                     aGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATg                                                                                                                           <  1:411369/46‑1 (MQ=255)                                                                                                     aGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGAGATGTTCATCCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTcc                              >  2:810048/1‑139 (MQ=255)                                                                                                          atgtACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGAc                         <  1:299314/136‑1 (MQ=255)                                                                                                                 atatCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCttttt                  >  1:83088/1‑139 (MQ=255)                                                                                                                     cGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACgg                                                             >  2:212030/1‑92 (MQ=255)                                                                                                                     cGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACgg                                                             <  1:212030/92‑1 (MQ=255)                                                                                                                     cGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTc              >  2:531171/1‑139 (MQ=255)                                                                                                                     cGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTc              >  2:538871/1‑139 (MQ=255)                                                                                                                     cGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTc              >  1:1264327/1‑139 (MQ=255)                                                                                                                     cGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTc              <  2:457778/139‑1 (MQ=255)                                                                                                                       tCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTaa                                                                                                      <  1:1034739/49‑1 (MQ=255)                                                                                                                       tCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTaa                                                                                                      >  2:1105068/1‑49 (MQ=255)                                                                                                                       tCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTaa                                                                                                      >  2:1034739/1‑49 (MQ=255)                                                                                                                       tCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTaa                                                                                                      <  1:1105068/49‑1 (MQ=255)                                                                                                                       tCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCaa            >  2:1245399/1‑139 (MQ=255)                                                                                                                        ccTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAag           >  2:792492/1‑139 (MQ=255)                                                                                                                          tGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAagag         <  2:1264327/139‑1 (MQ=255)                                                                                                                          tGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAagag         <  2:1084386/139‑1 (MQ=255)                                                                                                                          tGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAagag         <  2:83088/139‑1 (MQ=255)                                                                                                                          tGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAagag         <  1:332408/139‑1 (MQ=255)                                                                                                                          tGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAagag         <  2:990974/139‑1 (MQ=255)                                                                                                                              cGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGAt     >  2:1241664/1‑139 (MQ=255)                                                                                                                              cGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGAt     >  2:560823/1‑139 (MQ=255)                                                                                                                              cGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGAt     >  1:1123193/1‑139 (MQ=255)                                                                                                                                 aGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATg                                                                                          <  2:423136/51‑1 (MQ=255)                                                                                                                                 aGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATg                                                                                          >  1:423136/1‑51 (MQ=255)                                                                                                                                 aGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGtt  >  1:45040/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                   CATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTT  >  NC_000913/580104‑580370

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