Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	578,636	A→T	Y86* (TAT→TAA)	borD ←	DLP12 prophage; prophage lipoprotein BorD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	578,636	0	A	T	100.0%	132.5 / NA	36	Y86* (TAT→TAA)	borD	DLP12 prophage; prophage lipoprotein BorD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (20/16);  total (20/16)

AGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTG  >  NC_000913/578563‑578712                                                                          |                                                                             aGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTTTAAATGc                                                                       >  2:120965/1‑81 (MQ=255)    tAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTTTAAATGCCTa                                                                    <  2:120938/81‑1 (MQ=255)      gagCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTTTAAATGCCTa                                                                    >  1:120987/1‑79 (MQ=255)            aCCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTTTAAATGCCTAAAGtaat                                                             <  1:120939/80‑1 (MQ=255)            aCCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTTTAAATGCCTAAAGtaat                                                             <  2:120940/80‑1 (MQ=255)                tatCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTTTAAATGCCTAAAGTAATAAAAc                                                        >  2:120964/1‑81 (MQ=255)                    gATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTTTAAATGCCTAAAGTAATAAAACCGAg                                                    <  1:120941/81‑1 (MQ=255)                     aTGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTTTAAATGCCTAAAGTAATAAAACCGAg                                                    >  2:120966/1‑80 (MQ=255)                      tGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTTTAAATGCCTAAAGTAATAAAACCGAGCa                                                  <  1:120942/81‑1 (MQ=255)                       gggCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTTTAAATGCCTAAAGTAATAAAACCGAGCaa                                                 >  1:120974/1‑81 (MQ=255)                          cAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTTTAAATGCCTAAAGTAATAAAACCGAGCAATcc                                              <  1:120943/81‑1 (MQ=255)                          cAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTTTAAATGCCTAAAGTAATAAAACCGAGCAATc                                               >  2:120972/1‑80 (MQ=255)                           aaCTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTTTAAATGCCTAAAGTAATAAAACCGAGCAATCCa                                             >  1:120970/1‑81 (MQ=255)                           aaCTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTTTAAATGCCTAAAGTAATAAAACCGAGCAATCCa                                             >  2:120979/1‑81 (MQ=255)                           aaCTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTTTAAATGCCTAAAGTAATAAAACCGAGCAATCCa                                             >  2:120975/1‑81 (MQ=255)                             cTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTTTAAATGCCTAAAGTAATAAAACCGAGCAATCCAtt                                           >  1:120978/1‑81 (MQ=255)                                 tGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTTTAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACg                                       <  2:120944/81‑1 (MQ=255)                                   cAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTTTAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGaa                                     <  1:120945/81‑1 (MQ=255)                                      ttattTTGAGCAATACACACGCGCTTCCAGCGGAGTTTAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGt                                  <  1:120946/81‑1 (MQ=255)                                       tattTTGAGCAATACACACGCGCTTCCAGCGGAGTTTAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGtt                                 >  2:120980/1‑81 (MQ=255)                                              aGCAATACACACGCGCTTCCAGCGGAGTTTAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTggg                          >  2:120988/1‑81 (MQ=255)                                              aGCAATACACACGCGCTTCCAGCGGAGTTTAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTgg                           <  2:120947/80‑1 (MQ=255)                                               gCAATACACACGCGCTTCCAGCGGAGTTTAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGt                         >  1:120976/1‑81 (MQ=255)                                                cAATACACACGCGCTTCCAGCGGAGTTTAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGtt                        <  2:120948/81‑1 (MQ=255)                                                   tACACACGCGCTTCCAGCGGAGTTTAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCt                     <  2:120949/81‑1 (MQ=255)                                                    acacacGCGCTTCCAGCGGAGTTTAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTg                    >  1:120983/1‑81 (MQ=255)                                                    acacacGCGCTTCCAGCGGAGTTTAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTg                    >  2:120982/1‑81 (MQ=255)                                                      acacGCGCTTCCAGCGGAGTTTAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGtt                  >  2:120986/1‑81 (MQ=255)                                                       cacGCGCTTCCAGCGGAGTTTAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGttt                 <  1:120951/81‑1 (MQ=255)                                                       cacGCGCTTCCAGCGGAGTTTAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGttt                 <  1:120950/81‑1 (MQ=255)                                                         cgcgCTTCCAGCGGGGTTTAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGtttt                >  1:120984/1‑80 (MQ=255)                                                          gcgcTTCCAGCGGAGTTTAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTaa              >  2:120991/1‑81 (MQ=255)                                                            gcTTCCAGCGGAGTTTAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTaac             >  2:120977/1‑80 (MQ=255)                                                             cTTCCAGCGGAGTTTAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTaacaa           >  2:120981/1‑81 (MQ=255)                                                                ccAGCGGAGTTTAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAacaaca         <  2:120952/80‑1 (MQ=255)                                                                      gAGTTTAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTg  <  2:120953/81‑1 (MQ=255)                                                                          |                                                                             AGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTG  >  NC_000913/578563‑578712

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