Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	293,744	T→C	S59G (AGC→GGC)	yagL ←	CP4‑6 prophage; resolvase‑like catalytic domain‑containing protein YagL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	293,744	0	T	C	100.0%	142.7 / NA	40	S59G (AGC→GGC)	yagL	CP4‑6 prophage; resolvase‑like catalytic domain‑containing protein YagL
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/18);  total (22/18)

TGATGAATCAATTGTCTGAGTAAACGTTGCTCTGAAACAACCTCTTTAATATCTACCTGCTCGATCACAACATGTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCTTTCTTGCTGTTCGCTACTGCATATGGTGTAAAAAAATGAACGCATCTGATGCACTCCGGTGTCATAAATTATTTTTTGATACCTGAATAATGTCATTTTGCATAAAAAACAACCAAAAATAGAGAGT  >  NC_000913/293611‑293871                                                                                                                                      |                                                                                                                                tgatgaATCAATTGTCTGAGTAAACGTTGCTCTGAAACAACCTCTTTAATATCTACCTGCTCGATCACAACATGTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTg                                                                                                                           >  2:710906/1‑140 (MQ=255)      aaTCAATTGTCTGAGTAAACGTTGCTCTGAAACAACCTCTTTAATATCTACCTGCTCGATCACAACATGTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGtt                                                                                                                      <  2:482529/140‑1 (MQ=255)        tCAATTGTCTGAGTAAACGTTGCTCTGAAACAACCTCTTTAATATCTACCTGCTCGATCACAACATGTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCg                                                                                                                    <  2:67379/140‑1 (MQ=255)        tCAATTGTCTGAGTAAACGTTGCTCTGAAACAACCTCTTTAATATCTACCTGCTCGATCACAACATGTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCg                                                                                                                    <  1:557777/140‑1 (MQ=255)         cAATTGTCTGAGTAAACGTTGCTCTGAAACAACCTCTTTAATATCTACCTGCTCGATCACAACATGTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCGc                                                                                                                   <  1:371722/140‑1 (MQ=255)                cTGAGTAAACGTTGCTCTGAAGCAACCTCTTTAATATCTACCTGCTCGATCACAACATGTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCGCTa                                                                                                                 >  2:833630/1‑135 (MQ=255)                cTGAGTAAACGTTGCTCTGAAGCAACCTCTTTAATATCTACCTGCTCGATCACAACATGTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCGCTa                                                                                                                 <  1:833630/135‑1 (MQ=255)                 tGAGTAAACGTTGCTCTGAAACAACCTCTTTAATATCTACCTGCTCGATCACAACATGTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCGCTACTGCat                                                                                                           >  2:169396/1‑140 (MQ=255)                   aGTAAACGTTGCTCTGAAACAACCTCTTTAATATCTACCTGCTCGATCACAACATGTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCGCTACTGCatat                                                                                                         >  1:484715/1‑140 (MQ=255)                        aCGTTGCTCTGAAACAACCTCTTTAATATCTACCTGCTCGATCACAACATGTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCGCTACTGCATATGgtgt                                                                                                    >  1:610819/1‑140 (MQ=255)                             gCTCTGAAACAACCTCTTTAATATCTACCTGCTCGATCACAACATGTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCGCTACTGCATATGGTGTaaaaa                                                                                               >  1:490733/1‑140 (MQ=255)                              ctctGAAACAACCTCTTTAATATCTACCTGCTCGATCACAACATGTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCGCTACTGCATATGGTGTaaaaaa                                                                                              <  1:256946/140‑1 (MQ=255)                                   aaaCAACCTCTTTAATATCTACCTGCTCGATCACAACATGTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCGCTACTGCATATGGTGTAAAAAAATGaa                                                                                         <  1:710906/140‑1 (MQ=255)                                   aaaCAACCTCTTTAATATCTACCTGCTCGATCACAACATGTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCGCTACTGCATATGGTGTAAAAAAATGaa                                                                                         <  1:169396/140‑1 (MQ=255)                                      caacCTCTTTAATATCTACCTGCTCGATCACAACATGTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCGCTACTGCATATGGTGTAAAAAAATGAACGc                                                                                      <  2:490733/140‑1 (MQ=255)                                        acCTCTTTAATATCTACCTGCTCGATCACAACATGTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCGCTACTGCATATGGTGTAAAAAAATGAACGCAt                                                                                    <  2:610819/140‑1 (MQ=255)                                        acCTCTTTAATATCTACCTGCTCGATCACAACATGTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCGCTACTGCATATGGTGTAAAAAAATGAACGCAt                                                                                    >  2:428876/1‑140 (MQ=255)                                           tctTTAATATCTACCTGCTCGATCACAACATGTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCGCTACTGCATATGGTGTAAAAAAATGAACGCATCTg                                                                                 <  1:501179/140‑1 (MQ=255)                                             tttAATATCTACCTGCTCGATCACAACATGTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCGCTACTGCATATGGTGTAAAAAAATGAACGCATCTGAt                                                                               >  2:126186/1‑140 (MQ=255)                                                    tCTACCTGCTCGATCACAACATGTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCGCTACTGCATATGGTGTAAAAAAATGAACGCATCTGATGCACTcc                                                                        <  1:126186/140‑1 (MQ=255)                                                     cTACCTGCTCGATCACAACATGTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCGCTACTGCATATGGTGTAAAAAAATGAACGCATCTGATGCACTCCg                                                                       <  2:220506/140‑1 (MQ=255)                                                       aCCTGCTCGATCACAACATGTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCGCTACTGCATATGGTGTAAAAAAATGAACGCATCTGATGCACTCCGgt                                                                     >  1:104329/1‑140 (MQ=255)                                                       aCCTGCTCGATCACAACATGTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCGCTACTGCATATGGTGTAAAAAAATGAACGCATCTGATGCACTCCGgt                                                                     >  2:684886/1‑140 (MQ=255)                                                       aCCTGCTCGATCACAACATGTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCGCTACTGCATATGGTGTAAAAAAATGAACGCATCTGATGCACTCCGgt                                                                     >  2:19539/1‑140 (MQ=255)                                                             tCGATCACAACATGTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCGCTACTGCATATGGTGTAAAAAAATGAACGCATCTGATGCACTCCGGTGTCATa                                                               >  1:769033/1‑140 (MQ=255)                                                                 tCACAACATGTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCGCTACTGCATATGGTGTAAAAAAATGAACGCATCTGATGCACTCCGGTGTCATAAatt                                                           >  1:664413/1‑140 (MQ=255)                                                                  cacaACATGTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCGCTACTGCATATGGTGTAAAAAAATGAACGCATCTGATGCACTCCGGTGTCATAAatta                                                          >  2:71984/1‑140 (MQ=255)                                                                      acaTGTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCGCTACTGCATATGGTGTAAAAAAATGAACGCATCTGATGCACTCCGGTGTCATAAATTAtttt                                                      >  2:624413/1‑140 (MQ=255)                                                                        aTGTTAAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCGCTACTGCATATGGTGTAAAAAAATGAACGCATCTGATGCACTCCGGTGTCATAAATTAtttttt                                                    <  2:54964/140‑1 (MQ=255)                                                                          gTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCGCTACTGCATATGGTGTAAAAAAATGAACGCATCTGATGCACTCCGGTGTCATAAATTATTTTTTGa                                                  >  2:599119/1‑140 (MQ=255)                                                                          gTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCGCTACTGCATATGGTGTAAAAAAATGAACGCATCTGATGCACTCCGGTGTCATAAATTATTTTTTGa                                                  >  1:50601/1‑140 (MQ=255)                                                                                    aaTATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCGCTACTGCATATGGTGTAAAAAAATGAACGCATCTGATGCACTCCGGTGTCATAAATTATTTTTTGATACCTGaata                                        >  1:923326/1‑140 (MQ=255)                                                                                          aaaTTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCGCTACTGCATATGGTGTAAAAAAATGAACGCATCTGATGCACTCCGGTGTCATAAATTATTTTTTGATACCTGAATAATGTCa                                  >  2:629821/1‑140 (MQ=255)                                                                                                        tctcAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCGCTACTGCATATGGTGTAAAAAAATGAACGCATCTGATGCACTCCGGTGTCATAAATTATTTTTTGATACCTGAATAATGTCATTTTGCATaaaaaa                    <  1:134525/140‑1 (MQ=255)                                                                                                          tcAGCAAGTGAGTGATGATCTGTAATGCCTTCTTGCTGTTCGCTACTGCATATGGTGTAAAAAAATGAACGCATCTGATGCACTCCGGTGTCATAAATTATTTTTTGATACCTGAATAATGTCATTTTGCATAAAAaaca                  <  2:104329/140‑1 (MQ=255)                                                                                                                  tgagtgaTGATCTGTAATGCCTTCTTGCTGTTCGCTACTGCATATGGTGTAAAAAAATGAACGCATCTGATGCACTCCGGTGTCATAAATTATTTTTTGATACCTGAATAATGTCATTTTGCATAAAAAACAACCaaaaa          >  2:350644/1‑140 (MQ=255)                                                                                                                          gatCTGTAATGCCTTCTTGCTGTTCGCTACTGCATATGGTGTAAAAAAATGAACGCATCTGATGCACTCCGGTGTCATAAATTATTTTTTGATACCTGAATAATGTCATTTTGCATAAAAAACAACCAAAAATAGAGAGt  <  2:10411/140‑1 (MQ=255)                                                                                                                          gatCTGTAATGCCTTCTTGCTGTTCGCTACTGCATATGGTGTAAAAAAATGAACGCATCTGATGCACTCCGGTGTCATAAATTATTTTTTGATACCTGAATAATGTCATTTTGCATAAAAAACAACCAAAAATAGAGAGt  >  2:756553/1‑140 (MQ=255)                                                                                                                            tCTGTAATGCCTTCTTGCTGTTCGCTACTGCATATGGTGTAAAAAAATGAACGCATCTGATGCACTCCGGTGTCATAAATTATTTTTTGATACCTGAATAATGTCATTTTGCATAAAAAACAACCAAAAATa        <  1:222157/132‑1 (MQ=255)                                                                                                                            tCTGTAATGCCTTCTTGCTGTTCGCTACTGCATATGGTGTAAAAAAATGAACGCATCTGATGCACTCCGGTGTCATAAATTATTTTTTGATACCTGAATAATGTCATTTTGCATAAAAAACAACCAAAAATa        >  2:222157/1‑132 (MQ=255)                                                                                                                                      |                                                                                                                                TGATGAATCAATTGTCTGAGTAAACGTTGCTCTGAAACAACCTCTTTAATATCTACCTGCTCGATCACAACATGTTCAGGTAGAATATTAAATTTTTGGCATATCTCAGCAAGTGAGTGATGATCTGTAATGCTTTCTTGCTGTTCGCTACTGCATATGGTGTAAAAAAATGAACGCATCTGATGCACTCCGGTGTCATAAATTATTTTTTGATACCTGAATAATGTCATTTTGCATAAAAAACAACCAAAAATAGAGAGT  >  NC_000913/293611‑293871

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