Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,770,404	G→A	*234* (TAG→TAA)	yfjR →	CP4‑57 prophage, putative DNA‑binding transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,770,404	0	G	A	100.0%	156.2 / NA	49	*234* (TAG→TAA)	yfjR	CP4‑57 prophage, putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/31);  total (18/31)

TGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGAACAAAACCGGAGTGTCCGACAGCCTCACGGCCCTG  >  NC_000913/2770268‑2770528                                                                                                                                         |                                                                                                                             tGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAacg                                                                                                                            >  2:441281/1‑139 (MQ=255) tGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAacg                                                                                                                            >  1:735094/1‑139 (MQ=255) tGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAacg                                                                                                                            >  2:572422/1‑139 (MQ=255) tGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCGTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAacg                                                                                                                            <  2:735094/139‑1 (MQ=255) tGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGGCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAacg                                                                                                                            >  1:588020/1‑139 (MQ=255)    aGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAg                                                                                                                         <  1:971685/139‑1 (MQ=255)    aGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAg                                                                                                                         <  1:1064098/139‑1 (MQ=255)       agaCTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCg                                                                                                                      >  1:724420/1‑139 (MQ=255)        gaCTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCgt                                                                                                                     >  1:252306/1‑139 (MQ=255)           ttttCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCgttgt                                                                                                                  >  1:463361/1‑139 (MQ=255)                   gcAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAAcc                                                                                                          <  1:420843/139‑1 (MQ=255)                   gcAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAAcc                                                                                                          <  1:572422/139‑1 (MQ=255)                      aTGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCgg                                                                                                       <  1:570441/139‑1 (MQ=255)                      aTGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCgg                                                                                                       <  1:774170/139‑1 (MQ=255)                      aTGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCgg                                                                                                       <  2:463361/139‑1 (MQ=255)                      aTGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACGGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCgg                                                                                                       <  2:518809/139‑1 (MQ=255)                           aGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACa                                                                                                  <  2:539230/139‑1 (MQ=255)                                 tgtCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATggagg                                                                                            >  2:739481/1‑139 (MQ=255)                                 tgtCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATggag                                                                                             <  1:58374/138‑1 (MQ=255)                                 tgtCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATggag                                                                                             >  2:58374/1‑138 (MQ=255)                                     gCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGcc                                                                                        <  1:1126335/139‑1 (MQ=255)                                                     ggTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGc                                                                        >  1:870499/1‑139 (MQ=255)                                                     ggTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGc                                                                        >  2:611187/1‑139 (MQ=255)                                                     ggTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGc                                                                        >  1:870500/1‑139 (MQ=255)                                                        cTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCAt                                                                     <  1:203905/139‑1 (MQ=255)                                                        cTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCAt                                                                     <  2:50082/139‑1 (MQ=255)                                                             ttcaGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTGGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCg                                                                <  2:724420/139‑1 (MQ=255)                                                              tcaGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATTAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTAACGc                                                               <  2:252306/139‑1 (MQ=255)                                                              tcaGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGc                                                               <  2:870500/139‑1 (MQ=255)                                                              tcaGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGc                                                               <  2:870499/139‑1 (MQ=255)                                                               caGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGcc                                                              >  1:885430/1‑139 (MQ=255)                                                                aGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCa                                                             <  1:441281/139‑1 (MQ=255)                                                                 ggCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCgt                                                                                                                     <  1:216151/82‑1 (MQ=255)                                                                 ggCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCgt                                                                                                                     >  2:216151/1‑82 (MQ=255)                                                                 ggCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCgg                                                                                                       >  2:1022989/1‑96 (MQ=255)                                                                 ggCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCgg                                                                                                       <  1:1022989/96‑1 (MQ=255)                                                                     ttACCCCCTTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGtt                                                        <  2:925514/139‑1 (MQ=255)                                                                       aCCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCa                                                             >  2:685796/1‑132 (MQ=255)                                                                       aCCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCa                                                             <  1:685796/132‑1 (MQ=255)                                                                       aCCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTat                                                      >  1:275169/1‑139 (MQ=255)                                                                        ccccATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGcc                                                                       >  2:46725/1‑121 (MQ=255)                                                                        ccccATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGcc                                                                       <  1:46725/121‑1 (MQ=255)                                                                                       aTCCAGCATTCACTGCTTACGTTTGTTCCGGCCGTCAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGccc                                                                                       >  2:661758/1‑90 (MQ=255)                                                                                       aTCCAGCATTCACTGCTTACGTTTGTTCCGGCCGTCAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGccc                                                                                       <  1:661758/90‑1 (MQ=255)                                                                                         ccAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGaa                                    <  1:1057934/139‑1 (MQ=255)                                                                                                 cACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGAACAAAACCg                            <  2:320587/139‑1 (MQ=255)                                                                                                 cACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGAACAAAACCg                            <  1:131126/139‑1 (MQ=255)                                                                                                      cTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCg                                                                                                                      <  1:499426/44‑1 (MQ=255)                                                                                                      cTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCg                                                                                                                      >  2:499426/1‑44 (MQ=255)                                                                                                              tGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGAACAAAACCGGAGTGTCCGACAg               >  1:784516/1‑139 (MQ=255)                                                                                                                ttCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGAACAAAACCGGAGTGTCCGACAGcc             >  2:365935/1‑139 (MQ=255)                                                                                                                           aGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGAACAAAACCGGAGTGTCCGACAGCCTCACGGCCCTg  <  1:611187/139‑1 (MQ=255)                                                                                                                           aGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGAACAAAACCGGAGTGTCCGACAGCCTCACGGCCCTg  <  1:739481/139‑1 (MQ=255)                                                                                                                           aGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGAACAAAACCGGAGTGTCCGACAGCCTCACGGCCCTg  <  1:365935/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                             TGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGAACAAAACCGGAGTGTCCGACAGCCTCACGGCCCTG  >  NC_000913/2770268‑2770528

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