Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,436,065	(T)7→8	coding (501/1944 nt)	cpdB ←	2'3' cyclic nucleotide phosphodiesterase/3' nucleotidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,436,058	1	.	T	100.0%	184.0 / NA	53	Q170K (CAG→AAG)	cpdB	2'3' cyclic nucleotide phosphodiesterase/3' nucleotidase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (29/24);  total (29/24)

TTCTCGCGCATTTCAGGCACGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTG‑TTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAATGACGTTGGCATTTACATAAGGGAATTTCGCTCCTGCCAGCGCATTTTTCAG  >  NC_000913/4435922‑4436190                                                                                                                                          |                                                                                                                                     ttCTCGCGCATTTCAGGCACGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTG‑ttt                                                                                                                                   >  2:347944/1‑140 (MQ=255)      gcgcATTTCAGGCACGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTG‑tttttttt                                                                                                                              >  1:307034/1‑139 (MQ=255)         cATTTCAGGCACGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTG‑ttttttttccg                                                                                                                           >  2:717280/1‑136 (MQ=255)         cATTTCAGGCACGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTG‑ttttttttccg                                                                                                                           <  2:267775/140‑5 (MQ=255)          aTTTCAGGCACGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTG‑ttttttttccgt                                                                                                                          <  2:239264/140‑6 (MQ=255)            ttCAGGCACGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCt                                                                                                                         <  2:344246/140‑1 (MQ=255)            ttCAGGCACGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCt                                                                                                                         >  1:598321/1‑140 (MQ=255)                  cACGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCg                                                                                                                   <  2:598321/140‑1 (MQ=255)                    cGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGAc                                                                                                                 <  1:274623/140‑1 (MQ=255)                    cGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGAc                                                                                                                 <  1:237370/140‑1 (MQ=255)                    cGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGAc                                                                                                                 <  1:347944/140‑1 (MQ=255)                    cGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGAc                                                                                                                 <  1:528047/140‑1 (MQ=255)                     gTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGAcc                                                                                                                <  2:232483/140‑1 (MQ=255)                      tATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCa                                                                                                               >  2:395536/1‑140 (MQ=255)                      tATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCa                                                                                                               >  2:515208/1‑140 (MQ=255)                           gcgcACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCg                                                                                                          <  2:559473/140‑1 (MQ=255)                           gcgcACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCg                                                                                                          <  1:765109/140‑1 (MQ=255)                             gcACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCgga                                                                                                        >  1:634009/1‑139 (MQ=255)                                cGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTc                                                                                                     <  1:822034/140‑1 (MQ=255)                                 ggTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCt                                                                                                    >  1:108975/1‑140 (MQ=255)                                    ttCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTtta                                                                                                 <  2:565686/140‑1 (MQ=255)                                           atatCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAAt                                                                                          >  1:583073/1‑140 (MQ=255)                                           atatCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAAt                                                                                          <  2:752514/140‑1 (MQ=255)                                            tatCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATa                                                                                         <  2:511965/140‑1 (MQ=255)                                            tatCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATa                                                                                         <  2:209203/140‑1 (MQ=255)                                            tatCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATa                                                                                         <  1:395536/140‑1 (MQ=255)                                             atCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATAc                                                                                        >  1:55775/1‑140 (MQ=255)                                             atCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATAc                                                                                        >  1:170202/1‑140 (MQ=255)                                             atCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATAc                                                                                        <  2:55775/140‑1 (MQ=255)                                                 ttCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGgtg                                                                                    <  1:39865/140‑1 (MQ=255)                                                      cGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAAc                                                                               >  2:606847/1‑140 (MQ=255)                                                      cGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAAc                                                                               >  1:732718/1‑140 (MQ=255)                                                        tCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACAt                                                                             >  1:750324/1‑140 (MQ=255)                                                                ccGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGtt                                                                     >  1:207050/1‑140 (MQ=255)                                                                ccGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGtt                                                                     >  1:730042/1‑140 (MQ=255)                                                                   gATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTgg                                                                  >  1:489837/1‑140 (MQ=255)                                                                   gATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTgg                                                                  >  2:765551/1‑140 (MQ=255)                                                                       aaTTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCt                                                              <  2:108975/140‑1 (MQ=255)                                                                        aTTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTg                                                             >  1:521256/1‑140 (MQ=255)                                                                        aTTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTg                                                             >  1:115140/1‑140 (MQ=255)                                                                         ttAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTgg                                                            <  2:750324/140‑1 (MQ=255)                                                                               ttATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCaa                                                      <  2:285988/140‑1 (MQ=255)                                                                                tATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAAt                                                     >  2:32927/1‑140 (MQ=255)                                                                                  tCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGtt                                                                     <  2:317270/122‑1 (MQ=255)                                                                                  tCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGtt                                                                     >  1:317270/1‑122 (MQ=255)                                                                                    ccAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAATGACg                                                 <  1:765551/140‑1 (MQ=255)                                                                                      aGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAATGACGtt                                               <  2:739951/140‑1 (MQ=255)                                                                                          cATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAATGACGTTGGCa                                           >  1:236927/1‑140 (MQ=255)                                                                                            tGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAATGACGTTGGCAtt                                         >  2:287519/1‑140 (MQ=255)                                                                                            tGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAATGACGTTGGCAtt                                         >  2:781444/1‑140 (MQ=255)                                                                                                        gCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAATGACGTTGGCATTTACATAAGGGaa                             >  1:43753/1‑140 (MQ=255)                                                                                                             acgCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAATGACGTTGGCATTTACATAAGGGAATTTCg                        >  2:815271/1‑140 (MQ=255)                                                                                                                         ccAATCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAATGACGTTGGCATTTACATAAGGGAATTTCGCTCCTGCCAgcg            >  2:638019/1‑140 (MQ=255)                                                                                                                             tCTTCAGCGTCTGTTTTTTTTCTGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAATGACGTTGGCATTTACATAAGGGAATTTCGCTCCTGCCAGCGCAtt        >  2:774172/1‑140 (MQ=255)                                                                                                                             tCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAATGACGTTGGCATTTACATAAGGGAATTTCGCTCCTGCCAGCGCAtt        >  1:103128/1‑140 (MQ=255)                                                                                                                             tCTTCAGCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAATGACGTTGGCATTTACATAAGGGAATTTCGCTCCTGCCAGCGCAt         >  2:772723/1‑139 (MQ=255)                                                                                                                             tCTTCAGCGCCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAATGACGTTGGCATTTACATAAGGGAATTTCGCTCCTGCCAGCGCAt         <  1:772723/139‑1 (MQ=255)                                                                                                                                   gCGTCTGTTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAATGACGTTGGCATTTACATAAGGGAATTTCGCTCCTGCCAGCGCATTTTTCAg  >  1:336953/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                     TTCTCGCGCATTTCAGGCACGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTTTATCCCAGCCCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTG‑TTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAATGACGTTGGCATTTACATAAGGGAATTTCGCTCCTGCCAGCGCATTTTTCAG  >  NC_000913/4435922‑4436190

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