Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,771,295	A→G	T50A (ACC→GCC)	yfjT →	CP4‑57 prophage; uncharacterized protein YfjT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,771,295	0	A	G	100.0%	108.3 / NA	31	T50A (ACC→GCC)	yfjT	CP4‑57 prophage; uncharacterized protein YfjT
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/20);  total (11/20)

TCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGCATCACAGTGGTGATGGACCATA  >  NC_000913/2771160‑2771424                                                                                                                                        |                                                                                                                                  tCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTGCCAt                                                                                                                               <  2:100437/140‑1 (MQ=255)         ccGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTGCCATCAACTgga                                                                                                                       >  2:368999/1‑138 (MQ=255)             tttGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTGCCATCAACTGCACggg                                                                                                                   <  1:102255/140‑1 (MQ=255)                      gCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTGCCATCAACTGCACGGGAAATCCGGa                                                                                                          >  1:216438/1‑140 (MQ=255)                            aCAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTGCCATCAACTGCACGGGAAATCCGGACCAGaa                                                                                                    <  1:124278/140‑1 (MQ=255)                            aCAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTGCCATCAACTGCACGGGAAATCCGGACCAGaa                                                                                                    <  2:17345/140‑1 (MQ=255)                                tGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTGCCATCAACTGCACGGGAAATCCGGACCAGAACGGt                                                                                                >  2:293292/1‑140 (MQ=255)                                tGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTGCCATCAACTGCACGGGAAATCCGGACCAGAACGGt                                                                                                >  2:242724/1‑140 (MQ=255)                                       cAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTGCCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTAcc                                                                                         <  1:315506/140‑1 (MQ=255)                                             tACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTGCCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACt                                                                                   >  2:150381/1‑140 (MQ=255)                                                     ctGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTGCCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGttttt                                                                           >  1:501971/1‑140 (MQ=255)                                                                                      gTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTGCCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCt                                            <  2:326748/138‑1 (MQ=255)                                                                                      gTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTGCCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCt                                            >  1:326748/1‑138 (MQ=255)                                                                                           ggCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTGCCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCatga                                     >  2:520233/1‑140 (MQ=255)                                                                                            gCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTGCCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCatgat                                    <  2:115001/140‑1 (MQ=255)                                                                                            gCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTGCCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCatgat                                    <  2:482864/140‑1 (MQ=255)                                                                                            gCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTGCCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCatgat                                    <  1:242724/140‑1 (MQ=255)                                                                                            gCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTGCCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCatgat                                    <  2:231700/140‑1 (MQ=255)                                                                                            gCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTGCCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCatgat                                    <  1:25684/140‑1 (MQ=255)                                                                                                     aaTACAGTGTGACGGATGCCAGCGGAAAGATGTTTGCCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTagc                           >  2:503522/1‑138 (MQ=255)                                                                                                       tACAGTGTGACGGATGCCAGCGGAAAGATGTTTGCCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCgg                         <  1:293292/140‑1 (MQ=255)                                                                                                       tACAGTGTGACGGATGCCAGCGGAAAGATGTTTGCCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCgg                         <  1:410034/140‑1 (MQ=255)                                                                                                          aGTGTGACGGATGCCAGCGGAAAGATGTTTGCCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGCAt                      <  1:368999/140‑1 (MQ=255)                                                                                                           gtgtGACGGATGCCAGCGGAAAGATGTTTGCCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGCATc                     <  1:520233/140‑1 (MQ=255)                                                                                                            tgtgACGGATGCCAGCGGAAAGATGTTTGCCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGCATca                    <  2:216438/140‑1 (MQ=255)                                                                                                            tgtgACGGATGCCAGCGGAAAGATGTTTGCCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGCATca                    >  2:441637/1‑140 (MQ=255)                                                                                                             gtgACGGATGCCAGCGGAAAGATGTTTGCCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGCATcac                   <  1:124881/140‑1 (MQ=255)                                                                                                              tgACGGATGCCAGCGGAAAGATGTTTGCCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTAGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCGGCTCGCATGATGACAGTACCGGCATcaca                  >  2:292564/1‑140 (MQ=255)                                                                                                                      gCCAGCGGAAAGATGTTTGCCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGCATCACAGTGGTGAt          <  1:96607/140‑1 (MQ=255)                                                                                                                             gAAAGATGTTTGCCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGCATCACAGTGGTGATGGACCAt   <  2:501971/140‑1 (MQ=255)                                                                                                                              aaaGATGTTTGCCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGCATCACAGTGGTGATGGACCATa  <  2:438818/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                  TCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGCATCACAGTGGTGATGGACCATA  >  NC_000913/2771160‑2771424

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