Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,435,750	C→T	G272G (GGG→GGA)	cpdB ←	2',3'‑cyclic‑nucleotide 2'‑phosphodiesterase/3'‑nucleotidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,435,750	0	C	T	100.0%	116.4 / NA	35	G272G (GGG→GGA)	cpdB	2',3'‑cyclic‑nucleotide 2'‑phosphodiesterase/3'‑nucleotidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/19);  total (16/19)

TGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCCCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGAACGACAACAACATCGGCACCTTTCTCGCGCATTTCAGGCACGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCC  >  NC_000913/4435506‑4435985                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            tGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCg                                                                                                                                                                                                                                       >  1:158836/1‑251 (MQ=255)        gAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCAGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCtt                                                                                                                                                                                                                                         >  2:247091/1‑242 (MQ=255)          aCTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCa                                                                                                                                                                                                                              <  2:96973/251‑1 (MQ=255)                    ggcgtCGTGGTCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTa                                                                                                                                                                                                                     <  2:102161/246‑1 (MQ=255)                      cATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCt                                                                                                                                                                                                                  <  2:222066/251‑1 (MQ=255)                            gATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGaa                                                                                                                                                                                                             >  1:59984/1‑250 (MQ=255)                                      gAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTg                                                                                                                                                                                                    >  2:137720/1‑249 (MQ=255)                                      gAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGa                                                                                                                                                                                                   <  2:158836/250‑1 (MQ=255)                                       aGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTGATTACTGAGTTGTAAGTCGACCACCCCAAGATGATCGCCCCACATTCCTGGCATTACCGCCTGAACACCATTCAGCGTGCCTTTGTCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGc                                                                                                                                                                                                 <  1:154396/251‑1 (MQ=255)                                                  cAGCTTGCTGTCTTCCGCCGCGAGGGCTTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAAc                                                                                                                                                                                       >  2:48135/1‑250 (MQ=255)                                                        gCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATgg                                                                                                                                                                                <  1:109594/251‑1 (MQ=255)                                                             cgtCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGtt                                                                                                                                                                            <  2:59984/248‑1 (MQ=255)                                                                 cgccgcGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGc                                                                                                                                                                       >  1:302827/1‑251 (MQ=255)                                                                      cgAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCCCCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGat                                                                                                                                                                                                                                     <  2:220432/184‑1 (MQ=255)                                                                      cgAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGat                                                                                                                                                                                                                                     >  1:220432/1‑184 (MQ=255)                                                                                 ttattaGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTg                                                                                                                                                         >  1:122707/1‑249 (MQ=255)                                                                                  tattaGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAg                                                                                                                                                       <  2:160777/250‑1 (MQ=255)                                                                                              cGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGa                                                                                                                                           <  1:31170/250‑1 (MQ=255)                                                                                                              ttCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTAGACCACACCAAGATGATCGCCCCACATACCTGGCATTACCGCCGGAACACCATTCAGCGTGCATTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACtt                                                                                                                          >  2:147044/1‑251 (MQ=255)                                                                                                                                            aCCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAAtgtg                                                                                            >  2:105102/1‑251 (MQ=255)                                                                                                                                            aCCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACGGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAAtgt                                                                                             <  1:302576/250‑1 (MQ=255)                                                                                                                                                cTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGcc                                                                                          >  1:19278/1‑249 (MQ=255)                                                                                                                                                  gTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAg                                                                                        <  2:213339/249‑1 (MQ=255)                                                                                                                                                        aCTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGAACGacaa                                                                                <  2:19278/251‑1 (MQ=255)                                                                                                                                                                gTAAGTCGACCACATCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGAACGACAACAACATCg                                                                        >  2:125610/1‑251 (MQ=255)                                                                                                                                                                       gACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGAACGACAACAACATCGGCACCtt                                                                 >  2:117094/1‑251 (MQ=255)                                                                                                                                                                       gACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGAACGACAACAACATCGGCACCtt                                                                 <  1:137720/251‑1 (MQ=255)                                                                                                                                                                         ccACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGAACGACAACAACATCGGCACCTTtc                                                               <  1:102223/251‑1 (MQ=255)                                                                                                                                                                                      atCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGAACGACAACAACATCGGCACCTTTCTCGCGCATTTCAg                                                  <  1:12368/251‑1 (MQ=255)                                                                                                                                                                                               caTGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGAACGACAACAACATCGGCACCTTTCTCGCGCATTTCAGGCACGTAtt                                         >  2:143160/1‑251 (MQ=255)                                                                                                                                                                                               caTGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGAACGACAACAACATCGGCACCTTTCTCGCGCATTTCAGGCACGTAtt                                         <  1:117094/251‑1 (MQ=255)                                                                                                                                                                                                            aCCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGAACGACAACAACATCGGCACCTTTCTCGCGCATTTCAGGCACGTATTTGCGCACGGTTTc                            >  2:150743/1‑251 (MQ=255)                                                                                                                                                                                                                     acacCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAACCATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGAACGACAACAACATCGGCACCTTTCTCGCGCATTTCAGGCACGTATTTGCGCACGGTTTCGGTAATATc                   >  2:98330/1‑251 (MQ=255)                                                                                                                                                                                                                                   gCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGAACGACAACAACATCGGCACCTTTCTCGCGCATTTCAGGCACGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACt      <  2:187695/250‑1 (MQ=255)                                                                                                                                                                                                                                       ttGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGAACGACAACAACATCGGCACCTTTCTCGCGCATTTCAGGCACGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTcc  <  2:75381/250‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            TGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCCCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGAACGACAACAACATCGGCACCTTTCTCGCGCATTTCAGGCACGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCC  >  NC_000913/4435506‑4435985

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