Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,255,685	A→C	11.3%	D72A (GAC→GCC)	folM →	bifunctional dihydrofolate reductase/dihydromonapterin reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,255,685	0	A	C	11.3%	107.0 / 7.0	44	D72A (GAC→GCC)	folM	bifunctional dihydrofolate reductase/dihydromonapterin reductase
Reads supporting (aligned to +/- strand):  ref base A (23/16);  new base C (2/3);  total (25/19)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01

CCAGCCAGCGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGC  >  NC_002947/5255544‑5255814                                                                                                                                               |                                                                                                                                  ccagccagCGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAg                                                                                                                            >  2:73321/1‑149 (MQ=255)       agcGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTTGTCAGTGCCAGCCTCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGT‑GCAGTTTCAGGCCGACCTGTGCCTGGCGGTCGCCCGCCTGgcacag                                                                                                                      <  5:236632/149‑5 (MQ=255)       agcGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGT‑GCAGTTTCAGGCCGACCTGTGCCTGGCGGTCGCCCGCCTGgcacag                                                                                                                      <  5:295897/149‑5 (MQ=255)        gcgcGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTga                                                                                                                     <  3:47158/149‑1 (MQ=255)         cgcgTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTgat                                                                                                                    <  5:74178/149‑1 (MQ=255)                    ggccctggccctCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTg                                                                                                          >  8:176570/1‑148 (MQ=255)                     gccctggccctCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGc                                                                                                         <  7:176570/148‑1 (MQ=255)                             cctCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGtt                                                                                                >  6:301785/1‑149 (MQ=255)                              ctCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCa                                                                                                        >  6:113848/1‑140 (MQ=255)                              ctCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCa                                                                                                        <  5:113848/140‑1 (MQ=255)                                cGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTAt                                                                                              >  7:177707/1‑148 (MQ=255)                                     tGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGa                                                                                        >  8:107575/1‑149 (MQ=255)                                       gctcaGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGAc                                                                                       >  8:321949/4‑148 (MQ=255)                                         acaGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGT‑GCAGTTTCAGGCCGACCTGCGCCTGGCGGTCGCCCGCCTGgcaca                                                                                                                       <  6:223546/114‑4 (MQ=255)                                         acaGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGT‑GCAGTTTCAGGCCGACCTGCGCCTGGCGGTCGCCCGCCTGgcaca                                                                                                                       >  5:223546/1‑111 (MQ=255)                                             gccggccATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCCCCCCAACATCGT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGCCCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACgcgc                                                                                >  3:216473/1‑149 (MQ=255)                                             gccggccATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACgcgc                                                                                >  5:348990/1‑149 (MQ=255)                                                   cATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCcatcat                                                                          <  1:73321/149‑1 (MQ=255)                                                         gtggtCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGT‑GCAGTTTCAGGCCGACCTGTCCCAGGCTGGCGACCGCCAGCCGCTGATCGATCACCTGCACGATCGATAAGACGGACTACGCGCCATCATCCACaa                                                                    >  4:239738/1‑149 (MQ=255)                                                            gtCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGc                                                                 >  2:153087/1‑149 (MQ=255)                                                                  gCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGc                                                            >  5:294328/1‑148 (MQ=255)                                                                   ccAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGGGAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGc                                                            >  1:255887/1‑148 (MQ=255)                                                                      gCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATgg                                                       <  5:301785/149‑1 (MQ=255)                                                                      gCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATgg                                                       <  1:153087/149‑1 (MQ=255)                                                                           agaatacagcagcagCCAGCCCACCCCAACATCGT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTcga                                                  <  4:216473/143‑1 (MQ=255)                                                                               tCCAGCCGCAGGCAGCCCACCCCAACATCGT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCgacgac                                              >  4:346314/1‑149 (MQ=255)                                                                                 ctgCCGCTGGCTTCCCAGCCCATCATCGT‑GCAGTTTCAGGCCGACCTGTTCAAGGCGGGCGACCGCCAGGCGATGATCGATCACCTGCACGATCGTTATTACGGACTACGCGCTATCATCCACAACGCCTCGATATGGCTCGacgacg                                             <  3:239738/146‑1 (MQ=255)                                                                                 cAGCCGCAGGCAGCCCACCCCAACATCGT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGc                                                      >  7:73561/1‑139 (MQ=255)                                                                                 cAGCCGCAGGCAGCCCACCCCAACATCGT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGc                                                      <  8:73561/139‑1 (MQ=255)                                                                                         ggcagCCCACCCCAACATCGT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACCAGCGATATGACGAACAACCCGCCAACATCCACAACGCCTCGCTATGGCTCGCCGACGGCAACGAc                                     >  4:259774/1‑148 (MQ=255)                                                                                               ccACCCCAACATCGT‑GCAGTTTCAGGCCAACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTc                               >  1:203249/1‑148 (MQ=255)                                                                                                cACCCCAACATCGT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCga                             <  8:177707/149‑1 (MQ=255)                                                                                                 aCCCCAACATCGT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCcatc                                                                            <  4:106367/101‑1 (MQ=255)                                                                                                 aCCCCAACATCGT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCcatc                                                                            >  3:106367/1‑101 (MQ=255)                                                                                                   cccAACATCGT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGAc                          >  2:106954/1‑149 (MQ=255)                                                                                                       aCATCGT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATg                      <  7:107575/149‑1 (MQ=255)                                                                                                         aTCGT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCACGACAACCTCGAGACCATGtt                    >  7:31784/1‑149 (MQ=255)                                                                                                            gT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTcc                  <  7:321949/148‑1 (MQ=255)                                                                                                                 aGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTg             <  6:294328/149‑1 (MQ=255)                                                                                                                        ggCCGACCTGTGCCAGGCGGGCGACCGCCAGGCACTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGCCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTc       >  2:106181/1‑148 (MQ=255)                                                                                                                           cGACCTGTGTCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAg   >  1:258238/1‑149 (MQ=255)                                                                                                                            gACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAgc  >  1:303229/1‑149 (MQ=255)                                                                                                                                       ggcgggcgACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGc                <  5:5522/124‑1 (MQ=255)                                                                                                                                       ggcgggcgACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGc                >  6:5522/1‑124 (MQ=255)                                                                                                                                               |                                                                                                                                  CCAGCCAGCGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGT‑GCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGC  >  NC_002947/5255544‑5255814

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