Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,421,638	(C)7→8	coding (520/858 nt)	ydaU →	Rac prophage; DUF1376 domain‑containing protein YdaU

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,421,631	1	.	C	100.0%	114.5 / NA	32	D171D (GAT→GAC)	ydaU	Rac prophage; DUF1376 domain‑containing protein YdaU
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (20/12);  total (20/12)

GCAGAAAGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGAT‑‑CCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTG  >  NC_000913/1421495‑1421754                                                                                                                                          |                                                                                                                             gCAGAAAGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGAT‑‑ccc                                                                                                                          <  2:584811/140‑1 (MQ=255)     aaaGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGAT‑‑ccccccc                                                                                                                      >  1:139283/1‑140 (MQ=255)      aaGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGAT‑‑cccccccc                                                                                                                     >  1:885677/1‑139 (MQ=255)                   ttCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATC‑CCCCCCCTAAATCCCCCCCg                                                                                                         >  1:431552/1‑140 (MQ=255)                          aaaacaaGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATC‑CCCCCCCTAAATCCCCCCCGGGGGAAt                                                                                                  >  2:46300/1‑140 (MQ=255)                            aacaaGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATC‑CCCCCCCTAAAT‑CCCCCCGGGGGAATCGa                                                                                               >  1:302846/1‑140 (MQ=255)                              caaGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATC‑CCCCCCCTAAATCCCCCCCGGGGGAATCGAg                                                                                              <  1:274311/140‑1 (MQ=255)                                    agaAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATC‑CCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCa                                                                                        >  2:71827/1‑140 (MQ=255)                                             aCAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATC‑CCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTg                                                                               <  2:431552/140‑1 (MQ=255)                                             aCAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATC‑CCCCCCCTAAAT‑CCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGa                                                                              <  2:302846/140‑1 (MQ=255)                                              cAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATC‑CCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGa                                                                              <  2:90116/140‑1 (MQ=255)                                                ggTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATC‑CCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGAcc                                                                            >  2:307499/1‑140 (MQ=255)                                                  tgtTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATC‑CCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct                                                                          <  1:899477/140‑1 (MQ=255)                                                   gtTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct                                                                          <  2:978827/140‑1 (MQ=255)                                                            gATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATC‑CCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTa                                                                >  1:826284/1‑140 (MQ=255)                                                             aTGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATC‑CCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTAc                                                               >  1:687344/1‑140 (MQ=255)                                                                       ggtgTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATC‑CCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAAc                                                     >  1:231839/1‑140 (MQ=255)                                                                        gtgtTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATC‑CCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACt                                                    >  2:299060/1‑140 (MQ=255)                                                                            ttGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATC‑CCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGAt                                                <  2:343402/140‑1 (MQ=255)                                                                             tGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATC‑CCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGaa                                                      >  1:692158/1‑133 (MQ=255)                                                                             tGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATC‑CCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGaa                                                      <  2:692158/133‑1 (MQ=255)                                                                                    gatgTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATC‑CCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGtct                                        >  2:459742/1‑140 (MQ=255)                                                                                     atgTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATC‑CCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGtctc                                       >  2:957042/1‑140 (MQ=255)                                                                                     atgTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATC‑CCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGtctc                                       >  2:927141/1‑140 (MQ=255)                                                                                        tCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATC‑CCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCt                                    <  1:927141/140‑1 (MQ=255)                                                                                        tCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATC‑CCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCt                                    <  2:885677/140‑1 (MQ=255)                                                                                              caAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATC‑CCCCCCCTAAATCCCACCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCg                              >  1:486009/1‑140 (MQ=255)                                                                                                 aagCAACTAAAAAAGATAAAGATAAAAAGCTAAAATAAGATC‑CCCCCCCGAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCATAACTTGTTTTCTGTCTCGTTTTGGCGTGa                           <  2:139283/138‑1 (MQ=255)                                                                                                 aGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATC‑CCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGa                           >  1:590621/1‑140 (MQ=255)                                                                                                   gCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATC‑CCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGt                         <  2:826284/140‑1 (MQ=255)                                                                                                             aaGATACAGATAAAGATCTAAAAACAGATC‑CCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAAtt               >  2:663175/1‑140 (MQ=255)                                                                                                                  aCAGATAAAGATCTAAAAACAGATC‑CCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAAtt               >  1:535290/1‑135 (MQ=255)                                                                                                                  aCAGATAAAGATCTAAAAACAGATC‑CCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAAtt               <  2:535290/135‑1 (MQ=255)                                                                                                                  aCAGATAAAGATCTAAAAACAGATC‑CCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGcc          >  2:88895/1‑140 (MQ=255)                                                                                                                  aCAGATAAAGATCTAAAAACAGATC‑CCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGACGTGAGTGGGTTGAATTTCGcc          >  1:737056/1‑140 (MQ=255)                                                                                                                          aGATCTAAAAACAGATC‑CCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTg  >  2:393470/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                             GCAGAAAGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGAT‑‑CCCCCCCTAAATCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTG  >  NC_000913/1421495‑1421754

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