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%	151.9 / NA	47	*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 (23/24);  total (23/24)

TGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGAACAAAACCGGAGTGTCCGACAGCCTCACGGCC  >  NC_000913/2770268‑2770525                                                                                                                                         |                                                                                                                          tGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGATAGCCCACCGGAATAacg                                                                                                                         <  2:1010115/139‑1 (MQ=255) tGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAacg                                                                                                                         >  1:813247/1‑139 (MQ=255)       agaCTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCg                                                                                                                   <  2:1236963/139‑1 (MQ=255)           ttttCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCgttgt                                                                                                               >  1:1596431/1‑139 (MQ=255)                              cTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATgg                                                                                            >  2:1289422/1‑139 (MQ=255)                                     gCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGcc                                                                                     <  1:1289422/139‑1 (MQ=255)                                                  cccGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTAt                                                                        <  2:813247/139‑1 (MQ=255)                                                         ttcattcaGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCg                                                                            >  1:1102909/1‑128 (MQ=255)                                                         ttcattcaGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCg                                                                            <  2:1102909/128‑1 (MQ=255)                                                           cattcaGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTAt                                                               >  2:638518/1‑139 (MQ=255)                                                               caGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGcc                                                           <  1:259432/139‑1 (MQ=255)                                                               caGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGcc                                                           <  2:1596431/139‑1 (MQ=255)                                                                 ggCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGc                                                                               >  1:209029/1‑117 (MQ=255)                                                                 ggCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGc                                                                               <  2:209029/117‑1 (MQ=255)                                                                 ggCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCa                                                          <  1:360421/138‑1 (MQ=255)                                                                 ggCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCa                                                          >  2:360421/1‑138 (MQ=255)                                                                  gCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCa                                                          >  2:913583/1‑137 (MQ=255)                                                                  gCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCa                                                          <  1:913583/137‑1 (MQ=255)                                                                  gCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCAtt                                                        >  1:777545/1‑139 (MQ=255)                                                                       aCCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTAt                                                               >  1:765789/1‑127 (MQ=255)                                                                       aCCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTAt                                                               <  2:765789/127‑1 (MQ=255)                                                                       aCCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTa                                                    >  2:779216/1‑138 (MQ=255)                                                                       aCCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTa                                                    <  1:779216/138‑1 (MQ=255)                                                                        ccccATTGTCGCTTTATCCAGCATGCACTGCTTACGTTTGATCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTatc                                                  >  2:983468/1‑139 (MQ=255)                                                                          ccATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTatcat                                                <  2:777545/139‑1 (MQ=255)                                                                                       aTCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCa                                                                                   <  1:1471591/91‑1 (MQ=255)                                                                                       aTCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCa                                                                                   >  2:1471591/1‑91 (MQ=255)                                                                                       aTCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTg                                   <  2:835068/139‑1 (MQ=255)                                                                                             cATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCgg                                                                                                    >  1:227094/1‑68 (MQ=255)                                                                                             cATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCgg                                                                                                    <  2:227094/68‑1 (MQ=255)                                                                                                tCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGAACAAAAcc                          >  2:448845/1‑139 (MQ=255)                                                                                                  aCTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAg                                                                                                                      <  2:1432794/45‑1 (MQ=255)                                                                                                  aCTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAg                                                                                                                      >  1:1432794/1‑45 (MQ=255)                                                                                                  aCTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAAcc                                                                                                       <  1:825012/60‑1 (MQ=255)                                                                                                  aCTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAAcc                                                                                                       >  2:825012/1‑60 (MQ=255)                                                                                                       ttACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCa                                                                                                                       >  2:174477/1‑39 (MQ=255)                                                                                                       ttACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCa                                                                                                                       <  1:174477/39‑1 (MQ=255)                                                                                                       ttACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAg                                                                                  <  2:1023736/76‑1 (MQ=255)                                                                                                       ttACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAg                                                                                  >  1:1023736/1‑76 (MQ=255)                                                                                                             ttGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGAACAAAACCGGAGTGTCCGACa             <  1:638518/139‑1 (MQ=255)                                                                                                                ttCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGAACAAAACCGGAGTGTCCGACAGcc          >  2:1471492/1‑139 (MQ=255)                                                                                                                     gCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATgga                                                                                           >  2:685729/1‑53 (MQ=255)                                                                                                                     gCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATgga                                                                                           <  1:685729/53‑1 (MQ=255)                                                                                                                     gCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGcc                                                                                     <  2:1153776/59‑1 (MQ=255)                                                                                                                     gCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGcc                                                                                     >  1:1153776/1‑59 (MQ=255)                                                                                                                     gCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGAACAAAACCGGAGTGTCCGACAGCCTCACg     >  2:852420/1‑139 (MQ=255)                                                                                                                        gACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGAACAAAACCGGAGTGTCCGACAGCCTCACGGcc  <  1:852420/139‑1 (MQ=255)                                                                                                                         aCAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATgga                                                                                           <  2:88657/49‑1 (MQ=255)                                                                                                                         aCAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATgga                                                                                           >  1:88657/1‑49 (MQ=255)                                                                                                                                         |                                                                                                                          TGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGAACAAAACCGGAGTGTCCGACAGCCTCACGGCC  >  NC_000913/2770268‑2770525

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