Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,176,721	G→A	T38I (ACC→ATC)	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,721	0	G	A	100.0%	178.2 / NA	55	T38I (ACC→ATC)	cpdA	3',5' cAMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/29);  total (26/29)

CTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGAC  >  NC_000913/3176586‑3176856                                                                                                                                        |                                                                                                                                        cTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGatgt                                                                                                                                      <  2:864172/139‑1 (MQ=255)      gATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTAcc                                                                                                                                 <  2:657603/139‑1 (MQ=255)              cAGCGAAATGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAg                                                                                                                         <  1:32903/139‑1 (MQ=255)                     aTGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCg                                                                                                                  <  2:824218/139‑1 (MQ=255)                           aTAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCttt                                                                                                            >  2:746750/1‑139 (MQ=255)                                ccGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACa                                                                                                                          >  2:440733/1‑120 (MQ=255)                                ccGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACa                                                                                                                          <  1:440733/120‑1 (MQ=255)                                ccGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGc                                                                                                       >  2:112297/1‑139 (MQ=255)                                ccGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGc                                                                                                       >  1:274970/1‑139 (MQ=255)                                 cGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCa                                                                                                      >  1:1321492/1‑139 (MQ=255)                                  gCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCaa                                                                                                     >  1:791590/1‑139 (MQ=255)                                  gCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCaa                                                                                                     >  1:1407690/1‑139 (MQ=255)                                        gATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGt                                                                                               >  1:414285/1‑139 (MQ=255)                                         aTTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTg                                                                                              <  1:746750/139‑1 (MQ=255)                                                  gcgcTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAAc                                                                                                                           >  2:1198452/1‑101 (MQ=255)                                                  gcgcTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAAc                                                                                                                           <  1:1198452/101‑1 (MQ=255)                                                  gcgcTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCa                                                                                                      <  2:1348195/122‑1 (MQ=255)                                                  gcgcTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCa                                                                                                      >  1:1348195/1‑122 (MQ=255)                                                  gcgcTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGa                                                                                             <  1:325118/131‑1 (MQ=255)                                                  gcgcTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGa                                                                                             >  2:325118/1‑131 (MQ=255)                                                  gcgcTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGt                                                                                     <  2:370849/139‑1 (MQ=255)                                                  gcgcTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGt                                                                                     >  1:1091689/1‑139 (MQ=255)                                                   cgcTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTtgt                                                                                                         >  1:1426104/1‑118 (MQ=255)                                                   cgcTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTtgt                                                                                                         <  2:1426104/118‑1 (MQ=255)                                                   cgcTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAgg                                                                                                <  1:1511399/127‑1 (MQ=255)                                                   cgcTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAgg                                                                                                >  2:1511399/1‑127 (MQ=255)                                                   cgcTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTa                                                                                    <  1:563642/139‑1 (MQ=255)                                                     cTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAAt                                                                                  <  2:1321492/139‑1 (MQ=255)                                                          tCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGCTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTa                                                                             <  2:972018/139‑1 (MQ=255)                                                             ccTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTaaaa                                                                          >  1:489861/1‑139 (MQ=255)                                                                gTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATcc                                                                       <  2:918653/139‑1 (MQ=255)                                                                  cgcgACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTg                                                                                              <  2:277525/114‑1 (MQ=255)                                                                  cgcgACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTg                                                                                              >  1:277525/1‑114 (MQ=255)                                                                  cgcgACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTg                                                                     >  1:949030/1‑139 (MQ=255)                                                                      aCAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGActc                                                                 >  2:485945/1‑139 (MQ=255)                                                                           cAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGcc                                                            >  2:1245299/1‑139 (MQ=255)                                                                                cGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCAcc                                                       <  2:1091689/139‑1 (MQ=255)                                                                                cGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCAcc                                                       <  2:791590/139‑1 (MQ=255)                                                                                   aTTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCAccagc                                                    >  2:1012183/1‑139 (MQ=255)                                                                                   aTTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCAccagc                                                    <  1:1273202/139‑1 (MQ=255)                                                                                   aTTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCAccagc                                                    <  1:1245299/139‑1 (MQ=255)                                                                                        tGCGGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCATTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCagag                                               >  1:1035171/1‑139 (MQ=255)                                                                                            ggTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAg                                           >  2:1382385/1‑139 (MQ=255)                                                                                                       aTCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGc                                <  1:112297/139‑1 (MQ=255)                                                                                                          gCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCAc                                                        <  1:317536/112‑1 (MQ=255)                                                                                                          gCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCAc                                                        >  2:317536/1‑112 (MQ=255)                                                                                                                         tGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTtgt                                                                                                         >  2:1297807/1‑48 (MQ=255)                                                                                                                         tGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTtgt                                                                                                         <  1:1297807/48‑1 (MQ=255)                                                                                                                                 ctcCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCa      <  2:949030/139‑1 (MQ=255)                                                                                                                                 ctcCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCa      >  2:713596/1‑139 (MQ=255)                                                                                                                                  tcCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGgt                    >  2:247904/1‑124 (MQ=255)                                                                                                                                  tcCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGgt                    <  1:247904/124‑1 (MQ=255)                                                                                                                                  tcCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCac     <  2:489861/139‑1 (MQ=255)                                                                                                                                    ccAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCacga   <  2:414285/139‑1 (MQ=255)                                                                                                                                     cAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCacgac  <  2:1407690/139‑1 (MQ=255)                                                                                                                                     cAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCacgac  <  2:274970/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                        CTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGAC  >  NC_000913/3176586‑3176856

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