Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,387,191	C→G	12.1%	intergenic (‑88/‑50)	PP_2094 ← / → pyrD	nitrate‑binding protein NasS/dihydroorotate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,387,191	0	C	G	12.1%	83.1 / 4.7	33	intergenic (‑88/‑50)	PP_2094/pyrD	nitrate‑binding protein NasS/dihydroorotate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (12/17);  new base G (2/2);  total (14/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

CGCTGCCGTTGACCCAGGCCAGGGGCGTTGCTCGGGGTACAGTGTTCACAAGAGTCCTCGCCCTTACTGATAAGGCTGTGGGCTGAACGGGTGCAACCGATGTGCCAAGCGCTCCACGATTTGCTGTCCTGCCAGCCCCGCGCTGGCTATAATCGCCCCCTCACTCACCCTGAGCCTTGCCCGCCCATGTATACCCTGGCCCGCCAGCTGCTGTTCAAGCTTTCCCCGGAAACTTCCCACGACCTGTCCCTGGACCTGATCGGCGCCGGTGGCCGTTTGG  >  NC_002947/2387055‑2387334                                                                                                                                         |                                                                                                                                                cGCTGCCGTTGACCCAGGCCAGGGGCGTTGCTCGGGGTACAGTGTTCACAAGAGTCCTCGCCCTTACTGATAAGGCTGTGGGCTGAACGGGTGCAACCGATGTGCCAAGCGCTCCACGATTTGCTGTCCTGCCAGCGCGGGGCTGGcag                                                                                                                                     <  7:31214/149‑3 (MQ=255)      ccGTTGACCCAGGCCAGGGGCGTTGCTCGGGGTACAGTGTTCACAAGAGTCCTCGCCCTTACTGATAAGGCTGTGGGCTGAACGGGTGCAACCGATGTGCCAAGCGCTCCACGATTTGCTGTCCTGCCAGCCCCGCGCTGGCTATAATc                                                                                                                                <  6:162303/149‑1 (MQ=255)        gTTGACCCAGGCCAGGGGCGTTGCTCGGGGTACAGTGTTCACAAGAGTCCTCGCCCTTACTGATAAGGCTGTGGGCTGAACGGGTGCAACCGATGTGCCAAGCGCTCCACGATTTGCTGTCCTGCCAGCCCCGCGCTGGCTATAATCGc                                                                                                                              <  5:119685/149‑1 (MQ=255)            aCCCAGGCCAGGGGCGTTGCTCGGGGTACAGTGTTCACAAGAGTCCTCGCCCTTACTGATAAGGCTGTGGGCTGAACGGGTGCAACCGATGTGCCAAGCGCTCCACGATTTGCTGTCCTGCCAGCGCGGGGCTGGcagga                                                                                                                                   <  6:183457/140‑5 (MQ=255)            aCCCAGGCCAGGGGCGTTGCTCGGGGTACAGTGTTCACAAGAGTCCTCGCCCTTACTGATAAGGCTGTGGGCTGAACGGGTGCAACCGATGTGCCAAGCGCTCCACGATTTGCTGTCCTGCCAGCGCGGGGCTGGcagga                                                                                                                                   >  5:183457/1‑136 (MQ=255)             cccAGGCCAGGGGCGTTGCTCGGGGTACAGTGTTCACAAGAGTCCTCGCCCTTACTGATAAGGCTGTGGGCTGAACGGGTGCAACCGATGTGCCAAGCGCTCCACGATTTGCTGTCCTGCCAGCCCCGCGCTGGCTATAATCGCCCCct                                                                                                                         >  1:89535/1‑149 (MQ=255)               caggccaggGGCGTTGCTCGGTGTACAGTGTTCACAAGAGTCCTCGCCCTTACTGATAAGGCTGTTGGCTGAACGGTTGCAACCGATGTGCCAAGCGCTCCACGATTTGATGTCCTGCCAGCCCCGCGCTGGCTATAATCGCCCCctca                                                                                                                       <  5:186449/149‑1 (MQ=255)                aggccaggGGCGTTGCTCGGGGTACAGTGTTCACAAGAGTCCTCGCCCTTACTGATAAGGCTGTGGGCTGAACGGGTGCAACCGATGTGCCAAGCGCTCCACGATTTGCTGTCCTGCCAGCCCCGCGCTGGCTATAATCGCCCCctcac                                                                                                                      <  5:285927/149‑1 (MQ=255)                               cTCGGGGTACAGTGTTCACAAGAGTCCTCACCCTTACTGATAAGGCTGTGGGCTGAACGGGTGCAACCGATGTGCCAAGCGCTCCACGATTTGCTGTCCTGCCAGCCCCGCGCTGGCTATAATCGCCCCCTCACTCACCCTGAGCCTTg                                                                                                       <  8:190273/149‑1 (MQ=255)                                  ggggTACAGTGTTCACAAGAGTCCTCGCCCTTACTGATAAGGCTGTGGGCTGAACGGGTGCAACCGATGTGCCAAGCGCTCCACGATTTGCTGTCCTGCCAGCCCCGCGCTGGCTATAATCGCCCCCTCACTCACCCTGAGCCTTgccc                                                                                                    <  1:315479/149‑1 (MQ=255)                                  ggggTACAGTGTTCACAAGAGTCCTCGCCCTTACTGATAAGGCTGTGGGCTGAACGGGTGCAACCGATGTGCCAAGCGCTCCACGATTTGCTGTCCTGCCAGCCCCGCGCTGGCTATAATCGCCCCCTCACTCACCCTGAGCCTTgccc                                                                                                    >  1:99882/1‑149 (MQ=255)                                    ggTACAGTGTTCACAAGAGTCCTCGCCCTTACTGATAAGGCTGTGGGCTGAACGGGTGCAACCGATGTGCCAAGCGCTCCACGATTTGCTGTCCTGCCAGCCCcgcg                                                                                                                                            <  6:99557/107‑1 (MQ=255)                                    ggTACAGTGTTCACAAGAGTCCTCGCCCTTACTGATAAGGCTGTGGGCTGAACGGGTGCAACCGATGTGCCAAGCGCTCCACGATTTGCTGTCCTGCCAGCCCcgcg                                                                                                                                            >  5:99557/1‑107 (MQ=255)                                                   agagTCCTCGCCCTTACTGATAAGGCTGTGGGCTGAACGGGTGCAACCGATGTGCCAAGCGCTCCACGATTTGCTGTCCTGCCAGCCCCGCGCTGGCTATAATCGCCCCCTCACTCACCCTGAGCCTTGCCCGCCCATGTATACCCTgg                                                                                   <  1:126682/149‑1 (MQ=255)                                                                ttACTGATAAGGCTGTGGGCTGAACGGGTGCAACCGATGTGCCAAGCGCTCCACGATTTGCTGTCCTGCCAGCCCCGCGCTGGCTATAATCGCCCCCTCACTCACCCTGAGCCTTGCCCGCCCATGTATACCCTGGCCCGCCAgctgct                                                                      >  2:295865/1‑149 (MQ=255)                                                                               tgGGCTGAACGGGTGCAACCGATGTGCCAAGCGCTCCACGATTTGCTGTCCTGCCAGCCCCGCGCTGGCTATAATCGCCCCCTCACTCACCCTGAGCCTTGCCCGCCCATGTATACCCTGGCCCGCCAGCTGCTGTTCAAGCTTTcccc                                                       <  2:89535/149‑1 (MQ=255)                                                                                               aaCCGATGTGCCAAGCGCTCCACGATTTGCTGTCCTGCCAGCCCCGCGCTGGCTATAATCGCCCCCTCACTCACCCTGAGCCTTGCCCGCCCATGTATACCCTGGCCCGCCAGCTGCTGTTCAAGCTTTCCCCGGAAACTTCCCACGAc                                       >  2:12351/1‑149 (MQ=255)                                                                                                aCCGATGTGCCAAGCGCTCCACGATTTGCTGTCCTGCCAGCCCCGCGCTGGCTATAATCGCCCCCTCACTCACCCTGAGCCTTGCCCGCCCATGTATACCCTGGCCCGCCAGCTGCTGTTCAAGCTTTCCCCGGAAACTTCCCACGAcc                                      >  4:199570/1‑149 (MQ=255)                                                                                                   gATGTGCCAAGCGCTCCACGATTTGCTGTCCTGCCAGCCCCGCGCTGGCTATAATCGCCCCCTCACTCACCCTGAGCCTTGCCCGCCCATGTATACCCTGGCCCGCCAGCTGCTGTTCAAGCTTTCCCCGGAAAc                                                 >  5:310375/1‑135 (MQ=255)                                                                                                   gATGTGCCAAGCGCTCCACGATTTGCTGTCCTGCCAGCCCCGCGCTGGCTATAATCGCCCCCTCACTCACCCTGAGCCTTGCCCGCCCATGTATACCCTGGCCCGCCAGCTGCTGTTCAAGCTTTCCCCGGAAAc                                                 <  6:310375/135‑1 (MQ=255)                                                                                                       tgCCAAGCGCTCCACGTTTTGCAGTCCAGCCAGCCCCGCGCTGGCTATAATCGCCCCCTCACTCACGCTGAGCCTTGCCCGCGCATGAATACTCTGGCCCGCCAGCTGCTGTTCATGCTTTCCCCGGAAACTTCCCACGACCTGTCCCt                               <  4:219953/149‑1 (MQ=255)                                                                                                        gCCAAGCGCTCCACGATTTGCTGTCCTGCCAGCCCCGCGCTGGCTATAATCGCCCCCTCACTCCCCCTGAGCCTTGCCCGCCCATGTATACCCTGGCCCGCCAGCTGCTGTTCAAGCTTTCCCCGGAAACTTCCCACGACCTGTCCCTg                              >  6:10104/1‑149 (MQ=255)                                                                                                        gCCAAGCGCTCCACGATTTGCTGTCCTGCCAGCCCCGCGCTGGCTATAATCGCCCCCTCACTCACCCTGAGCCTTGCCCGCCCATGTATACCCTGGCCCGCCAGCTGCTGTTCAAGCTTTCCCCGGAAACTTCCCACGACCTGTCCCTg                              <  6:102607/149‑1 (MQ=255)                                                                                                            aGCGCTCCACGATTTGCTGTCCTGCCAGCCCCGCGCTGGCTATAATCGCCCCCTCACTCACCCTGAGCCTTGCCCGCCCATGTATACCCTGGCCCGCCAGCTGCTGTTCAAGCTTTCCCCGGAAACTTCCCACGACCTGTCCCTGGAcc                          <  3:199570/149‑1 (MQ=255)                                                                                                                    aCGATTTGCTGTCCTGCCAGCCCCGCGCTGGCTATAATCGCCCCCTCACTCACCCTGAGCCTTGCCCGCCCATGTATACCCTGGCCCGCCAGCTGCTGTTCAAGCTTTCCCCGGAAACTTCCCACGACCTGTCCCTGGACCTGATCGgc                  >  1:11541/1‑149 (MQ=255)                                                                                                                        tttGCTGTCCTGCCAGCCCCGCGCTGGCTATAATCGCCCCCTCACTCACCCTGAGCCTTGCCCGCCCATGTATAc                                                                                        >  1:292318/1‑75 (MQ=255)                                                                                                                        tttGCTGTCCTGCCAGCCCCGCGCTGGCTATAATCGCCCCCTCACTCACCCTGAGCCTTGCCCGCCCATGTATAc                                                                                        <  2:292318/75‑1 (MQ=255)                                                                                                                                 tagCCAGCGCGGGGCTGGCTATAATCGCCCCCTCACTCACCCTGAGCCTTGCCCGCCCATGTATACCCTGGCCCGCCAGCTGCTGTTCAAGCTTTCCCCGGAAACTTCCCACGACCTGTCCCTGGACCTGATCGGCGCCGGTGGCCGtt     >  1:87708/3‑149 (MQ=255)                                                                                                                                    ccAGCCCCGCGCTGGCTATAATCGCCCCCTCACTCCCACAGAGGCTTGCCCGCCCAAGTATACCCGGGCCCGCCAGCGGCTGTTCAAGCTTTCCCCGGAAACTACCAACGACCTGTCCCAGGACCTGATCGGAGCCGGTGGCCGTTTgg  >  8:102802/1‑149 (MQ=255)                                                                                                                                    ccAGCCCCGCGCTGGCTATAATCGCCCCCTCACTCACCCTGAGCCTTGCCCGCCCATGTATACCCTGGCCCGCCAGCTGCTGTTCAAGCTTTCCCCGGAAACTTCCCACGACCTGTCCCTGGACCTGATCGGCGCCGGTGGCCGTTTg   <  1:295865/148‑1 (MQ=255)                                                                                                                                    ccAGCCCCGCGCTGGCTATAATCGCCCCCTCACTCACCCTGAGCCTTGCCCGCCCATGTATACCCTGGCCCGCCAGCTGCTGTTCAAGCTTTCCCCGGAAACTTCCCACGACCTGTCCCTGGACCTGATCGGCGCCGGGGGCCGTTTgg  <  7:102802/149‑1 (MQ=255)                                                                                                                                       gCCCCGCGCTGGCTATAATCGCCCCCTCACTCACCCTGAGCCTTGCCCGCCCATGTATACCCTGGCCCGCCAGCTGCTGTTCAAg                                                               >  6:263968/1‑85 (MQ=255)                                                                                                                                       gCCCCGCGCTGGCTATAATCGCCCCCTCACTCACCCTGAGCCTTGCCCGCCCATGTATACCCTGGCCCGCCAGCTGCTGTTCAAg                                                               <  5:263968/85‑1 (MQ=255)                                                                                                                                         |                                                                                                                                                CGCTGCCGTTGACCCAGGCCAGGGGCGTTGCTCGGGGTACAGTGTTCACAAGAGTCCTCGCCCTTACTGATAAGGCTGTGGGCTGAACGGGTGCAACCGATGTGCCAAGCGCTCCACGATTTGCTGTCCTGCCAGCCCCGCGCTGGCTATAATCGCCCCCTCACTCACCCTGAGCCTTGCCCGCCCATGTATACCCTGGCCCGCCAGCTGCTGTTCAAGCTTTCCCCGGAAACTTCCCACGACCTGTCCCTGGACCTGATCGGCGCCGGTGGCCGTTTGG  >  NC_002947/2387055‑2387334

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