Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	285,115	A→G	T639A (ACC→GCC)	yagF →	CP4‑6 prophage; D‑xylonate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	285,115	0	A	G	100.0%	120.7 / NA	34	T639A (ACC→GCC)	yagF	CP4‑6 prophage; D‑xylonate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/13);  total (21/13)

GGCGCGCGCGAGCTGGCGCGGCGGCAGACGCACCCGGACCTGCACGCCCACGACTTTTTGCCGGACGACACCCGGCTGTGGGCGGCACTGCAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATT  >  NC_000913/284983‑285250                                                                                                                                     |                                                                                                                                        ggCGCGCGCGAGCTGGCGCGGCGGCAGACGCACCCGGACCTGCACGCCCACGACTTTTTGCCGGACGACACCCGGCTGTGGGCGGCACTGCAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACGCCGATaa                                                                                                                                  >  1:186449/1‑140 (MQ=255)  gcgcgcgcGAGCTGGCGCGGCGGCAGACGCACCCGGACCTGCACGCCCACGACTTTTTGCCGGACGACACCCGGCTGTGGGCGGCACTGCAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACGCCGATaaa                                                                                                                                 >  1:197284/1‑140 (MQ=255)      gcgcgAGCTGGCGCGGCGGCAGACGCACCCGGACCTGCACGCCCACGACTTTTTGCCGGACGACACCCGGCTGTGGGCGGCACTGCAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACGCCGATAAAATTa                                                                                                                             >  2:974439/1‑140 (MQ=255)          gAGCTGGCGCGGCGGCAGACGCACCCGGACCTGCACGCCCACGACTTTTTGCCGGACGACACCCGGCTGTGGGCGGCACTGCAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACGCCGATAAAATTATCGa                                                                                                                         <  1:798810/140‑1 (MQ=255)           aGCTGGCGCGGCGGCAGACGCACCCGGACCTGCACGCCCACGACTTTTTGCCGGACGACACCCGGCTGTGGGCGGCACTGCAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACGCCGATAAAATTATCGAg                                                                                                                        <  2:1149117/140‑1 (MQ=255)            gCTGGCGCGGCGGCAGACGCACCCGGACCTGCACGCCCACGACTTTTTGCCGGACGACACCCGGCTGTGGGCGGCACTGCAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACGCCGATAAAATTATCGAgg                                                                                                                       <  1:909314/140‑1 (MQ=255)                      cggcAGACGCACCCGGACCTGCACGCCCACGACTTTTTGCCGGACGACACCCGGCTGTGGGCGGCACTGCAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACGCCGATAAAATTATCGAGGTAATTAACGc                                                                                                             <  2:162062/140‑1 (MQ=255)                                 cccGGACCTGCACGCCCACGACTTTTTGCCGGACGACACCCGGCTGTGGGCGGCACTGCAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACGCCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAg                                                                                                  <  1:490444/140‑1 (MQ=255)                                  ccGGACCTGCACGCCCACGACTTTTTGCCGGACGACACCCGGCTGTGGGCGGCACTGCAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACGCCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAgc                                                                                                 <  1:974439/140‑1 (MQ=255)                                          gCACGCCCACGACTTTTTGCCGGACGACACCCGGCTGTGGGCGGCACTGCAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACGCCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGaa                                                                                         >  1:493175/1‑140 (MQ=255)                                            aCGCCCACGACTTTTTGCCGGACGACACCCGGCTGTGGGCGGCACTGCAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACGCCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAAtt                                                                                       >  1:78666/1‑140 (MQ=255)                                            aCGCCCACGACTTTTTGCCGGACGACACCCGGCTGTGGGCGGCACTGCAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACGCCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAAtt                                                                                       >  1:227510/1‑140 (MQ=255)                                              gCCCACGACTTTTTGCCGGACGACACCCGGCTGTGGGCGGCACTGCAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACGCCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTa                                                                                     <  2:493175/140‑1 (MQ=255)                                                ccACGACTTTTTGCCGGACGACACCCGGCTGTGGGCGGCACTGCAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACGCCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAat                                                                                   >  2:346373/1‑140 (MQ=255)                                                       tttttGCCGGACGACACCCGGCTGTGGGCGGCACTGCAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACGCCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTa                                                                            <  2:186449/140‑1 (MQ=255)                                                       tttttGCCGGACGACACCCGGCTGTGGGCGGCACTGCAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACGCCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTa                                                                            <  1:346373/140‑1 (MQ=255)                                                                          ggCTGTGGGCGGCACTGCAGTCGGTGAGCGGCGGCACCTGGAAAGGTTGTATTTATGACGCCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCg                                                         >  2:606098/1‑140 (MQ=255)                                                                          ggCTGTGGGCGGCACTGCAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACGCCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCg                                                         >  2:376206/1‑140 (MQ=255)                                                                            cTGTGGGCGGCACTGCAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACGCCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGa                                                       >  2:461712/1‑140 (MQ=255)                                                                               tgGGCGGCACTGCAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACGCCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAata                                                    >  2:303092/1‑140 (MQ=255)                                                                                     gcACTGCAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACGCCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTccc                                              <  2:197284/140‑1 (MQ=255)                                                                                      cACTGCAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACGCCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCg                                             >  2:1047105/1‑140 (MQ=255)                                                                                         tGCAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACGCCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAg                                          >  2:385645/1‑140 (MQ=255)                                                                                           cAGTCGGCGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACGCCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGAc                                        >  2:971374/1‑140 (MQ=255)                                                                                             gTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACGCCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACgg                                      <  1:1009319/140‑1 (MQ=255)                                                                                                    aGCGGCGGCACCTGGAAAGGCTGTATTTATGACGCCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTttgtt                               <  1:971374/140‑1 (MQ=255)                                                                                                     gcggcgGCACCTGGAAAGGCTGTATTTATGACGCCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTttgttg                              >  1:958377/1‑140 (MQ=255)                                                                                                        gcggcACCTGGAAAGGCTGTATTTATGACGCCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGttt                           >  1:551985/1‑140 (MQ=255)                                                                                                                   aaaGGCTGTATTTATGACGCCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATttat                >  2:264486/1‑140 (MQ=255)                                                                                                                       gCTGTATTTATGACGCCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAAt            >  2:431807/1‑140 (MQ=255)                                                                                                                          gTATTTATGACGCCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTa         >  1:1034356/1‑140 (MQ=255)                                                                                                                           tATTTATGACGCCGATAAAATTATCAAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTaa        <  1:785186/140‑1 (MQ=255)                                                                                                                            aTTTATGACGCCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTaaa       >  1:839630/1‑140 (MQ=255)                                                                                                                                 tGACGCCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAAtt  >  1:794731/1‑140 (MQ=255)                                                                                                                                     |                                                                                                                                        GGCGCGCGCGAGCTGGCGCGGCGGCAGACGCACCCGGACCTGCACGCCCACGACTTTTTGCCGGACGACACCCGGCTGTGGGCGGCACTGCAGTCGGTGAGCGGCGGCACCTGGAAAGGCTGTATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATT  >  NC_000913/284983‑285250

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: ‑