Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,359,353	C→A	12.8%	T106T (ACC→ACA)	modA →	molybdate‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,359,353	0	C	A	12.8%	162.1 / 5.6	52	T106T (ACC→ACA)	modA	molybdate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (23/22);  new base A (3/3);  total (26/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.77e-01

CTACGCGCAGATCAAGAACGGTGCACCGTTCGAAGTGTTCCTCGCCGCTGACGACTCCACCCCCAAGAAGCTGGAAGCAGAAAGGGAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAGGGTGAAGTGCTGAAAAAGAACGAATACAAACACCTGTCCATCGCCAACCCGAAAGCGGCACCTTACGGCCTGGCCGCCACCCAGGTACTG  >  NC_002947/4359227‑4359485                                                                                                                               |                                                                                                                                     cTACGCGCAGATCAAGAACGGTGCACCGTTCGAAGTGTTCCTCGCCGCTGACGACTCCACCCCCAAGAAGCTGGAAGCAGAAAGGGAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCTtg                                                                                                                             <  1:178115/136‑1 (MQ=255)       gcAGATCAAGAACGGTGCACCGTTCGAAGTGTTCCTCGCCGCTGACGACTCCACCCCCAAGAAGCTGGAAGCAGAAAGGGAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCATGTGGTCt                                                                                                                       >  1:59506/1‑136 (MQ=255)            tCAAGAACGGTGCACCGTTCGAAGTGTTCCTCGCCGCTGACGACTCCACCCCCAAGAAGCTGGAAGCAGAAAGGGAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCgaa                                                                                                                  <  2:457066/136‑1 (MQ=255)               aGAACTGTGCAACGTTCGAAGTGTTCCTCGCCGCTGACTACTCCACCACCAAGAAGCTGGAATCATAAAGGGAAACCGTACCTGGATCACTCTTTACCTACGCCATCTGCATTATTTCCTTGTTGTCTTCgaagga                                                                                                               <  2:523104/136‑1 (MQ=255)               aGAACGGTGCACCGTTCGAAGTGTTCCTCGCCGCTGACGACTCCACCCCCAAGAAGCTGGAAGCAGAAAGGGAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCgaagga                                                                                                               <  2:529589/136‑1 (MQ=255)                        cACCGTTCGAAGTGTTCCTCGCCGCTGACGACTCCACCCCCAAGAAGCTGGAAGCAGAAAGGGAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGt                                                                                                      >  2:188940/1‑136 (MQ=255)                        cACCGTTCGAAGTGTTCCTCGCCGCTGACGACTCCACCCCCAAGAAGCTGGAAGCAGAAAGGGAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACACTGGCCTTGTGGTCTGCGAAGGAAGGCTACGt                                                                                                      >  1:320814/1‑136 (MQ=255)                        cACCGTTCGAAGTGTTCCTCGCCGCTGACGACTCCACCCCCAAGAAGCTGGAAGCAGAAAGGGAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACACTGGCCTTGTGGTCTGCGAAGGAAGGCTACGt                                                                                                      >  1:353285/1‑136 (MQ=255)                        cACCGTTCGAAGTGTTCCTCGCCGCTGACGACTCCACCCCCAAGAAGCTGGAAGCAGAAAGGGAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACACTGGCCTTGTGGTCTGCGAAGGAAGGCTACGt                                                                                                      >  1:497663/1‑136 (MQ=255)                           cGTTGGAAGTTTTCCTCTCCTCTGACGACTCCACCCCCAAGAAGCTGGAAGCAGAAAGGGAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGa                                                                                                   <  2:59506/136‑1 (MQ=255)                                gAAGTGTTCCTCGCCGCTGACGACTCCACCCCCAAGAAGCTGGAAGCAGAAAGGGAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTAc                                                                                                        <  2:426464/126‑1 (MQ=255)                                gAAGTGTTCCTCGCCGCTGACGACTCCACCCCCAAGAAGCTGGAAGCAGAAAGGGAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTAc                                                                                                        >  1:426464/1‑126 (MQ=255)                                gAAGTGTTCCTCGCCGCTGACGACTCCACCACCAAGAAGCTGGAAGCAGAAAGGGAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCa                                                                                              >  2:521082/1‑136 (MQ=255)                                gAAGTGTTCCTCGCCGCTGACGACTCCACCACCAAGAAGCTGGAAGCAGAAAGGGAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCa                                                                                              >  2:214549/1‑136 (MQ=255)                                   gtgtTCGTCGCCGCTGACGACTCCACCCCCAAGAAGCTGGAAGCAGAAAGGGAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAgg                                                                                           <  2:545332/136‑1 (MQ=255)                                   gtgtTCCTCGCCGCTGACGACTCCACCCCCAAGAAGCTGGAAGCAGAAAGGGAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAgg                                                                                           >  1:279533/1‑136 (MQ=255)                                   gtgtTCCTCGCCGCTGACGACTCCACCCCCAAGAAGCTGGAAGCAGAAAGGGAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAgg                                                                                           <  2:485189/136‑1 (MQ=255)                                   atgtTCCTCGCCGCTGACGACTCCACCCCCAAGAAGCTGGAAGCAGAAAGGGAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAgg                                                                                           >  1:51861/2‑136 (MQ=255)                                     gtTCCTCGCCGCTGACGACTCCACCCCCAAGAAGCTGGAAGCAGAAAGGGAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAGGGt                                                                                         >  2:347845/1‑136 (MQ=255)                                       cccTCGCCGCTGACGACTCCACCCCCAAGAAGCTGGAAGCAGAAAGGGAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAGGGTGa                                                                                       <  2:87183/135‑1 (MQ=255)                                        ccTCGCCGCTGACGACTCCACCCCCAAGAAGCTGGAAGCAGAAAGGGAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAGGGTGaa                                                                                      <  2:543507/136‑1 (MQ=255)                                        ccTCGCCGCTGACGACTCCACCCCCAAGAAGCTGGAAGCAGAAAGGGAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAGGGTGaa                                                                                      >  2:108334/1‑136 (MQ=255)                                                   acgacTCCACCCCCAAGAAGCTGGAAGCAGAAAGGGAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAGGGTGAAGTGCTGaaaaa                                                                           <  1:347845/136‑1 (MQ=255)                                                           aCCCCCAAGAAGCTGGAAGCAGAAAGGGAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAGGGTGAAGTGCTGAAAAAGAACGAAt                                                                   <  1:87154/136‑1 (MQ=255)                                                              cccAAGAAGCTGGAAGCAGAAAGGGAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGACGGCTACGTAGATGCCAAGGGTGAAGTGCTGAAAAAGAACGAATACa                                                                >  2:293893/1‑136 (MQ=255)                                                              cccAAGAAGCTGGAAGCAGAAAGGGAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCTTGTGGTCGGCGAAGGAAGGCTACGTAGATGCCAAGGGTGAAGTGCTGAAAAAGAACGAATACa                                                                >  2:4106/1‑136 (MQ=255)                                                                            aGCAGAAAGGGAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAGGGTGAAGTGCTGAAAAAGAACGAATACAAACACCTGTCCATc                                                  <  1:188940/136‑1 (MQ=255)                                                                             gCAGAAAGGGAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGCCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAGGGTGAAGTGCTGAAAAAGAACGAATACAAACACCTGTCCATCg                                                 >  1:138254/1‑136 (MQ=255)                                                                                    gggAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCTtgt                                                                                                                            <  1:8737/54‑1 (MQ=255)                                                                                    gggAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCTtgt                                                                                                                            >  2:8737/1‑54 (MQ=255)                                                                                    gggAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAGGGTGAAGTGCTGAAAAAGAACGAATACAAACACCTGTCCATCGCCAAccc                                          >  2:316587/1‑136 (MQ=255)                                                                                         cccGTACCTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAGGGTGAAGTGCTGAAAAAGAACGAATACAAACACCTGTCCATCGCCAACCCGAAAg                                     <  1:82084/135‑1 (MQ=255)                                                                                         aCCGTACCTGTCTCAATCTTTACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGTCTACGTATATTCCAAGTGTGAAGTGCTGAAAAAGAACGAATACAAACACCTGTCCATCGCCAACCCGAAAg                                     <  1:368783/136‑1 (MQ=255)                                                                                         aCCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACACTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAGGGTGAAGTGCTGAAAAAGAACGAATACAAACACCTGTCCATCGCCAACCCGAAAg                                     <  2:320814/136‑1 (MQ=255)                                                                                         aCCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACACTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAGGGTGAAGTGCTGAAAAAGAACGAATACAAACACCTGTCCATCGCCAACCCGAAAg                                     <  2:353285/136‑1 (MQ=255)                                                                                         aCCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACACTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAGGGTGAAGTGCTGAAAAAGAACGAATACAAACACCTGTCCATCGCCAACCCGAAAg                                     <  2:497663/136‑1 (MQ=255)                                                                                          ccGTACCTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAGGGTGAAGTGCTGAAAAAGAACGAATACAAACACCTGTCCATCGCCAACCCGAAAGc                                    <  1:325367/136‑1 (MQ=255)                                                                                                cTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAGGGTGAAGTGCTGAAAAAGAACGAATACAAACACCTGTCCATCGCCAACCCGAAAGCGGCAcc                              >  1:408186/1‑136 (MQ=255)                                                                                                cTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAGGGTGAAGTGCTGAAAAAGAACGAATACAAACACCTGTCCATCGCCAACCCGAAAGCGGCAcc                              >  1:382450/1‑136 (MQ=255)                                                                                                        cGCTTTACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAGGGTGAAGTGCTGAAAAAGAACGAATACAAACACCTGTCCATc                                                  <  1:42601/108‑1 (MQ=255)                                                                                                        cGCTTTACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAGGGTGAAGTGCTGAAAAAGAACGAATACAAACACCTGTCCATc                                                  >  2:42601/1‑108 (MQ=255)                                                                                                            ttACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAGGGTGAAGTGCTGAAAAAGAACGAATACAAACACCTGTCCATCGCCAACCCGAAAGCGGCACCTTACGGCCTGgc                  <  1:48407/136‑1 (MQ=255)                                                                                                            ttACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAGGGTGAAGTGCTGAAAAAGAACGAATACAAACACCTGTCCATCGCCAACCCGAAAGCGGCACCTTACGGCCTGgc                  <  1:156931/136‑1 (MQ=255)                                                                                                            ttACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAGGGTGAAGTGCTGAAAAAGAACGAATACAAACACCTGTCCATCGCCAACCCGAAAGCGGCACCTTACGGCCTGgc                  <  1:316587/136‑1 (MQ=255)                                                                                                                 tACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAGGGTGAAGTGCTGAAAAAGAACGAATACAAACACCTGTCCATCGCCAACCCGAAAGCGGCACCTTACGGCCTGGCCGCCa             <  1:318402/136‑1 (MQ=255)                                                                                                                      cATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAGGGTGAAGTGCTGAAAAAGAACGAATACAAACACCTGTCCATCGCCAACCCGAAAGCGGCACCTTACGGCCTGGCCGCCACCCAg        >  1:506902/1‑136 (MQ=255)                                                                                                                      cATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAGGGTGAAGTGCTGAAAAAGAACGAATACAAACACCTGTCCATCGCCAACCCGAAAGCGGCACCTTACGGCCTGGCCGCCACCCAg        >  1:344338/1‑136 (MQ=255)                                                                                                                      cATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAGGGTGAAGTGCTGAAAAAGAACGAATACAAACACCTGTCCATCGCCAACCCGAAAGCGGCACCTTACGGCCTGGCCGCCACCCAg        >  1:376889/1‑136 (MQ=255)                                                                                                                      cATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAGGGTGAAGTGCTGAAAAAGAACGAATACAAACACCTGTCCATCGCCAACCCGAAAGCGGCACCTTACGGCCTGGCCGCCACCCAg        >  2:258993/1‑136 (MQ=255)                                                                                                                            cACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAGGGTGAAGTGCTGAAAAAGAACGAATACAAACACCTGTCCATCGCCAACCCGAAAGCGGCACCTTACGGCCTGGCCGCCACCCAGGTACTg  >  1:105366/1‑136 (MQ=255)                                                                                                                             aCCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAGGGTGAAGTGCTGAAAAAGAACGAATACAAACACCTGTCCATCGCCAACCCGAAAg                                     >  1:32237/1‑100 (MQ=255)                                                                                                                             aCCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAGGGTGAAGTGCTGAAAAAGAACGAATACAAACACCTGTCCATCGCCAACCCGAAAg                                     <  2:32237/100‑1 (MQ=255)                                                                                                                               |                                                                                                                                     CTACGCGCAGATCAAGAACGGTGCACCGTTCGAAGTGTTCCTCGCCGCTGACGACTCCACCCCCAAGAAGCTGGAAGCAGAAAGGGAAACCGTACCTGGCTCACGCTTTACCTACGCCATCGGCACCCTGGCCTTGTGGTCTGCGAAGGAAGGCTACGTAGATGCCAAGGGTGAAGTGCTGAAAAAGAACGAATACAAACACCTGTCCATCGCCAACCCGAAAGCGGCACCTTACGGCCTGGCCGCCACCCAGGTACTG  >  NC_002947/4359227‑4359485

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