Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,176,429	C→T	8.0%	P135P (CCG→CCA)	cpdA ←	cAMP phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,176,429	0	C	T	8.0%	149.2 / 2.9	50	P135P (CCG→CCA)	cpdA	cAMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (22/24);  new base T (2/2);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

CGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGC  >  NC_000913/3176281‑3176578                                                                                                                                                     |                                                                                                                                                      cGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCg                                                                                                                                                      <  2:341124/150‑1 (MQ=255)    ccGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCa                                                                                                                                                   <  1:1161374/150‑1 (MQ=255)        gTTACACAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGcc                                                                                                                                               <  2:1211511/150‑1 (MQ=255)          tACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCATCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCaa                                                                                                                                             >  1:993363/1‑150 (MQ=255)                          aTCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAg                                                                                                                             <  2:903888/150‑1 (MQ=255)                           tCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGc                                                                                                                            <  2:351129/150‑1 (MQ=255)                            cGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCa                                                                                                                           >  1:765264/1‑150 (MQ=255)                                ccAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAg                                                                                                                       <  1:348176/150‑1 (MQ=255)                                     tACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACg                                                                                                                                                 >  2:795406/1‑119 (MQ=255)                                     tACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACg                                                                                                                                                 <  1:795406/119‑1 (MQ=255)                                     tACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGAttt                                                                                                                  >  1:31198/1‑150 (MQ=255)                                      aCAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTg                                                                                                                 >  1:281425/1‑150 (MQ=255)                                             gCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCTGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGc                                                                                                          <  1:1163082/150‑1 (MQ=255)                                              cAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCTGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGTGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCt                                                                                                         >  2:875462/1‑150 (MQ=255)                                              cAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCTGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGTGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCt                                                                                                         >  2:880689/1‑150 (MQ=255)                                                ggTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCa                                                                                                       <  1:1248162/150‑1 (MQ=255)                                                ggTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCGGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCa                                                                                                       <  2:991585/150‑1 (MQ=255)                                                         tgatgatgCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGGACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAac                                                                                              <  1:99351/150‑1 (MQ=255)                                                          gatgatgCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAaca                                                                                             <  2:31198/150‑1 (MQ=255)                                                             gatgCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACAcgc                                                                                          >  2:260679/1‑150 (MQ=255)                                                              atgCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACAcgcg                                                                                         <  1:699469/150‑1 (MQ=255)                                                                 cagcagcagcaACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCtt                                                                                      >  2:1019234/1‑150 (MQ=255)                                                                  agcagcagcaACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTc                                                                                     >  2:104115/1‑150 (MQ=255)                                                                       cagcaACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCgg                                                                                <  1:1358661/150‑1 (MQ=255)                                                                         gcaACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCgggg                                                                              <  2:281425/150‑1 (MQ=255)                                                                           aaCGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGgag                                                                            >  2:806695/1‑150 (MQ=255)                                                                                aTGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATAcc                                                                       <  1:474355/150‑1 (MQ=255)                                                                                 tGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATAccc                                                                      >  1:1392712/1‑150 (MQ=255)                                                                                          ggCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGt                                                             >  1:99162/1‑150 (MQ=255)                                                                                           gcgcATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTa                                                            <  1:139548/150‑1 (MQ=255)                                                                                                 cagacAGTTTACGTTCCAAACACTCAAGCTGAAACTAGCTCAGCTAACAGTGCGTAACGAAAAAAACCTCGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACgcgc                                                      >  2:517023/5‑150 (MQ=255)                                                                                                   gCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTg                                                    <  1:260679/150‑1 (MQ=255)                                                                                                   gCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTg                                                    >  1:228444/1‑150 (MQ=255)                                                                                                    ccAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGt                                                   <  2:1392712/150‑1 (MQ=255)                                                                                                    ccAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGt                                                   >  2:296216/1‑150 (MQ=255)                                                                                                     cAGTTTACGTTCCTGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGTGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTa                                                  <  1:880689/150‑1 (MQ=255)                                                                                                     cAGTTTACGTTCCTGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGTGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTa                                                  <  1:875462/150‑1 (MQ=255)                                                                                                     cAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTa                                                  >  2:486239/1‑150 (MQ=255)                                                                                                         ttACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATc                                              >  2:275718/1‑150 (MQ=255)                                                                                                                     cACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGga                                                                             <  2:1418969/107‑1 (MQ=255)                                                                                                                     cACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGga                                                                             >  1:1418969/1‑107 (MQ=255)                                                                                                                       cTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATc                                >  1:1377582/1‑150 (MQ=255)                                                                                                                          aaGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTg                             <  2:993363/150‑1 (MQ=255)                                                                                                                           aGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTgg                            <  1:963695/150‑1 (MQ=255)                                                                                                                                 aaCTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGccc                      >  2:845602/1‑150 (MQ=255)                                                                                                                                   cTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCgg                    >  1:994057/1‑150 (MQ=255)                                                                                                                                         cagctcaCCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCTGcgtctg              >  1:1195548/1‑145 (MQ=255)                                                                                                                                         cagctcaCCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCa              >  1:1088723/1‑150 (MQ=255)                                                                                                                                             tcaCCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACg          <  1:104115/150‑1 (MQ=255)                                                                                                                                                  gTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAggg     <  1:275718/150‑1 (MQ=255)                                                                                                                                                     cGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGgcgc  >  1:965494/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                      CGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGC  >  NC_000913/3176281‑3176578

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