Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	572,454	C→T	C29C (TGC→TGT)	ylcH →	uncharacterized protein, DLP12 prophage

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	572,454	0	C	T	100.0%	144.5 / NA	45	C29C (TGC→TGT)	ylcH	uncharacterized protein, DLP12 prophage
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/20);  total (25/20)

GGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAACC  >  NC_000913/572322‑572578                                                                                                                                     |                                                                                                                             ggAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAggg                                                                                                                        <  2:1043807/139‑1 (MQ=255)  gAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAgggg                                                                                                                       >  2:1073180/1‑139 (MQ=255)        gAACAACTCCGGTGCCTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTACTCCCTTCACAGATGCAGAAACGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGt                                                                                                                 <  1:806413/139‑1 (MQ=255)           caacTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGaa                                                                                                              <  1:21903/139‑1 (MQ=255)           caacTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGaa                                                                                                              <  1:602774/139‑1 (MQ=255)                   gTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTcac                                                                                                      >  1:206123/1‑139 (MQ=255)                        ttCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTACAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGAc                                                                                                 >  2:804540/1‑139 (MQ=255)                            gcgcTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa                                                                                             <  1:1243373/139‑1 (MQ=255)                                  aaCGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAAt                                                                                       >  1:230023/1‑139 (MQ=255)                                  aaCGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAAt                                                                                       >  1:924490/1‑139 (MQ=255)                                       ggACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGt                                                                                                                 >  2:801399/1‑108 (MQ=255)                                       ggACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGt                                                                                                                 <  1:801399/108‑1 (MQ=255)                                                   gcACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTc                                                                      <  2:230023/139‑1 (MQ=255)                                                      caTACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCAcc                                                                   >  1:505170/1‑139 (MQ=255)                                                               cAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa                                                                                             >  1:587084/1‑104 (MQ=255)                                                               cAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa                                                                                             <  2:204410/104‑1 (MQ=255)                                                               cAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa                                                                                             <  2:587084/104‑1 (MQ=255)                                                               cAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa                                                                                             >  1:204410/1‑104 (MQ=255)                                                               cAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTAtcggt                                                          >  2:50012/1‑139 (MQ=255)                                                                 aaCTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGaa                                                                                                              <  2:386761/85‑1 (MQ=255)                                                                 aaCTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGaa                                                                                                              >  1:386761/1‑85 (MQ=255)                                                                 aaCTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCGAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTAtcggtcg                                                        >  1:888547/1‑139 (MQ=255)                                                                        ttctCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGAta                                                 <  2:825066/139‑2 (MQ=255)                                                                             ccTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGcc                                            <  2:924490/139‑1 (MQ=255)                                                                             ccTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGcc                                            <  2:206123/139‑1 (MQ=255)                                                                             ccTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGcc                                            <  2:888547/139‑1 (MQ=255)                                                                                tCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTa                                         >  2:90891/1‑139 (MQ=255)                                                                                tCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTa                                         >  1:120964/1‑139 (MQ=255)                                                                                  acaGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATc                                       >  2:553918/1‑139 (MQ=255)                                                                                     gATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTgg                                    >  2:18507/1‑139 (MQ=255)                                                                                       tGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACgg                                                                                               >  1:518182/1‑78 (MQ=255)                                                                                       tGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACgg                                                                                               <  2:518182/78‑1 (MQ=255)                                                                                       tGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTggag                                  >  2:1266482/1‑139 (MQ=255)                                                                                                gCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGc                                             <  2:1098846/119‑1 (MQ=255)                                                                                                gCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGc                                             >  1:1098846/1‑119 (MQ=255)                                                                                                 cACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAAt                        >  2:695182/1‑139 (MQ=255)                                                                                                 cACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAAt                        >  1:1280555/1‑139 (MQ=255)                                                                                                    cTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCt                     >  1:657734/1‑139 (MQ=255)                                                                                                           tCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCAt              <  2:505170/139‑1 (MQ=255)                                                                                                           tCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCAt              <  1:90891/139‑1 (MQ=255)                                                                                                            cTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATg             <  2:1106748/139‑1 (MQ=255)                                                                                                              ttGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGt           >  2:23754/1‑139 (MQ=255)                                                                                                                gCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGc         >  1:119648/1‑139 (MQ=255)                                                                                                                       gCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAAcc  >  2:516840/1‑139 (MQ=255)                                                                                                                       gCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAAcc  <  1:50012/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                             GGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAACC  >  NC_000913/572322‑572578

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