Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,192,484	A→G	S50P (TCG→CCG)	hflD ←	lysogenization regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,192,484	0	A	G	100.0%	88.8 / NA	26	S50P (TCG→CCG)	hflD	lysogenization regulator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/16);  total (10/16)

AGCACCATCAAGCTGAGTGTGTAGCGGGTTAATTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGG  >  NC_000913/1192348‑1192606                                                                                                                                         |                                                                                                                           aGCACCATCAAGCTGAGTGTGTAGCGGGTTAATTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGGGCt                                                                                                                         >  1:650975/1‑140 (MQ=255)     ccATCAAGCTGAGTGTGTAGCGGGTTAATTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGGGCTgggg                                                                                                                     <  1:14283/140‑1 (MQ=255)          aaGCTGAGTGTGTAGCGGGTTAATTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGGGCTGGGGTTCat                                                                                                                >  2:230182/1‑140 (MQ=255)            gCTGAGTGTGTAGCGGGTTAATTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGGGCTGGGGTTCatat                                                                                                              <  1:443246/140‑1 (MQ=255)                  tgtgtAGCGGGTTAATTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGGGCTGGGGTTCATATCaataa                                                                                                        >  1:298021/1‑140 (MQ=255)                      tAGCGGGTTAATTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGGGCTGGGGTTCATATCAATAATACt                                                                                                    <  1:230182/140‑1 (MQ=255)                            gTTAATTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGGGCTGGGGTTCATATCAATAATACTGTtgag                                                                                              >  1:159804/1‑140 (MQ=255)                            gTTAATTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGGGCTGGGGTTCATATCAATAATACTGTtgag                                                                                              >  2:889116/1‑140 (MQ=255)                            gTTAATTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGGGCTGGGGTTCATATCAATAATACTGTtgag                                                                                              >  1:287096/1‑140 (MQ=255)                             ttAATTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCCCGCGGAGCAGGGTTCCCAGCCCGACGCGCAGTTTGGTTTCGCTACCGCCAAAAACCGCCAGCGTCGGGCTGGGGTTCATATCAATAATACTGTtgagt                                                                                             <  2:476401/140‑1 (MQ=255)                             ttAATTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGGGCTGGGGTTCATATCAATAATACTGTtgagt                                                                                             >  1:107489/1‑140 (MQ=255)                                aTTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGGGCTGGGGTTCATATCAATAATACTGTtgagtgag                                                                                          <  1:889116/140‑1 (MQ=255)                                      cGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGGGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACgtgt                                                                                    <  2:298021/140‑1 (MQ=255)                                        ttaaggCCCTGGCGCCGGCTGCCATTGGCCCCGGCGTGTAGGGTTTCCAGCCGGGCGCGCGGGTTGGGTTCGCTCCCCCCAAAAACCGCCAGCGTCGGGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAg                                                                                  <  2:892599/135‑1 (MQ=255)                                           aaGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGGGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAgcgc                                                                               <  2:107489/140‑1 (MQ=255)                                                  ggCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGGGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCa                                                                        >  1:23606/1‑140 (MQ=255)                                                       ctgctgGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGGGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATcaca                                                                   >  1:103415/1‑140 (MQ=255)                                                                       cGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGGGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCg                                                   <  2:287096/140‑1 (MQ=255)                                                                         ccGGGCAGGGTTTCCAGCCAGTCGCGCGGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGGGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAg                                                 <  2:650975/140‑1 (MQ=255)                                                                          cGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGGGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAgt                                                >  2:177928/1‑140 (MQ=255)                                                                            aGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGGGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAgttg                                              <  2:39138/140‑1 (MQ=255)                                                                                              aCGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGGGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGa                            <  2:775604/140‑1 (MQ=255)                                                                                                   cAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGGGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGAc                       <  1:49811/140‑1 (MQ=255)                                                                                                             tCGCTACCGCCAAAAACCGCCAGCGTCGGGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGc             <  1:577788/140‑1 (MQ=255)                                                                                                             tCGCTACCGCCAAAAACCGCCAGCGTCGGGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGc             <  2:325717/140‑1 (MQ=255)                                                                                                                        aaaaaCCGCCAGCGTCGGGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAggg  <  2:23606/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                           AGCACCATCAAGCTGAGTGTGTAGCGGGTTAATTCGGCGTTTAAGCCCTGGCGACTGCTGGCATTGAGCACGCCGAGCAGGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGG  >  NC_000913/1192348‑1192606

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