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%	194.3 / NA	55	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 (26/29);  total (26/29)

ATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTT  >  NC_000913/580105‑580375                                                                                                                                        |                                                                                                                                        aTTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAgcgc                                                                                                                                      <  2:681128/139‑1 (MQ=255)  ttCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGcgcg                                                                                                                                     <  1:1269524/139‑1 (MQ=255)   tCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCgg                                                                                                                                    <  1:522105/139‑1 (MQ=255)    cAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCggg                                                                                                                                   <  1:725587/139‑1 (MQ=255)    ataagagacagaCCTGTGACAGAGCCAACACGCAGTCTGTCACTGGCAGGAAAGTGGTAAATCTACAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAAGGAAGAGCGCggg                                                                                                                                   <  2:1125694/128‑1 (MQ=255)           ctTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCa                                                                                                                            >  1:931795/1‑139 (MQ=255)                   gtgACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTtcat                                                                                                                    >  2:1316292/1‑139 (MQ=255)                          agCCAACATGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACtttt                                                                                                             >  2:458846/1‑139 (MQ=255)                          agCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCAttctctgtctcttat                                                                                                             >  1:1125694/1‑128 (MQ=255)                          agCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACtttt                                                                                                             >  1:499491/1‑139 (MQ=255)                                  cGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATg                                                                                                     <  1:1344137/139‑1 (MQ=255)                                  cGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATg                                                                                                     <  2:1049046/139‑1 (MQ=255)                                      gTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGtata                                                                                                 >  1:171730/1‑139 (MQ=255)                                          gtcactgtcaGGAAAGTGGTAAAACTGCAACTCAAGTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGct                                                                                             <  2:1359735/139‑1 (MQ=255)                                           tcactgtcaGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAgcgc                                                                                                                                      <  2:829990/97‑1 (MQ=255)                                           tcactgtcaGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAgcgc                                                                                                                                      >  1:829990/1‑97 (MQ=255)                                            cactgtcaGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGctct                                                                                           <  1:1316292/139‑1 (MQ=255)                                                 tcaGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCtttt                                                                                      >  1:283449/1‑139 (MQ=255)                                                  caGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTc                                                                                     <  2:1270227/139‑1 (MQ=255)                                                  caGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTc                                                                                     <  1:458846/139‑1 (MQ=255)                                                  caGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTc                                                                                     <  2:499491/139‑1 (MQ=255)                                                         gTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGtt                                                                              >  1:529034/1‑139 (MQ=255)                                                            gTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGCAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGc                                                                           >  1:569888/1‑139 (MQ=255)                                                               aaaCTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCggg                                                                                                                                   <  1:62367/80‑1 (MQ=255)                                                               aaaCTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCggg                                                                                                                                   >  2:62367/1‑80 (MQ=255)                                                                  cTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGa                                                                     <  2:931795/139‑1 (MQ=255)                                                                     cAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAAt                                                                                                          <  1:686851/99‑1 (MQ=255)                                                                     cAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAAt                                                                                                          >  2:686851/1‑99 (MQ=255)                                                                                   cAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGAc                                                    <  2:283449/139‑1 (MQ=255)                                                                                     aTGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAgcgc                                                                                                                                      >  1:417789/1‑55 (MQ=255)                                                                                     aTGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAgcgc                                                                                                                                      <  2:417789/55‑1 (MQ=255)                                                                                     aTGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGct                                                                                             >  1:1193113/1‑96 (MQ=255)                                                                                     aTGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGct                                                                                             <  2:1193113/96‑1 (MQ=255)                                                                                     aTGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGt                                                                               >  2:332439/1‑110 (MQ=255)                                                                                     aTGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGt                                                                               <  1:332439/110‑1 (MQ=255)                                                                                      tGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAAt                                                                                                          >  1:374285/1‑82 (MQ=255)                                                                                      tGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAAt                                                                                                          <  2:374285/82‑1 (MQ=255)                                                                                       gCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGc                                                             >  2:1420639/1‑126 (MQ=255)                                                                                       gCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGc                                                             <  1:1420639/126‑1 (MQ=255)                                                                                          cTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCt                                             >  2:685497/1‑139 (MQ=255)                                                                                           tCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACt                                                                                                                >  2:1388106/1‑71 (MQ=255)                                                                                           tCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACt                                                                                                                <  1:1388106/71‑1 (MQ=255)                                                                                                  tAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAAtt                                     >  2:684115/1‑139 (MQ=255)                                                                                                       gAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCgg                                >  1:1482656/1‑139 (MQ=255)                                                                                                       gAATTTACAATATCGTCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCgg                                <  1:684115/139‑1 (MQ=255)                                                                                                                   tCGTCCTGTTCGGAGTGAAGAGCGCGTGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCt                    <  2:569888/139‑1 (MQ=255)                                                                                                                      tCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACtttt                                                                                                             <  1:1016293/47‑1 (MQ=255)                                                                                                                      tCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACtttt                                                                                                             >  2:1016293/1‑47 (MQ=255)                                                                                                                      tCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTaa                                                                                                           <  2:658003/49‑1 (MQ=255)                                                                                                                      tCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTaa                                                                                                           >  1:658003/1‑49 (MQ=255)                                                                                                                      tCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGct                                                                                             <  2:1446833/63‑1 (MQ=255)                                                                                                                      tCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGct                                                                                             >  1:1446833/1‑63 (MQ=255)                                                                                                                      tCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGAc                              <  1:1298476/126‑1 (MQ=255)                                                                                                                      tCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGAc                              >  2:1298476/1‑126 (MQ=255)                                                                                                                      tCCTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCaa                 >  2:65682/1‑139 (MQ=255)                                                                                                                        cTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACgg                                                                  <  2:356021/88‑1 (MQ=255)                                                                                                                        cTGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACgg                                                                  >  1:356021/1‑88 (MQ=255)                                                                                                                         tGTTCGGAGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAagag              <  1:685497/139‑1 (MQ=255)                                                                                                                                aGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGa                                                                                                       <  1:101075/43‑1 (MQ=255)                                                                                                                                aGGGAAGAGCGCGGGATGTTCATTCTTCATCACTTTTAATTGa                                                                                                       >  2:101075/1‑43 (MQ=255)                                                                                                                                     agagCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAAttt  >  2:337399/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                        ATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTT  >  NC_000913/580105‑580375

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