Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,196,817	A→G	intergenic (+444/+50)	icd → / ← ymfD	isocitrate dehydrogenase/e14 prophage; putative SAM‑dependent methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,196,817	0	A	G	100.0%	128.1 / NA	37	intergenic (+444/+50)	icd/ymfD	isocitrate dehydrogenase/e14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (6/31);  total (6/31)

AAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGAC  >  NC_000913/1196689‑1196924                                                                                                                                 |                                                                                                            aaaTTGCCGGTTTTGTATCAACTACTCA‑CCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGTAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCgcatac                                                                                                        <  1:467176/133‑1 (MQ=255)  aagtGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCATCgcatac                                                                                                        <  2:421656/130‑1 (MQ=255)  aaTTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCgcatac                                                                                                        <  1:6255/133‑1 (MQ=255)  aaTTGCCGATTTTGTCTCGACTCCTCACCCCGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCgcatac                                                                                                        <  1:172138/133‑1 (MQ=255)    ttGCCGGTCTTGTATCAAATACTCACCCGGCACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCgcatacat                                                                                                      <  1:200205/133‑1 (MQ=255)       ccGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTcgca                                                                                                           >  1:283211/1‑125 (MQ=255)        cGGTTTTGTACCAACTACTCACCCTGGACTCGCCAGGGGACTGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCttt                                                                                                  <  2:215183/133‑1 (MQ=255)              tGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACt                                                                                            <  2:473762/133‑1 (MQ=255)                 aTCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTgag                                                                                         <  2:92767/133‑1 (MQ=255)                 aTCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTgag                                                                                         <  1:428826/133‑1 (MQ=255)                  tCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTgaga                                                                                        <  1:308496/133‑1 (MQ=255)                     actactCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTGAGACTg                                                                                     >  1:413277/1‑133 (MQ=255)                     actactCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTGAGACTg                                                                                     >  1:285340/1‑133 (MQ=255)                       tactCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTGAGACTGtt                                                                                   <  1:438176/133‑1 (MQ=255)                               cGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTGAGACTGTTTAACCTAg                                                                           <  1:347657/133‑1 (MQ=255)                                gggACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTGAGACTGTTTAACCTAgg                                                                          <  2:413277/133‑1 (MQ=255)                                 ggACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACCTCTTTAGCACTTAGTCTGTTTAACCTAggg                                                                         >  1:190169/1‑133 (MQ=255)                                   aCTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTGAGACTGTTTAACCTAGGGaa                                                                       <  1:358233/133‑1 (MQ=255)                                      cGCCAGGGGACAGCCCACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAatta                                                                    <  1:191792/133‑1 (MQ=255)                                      cGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAatta                                                                    <  1:95491/133‑1 (MQ=255)                                       gCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAattat                                                                   <  2:285340/133‑1 (MQ=255)                                         cAGGGGACAGCCAACCGGCATTGGTTGCAATCGCCTTAGCGGTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAAttatta                                                                 <  2:167312/133‑1 (MQ=255)                                         cAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAAttatta                                                                 <  2:293026/133‑1 (MQ=255)                                            gggACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCg                                                              >  1:406269/1‑133 (MQ=255)                                              gACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGta                                                            <  1:58270/133‑1 (MQ=255)                                              gACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGta                                                            <  2:404767/133‑1 (MQ=255)                                                 aGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGtatta                                                         <  2:87745/133‑1 (MQ=255)                                                   ccAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATtata                                                       <  1:382576/133‑1 (MQ=255)                                                   ccAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATtata                                                       <  1:356968/133‑1 (MQ=255)                                                    cAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTatat                                                      <  2:406269/133‑1 (MQ=255)                                                                ggTGCAATCATCTTAGCGTTCAGGTACAAGCGAAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACt                                          <  2:190169/133‑1 (MQ=255)                                                                         aCCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGt                                 <  2:351960/133‑1 (MQ=255)                                                                         aCCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGt                                 <  1:506999/133‑1 (MQ=255)                                                                         aCCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGt                                 <  2:421565/133‑1 (MQ=255)                                                                         aCCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGt                                 <  2:469397/133‑1 (MQ=255)                                                                         aCCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGt                                 <  2:47467/133‑1 (MQ=255)                                                                         aCCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGt                                 <  1:415618/133‑1 (MQ=255)                                                                                     aGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACt                     <  1:475136/133‑1 (MQ=255)                                                                                     aGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACt                     >  2:475136/1‑133 (MQ=255)                                                                                             gCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACTACTTCCAACAGTGACAAATGAACTTCCTGCtt             >  2:216162/1‑133 (MQ=255)                                                                                             gCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCtt             <  2:299289/133‑1 (MQ=255)                                                                                                    gTAAAAAAGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCa      <  1:129532/133‑1 (MQ=255)                                                                                                        aaaaGGCCGCGAGCGCGGCCCCTTCGCATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGAc  <  2:394898/133‑1 (MQ=255)                                                                                                                                 |                                                                                                            AAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGAC  >  NC_000913/1196689‑1196924

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