Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,176,394	A→T	10.4%	L147Q (CTG→CAG)	cpdA ←	cAMP phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,176,394	0	A	T	10.4%	150.4 / 5.1	55	L147Q (CTG→CAG)	cpdA	cAMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (23/26);  new base T (3/3);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATC  >  NC_000913/3176248‑3176534                                                                                                                                                   |                                                                                                                                             cGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGcc                                                                                                                                           <  1:874219/150‑1 (MQ=255)    gCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACt                                                                                                                                        >  2:1358661/1‑150 (MQ=255)      ggAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCa                                                                                                                                      >  1:508370/1‑150 (MQ=255)      ggAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCAGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCa                                                                                                                                      >  2:139548/1‑150 (MQ=255)         aaCTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGc                                                                                                                                   <  2:1358910/150‑1 (MQ=255)            ttCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGa                                                                                                                                <  2:508370/150‑1 (MQ=255)            ttCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGa                                                                                                                                >  1:517200/1‑150 (MQ=255)             tCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGaa                                                                                                                               >  1:351129/1‑150 (MQ=255)              cGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGaaa                                                                                                                              <  1:1401485/150‑1 (MQ=255)                   gCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCgc                                                                                                                         <  2:1382345/150‑1 (MQ=255)                      cGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCTGCCACTCAAGCTGAAACTCgctca                                                                                                                      >  2:1163082/1‑150 (MQ=255)                       ggTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCgctcag                                                                                                                     <  2:517200/150‑1 (MQ=255)                        gTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCgctcagc                                                                                                                    <  2:888040/150‑1 (MQ=255)                        gTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCgctcagc                                                                                                                    <  2:49219/150‑1 (MQ=255)                          aTCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCgctcagctc                                                                                                                  >  1:486930/1‑150 (MQ=255)                          aTCCAGTTCGCCCGCGTTACGCAGACTGTGTAGATCTAGCAAAAAACAAAAACAATAAAAAAAATAATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCgctcagctc                                                                                                                  >  2:43691/1‑150 (MQ=255)                                  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                                                                                                     <  1:795406/119‑1 (MQ=255)                                                                      tACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACg                                                                                                     >  2:795406/1‑119 (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:880689/1‑150 (MQ=255)                                                                               cAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCTGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGTGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCt                                                             >  2:875462/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:296216/1‑150 (MQ=255)                                                                                                                                     ccAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGt       <  2:1392712/150‑1 (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)                                                                                                                                                   |                                                                                                                                             CGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATC  >  NC_000913/3176248‑3176534

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