Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,176,112	G→A	P241L (CCC→CTC)	cpdA ←	3',5' cAMP phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,176,112	0	G	A	100.0%	137.6 / NA	42	P241L (CCC→CTC)	cpdA	3',5' cAMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/21);  total (21/21)

CTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGG  >  NC_000913/3175972‑3176254                                                                                                                                             |                                                                                                                                               cTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATg                                                                                                                                        <  5:74403/149‑1 (MQ=255)  tGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATgg                                                                                                                                       >  1:29223/1‑149 (MQ=255)               taaataaaGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTaaa                                                                                                                          <  3:3650/149‑1 (MQ=255)                 aataaaGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGt                                                                                                                        >  1:176155/1‑149 (MQ=255)                 aataaaGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGt                                                                                                                        >  2:128257/1‑149 (MQ=255)                  ataaaGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGtt                                                                                                                       >  2:178448/1‑149 (MQ=255)                    aaaGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGAt                                                                                                                                         <  3:123861/129‑1 (MQ=255)                    aaaGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGAt                                                                                                                                         >  4:123861/1‑129 (MQ=255)                      agaagCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGaa                                                                                                                   >  6:18393/1‑149 (MQ=255)                            gTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGc                                                                                                             >  8:32900/1‑149 (MQ=255)                                     aGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTaaa                                                                                                     <  7:101719/148‑1 (MQ=255)                                     aGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTaaa                                                                                                     <  1:128257/148‑1 (MQ=255)                                      gTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACt                                                                                                   >  7:41469/1‑149 (MQ=255)                                         gCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGca                                                                                                <  2:29223/149‑1 (MQ=255)                                            ttCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGcaca                                                                                              >  7:60278/1‑148 (MQ=255)                                               tGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAgg                                                                                          >  8:110811/1‑149 (MQ=255)                                                 caGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTc                                                                                        >  7:33016/2‑149 (MQ=255)                                                                  aaaCGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCCGCCGAGCGCGATGGTATCCAGCGTAAAGTGGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCcagcag                                                                       <  3:19599/149‑1 (MQ=255)                                                                                cAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTc                                                          <  1:95439/148‑1 (MQ=255)                                                                                cAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTc                                                          >  2:95439/1‑148 (MQ=255)                                                                                      aTGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCga                                                   <  7:32900/149‑1 (MQ=255)                                                                                      aTGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCga                                                   <  8:41469/149‑1 (MQ=255)                                                                                         cACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTc                                                     >  4:24550/1‑144 (MQ=255)                                                                                           ccTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATcc                                              <  3:24550/149‑1 (MQ=255)                                                                                               ggtggtCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAAc                                                                                                    <  7:132353/91‑1 (MQ=255)                                                                                               ggtggtCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAAc                                                                                                    >  8:132353/1‑91 (MQ=255)                                                                                                 tggtCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTcc                                        <  2:176155/149‑1 (MQ=255)                                                                                                      aGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCtgat                                    <  1:178448/148‑1 (MQ=255)                                                                                                              aTCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGAc                           >  7:8690/1‑149 (MQ=255)                                                                                                              aTCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGAc                           <  5:18393/149‑1 (MQ=255)                                                                                                              aTCAGCATGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGAc                           <  8:60278/149‑1 (MQ=255)                                                                                                                     tGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCtgatga                                  <  2:61346/135‑1 (MQ=255)                                                                                                                     tGTAACTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCtgatga                                  >  1:61346/1‑135 (MQ=255)                                                                                                                         aCTCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCaa                 <  8:8690/148‑1 (MQ=255)                                                                                                                           tCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTa              >  7:85940/1‑149 (MQ=255)                                                                                                                           tCGAGAGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTa              >  3:48366/1‑149 (MQ=255)                                                                                                                             gagaGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTAtt            >  2:64927/1‑149 (MQ=255)                                                                                                                             gagaGTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTAt             <  8:85940/148‑1 (MQ=255)                                                                                                                              agagTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTAttt           <  4:48366/149‑1 (MQ=255)                                                                                                                               gagTACGCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTAttt           >  7:2678/1‑148 (MQ=255)                                                                                                                                     gCCAGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGc    <  8:28164/149‑1 (MQ=255)                                                                                                                                        aGCCGAGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCgg  >  3:1973/1‑148 (MQ=255)                                                                                                                                             |                                                                                                                                               CTGTTGAAACCGTGTAAATAAAGAAGCGTAGACATCAGTAGCCTTCTGAAGCGGTATCAGGTTGGAAACGTGTGTCCGCCAGGCGATGCACCTCGGTGGTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGG  >  NC_000913/3175972‑3176254

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