Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	277,597	T→C	C295R (TGC→CGC)	mmuM →	CP4‑6 prophage; S‑methylmethionine:homocysteine methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	277,597	0	T	C	100.0%	118.4 / NA	39	C295R (TGC→CGC)	mmuM	CP4‑6 prophage; S‑methylmethionine:homocysteine methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/25);  total (14/25)

GCTGCCGCTGGTGGTGTATCCGAACTCGGGCGAGCATTACGATGCCGTGAGCAAAACCTGGCATCATCACGGTGAACATTGCGCGCAGCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGTGCTGTCGCACCACGCCTGCGGATATCG‑CCGCGTTAAAAGCGCGAAGCTGAGGGTTTATCGGGTCTACATCGTTCATTGTAGGCCTGATAAGCGCAGCGCATCAGGCATTGCCGGATGGCGGCACAAGCGCCTTATCCGG  >  NC_000913/277455‑277736                                                                                                                                               |                                                                                                                                             gctgcCGCTGGTGGTGTATCCGAACTCGGGCGAGCATTACGATGCCGTGAGCAAAACCTGGCATCATCACGGTGAACATTGCGCGCAGCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCg                                                                                                                                       <  1:285170/150‑1 (MQ=255)                gtATCCGAACTCGGGCGAGCATTACGATGCCGTGAGCAAAACCTGGCATCATCACGGTGAACATTGCGCGCAGCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCACCACGCCTGCGGa                                                                                                                        <  1:102544/150‑1 (MQ=255)                  aTCCGAACTCGGGCGAGCATTACGATGCCGTGAGCAAAACCTGGCATCATCACGGTGAACATTGCGCGCAGCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCACCACGCCTGCGGatat                                                                                                                     >  2:582354/1‑151 (MQ=255)                    ccGAACTCGGGCGAGCATTACGATGCCGTGAGCAAAACCTGGCATCATCACGGTGAACATTGCGCGCAGCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCACCACGCCTGCGGATATc                                                                                                                    <  1:246072/150‑1 (MQ=255)                      gAACTCGGGCGAGCATTACGATGCCGTGAGCAAAACCTGGCATCATCACGGTGAACATTGCGCGCAGCTGGCGGAGTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCACCACGCCTGCGGatat                                                                                                                     <  1:460728/147‑1 (MQ=255)                          tCGGGCGAGCATTACGATGCCGTGAGCAAAACCTGGCATCATCACGGTGAACATTGCGCGCAGCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCGGTCGCACCACCCCCGCGGGGATcg‑cc                                                                                                                >  1:407906/1‑147 (MQ=255)                           cGGGCGAGCATTACGATGCCGTGAGCAAAACCTGGCATCATCACGGTGAACATTGCGCGCAGCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCACCACGCCTGCGGATATCG‑Ccgcg                                                                                                             <  1:450549/149‑1 (MQ=255)                           cGGGCGAGCATTACGATGCCGTGAGCAAAACCTGGCATCATCACGGTGAACATTGCGCGCAGCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCACCACGCCTGCGGATATCG‑CCGCGtt                                                                                                           <  1:285938/151‑1 (MQ=255)                              gCGAGCATTACGATGCCGTGAGCAAAACCTGGCATCATCACGGTGAACATTGCGCGCAGCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCACCACGCCTGCGGATATCG‑CCGCGTTaaa                                                                                                        <  2:607613/151‑1 (MQ=255)                               cGAGCATTACGATGCCGTGAGCAAAACCTGGCATCATCACGGTGAACATTGCGCGCAGCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCACCACGCCTGCGGATATCG‑CCGCGTTaaa                                                                                                        >  2:374694/1‑150 (MQ=255)                                gAGCATTACGATGCCGTGAGCAAAACCTGGCATCATCACGGTGAACATTGCGCGCAGCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCACCACGCCTGCGGATATCG‑CCGCGTTAAAAg                                                                                                      <  1:94889/151‑1 (MQ=255)                                     ttACGATGCCGTGAGCAAAACCTGGCATCATCACGGTGAACATTGCGCGCAGCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCACCACGCCTGCGGATATCG‑CCGCGTTAAAAGcgcg                                                                                                  >  2:236384/1‑150 (MQ=255)                                       aCGATGCCGTGAGCAAAACCTGGCATCATCACGGTGAACATTGTGCGCAGCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCCCCGTCGCACCCCCCCCGCGGATATCGCCCGCGTaaaa                                                                                                        >  1:434703/1‑143 (MQ=255)                                               gTGAGCAAAACCTGGCATCATCACGGTGAACATTGCGCGCAGCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCACCACGCCTGCGGATATCG‑CCGCGTTAAAAGCGCGAAGCTGAGGGt                                                                                       <  2:282154/151‑1 (MQ=255)                                                 gAGCAAAACCTGGCATCATCACGGTGAACATTGCGCGCAGCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCACCACGCCTGCGGATATCG‑CCGCGTTAAAAGCGCGAAGCTGAGGGtt                                                                                      <  1:535320/150‑1 (MQ=255)                                                    cAAAACCTGGCATCATCACGGTGAACATTGCGCGCAGCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCACCACGCCTGCGGATATCG‑CCGCGTTAAAAGCGCGAAGCTGAGGGTTTATc                                                                                  <  2:539013/151‑1 (MQ=255)                                                               atcatcACGGTGAACATTGCGCGCAGCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCACCACGCCTGCGGATATCG‑CCGCGTTAAAAGCGCGAAGCTGAGGGTTTATCGGGTCTACATc                                                                       <  1:376565/151‑1 (MQ=255)                                                                  atcaCGGTGAACATTGCGCGCAGCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCACCACGCCTGCGGATATCG‑CCGCGTTAAAAGCGCGAAGCTGAGGGTTTATCGGGTCTACATCGtt                                                                    <  2:279522/151‑1 (MQ=255)                                                                     aCGGTGAACATTGCGCGCAGCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCCCCCAGCCTGCGGGTATCC‑CCGCGTTAAAAGCGCGGAGCCGGGGGGTTATCGGGTCTACATCGTTCAt                                                                 >  2:549091/1‑151 (MQ=255)                                                                         tGAACATTGCGCGCAGCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCACCACGCCTGCGGATATCG‑CCGCGTTAAAAGCGCGAAGCTGAGGGGTTATCGGGTCTACATCGTTCATTGt                                                              >  2:303119/1‑150 (MQ=255)                                                                                 gcgcgcAGCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCACCACGCCTGCGGATATCG‑CCGCGTTAAAAGCGCGAAGCTGAGGGTTTATCGGGTCTACATCGTTCATTGTAGGCCTGAt                                                     <  2:232621/151‑1 (MQ=255)                                                                                 gcgcgcAGCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCACCACGCCTGCGGATATCG‑CCGCGTTAAAAGCGCGAAGCTGAGGGTTTATCGGGTCTACATCGTTCATTGTAGGCCTGAt                                                     <  2:633548/151‑1 (MQ=255)                                                                                  cgcgcAGCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCACCACGCCTGCGGATATCG‑CCGCGTTAAAAGCGCGAAGCTGAGGGGTTATCGGGTCTACATCGTTCATTGTAGGcc                                                         >  2:566567/1‑146 (MQ=255)                                                                                    cgcAGCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCACCACGCCTGCGGATATCG‑CCGCGTTAAAAGCGCGAAGCTGAGGGTTTATCGGGTCTACATCGTTCATTGTAGGCCTGATAag                                                  <  2:554817/151‑1 (MQ=255)                                                                                    cgcAGCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCACCACCACTGGGGAAACCG‑CCGCGTTAAAAGCGCGAAACTGAGGGGCTATTGGGTCCACATCCTGCATTGTTGGCCTgct                                                     >  1:294198/1‑146 (MQ=255)                                                                                      cAGCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCACCACGCCTGCGGATATCG‑CCGCGTTAAAAGCGCGAAGCTGAGGGTTTATCGGGTCTACATCGTTCATTGTAGGCCTGATAagcg                                                <  2:89435/151‑1 (MQ=255)                                                                                       aGCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCCGTCGCACCACGCCTGCGGATATCG‑CCGCGTTAAAAGCGCGAAGCTGAGGGGTTTTCGGGGCTACATCGTTCATTTTAGGCCTGATAagcgc                                               >  2:196932/1‑151 (MQ=255)                                                                                        gCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCACCACGCCTGCGGATATCG‑CCGCGTTAAAAGCGCGAAGGTGAGGGGTTATCGGGTCTACATCGTTTATTGTAGGCCTGATAagcgc                                               >  2:661993/1‑150 (MQ=255)                                                                                            gCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCACCACGCCCGCGGATATCG‑CCGCGTTAAAAGGGCGAAGCTGGGGGGTTATCGGGTCTACATCGGTCTTTGTAGGGCTGATAaacgcaacg                                           >  2:339852/1‑150 (MQ=255)                                                                                                       gCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCACCACGCCTGCGGATATCG‑CCGCGTTAAAAGCGCGAAGCTGAGGGTTTATCGGGTCTACATCGTTCATTGTAGGCCTGATAAGCGCAGCGCATCAGGCa                                  <  2:449002/148‑1 (MQ=255)                                                                                                       gCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCACCACGCCTGCGGATATCG‑CCGCGTTAAAAGCGCGAAGCTGAGGGTTTATCGGGTCTACATCGTTCATTGTAGGCCTGATAAGCGCAGCGCATCAGGCAt                                 <  2:395149/149‑1 (MQ=255)                                                                                                       gCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCACCACGCCTGCGGATATCG‑CCGCGTTAAAAGCGCGAAGCTGAGGGTTTATCGGGTCTACATCGTTCATTGTAGGCCTGATAAGCGCAGCGCATCAGGCATTg                               <  1:446984/151‑1 (MQ=255)                                                                                                             gTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCACCACGCCTGCGGATATCG‑CCGCGTTAAAAGCGCGAAGCTGAGGGTTTATCGGGTCTACATCGTTCATTGTAGGCCTGATAAGCGCAGCGCATCAGGCATTGCCGGAt                         <  1:566539/151‑1 (MQ=255)                                                                                                             gTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCACCACGCCTGCGGATATCG‑CCGCGTTAAAAGCGCGAAGCTGAGGGTGTATCGGGTCTACATCGTTCATTGTAGGCCTGATAAGCGCAGCGCATCAGGCATTGCCGGa                          <  1:396867/150‑1 (MQ=255)                                                                                                                  aGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCACCACGCCTGCGGATATCG‑CCGCGTTAAAAGCGCGAAGCTGAGGGTTTATCGGGTCTACATCGTTCATTGTAGGCCTGATAAGCGCAGCGCATCAGGCATTGCCGGATggcgg                    >  1:190593/1‑151 (MQ=255)                                                                                                                  aGGCCGCTGGCGCACGGTTGATTGGCGGGCGCTGTCGCACCACGCCTGCGGATATCG‑CCGCGTTAAAAGCGCGAAGCTGAGGGGTTATCGGGTCTACATCGTTCATTGTAGGCCTGATAAGCGCAGCGCATCAGGGATTGCCGGATggcgg                    >  1:404028/1‑151 (MQ=255)                                                                                                                          ggCGCACGGTTGATTGGCGGGCGCTGTCGCACCACGCCTGCGGATATCG‑CCGCGTTAAAAGCGCGAAGCTGAGGGTTTATCGGGTCTACATCGTTCATTGTAGGCCTGATAAGCGCAGCGCATCAGGCATTGCCGGATGGCGGCACAAgc             <  1:607650/150‑1 (MQ=255)                                                                                                                                cGGTTGATTGGCGGGCGCTGTCGCACCACGCCTGCGGATATCG‑CCGCGTTAAAAGCGCGAAGCTGAGGGTTTATCGGGTCTACATCGTTCATTGTAGGCCTGATAATCGCAGCGCATCAGGCATTGCCGGATGGCGGCACAAGCGCCTTAt      <  2:321741/151‑1 (MQ=255)                                                                                                                                     gATTGGCGGGCGCTGTCGCACCACGCCTGCGGATATCG‑CCGCGTTAAAAGCGCGAAGCTGAGGGTTTATCGGGTCTACATCGTTCATTGTAGGCCTGATAAGCGCAGCGCATCAGGCATTGCCGGATGGCGGCACAAGCGCCTTATCCgg  <  2:585309/150‑1 (MQ=255)                                                                                                                                               |                                                                                                                                             GCTGCCGCTGGTGGTGTATCCGAACTCGGGCGAGCATTACGATGCCGTGAGCAAAACCTGGCATCATCACGGTGAACATTGCGCGCAGCTGGCGGATTATCTGCCGCAGTGGCAGGCCGCTGGCGCACGGTTGATTGGCGGGTGCTGTCGCACCACGCCTGCGGATATCG‑CCGCGTTAAAAGCGCGAAGCTGAGGGTTTATCGGGTCTACATCGTTCATTGTAGGCCTGATAAGCGCAGCGCATCAGGCATTGCCGGATGGCGGCACAAGCGCCTTATCCGG  >  NC_000913/277455‑277736

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