Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	571,663	T→C	Y70Y (TAT→TAC)	ybcM →	DLP12 prophage; putative DNA‑binding transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	571,663	0	T	C	100.0%	105.8 / NA	29	Y70Y (TAT→TAC)	ybcM	DLP12 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/14);  total (15/14)

TTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGTCAGGGAAAGAAAGTTATCAATTGGGAGACGCATACTTTAACTGTTGATAATAATCAAATAATATTATTTCCTTGTGGTTATGAATTTTATATTGAGAATTATCCTGAAGCAGGGCTTTATCTTGCAGAAATGCTTTACTTACCCATTGATTTAATTGAGAGTTTCCAAAAACTTTATACGGTAACTGATCAAATACGTAACAAAACAAGTTTCTTTTTACCTCAGAATCCTGAGTTA  >  NC_000913/571527‑571798                                                                                                                                         |                                                                                                                                        ttAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGTCAGGGAAAGAAAGTTATCAATTGGGAGACGCATACTTTAACTGTTGATAATAATCAAATAATATTATTTCCTTGTGGTTATGAATTTTATATTGAGAATTACCCt                                                                                                                                      <  2:2001371/140‑1 (MQ=255)     gCCTCTTATGCCAGCTTTGTGTCGTATCCGTCAGGGAAAGAAAGTTATCAATTGGGAGACGCATACTTTAACTGTTGATAATAATCAAATAATATTATTTCCTTGTGGTTATGAATTTTATATTGAGAATTACCCTGAAg                                                                                                                                  <  2:2118949/140‑1 (MQ=255)                    tttGTGTCGTATCCGTCAGGGAAAGAAAGTTATCAATTGGGAGACGCATACTTTAACTGTTGATAATAATCAAATAATATTATTTCCTTGTGGTTATGAATTTTATATTGAGAATTACCCTGAAGCAGGGCTTTATCTTg                                                                                                                   >  1:1730660/1‑140 (MQ=255)                            gTATCCGTCAGGGAAAGAAAGTTATCAATTGGGAGACGCATACTTTAACTGTTGATAATAATCAAATAATATTATTTCCTTGTGGTTATGAATTTTATATTGAGAATTACCCTGAAGCAGGGCTTTATCTTGCAGAAATg                                                                                                           <  2:2251162/140‑1 (MQ=255)                             tATCCGTCAGGGAAAGAAAGTTATCAATTGGGAGACGCATACTTTAACTGTTGATAATAATCAAATAATATTATTTCCTTGTGGTTATGAATTTTATATTGAGAATTACCCTGAAGCAGGGCTTTATCTTGCAGAAATGc                                                                                                          >  2:1692004/1‑140 (MQ=255)                                      gggAAAGAAAGTTATCAATTGGGAGACGCATACTTTAACTGTTGATAATAATCAAATAATATTATTTCCTTGTGGTTATGAATTTTATATTGAGAATTACCCTGAAGCAGGGCTTTATCTTGCAGAAATGCTttacttac                                                                                                 <  2:1977158/140‑1 (MQ=255)                                       ggAAAGAAAGTTATCAATTGGGAGACGCATACTTTAACTGTTGATAATAATCAAATAATATTATTTCCTTGTGGTTATGAATTTTATATTGAGAATTACCCTGAAGCAGGGCTTTATCTTGCAGAAATGCTTTACTTAcc                                                                                                >  2:1997095/1‑140 (MQ=255)                                                      aaTTGGGAGACGCATACTTTAACTGTTGATAATAATCAAATAATATTATTTCCTTGTGGTTATGAATTTTATATTGAGAATTACCCTGAAGCAGGGCTTTATCTTGCAGAAATGCTTTACTTACCCATTGATTTAATTga                                                                                 <  2:2253973/140‑1 (MQ=255)                                                                cGCATACTTTAACTGTTGATAATAATCAAATAATATTATTTCCTTGTGGTTATGAATTTTATATTGAGAATTACCCTGAAGCAGGGCTTTATCTTGCAGAAATGCTTTACTTACCCATTGATTTAATTGAGAGTTTCCaa                                                                       >  2:2244228/1‑140 (MQ=255)                                                                           aCTGTTGATAATAATCAAATAATATTATTTCCTTGTGGTTAGGAATTTTATATTGAGAATTACCCTGAAGCAGGGCTTTATCTTGCAGAAATGCTTTACTTACCCATTGATTTAATTGAGAGTTTCCAAAAACTTTATAc                                                            >  1:2095510/1‑140 (MQ=255)                                                                              gTTGATAATAATCAAATAATATTATTTCCTTGTGGTTATGAATTTTATATTGAGAATTACCCTGAAGCAGGGCTTTATCTTGCAGAAATGCTTTACTTACCCATTGATTTAATTGAGAGTTTCCAAAAACTTTATACGGt                                                         <  1:1997095/140‑1 (MQ=255)                                                                              gTTGATAATAATCAAATAATATTATTTCCTTGTGGTTATGAATTTTATATTGAGAATTACCCTGAAGCAGGGCTTTATCTTGCAGAAATGCTTTACTTACCCATTGATTTAATTGAGAGTTTCCAAAAACTTTATACGGt                                                         <  1:2228879/140‑1 (MQ=255)                                                                               ttGATAATAATCAAATAATATTATTTCCTTGTGGTTATGAATTTTATATTGAGAATTACCCTGAAGCAGGGCTTTATCTTGCAGAAATGCTTTACTTACCCATTGATTTAATTGAGAGTTTCCAAAAACTTTATACGGTa                                                        <  2:1829383/140‑1 (MQ=255)                                                                               ttGATAATAATCAAATAATATTATTTCCTTGTGGTTATGAATTTTATATTGAGAATTACCCTGAAGCAGGGCTTTATCTTGCAGAAATGCTTTACTTACCCATTGATTTAATTGAGAGTTTCCAAAAACTTTATACGGTa                                                        >  2:1746951/1‑140 (MQ=255)                                                                                 gATAATAATCAAATAATATTATTTCCTTGTGGTTATGAATTTTATATTGAGAATTACCCTGAAGCAGGGCTTTATCTTGCAGAAATGCTTTACTTACCCATTGATTTAATTGAGAGTTTCCAAAAACTTTATACGGTAAc                                                      >  1:1559787/1‑140 (MQ=255)                                                                                 gATAATAATCAAATAATATTATTTCCTTGTGGTTATGAATTTTATATTGAGAATTACCCTGAAGCAGGGCTTTATCTTGCAGAAATGCTTTACTTACCCATTGATTTAATTGAGAGTTTCCAAAAACTTTATACGGTAAc                                                      >  2:1786105/1‑140 (MQ=255)                                                                                         tCAAATAATATTATTTCCTTGTGGTTATGAATTTTATATTGAGAATTACCCTGAAGCAGGGCTTTATCTTGCAGAAATGCTTTACTTACCCATTGATTTAATTGAGAGTTTCCAAAAACTTTATACGGTAACTGATCaaa                                              >  1:1991021/1‑140 (MQ=255)                                                                                           aaaTAATATTATTTCCTTGTGGTTATGAATTTTATATTGAGAATTACCCTGAAGCAGGGCTTTATCTTGCAGAAATGCTTTACTTACCCATTGATTTAATTGAGAGTTTCCAAAAACTTTATACGGTAACTGATCAAATa                                            <  1:1714267/140‑1 (MQ=255)                                                                                                   ttattTCCTTGTGGTTATGAATTTTATATTGAGAATTACCCTGAAGCAGGGCTTTATCTTGCAGAAATGCTTTACTTACCCATTGATTTAATTGAGAGTTTCCAAAAACTTTATACGGTAACTGATCAAATACGTaacaa                                    >  2:1838386/1‑140 (MQ=255)                                                                                                       ttCCTTGTGGTTATGAATTTTATATTGAGAATTACCCTGAAGCAGGGCTTTATCTTGCAGAAATGCTTTACTTACCCATTGATTTAATTGAGAGTTTCCAAAAACTTTATACGGTAACTGATCAAATACGTaacaaaaca                                >  1:1781695/1‑140 (MQ=255)                                                                                                             gtgGTTATGAATTTTATATTGAGAATTACCCTGAAGCAGGGCTTTATCTTGCAGAAATGCTTTACTTACCCATTGATTTAATTGAGAGTTTCCAAAAACTTTATACGGTAACTGATCAAATACGTAACAAAACAAGTTTc                          >  1:1949517/1‑140 (MQ=255)                                                                                                                   aTGAATTTTATATTGAGAATTACCCTGAAGCAGGGCTTTATCTTGCAGAAATGCTTTACTTACCCATTGATTTAATTGAGAGTTTCCAAAAACTTTATACGGTAACTGATCAAATACGTAACAAAACAAGTTTCTTTTTa                    >  1:1886489/1‑140 (MQ=255)                                                                                                                        ttttATATTGAGAATTACCCTGAAGCAGGGCTTTATCTTGCAGAAATGCTTTACTTACCCATTGATTTAATTGAGAGTTTCCAAAAACTTTATACGGTAACTGATCAAATACGTAACAAAACAAGTTTCTTTTTACCTCa               <  1:2326835/140‑1 (MQ=255)                                                                                                                          ttATATTGAGAATTACCCTGAAGCAGGGCTTTATCTTGCAGAAATGCTTTACTTACCCATTGATTTAATTGAGAGTTTCCAAAAACTTTATACGGTAACTGATCAAATACGTAACAAAACAAGTTTCTTTTTACCTCAGa             >  2:2027399/1‑140 (MQ=255)                                                                                                                             tatTGAGAATTACCCTGAAGCAGGGCTTTATCTTGCAGAAATGCTTTACTTACCCATTGATTTAATTGAGAGTTTCCAAAAACTTTATACGGTAACTGATCAAATACGTAACAAAACAAGTTTCTTTTTACCTCAGAATc          <  1:1746951/140‑1 (MQ=255)                                                                                                                              atTGAGAATTACCCTGAAGCAGGGCTTTATCTTGCAGAAATGCTTTACTTACCCATTGATTTAATTGAGAGTTTCCAAAAACTTTATACGGTAACTGATCAAATACGTAACAAAACAAGTTTCTTTTTACCTCAGAATcc         <  1:2244228/140‑1 (MQ=255)                                                                                                                                    aaTTACCCTGAAGCAGGGCTTTATCTTGCAGAAATGCTTTACTTACCCATTGATTTAATTGAGAGTTTCCAAAAACTTTATACGGTAACTGATCAAATACGTAACAAAACAAGTTTCTTTTTa                    <  2:2156262/123‑1 (MQ=255)                                                                                                                                    aaTTACCCTGAAGCAGGGCTTTATCTTGCAGAAATGCTTTACTTACCCATTGATTTAATTGAGAGTTTCCAAAAACTTTATACGGTAACTGATCAAATACGTAACAAAACAAGTTTCTTTTTa                    >  1:2156262/1‑123 (MQ=255)                                                                                                                                     aTTACCCTGAAGCAGGGCTTTATCTTGCAGAAATGCTTTACTTACCCATTGATTTAATTGAGAGTTTCCAAAAACTTTATACGGTAACTGATCAAATACGTAACAAAACAAGTTTCTTTTTACCTCAGAATCCTGAGTTa  <  2:1897587/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                        TTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGTCAGGGAAAGAAAGTTATCAATTGGGAGACGCATACTTTAACTGTTGATAATAATCAAATAATATTATTTCCTTGTGGTTATGAATTTTATATTGAGAATTATCCTGAAGCAGGGCTTTATCTTGCAGAAATGCTTTACTTACCCATTGATTTAATTGAGAGTTTCCAAAAACTTTATACGGTAACTGATCAAATACGTAACAAAACAAGTTTCTTTTTACCTCAGAATCCTGAGTTA  >  NC_000913/571527‑571798

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