Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,176,209	C→A	D209Y (GAT→TAT)	cpdA ←	cAMP phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,176,209	0	C	A	100.0%	116.1 / NA	40	D209Y (GAT→TAT)	cpdA	cAMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (21/19);  total (21/19)

GGTCAGCGTGCCATCAGCA‑TGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCA  >  NC_000913/3176069‑3176346                                                                                                                                              |                                                                                                                                          ggtCAGCGTGCCATCAGCA‑TGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACa                                                                                                                                         >  2:259814/1‑143 (MQ=255) cgtCAGGGCGCCATCAGCACTGT‑ACGCGAGAGGACGCCAGCCTGGCGCGATGGTGTCCAGCGTAAAGTTGGAACTGTCATGGTTAAGTTGCACACAGGTGGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACa                                                                                                                                         <  2:361970/142‑1 (MQ=255)  gtCAGCGTGCCATCAGCA‑TGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAg                                                                                                                                        >  1:256076/1‑143 (MQ=255)   tCAGCGTGCCATCAGCA‑TGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGc                                                                                                                                       >  1:566750/1‑143 (MQ=255)     aGCGTGCCGTCTGCA‑TATAACTCGAGAGTACGCCATCCGGGCGCGATGGTGTCAAGCGATAGGTTGGAACAGGGCGGCTTAAACTGCACACAGGTCGACGGCGTGGCCAGCCGGCGGCGACCACTCCAGTCAAGATACAGCTc                                                                                                                                     >  2:626167/1‑143 (MQ=255)       cGTGCCATCAGCA‑TGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCt                                                                                                                                   >  1:616495/1‑143 (MQ=255)            cATCAGCA‑TGTAACTCGGGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCtgatga                                                                                                                              <  2:611717/143‑1 (MQ=255)              tCAGCA‑TGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCGAGCCGAGATACAGCTGCTGCTGaag                                                                                                                            >  1:346108/1‑142 (MQ=255)                      gTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGGCGAGATACAGCTCCTGATGAATATGACCg                                                                                                                     >  1:889996/1‑143 (MQ=255)                                   gcGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTGAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTAtt                                                                                                        <  2:21268/142‑1 (MQ=255)                                       cAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGGAAGTATTTGAc                                                                                                    >  1:539597/1‑143 (MQ=255)                                             ggCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCgg                                                                                              >  1:550537/1‑143 (MQ=255)                                             ggCGCGACGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCgg                                                                                              <  1:883436/143‑1 (MQ=255)                                               cgcgATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGaa                                                                                            >  1:590048/1‑143 (MQ=255)                                                gcgATGGTATCCAGCGTAACGTTGGCACAGTGCGGCTCAAACTCCACACATGTCCACGGCGTTGCCACCAGCCTCCAACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGaaa                                                                                           <  1:230701/143‑1 (MQ=255)                                                  gATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAAcc                                                                                         >  2:635737/1‑142 (MQ=255)                                                    tcgTATCCAGCGTAAAGCTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTc                                                                                       <  1:768560/141‑1 (MQ=255)                                                    agctATCCAGCGTAAAGTGCGAACAGTGCCTCGTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTc                                                                                       <  1:326960/140‑1 (MQ=255)                                                            aGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCAc                                                                               <  1:407518/143‑1 (MQ=255)                                                                gacaGTTGGAACAGTGCGGCTTAAAGTGCACACAGGTCGACGGCGTTGCGAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTa                                                                           <  2:489226/140‑1 (MQ=255)                                                                  aaGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACCGCTCCCGACGAATGTGACCGCACCGCACGCACCTGACGCGCGGAAACTCCGCCAGCAccgcaca                                                                         >  2:143878/1‑137 (MQ=255)                                                                    tttGGAACAGTGCGGCTTAAACTGCCCACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCTTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCa                                                                       <  1:113616/142‑1 (MQ=255)                                                                       ggAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGtt                                                                    <  2:256076/143‑1 (MQ=255)                                                                        gAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACAGCACAGCACGTATTCGACGTGCGGAAACCTCGACAGCACGGTATCCAGtac                                                                   >  2:379293/1‑141 (MQ=255)                                                                              tGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGCTACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCcgc                                                             >  1:228826/1‑143 (MQ=255)                                                                              tGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGACACTTCGCCCGCCCGGTATCCCGTTCGCCcgc                                                             >  1:234182/1‑143 (MQ=255)                                                                                cGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGt                                                           >  1:604650/1‑143 (MQ=255)                                                                                     tAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCAAGATACAGCCCCTGATGAATATGACCGCACAGCAACTATTTGACGTACGGAAACTTCGCCAGCACGATAACCAGTTCGCCCGCCTAACGc                                                      >  2:316611/1‑143 (MQ=255)                                                                                     tAAACGGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGc                                                      <  1:245471/143‑1 (MQ=255)                                                                                        aCTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGa                                                   <  1:16138/143‑1 (MQ=255)                                                                                             acacgGGTCGTCGGCGTGGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGGCGTGGGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACtgtg                                              <  2:585068/143‑1 (MQ=255)                                                                                                             ttGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTa                              >  1:511991/1‑143 (MQ=255)                                                                                                               gCCAGCAGGCGGCAACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGGTCGCCCGCGTTACGCGGACTGTGTTGATCGCGCCACCTAcc                            >  1:187501/1‑142 (MQ=255)                                                                                                                 cagcagGCGGCTACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAAc                          <  2:159503/143‑1 (MQ=255)                                                                                                                  agcagGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGAATCCAGTTCGCCCGCGTCACGCAGACAGGGCTGCACGAGCCACCTACAAcc                         >  2:862249/1‑143 (MQ=255)                                                                                                                  agcagGCGGCGACCATGCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAAcc                         <  1:448579/143‑1 (MQ=255)                                                                                                                       gcggcgACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAg                    <  2:445809/143‑1 (MQ=255)                                                                                                                              ccATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCgg             <  2:616495/143‑1 (MQ=255)                                                                                                                              ccATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCgg             <  1:292395/143‑1 (MQ=255)                                                                                                                                         gagaTACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGCGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgca  >  1:817536/1‑143 (MQ=255)                                                                                                                                              |                                                                                                                                          GGTCAGCGTGCCATCAGCA‑TGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCA  >  NC_000913/3176069‑3176346

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑