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%	105.8 / NA	34	C29C (TGC→TGT)	ylcH	uncharacterized protein, DLP12 prophage
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/17);  total (17/17)

GAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAACCGTGGC  >  NC_000913/572323‑572583                                                                                                                                    |                                                                                                                                  gAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAgggg                                                                                                                            <  1:928185/139‑1 (MQ=255)            acTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAAcc                                                                                                                 >  2:789937/1‑139 (MQ=255)             cTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACct                                                                                                                >  1:776378/1‑139 (MQ=255)                       ttCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTAACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGAc                                                                                                      >  2:294822/1‑139 (MQ=255)                           gcgcTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa                                                                                                  <  2:776378/139‑1 (MQ=255)                               tAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTg                                                                                                                     >  1:254577/1‑116 (MQ=255)                               tAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTg                                                                                                                     <  2:254577/116‑1 (MQ=255)                                    ggggACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTAc                                                                                         >  1:984390/1‑139 (MQ=255)                                               tGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAAc                                                                              <  2:984390/139‑1 (MQ=255)                                                    acaTACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCAc                                                                         <  1:1013607/139‑1 (MQ=255)                                                             ccAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACgg                                                                                                    >  1:760720/1‑103 (MQ=255)                                                             ccAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTAtcgg                                                                <  1:789937/139‑1 (MQ=255)                                                             ccAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGAGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACgg                                                                                                    <  2:760720/103‑1 (MQ=255)                                                                  ctctctTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGc                                                           <  1:294822/139‑1 (MQ=255)                                                                       ttctCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATc                                                      >  2:665401/1‑139 (MQ=255)                                                                           cccTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGc                                                  <  1:372254/139‑1 (MQ=255)                                                                            ccTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGcc                                                 <  2:33939/139‑1 (MQ=255)                                                                               tCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTa                                              >  2:295682/1‑139 (MQ=255)                                                                                      tGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAAGCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGcc                                                 >  2:742460/1‑129 (MQ=255)                                                                                      tGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCAc                                                                         <  2:845616/105‑1 (MQ=255)                                                                                      tGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCAc                                                                         >  1:845616/1‑105 (MQ=255)                                                                                      tGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGcc                                                 <  1:742460/129‑1 (MQ=255)                                                                                                cACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGc                                                  <  2:122307/118‑1 (MQ=255)                                                                                                cACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGc                                                  >  1:122307/1‑118 (MQ=255)                                                                                                cACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAAt                             >  1:824738/1‑139 (MQ=255)                                                                                                     gCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa                                                                                                  <  2:387971/65‑1 (MQ=255)                                                                                                     gCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa                                                                                                  >  1:387971/1‑65 (MQ=255)                                                                                                        aTTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGc                     >  1:363220/1‑139 (MQ=255)                                                                                                          tCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCAt                   <  1:1038259/139‑1 (MQ=255)                                                                                                                     aGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCTTCAAATTACATCGCGCTGACTTCTACGCTAACGGTCGGCAGATCAAGCCTTATGTGGAGGAGGGCAAGTGCTTTCGCATGGTGCTTAAAc        >  2:3887/1‑139 (MQ=255)                                                                                                                     aGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAAc        >  1:834458/1‑139 (MQ=255)                                                                                                                           tGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAACCGTGGc  <  2:824738/139‑1 (MQ=255)                                                                                                                             aCCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCAt                   >  2:488090/1‑120 (MQ=255)                                                                                                                             aCCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCAt                   <  1:488090/120‑1 (MQ=255)                                                                                                                                    |                                                                                                                                  GAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAACCGTGGC  >  NC_000913/572323‑572583

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