Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,776,420	T→C	V12A (GTA→GCA)	yfjY →	CP4‑57 prophage; RadC‑like JAB domain‑containing protein YfjY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,776,420	0	T	C	100.0%	142.9 / NA	40	V12A (GTA→GCA)	yfjY	CP4‑57 prophage; RadC‑like JAB domain‑containing protein YfjY
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/15);  total (25/15)

TGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGG  >  NC_000913/2776298‑2776554                                                                                                                           |                                                                                                                                       tGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGCACTTCCACTGACCGcac                                                                                                                       >  1:505391/1‑140 (MQ=255)     ggAGCACTATTACCGCCCGCGGAATTACGCTCTGTTGCAACCTGAATCTAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGCACTTCCACTGACCGCACAACg                                                                                                                   <  2:206155/140‑1 (MQ=255)        gCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGCACTTCCACTGACCGCACAACGCAc                                                                                                                >  2:689450/1‑140 (MQ=255)                 ccGCCTGCGGGATTACGCTCTGTAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGCACTTCCAGTGGCCGCACTACGCACGGTAAAAcc                                                                                                       >  1:206155/1‑139 (MQ=255)                 ccGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAAcg                                                                                                       <  1:460853/140‑1 (MQ=255)                        cGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGCACTTCCACTGACCGTACTACGCACTGTAAAACGCGCCTTa                                                                                                >  1:347014/1‑140 (MQ=255)                        cGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTa                                                                                                >  1:156634/1‑140 (MQ=255)                        cGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTa                                                                                                >  1:524140/1‑140 (MQ=255)                        cGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTa                                                                                                >  2:405866/1‑140 (MQ=255)                         gggATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTaa                                                                                               <  2:325976/140‑1 (MQ=255)                                cGCTCTGCCGCACCCTGAATCCCGCGCCATTATTCTCATTATTGATTATGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACgctgct                                                                                        <  2:347014/140‑1 (MQ=255)                                cGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACgctgct                                                                                        <  2:505391/140‑1 (MQ=255)                                cGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACgctgct                                                                                        <  2:156634/140‑1 (MQ=255)                                cGCTCTGCAGCACCCCGAATCCAGCCCCAATATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACgctgct                                                                                        <  1:242961/140‑1 (MQ=255)                                     tGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGAcc                                                                                   >  2:601972/1‑140 (MQ=255)                                            cccTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGa                                                                                 >  2:642579/1‑135 (MQ=255)                                            cccTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGa                                                                                 <  1:642579/135‑1 (MQ=255)                                             ccTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGCACTTCCACTGACTGCACAACGAACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTg                                                                           >  2:287625/1‑140 (MQ=255)                                                    ccAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGaaa                                                                    <  1:405866/140‑1 (MQ=255)                                                      aGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACa                                                                  >  1:400848/1‑140 (MQ=255)                                                          ccATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCg                                                              <  1:554456/140‑1 (MQ=255)                                                             ttATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTgg                                                           <  1:689450/140‑1 (MQ=255)                                                                    cATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGAGCGAAAGAGGCGTGGCATTCAc                                                    >  1:240833/1‑140 (MQ=255)                                                                          tGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCAcc                                                   <  2:14516/135‑1 (MQ=255)                                                                          tGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCAcc                                                   >  1:14516/1‑135 (MQ=255)                                                                                   aTACTCATGATGGAACAGTCACTCATCCCACAGACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGgctgc                                     <  1:449101/140‑1 (MQ=255)                                                                                   aTACTCATGATGGAACAGTCACTCATCCCACAGACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGgctgc                                     >  1:125350/1‑140 (MQ=255)                                                                                         atgatgGAACAGTCACTCATCCCACAGACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGa                               >  1:611422/1‑140 (MQ=255)                                                                                         atgatgGAACAGTCACTCATCCCACAGACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGa                               >  1:147249/1‑140 (MQ=255)                                                                                            atgGAACAGTCACTCATCCCACAGACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTg                            >  2:424038/1‑140 (MQ=255)                                                                                                    gTCACTCATCCCACAGACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGgctgaagc                    >  1:424678/1‑140 (MQ=255)                                                                                                     tCACTCATCCCACAGACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGgttgaagct                   >  2:596818/1‑140 (MQ=255)                                                                                                     tCACTCATCCCACAGACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGgctgaagct                   >  1:479170/1‑140 (MQ=255)                                                                                                       aCTCATCCCACAGACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGgctgaagctga                 >  1:534213/1‑140 (MQ=255)                                                                                                          cATCCCACAGACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGaaaa              >  2:508124/1‑140 (MQ=255)                                                                                                          cATCCCACAGACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGaaaa              >  2:97773/1‑140 (MQ=255)                                                                                                             cccACAGACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATgg           <  1:508124/140‑1 (MQ=255)                                                                                                                   gACACCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCggggg     >  2:394716/1‑139 (MQ=255)                                                                                                                    acacCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCt    <  1:45115/140‑1 (MQ=255)                                                                                                                      acCGGCACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTgg  >  2:518325/1‑140 (MQ=255)                                                                                                                           |                                                                                                                                       TGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGGCGGGGCTGG  >  NC_000913/2776298‑2776554

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