Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	283,339	A→G	T47A (ACG→GCG)	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	283,339	0	A	G	100.0%	119.6 / NA	33	T47A (ACG→GCG)	yagF	CP4‑6 prophage; D‑xylonate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/16);  total (17/16)

ACCATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCATTG  >  NC_000913/283204‑283469                                                                                                                                        |                                                                                                                                   aCCATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGGCGCa                                                                                                                                <  1:81578/140‑1 (MQ=255)           aaaTTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGGCGCAGAACGCCggg                                                                                                                      <  2:641648/140‑1 (MQ=255)           aaaTTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGGCGCAGAACGCCggg                                                                                                                      <  1:455409/140‑1 (MQ=255)                  cACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGGCGCAGAACGCCGGGATGGGCt                                                                                                               <  2:25986/140‑1 (MQ=255)                      ccGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGGCGCAGAACGCCGGGATGGGCTGGGa                                                                                                           <  2:589363/140‑1 (MQ=255)                         cAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCAGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGGCGCAGAACGCCGGGATGGGCTGGGACGc                                                                                                        >  2:38534/1‑140 (MQ=255)                                acgCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGGCGCAGAACGCCGGGATGGGCTGGGACGCCAACAAg                                                                                                 <  1:447985/140‑1 (MQ=255)                                            cGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGGCGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCaaa                                                                                     <  2:58112/140‑1 (MQ=255)                                             ggTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGGCGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAag                                                                                    <  2:704709/140‑1 (MQ=255)                                               tGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGGCGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAagag                                                                                  >  1:435732/1‑140 (MQ=255)                                                gATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGGCGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAgg                                                                                 >  2:223161/1‑140 (MQ=255)                                                   cACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGGCGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGc                                                                              <  1:169225/140‑1 (MQ=255)                                                         cgcgGTGGGGCCGCAGGGCGCTCTGTCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGGCGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTa                                                                        >  1:257362/1‑140 (MQ=255)                                                         cgcgGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGGCGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTa                                                                        >  1:92417/1‑140 (MQ=255)                                                         cgcgGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCAGCAACCTGTTCGGCATGGCGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTa                                                                        >  2:275987/1‑140 (MQ=255)                                                         cgcgGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCCGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGGCGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTa                                                                        >  2:293582/1‑140 (MQ=255)                                                           cggcggGGCCGCAGGGCGCTCTGCCGCTAACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGGCGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATc                                                                      >  2:211840/1‑140 (MQ=255)                                                                             ctctGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGGCGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATc                                                    >  1:179947/1‑140 (MQ=255)                                                                             ctctGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGGCGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATc                                                    >  2:123596/1‑140 (MQ=255)                                                                                     cTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGGCGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCgg                                            >  1:43838/1‑140 (MQ=255)                                                                                       gACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGGCGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGgag                                          <  1:211840/140‑1 (MQ=255)                                                                                                   gCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGGCGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGGGACGGACGCCCaa                              <  1:123596/140‑1 (MQ=255)                                                                                                       aTGGAATCTCCCAGCGGCAACCTGTTCGGCATGGCGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATcgc                          <  1:223161/140‑1 (MQ=255)                                                                                                        tGGAATCTCCCAGCGGCAACCTGTTCGGCATGGCGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATcgcg                         >  1:672522/1‑140 (MQ=255)                                                                                                           aaTCTCCCAGCGGCAACCTGTTCGGCATGGCGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGcc                                 <  2:158002/129‑1 (MQ=255)                                                                                                           aaTCTCCCAGCGGCAACCTGTTCGGCATGGCGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGcc                                 >  1:158002/1‑129 (MQ=255)                                                                                                           aaTCTCCCAGCGGCAACCTGTTCGGCATGGCGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTg                      >  2:665604/1‑140 (MQ=255)                                                                                                             tctcCCAGCGGCAACCTGTTCGGCATGGCGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTggg                    >  1:235951/1‑140 (MQ=255)                                                                                                                  cAGCGGCAACCTGTTCGGCATGGCGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTAcc               >  1:701537/1‑140 (MQ=255)                                                                                                                         aaCCTGTTCGGCATGGCGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCgg        <  2:672522/140‑1 (MQ=255)                                                                                                                          aCCTGTTCGGCATGGCGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCggg       <  1:293582/140‑1 (MQ=255)                                                                                                                              gTTCGGCATGGCGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCAtt   <  2:435732/140‑1 (MQ=255)                                                                                                                               ttCGGCATGGCGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACTGACGCCCAATCGCGCTGGGCTACCACACCGGGCATTg  >  1:162166/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                   ACCATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCATTG  >  NC_000913/283204‑283469

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