Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,422,408	T→C	S142S (AGT→AGC)	ydaV →	Rac prophage; putative ATP‑binding protein YdaV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,422,408	0	T	C	100.0%	104.6 / NA	31	S142S (AGT→AGC)	ydaV	Rac prophage; putative ATP‑binding protein YdaV
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/15);  total (16/15)

CGGGTGTACAAATTTTGCGTTCAGCGGAAAACCTGGTACCGGGAAGAACCACTTAGCGGCAGCTATCGGGAATCGCCTGCTGAAAGACGGTCAGACAGTGATTGTGGTTACCGTGGCTGATGTTATGAGTGCCCTGCACGCCAGCTATGACGATGGGCAGTCAGGCGAAAAATTTTTGCGGGAACTGTGCGAAGTGGATCTGCTGGTTCTTGATGAAATTGGCATTCAGCGCGAGACGAAAAACGAGCAGG  >  NC_000913/1422279‑1422529                                                                                                                                  |                                                                                                                          cGGGTGTACAAATTTTGCGTTCAGCGGAAAACCTGGTACCGGGAAGAACCACTTAGCGGCAGCTATCGGGAATCGCCTGCTGCAAGACGGTCAGACAGTGATTGTGGTTACCGTGGCTGATGTTATGAGCGcc                                                                                                                        >  2:13849/1‑133 (MQ=255)   ggTGTACAAATTTTGCGTTCAGCGGAAAACCTGGTACCGGGAAGAACCACTTAGCGGCAGCTATCGGGAATCGCCTGCTGAAAGACGGTCAGACAGTGATTGTGGTTACCGTGGCTGATGTTATGAGCGCCCt                                                                                                                      >  2:357724/1‑133 (MQ=255)   ggTGTACAAATTTTGCGTTCAGCGGAAAACCTGGTACCGGGAAGAACCACTTAGCGGCAGCTATCGGGAATCGCCTGCTGAAAGACGGTCAGACAGTGATTGTGGTTACCGTGGCTGATGTTATGAGCGCCCt                                                                                                                      >  1:33970/1‑133 (MQ=255)       tACAAATTTTGCGTTCAGCGGAAAACCTGGTACCGGGAAGAACCACTTAGCGGCAGCTATCGGGAATCGCCTGCTGAAAGACGGTCAGACAGTGATTGTGGTTACCGTGGCTGATGTTATGAGCGCCCTGCAc                                                                                                                  <  2:365702/133‑1 (MQ=255)        aCAAATTTTGCGTTCAGCGGAAAACCTGGTACCGGGAAGAACCACTTAGCGGCAGCTATCGGGAATCGCCTGCTGAAAGACGGTCAGACAGTGATTGTGGTTACCGTGGCTGATGTTATGAGCGCCCTGCACg                                                                                                                 <  2:33970/133‑1 (MQ=255)          aaaTTTTGCGTTCAGCGGAAAACCTGGTACCGGGAAGAACCACTTAGCGGCAGCTATCGGGAATCGCCTGCTGAAAGACGGTCAGACAGTGATTGTGGTTACCGTGGCTGATGTTATGAGCGCCCTGCACGcc                                                                                                               >  2:328859/1‑133 (MQ=255)              tttGCGTTCAGCGGAAAACCTGGTACCGGGAAGAACCACTTAGCGGCAGCTATCGGGAATCGCCTGCTGAAAGACGGTCAGACAGTGATTGTGGTTACCGTGGCTGATGTTATGAGCGCCCTGCACGCCAGCt                                                                                                           >  2:361804/1‑133 (MQ=255)                              aaCCTGGTACCGGGAAGAACCACTTAGCGGCAGCTATCGGGAATCGCCTGCTGAAAGACGGTCAGACAGTGATTGTGGTTACCGTGGCTGATGTTATGAGCGCCCTGCACGCCAGCTAAGACGATGGGCAGTc                                                                                           >  1:332741/1‑133 (MQ=255)                                          ggAAGAACCACTTAGCGGCAGCTATCGGGAATCGCCTGCTGAAAGACGGTCAGACAGTGATTGTGGTTACCGTGGCTGATGTTATGAGCGCCCTGCACGCCAGCTATGACGATGGGCGGTCAGGAGAAAAAtt                                                                               >  1:134075/1‑133 (MQ=255)                                               aaCCACTTAGCGGCAGCTATCGGGAATCGCCTGCTGAAAGACGGTCAGACAGTGATTGTGGTTACCGTGGCTGATGTTATGAGCGCCCTGCACGCCAGCTATGACGATGGGCAGTCAGGCGAAAAATTTTTGc                                                                          >  2:279198/1‑133 (MQ=255)                                                  cACTTAGCGGCAGCTATCGGGAATCGCCTGCTGAAAGACGGTCAGACAGTGATTGTGGTTACCGTGGCTGATGTTATGAGCGCCCTGCACGCCAGCTATGACGATGGGCAGTCAGGCGAAAAATTTTTGCggg                                                                       >  1:67419/1‑133 (MQ=255)                                                    cTTAGCGGCAGCTATCGGGAATCGCCTGCTGAAAGACGGTCAGACAGTGATTGTGGTTACCGTGGCTGATGTTATGAGCGCCCTGCACGCCAGCTATAACGATGGGCAGTCAGGCGAAAAATTTTTGCGGGaa                                                                     <  1:361804/133‑1 (MQ=255)                                                       aGCGGCAGCTATCGGGAATCGCCTGCTGAAAGACGGTCAGACAGTGATTGTGGTTACCGTGGCTGATGTTATGAGCGCCCTGCACGCCAGCTATGACGATGGGCAGTCAGGCGAAAAATTTTTGCGGGAACtg                                                                  >  2:290704/1‑133 (MQ=255)                                                         cGGCAGCTATCGGGAATCGCCTGCTGAAAGACGGTCAGACAGTGATTGTGGTTACCGTGGCTGATGTTATGAGCGCCCTGCACGCCAGCTATGACGATGGGCAGTCAGGCGAAAAATTTTTGCGGGAACtgtg                                                                <  1:290704/133‑1 (MQ=255)                                                         cGGCAGCTATCGGGAATCGCCTGCTGAAAGACGGTCAGACAGTGATTGTGGTTACCGTGGCTGATGTTATGAGCGCCCTGCACGCCAGCTATGACGATGGGCAGTCAGGCGAAAAATTTTTGCGGGAACtgtg                                                                <  1:328859/133‑1 (MQ=255)                                                            cAGCTATCGGGAATCGCCTGCTGAAAGACGGTCAGACAGTGATTGTGGTTACCGTGGCTGATGTTATGAGCGCCCTGCACGCCAGCTATGACGATGGGCAGTCAGGCGAAAAATTTTTGCGGGAACTGTGCGa                                                             <  1:82270/133‑1 (MQ=255)                                                            cAGCTATCGGGAATCGCCTGCTGAAAGACGGTCAGACAGTGATTGTGGTTACCGTGGCTGATGTTATGAGCGCCCTGCACGCCAGCTATGACCATGGGCAGTCAGGCGAAAAATTTTTGCGGGAACTGTGCGa                                                             <  1:13849/133‑1 (MQ=255)                                                               cTATCGGGAATCGCCTGCTGAAAGACGGTCAGACAGTGATTGTGGTTACCGTGGCTGATGTTATGAGCGCCCTGCACGCCAGCTATGACGATGGGCAGTCAGGCGAAAAATTTTTGCGGGAACTGTGCGAAGt                                                          >  2:37096/1‑133 (MQ=255)                                                               cTATCGGGAATCGCCTGCTGAAAGACGGTCAGACAGTGATTGTGGTTACCGTGGCTGATGTTATGAGCGCCCTGCACGCCAGCTAAGACGATGGGCAGTCAGGCGAAAAATTTTTGCGGGAACTGTGCGAAGt                                                          <  2:332741/133‑1 (MQ=255)                                                                      gAATCGCCTGCTGAAAGACGGTCAGACAGTGATTGTGGTTACCGTGGCTGATGTTATGAGCGCCCTGCACGCCAGCTATGACGATGGGCAGTCAGGCGAAAAATTTTTGCGGGAACTGTGCGAAGTGGATctg                                                   >  2:489390/1‑133 (MQ=255)                                                                            ccTGCTGAAAGACGGTCAGACAGTGATTGTGGTTACCGTGGCTGATGTTATGAGCGCCCTGCACGCCAGCTATGACGATGGGCAGTCAGGCTAAAAATTTGTGCAGGAACTGTGCGAGGTGGCTCTGCTGGtt                                             >  2:152123/1‑133 (MQ=255)                                                                                ctgAAAGACGGTCAGACAGTGATTGTGGTTACCGTGGCTGATGTTATGAGCGCCCTGCACGCCAGCTATGACGATGGGCAGTCAGGCGAAAAATTTTTGCGGGAACTGTGCGAAGTGGATCTGCTGGTTCTtg                                         >  2:491244/1‑133 (MQ=255)                                                                                         ggTCAGACAGTGATTGTGGTTACCGTGGCTGATGTTATGAGCGCCCTGCACGCCAGCTATGACGATGGGCAGTCAGGCGAAAAATTTTTGCGGGAACTGTGCGAAGTGGATCTGCTGGTTCTTGATGAAATTg                                <  1:279198/133‑1 (MQ=255)                                                                                                     aTTGTGGTTACCGTGGCTGATGTTATGAGCGCCCTGCACGCCAGCTATGACGATGGGCAGTCAGGCGAAAAATTTTTGCGGGAACTGTGCGAAGTGGATCTGCTGGTTCTTGATGAAATTGGCATTCAGcgcg                    >  1:108184/1‑133 (MQ=255)                                                                                                           gTTACCGTGGCTGATGTTATGAGCGCCCTGCACGCCAGCTATGACGATGGGCAGTCAGGCGAAAAATTTTTGCGGGAACTGTGCGAAGTGGATCTGCTGGTTCTTGATGAAATTGGCATTCAGCGCGAGACGa              >  1:106961/1‑133 (MQ=255)                                                                                                             tACCGTGGCTGATGTTATGAGCGCCCTGCACGCCAGCTATGACGATGGGCAGTCAGGCGAAAAATTTTTGCGGGAACTGTGCGAAGTGGATCTGCTGGTTCTTGATGAAATTGGCATTCAGCGCGAGACGaaa            >  1:34462/1‑133 (MQ=255)                                                                                                                  tGGCTGATGTTATGAGCGCCCTGCACGCCAGCTATGACGATGGGCAGTCAGGCGAAAAATTTTTGCGGGAACTGTGCGAAGTGGATCTGCTGGTTCTTGATGAAATTGGCATTCAGCGCGAGACGAAAAACGa       <  1:357724/133‑1 (MQ=255)                                                                                                                  tGGCTGATGTTATGAGCGCCCTGCACGCCAGCTATGACGATGGGCAGTCAGGCGAAAAATTTTTGCGGGAACTGTGCGAAGTGGATCTGCTGGTTCTTGATGAAATTGGCATTCAGCGCGAGACGAAAAACGa       <  2:67419/133‑1 (MQ=255)                                                                                                                  tGGCTGATCTTATGACCGCCCTCCACGCCACCTATGACGATGGGCAGTCAGGCGAAAAATTTTTGCGGGAACTGTGCGAAGTGGATCTGCTGGTTCTTGATGAAATTGGCATTCAGCGCGAGACGAAAAACGa       <  1:145812/133‑1 (MQ=255)                                                                                                                  tGGCTGATCTTATGAACGTCCTACACGCCACCTATGACGATGGGCAGTCAGGCGAAAAATTTTTGCGGGAACTGTGCGAAGTGGATCTGCTGGTTCTTGATGAAATTGGCATTCAGCGCGAGACGAAAAACGa       <  1:152123/133‑1 (MQ=255)                                                                                                                     cTGATGTTATGAGCGCCCTGCACGCCAGCTATGACGATGGGCAGTCAGGCGAAAAATTTTTGCGGGAACTGTGCGAAGTGGATCTGCTGGTTCTTGATGAAATTGGCATTCAGCGCGAGACGAAAAACGAGCa    <  1:491244/133‑1 (MQ=255)                                                                                                                       gATGTTATGAGCGCCCTGCACGCCAGCTATGACGATGGGCAGTCAGGCGAAAAATTTTTGCGGGAACTGTGCGAAGTGGATCTGCTGGTTCTTGATGAAATTGGCATTCAGCGCGAGACGAAAAACGAGCAgg  <  1:37096/133‑1 (MQ=255)                                                                                                                                  |                                                                                                                          CGGGTGTACAAATTTTGCGTTCAGCGGAAAACCTGGTACCGGGAAGAACCACTTAGCGGCAGCTATCGGGAATCGCCTGCTGAAAGACGGTCAGACAGTGATTGTGGTTACCGTGGCTGATGTTATGAGTGCCCTGCACGCCAGCTATGACGATGGGCAGTCAGGCGAAAAATTTTTGCGGGAACTGTGCGAAGTGGATCTGCTGGTTCTTGATGAAATTGGCATTCAGCGCGAGACGAAAAACGAGCAGG  >  NC_000913/1422279‑1422529

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