Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	285,648	C→G	A85G (GCC→GGC)	yagG →	CP4‑6 prophage; putative sugar transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	285,648	0	C	G	100.0%	136.2 / NA	38	A85G (GCC→GGC)	yagG	CP4‑6 prophage; putative sugar transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (20/18);  total (20/18)

ACGTCTTCGGCCTGTCGGCGGGGATTATGGGCACGCTGTTTTTGGTCTCCCGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCAT  >  NC_000913/285501‑285786                                                                                                                                                    |                                                                                                                                           aCGTCTTCGGCCTGTCGGCGGGGATTATGGGCACGCTGTTTTTGGTCTCCCGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCAt                                                                                                                                         <  2:160533/151‑1 (MQ=255)  cGTCTTCGGCCTGTCGGCGGGGATTATGGGCACGCTGTTTTTGGTCTCCCGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATc                                                                                                                                        >  2:525211/1‑151 (MQ=255)      ttCGGCCTGTCGGCGGGGATTATGGGCACGCTGTTTTTGGTCTCCCGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCg                                                                                                                                     >  2:81261/1‑150 (MQ=255)                 ggcggGGATTATGGGCACGCTGTTTTTGGTCTCCCGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTc                                                                                                                         >  2:264643/1‑151 (MQ=255)                    ggggATTATGGGCACGCTGTTTTTGGTCTCCCGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGc                                                                                                                      >  2:413603/1‑151 (MQ=255)                               gCACGCTGTTTTTGGTCTCCCGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCtt                                                                                                           >  2:328575/1‑151 (MQ=255)                                  cGCTGTTTTTGGTCTCCCGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTa                                                                                                        <  2:290062/151‑1 (MQ=255)                                          ttGGTCTCCCGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCgg                                                                                                >  1:402808/1‑151 (MQ=255)                                               ctcCCGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTtc                                                                                           >  2:313685/1‑151 (MQ=255)                                                 cccGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACtt                                                                                            <  2:573740/148‑1 (MQ=255)                                                    gcgTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGc                                                                                      <  2:402801/151‑1 (MQ=255)                                                              cGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCa                                                                              >  1:33329/1‑149 (MQ=255)                                                              cGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAg                                                                            <  2:132615/151‑1 (MQ=255)                                                                  gTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGatca                                                                        <  1:317610/151‑1 (MQ=255)                                                                  gTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGatca                                                                        <  1:598715/151‑1 (MQ=255)                                                                      ccGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTa                                                                    <  2:318914/151‑1 (MQ=255)                                                                       cGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTAc                                                                   <  2:139872/151‑1 (MQ=255)                                                                       cGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTAc                                                                   <  1:234870/151‑1 (MQ=255)                                                                         aCCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGc                                                                 >  1:7595/1‑151 (MQ=255)                                                                          cccGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGcc                                                                >  1:623598/1‑151 (MQ=255)                                                                                 aTGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTg                                                          >  2:324857/1‑150 (MQ=255)                                                                                   ggggCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGAcc                                                       >  2:421173/1‑151 (MQ=255)                                                                                           tgGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATtct                                               >  2:351868/1‑151 (MQ=255)                                                                                               aGACCGCGCCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTg                                           <  2:279374/151‑1 (MQ=255)                                                                                                gACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTg                                           >  1:216352/1‑150 (MQ=255)                                                                                                 aCCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGAc                                         >  1:290175/1‑151 (MQ=255)                                                                                                    gcacccgcacGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCt                                      <  2:16568/151‑1 (MQ=255)                                                                                                     cacccgcacGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTg                                     >  1:472278/1‑151 (MQ=255)                                                                                                         cgcacGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCt                                 <  1:573666/151‑1 (MQ=255)                                                                                                                ggCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCtt                          <  2:597029/151‑1 (MQ=255)                                                                                                                 gCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTc                         <  1:301885/151‑1 (MQ=255)                                                                                                                    cGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGtt                      >  2:319501/1‑151 (MQ=255)                                                                                                                        cAGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCTTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAAc                   <  2:401065/150‑1 (MQ=255)                                                                                                                         aGTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGt                 >  1:356796/1‑151 (MQ=255)                                                                                                                          gTTCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTg                >  1:412392/1‑151 (MQ=255)                                                                                                                            tCCGCCCGTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGc               <  1:328582/150‑1 (MQ=255)                                                                                                                                   gTTCCTGCTGTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTgc       >  2:504891/1‑151 (MQ=255)                                                                                                                                         gctgTGGGGGGGCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCAt  <  1:183407/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                           ACGTCTTCGGCCTGTCGGCGGGGATTATGGGCACGCTGTTTTTGGTCTCCCGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCTGTGGGGGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGCGTGACCTACATTCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCAT  >  NC_000913/285501‑285786

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