Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,436,450	A→G	F39S (TTC→TCC)	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,450	0	A	G	100.0%	189.5 / NA	53	F39S (TTC→TCC)	cpdB	2'3' cyclic nucleotide phosphodiesterase/3' nucleotidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/29);  total (24/29)

GGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTGTTTTTCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGAAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCAGATCGACCGTCGCTGCATTCACACTGGCGGCAATCAGCGTGGCCAGGAGCGTTGCGCTAAACTTAATCATCAGGGACATCCTTTTATCAT  >  NC_000913/4436314‑4436585                                                                                                                                         |                                                                                                                                        ggCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTGTTTTTCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAAt                                                                                                                                      >  1:107507/1‑140 (MQ=255)  gCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTGTTTTTCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATc                                                                                                                                     >  2:438676/1‑140 (MQ=255)  gCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTGTTTTTCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATc                                                                                                                                     <  2:88493/140‑1 (MQ=255)  gCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTGTTTTTCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATc                                                                                                                                     <  2:444230/140‑1 (MQ=255)   ccAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTGTTTTTCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATcc                                                                                                                                    >  2:535162/1‑140 (MQ=255)   ccAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTGTTTTTCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATcc                                                                                                                                    >  1:63148/1‑140 (MQ=255)      gCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTGTTTTTCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCcatc                                                                                                                                 >  2:163669/1‑140 (MQ=255)       cGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTGTTTTTCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCcatca                                                                                                                                <  1:361734/140‑1 (MQ=255)       cGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTGTTTTTCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCcatca                                                                                                                                >  2:199404/1‑140 (MQ=255)       cGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTGTTTTTCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCcatca                                                                                                                                <  2:491262/140‑1 (MQ=255)         gACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTGTTTTTCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATg                                                                                                                              <  2:33807/140‑1 (MQ=255)                 tGAATCAAATCGCCGTTATCAACCAGTACGCTGTTTTTCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGcatag                                                                                                                      <  1:325167/140‑5 (MQ=255)                 tGAATCAAATCGCCGTTATCAACCAGTACGCTGTTTTTCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATg                                                                                                                      <  2:28671/140‑1 (MQ=255)                   aaTCAAATCGCCGTTATCAACCAGTACGCTGTTTTTCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCa                                                                                                                    <  2:63148/140‑1 (MQ=255)                      cAAATCGCCGTTATCAACCAGTACGCTGTTTTTCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGAt                                                                                                                 <  1:535162/140‑1 (MQ=255)                      cAAATCGCCGTTATCAACCAGTACGCTGTTTTTCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGAt                                                                                                                 <  1:483554/140‑1 (MQ=255)                        aaTCGCCGTTATCAACCAGTACGCTGTTTTTCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCa                                                                                                               <  2:441588/140‑1 (MQ=255)                          tCGCCGTTATCAACCAGTACGCTGTTTTTCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCAGt                                                                                                             <  1:199404/140‑1 (MQ=255)                             ccGTTATCAACCAGTACGCTGTTTTTCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCAGTGGt                                                                                                          >  2:321371/1‑140 (MQ=255)                                 tATCAACCAGTACGCTGTTTTTCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTcc                                                                                                      >  1:43608/1‑140 (MQ=255)                                   tCAACCAGTACGCTGTTTTTCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCAt                                                                                                    <  1:321371/140‑1 (MQ=255)                                          gTACGCTGTTTTTCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTcc                                                                                                      <  2:387916/131‑1 (MQ=255)                                          gTACGCTGTTTTTCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTcc                                                                                                      >  1:387916/1‑131 (MQ=255)                                          gTACGCTGTTTTTCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCATGATACGc                                                                                             >  1:466547/1‑140 (MQ=255)                                           tACGCTGTTTTTCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCa                                                                                            >  2:30669/1‑140 (MQ=255)                                           tACGCTGTTTTTCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCa                                                                                            >  2:496935/1‑140 (MQ=255)                                           tACGCTGTTTTTCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCa                                                                                            >  2:8825/1‑140 (MQ=255)                                           tACGCTGTTTTTCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCa                                                                                            >  1:355765/1‑140 (MQ=255)                                               cTGTTTTTCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGGAATAATCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCAGATc                                                                                        >  1:458358/1‑140 (MQ=255)                                                  tttttCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATTATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCAGATCGAc                                                                                     >  2:441131/1‑140 (MQ=255)                                                  tttttCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCAGATCGAc                                                                                     <  1:30669/140‑1 (MQ=255)                                                        aCTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCAGATCGACCGTCGc                                                                               <  1:163669/140‑1 (MQ=255)                                                          ttCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCAGATCGACCGTCGCTg                                                                             >  2:281656/1‑140 (MQ=255)                                                            cATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCAGATCGACCGTCGCTGCa                                                                           >  2:297987/1‑140 (MQ=255)                                                                   ggCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCAGATCGACCGTCGCTGCATTCacac                                                                    <  2:257879/140‑1 (MQ=255)                                                                   ggCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCAGATCGACCGTCGCTGCATTCacac                                                                    <  1:530680/140‑1 (MQ=255)                                                                      aTCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCAGATCGACCGTCGCTGCATTCACACTgg                                                                 >  2:8487/1‑140 (MQ=255)                                                                                 ggCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCAGATCGACCGTCGCTGCATTCACACTGGCGGCAATCAGc                                                      <  2:107507/140‑1 (MQ=255)                                                                                 ggCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCAGATCGACCGTCGCTGCATTCACACTGGCGGCAATCAGc                                                      <  1:8487/140‑1 (MQ=255)                                                                                 ggCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCAGATCGACCGTCGCTGCATTCACACTGGCGGCAATCAGc                                                      <  2:450571/140‑1 (MQ=255)                                                                                            cgtacCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCAGATCGACCGTCGCTGCATTCACACTGGCGGCAATCAGc                                                      <  2:86725/129‑1 (MQ=255)                                                                                            cgtacCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCAGATCGACCGTCGCTGCATTCACACTGGCGGCAATCAGc                                                      >  1:86725/1‑129 (MQ=255)                                                                                                   gTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCAGATCGACCGTCGCTGCATTCACACTGGCGGCAATCAGCGTGGCCAGGAGCGTTgcg                                    <  1:297987/140‑1 (MQ=255)                                                                                                    tCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATTATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCAGATCGACCGTCGCTGCATTCACACTGGCGGCAATCAGCGTGGCCAGGAGCGTTgcgc                                   <  1:441131/140‑1 (MQ=255)                                                                                                    tCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCAGATCGACCGTCGCTGCATTCACACTGGCGGCAATCAGCGTGGCCAGGAGCGTTgcgc                                   <  2:325282/140‑1 (MQ=255)                                                                                                    tCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCAGATCGACCGTCGCTGCATTCACACTGGCGGCAATCAGCGTGGCCAGGAGCGTTgcgc                                   >  1:506297/1‑140 (MQ=255)                                                                                                    tCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCAGATCGACCGTCGCTGCATTCACACTGGCGGCAATCAGCGTGGCCAGGAGCGTTgcgc                                   >  2:27860/1‑140 (MQ=255)                                                                                                    tCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCAGATCGACCGTCGCTGCATTCACACTGGCGGCAATCAGCGTGGCCAGGAGCGTTgcgc                                   >  2:348710/1‑140 (MQ=255)                                                                                                    tCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCAGATCGACCGTCGCTGCATTCACACTGGCGGCAATCAGCGTGGCCAGGAGCGATgcgc                                   >  2:481580/1‑140 (MQ=255)                                                                                                     ccGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCAGATCGACCGTCGCTGCATTCACACTGGCGGCAATCAGCGTGGCCAGGAGCGTTGCGCt                                  <  1:27860/140‑1 (MQ=255)                                                                                                            tttCCGTGGCGGTGTCTTTGTAATAATCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCAGATCGACCGTCGCTGCATTCACACTGGCGGCAATCAGCGTGGCCAGGAGCGTTGCGCTAAACTTa                           <  2:43608/140‑1 (MQ=255)                                                                                                                                    aatCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCAGATCGACCGTCGCTGCATTCACACTGGCGGCAATCAGCGTGGCCAGGAGCGTTGCGCTAAACTTAATCATCAGGGACATCCTTTTatca   <  2:355765/140‑1 (MQ=255)                                                                                                                                     atCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCAGATCGACCGTCGCTGCATTCACACTGGCGGCAATCAGCGTGGCCAGGAGCGTTGCGCTAAACTTAATCATCAGGGACATCCTTTTatcat  >  1:506920/1‑140 (MQ=255)                                                                                                                                      tCGGAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCAGATCGACCGTCGCTGCATTCACACTGGCGGCAATCAGCGTGGCCAGGAGCGTTGCGCTAAACTTAATCATCAGGGAc               <  2:406164/126‑1 (MQ=255)                                                                                                                                         |                                                                                                                                        GGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTGTTTTTCACTTCATTGCGGGCATCGTTAATCAGGCTTGCCGTACGTACCAGTCCGAATTTTTCCGTGGCGGTGTCTTTGTAATAATCGAAATCCATCATGTTGCTATGCAGATCAGTGGTTTCCATGATACGCAGATCGACCGTCGCTGCATTCACACTGGCGGCAATCAGCGTGGCCAGGAGCGTTGCGCTAAACTTAATCATCAGGGACATCCTTTTATCAT  >  NC_000913/4436314‑4436585

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