Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,434,695	C→T	G624E (GGG→GAG)	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,434,695	0	C	T	100.0%	112.0 / NA	33	G624E (GGG→GAG)	cpdB	2'3' cyclic nucleotide phosphodiesterase/3' nucleotidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/15);  total (18/15)

AATAAATGCCGGATGCGGTTGATGATTTATAGGCCGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCCCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTATCGCCAGCTATCGGTGCTAAACGCCAGTTGTTATCTGCCGCCGGGTGAATTTCCCCCGCACGTT  >  NC_000913/4434573‑4434827                                                                                                                           |                                                                                                                                     aaTAAATGCCGGATGCGGTTGATGATTTATAGGCCGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCTCTTTCTCTTTAATAGAc                                                                                                                     >  1:825011/1‑140 (MQ=255)   tAAATGCCGGATGCGGTTGATGATTTATAGGCCGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCTCTTTctct                                                                                                                              <  2:190082/129‑1 (MQ=255)   tAAATGCCGGATGCGGTTGATGATTTATAGGCCGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCCGATATCATCGGACGCGACTTTATTCATCGGATACTGCTCTTTctct                                                                                                                              >  1:190082/1‑129 (MQ=255)     aaTGCCGGATGCGGTTGATGATTTATAGGCCGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCTCTTTCTCTTTAATAAACGCTg                                                                                                                 <  2:263048/140‑1 (MQ=255)         ccGGATGCGGTTGATGATTTATAGGCCGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCTCTTTCTCTTTAATAAACGCTGCGGc                                                                                                             >  1:817215/1‑140 (MQ=255)            gATGCGGTTGATGATTTATAGGCCGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCTCTTTCTCTTTAATAAACGCTGCGGCttt                                                                                                          >  1:458682/1‑140 (MQ=255)               gCGGTTGATGATTTATAGGCCGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCTCTTTCTCTTTAATAGACGCTGCGGCTTTATc                                                                                                       <  2:825011/140‑1 (MQ=255)               gCGGTTGATGATTTATAGGCCGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCTCTTTCTCTTTAATAAACGCTGCGGCTTTATc                                                                                                       <  2:756711/140‑1 (MQ=255)                      atgatTTATAGGCCGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCTCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGga                                                                                                >  1:254234/1‑140 (MQ=255)                             ataGGCCGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCTCTTTCTCTTTTATAAACGCTGCGGCTTTATCTGACGGAGAGGttt                                                                                         >  2:41140/1‑140 (MQ=255)                               aGGCCGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCTCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCg                                                                                       >  2:388355/1‑140 (MQ=255)                                  ccGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCTCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGaaa                                                                                    <  1:871680/140‑1 (MQ=255)                                     aaaaaGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCTCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACgg                                                                                 <  2:334231/140‑1 (MQ=255)                                              tgtTTTACTTACTCAAATCCACCTGATAAATCGCAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCTCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGt                                                                        <  1:671823/140‑1 (MQ=255)                                               gtTTTACTTACTCAAATCCACCTGATAAATCGCAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCTCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGtt                                                                       >  1:441282/1‑140 (MQ=255)                                                     cttactCAAATCCACCTGATAAATCGCAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCTCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTAt                                                                 >  2:1003248/1‑140 (MQ=255)                                                       tactCAAATCCACCTGATAAATCGCAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCTCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTATCg                                                               <  2:868597/140‑1 (MQ=255)                                                              aTCCACCTGATAAATCGCAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCTCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTATCGCCAGCTa                                                        >  2:488070/1‑140 (MQ=255)                                                               tCCACCTGATAAATCGCAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCTCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTATCGCCAGCTAt                                                       >  2:43905/1‑140 (MQ=255)                                                                   ccTGATAAATCGCAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCTCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTATCGCCAGCTATCGGt                                                   >  1:628932/1‑140 (MQ=255)                                                                   ccTGATAAATCGCAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCTCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTATCGCCAGCTATCGGt                                                   >  2:84128/1‑140 (MQ=255)                                                                      gATAAATCGCAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCTCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTATCGCCAGCTATCGGTGCt                                                >  1:140078/1‑140 (MQ=255)                                                                            tCGCAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCTCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTATCGCCAGCTATCGGTGCTAAACGc                                          <  1:84128/140‑1 (MQ=255)                                                                             cGCAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCTCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTATCGCCAGCTATCGGTGCTAAACGcc                                         <  2:384408/140‑1 (MQ=255)                                                                             cGCAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCTCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTATCGCCAGCTATCGGTGCTAAACGcc                                         <  2:254234/140‑1 (MQ=255)                                                                             cGCAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCTCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTATCGCCAGCTATCGGTGCTAAACGc                                          >  1:143237/1‑139 (MQ=255)                                                                             cGCAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCTCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTATCGCCAGCTATCGGTGCTAAACGc                                          <  2:143237/139‑1 (MQ=255)                                                                               cAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCTCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTATCGCCAGCTATCGGTGCTAAACGCCAg                                       <  1:168384/140‑1 (MQ=255)                                                                                 aaCCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCTCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTATCGCCAGCTATCGGTGCTAAACGCCAgtt                                     >  2:749804/1‑140 (MQ=255)                                                                                             tcGGTCGCGACTTTATTCATCGGATACTGCTCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTGTCTTATCGCCAGCTATCGGTGCTAAACGCCAGTTGTTATCTgccgt                         <  2:441282/140‑2 (MQ=255)                                                                                                           aTTCATCGGATACTGCTCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTATCGCCAGCTATCGGTGCTAAACGCCAGTTGTTATCTGCCGCCGGGTGAATTTccc           >  1:829147/1‑140 (MQ=255)                                                                                                           aTTCATCGGATACTGCTCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTATCGCCAGCTATCGGTGCTAAACGCCAGTTGTTATCTGCCGCCGGGTGAATTTccc           >  2:984379/1‑140 (MQ=255)                                                                                                                tCGGATACTGCTCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTATCGCCAGCTATCGGTGCTAAACGCCAGTTGTTATCTGCCGCCGGGTGAATTTCCCCCGCa      <  1:43905/140‑1 (MQ=255)                                                                                                                    aTACTGCTCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTATCGCCAGCTATCGGTGCTAAACGCCAGTTGTTATCTGCCGCCGGGTGAATTTCCCCCGCACGtt  >  2:770530/1‑140 (MQ=255)                                                                                                                           |                                                                                                                                     AATAAATGCCGGATGCGGTTGATGATTTATAGGCCGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCCCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTATCGCCAGCTATCGGTGCTAAACGCCAGTTGTTATCTGCCGCCGGGTGAATTTCCCCCGCACGTT  >  NC_000913/4434573‑4434827

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