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, conserved protein

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%	182.1 / NA	51	D171D (GAT→GAC)	ydaU	Rac prophage, conserved protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (26/25);  total (26/25)

AAGCAGAAAGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGAT‑CCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAA  >  NC_000913/1421493‑1421758                                                                                                                                            |                                                                                                                                 aagCAGAAAGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGAt                                                                                                                                   <  1:1513387/139‑1 (MQ=255)                 gATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGAT‑c                                                                                                                                 >  2:1527885/1‑124 (MQ=255)                 gATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGAT‑c                                                                                                                                 <  1:1527885/124‑1 (MQ=255)                      tCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCg                                                                                                             >  2:966944/1‑139 (MQ=255)                                    agagaAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGgtg                                                                                               >  2:906977/1‑139 (MQ=255)                                    agagaAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGgtg                                                                                               >  2:1164338/1‑139 (MQ=255)                                         aaTTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCaaa                                                                                          <  1:906977/139‑1 (MQ=255)                                                   gtgtTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGAT‑cc                                                                                                                                >  2:899601/1‑91 (MQ=255)                                                   gtgtTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGAT‑cc                                                                                                                                <  1:899601/91‑1 (MQ=255)                                                     gtTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct                                                                              <  2:246246/139‑1 (MQ=255)                                                     gtTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct                                                                              <  1:966944/139‑1 (MQ=255)                                                       tCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCctct                                                                            <  1:1164338/139‑1 (MQ=255)                                                               aTGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAgg                                                                                                 >  2:94920/1‑110 (MQ=255)                                                               aTGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAgg                                                                                                 <  1:94920/110‑1 (MQ=255)                                                                     aaGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct                                                                              <  1:180558/123‑1 (MQ=255)                                                                     aaGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct                                                                              >  2:180558/1‑123 (MQ=255)                                                                         ggtgTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTa                                                                    >  2:286458/1‑129 (MQ=255)                                                                         ggtgTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTa                                                                    <  1:286458/129‑1 (MQ=255)                                                                                     tgatgTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAAt                                                                                                                      <  2:698932/68‑1 (MQ=255)                                                                                     tgatgTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAAt                                                                                                                      >  1:698932/1‑68 (MQ=255)                                                                                     tgatgTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGAAATCGAGGTGTCAAAAAGTTTGACCct                                                                              >  2:870526/1‑107 (MQ=255)                                                                                     tgatgTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGAAATCGAGGTGTCAAAAAGTTTGACCct                                                                              <  1:870526/107‑1 (MQ=255)                                                                                      gatgTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGAT‑ccccc                                                                                                                             <  2:886183/59‑1 (MQ=255)                                                                                      gatgTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGAT‑ccccc                                                                                                                             >  1:886183/1‑59 (MQ=255)                                                                                                   aGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGa                                                                                  <  1:145851/89‑1 (MQ=255)                                                                                                   aGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGa                                                                                  >  2:145851/1‑89 (MQ=255)                                                                                                      cAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct                                                                              <  1:701135/90‑1 (MQ=255)                                                                                                      cAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct                                                                              >  2:701135/1‑90 (MQ=255)                                                                                                        aCTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCgg                                                                                                            <  1:1546211/58‑1 (MQ=255)                                                                                                        aCTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCgg                                                                                                            >  2:1546211/1‑58 (MQ=255)                                                                                                               aaGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct                                                                              <  1:838108/81‑1 (MQ=255)                                                                                                               aaGATACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct                                                                              >  2:838108/1‑81 (MQ=255)                                                                                                                   tACAGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCg                >  1:1088570/1‑139 (MQ=255)                                                                                                                    aCAGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGAc                                                                                 <  2:1201155/73‑1 (MQ=255)                                                                                                                    aCAGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGAc                                                                                 >  1:1201155/1‑73 (MQ=255)                                                                                                                     cAGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGc                                   >  2:1432641/1‑118 (MQ=255)                                                                                                                     cAGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGc                                   <  1:1432641/118‑1 (MQ=255)                                                                                                                     cAGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGcc              <  2:1088570/139‑1 (MQ=255)                                                                                                                      aGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTg                                                                                   >  2:1658717/1‑69 (MQ=255)                                                                                                                      aGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTg                                                                                   <  1:1658717/69‑1 (MQ=255)                                                                                                                      aGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct                                                                              >  2:569307/1‑74 (MQ=255)                                                                                                                      aGATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct                                                                              <  1:569307/74‑1 (MQ=255)                                                                                                                       gATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGcc                                                             <  2:147351/90‑1 (MQ=255)                                                                                                                       gATAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGcc                                                             >  1:147351/1‑90 (MQ=255)                                                                                                                         tAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCaaa                                                                                          <  2:210927/59‑1 (MQ=255)                                                                                                                         tAAAGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCaaa                                                                                          >  1:210927/1‑59 (MQ=255)                                                                                                                          aaaGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCa         >  2:629239/1‑139 (MQ=255)                                                                                                                            aGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCctc                                                                             >  2:1412175/1‑69 (MQ=255)                                                                                                                            aGATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCctc                                                                             <  1:1412175/69‑1 (MQ=255)                                                                                                                             gATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct                                                                              <  2:110604/67‑1 (MQ=255)                                                                                                                             gATCTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct                                                                              >  1:110604/1‑67 (MQ=255)                                                                                                                                cTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCctct                                                                            <  2:1210422/66‑1 (MQ=255)                                                                                                                                cTAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCccct                                                                            >  1:1210422/1‑66 (MQ=255)                                                                                                                                 tAAAAACAGATCCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGaa  >  1:1126660/1‑139 (MQ=255)                                                                                                                                  aaaaaCAGATCCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct                                                                              >  1:216621/1‑63 (MQ=255)                                                                                                                                  aaaaaCAGATCCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct                                                                              >  1:925916/1‑63 (MQ=255)                                                                                                                                  aaaaaCAGATCCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct                                                                              <  2:216621/63‑1 (MQ=255)                                                                                                                                  aaaaaCAGATCCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct                                                                              <  2:925916/63‑1 (MQ=255)                                                                                                                                            |                                                                                                                                 AAGCAGAAAGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGAT‑CCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAA  >  NC_000913/1421493‑1421758

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