Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,255,730	G→T	7.4%	G71V (GGC→GTC)	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,730	0	G	T	7.4%	160.3 / 2.7	54	G71V (GGC→GTC)	folM	bifunctional dihydrofolate reductase/dihydromonapterin reductase
Reads supporting (aligned to +/- strand):  ref base G (23/27);  new base T (2/2);  total (25/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 = 6.99e-01

GTTACCGGAGCCAGCCAGCGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAA  >  NC_002947/5255582‑5255877                                                                                                                                                      |                                                                                                                                                    gTTACCGGAGCCAGCCAGCGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAggcgggc                                                                                                                                                    >  2:1101190/1‑150 (MQ=255)       ggAGCCAGCCAGCGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGc                                                                                                                                              >  2:509476/1‑150 (MQ=255)       ggAGCCAGCCAGCGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGc                                                                                                                                              <  1:53160/150‑1 (MQ=255)         agccagccAGCGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCa                                                                                                                                            <  1:718372/150‑1 (MQ=255)          gccagccaGCGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCTGGCGGTcgcccgcctg                                                                                                                                           >  2:488084/1‑148 (MQ=255)          gccagccaGCGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAAGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAg                                                                                                                                           <  2:709288/150‑1 (MQ=255)             agccagcGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAACCGCAACATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGgcg                                                                                                                                        >  1:319970/1‑150 (MQ=255)                cagcGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCGG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTg                                                                                                                                     >  1:119789/1‑150 (MQ=255)                cagcGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCTGGCGGTCGCCCGCCTGgcacag                                                                                                                                     <  1:488084/150‑5 (MQ=255)                 agcGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTga                                                                                                                                    <  2:59406/150‑1 (MQ=255)                  gcgcGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTgat                                                                                                                                   <  1:954818/150‑1 (MQ=255)                   cgcgTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGCCCTCGGCGGGCGACCGCCAGGCGCTgatc                                                                                                                                  <  1:1069034/150‑1 (MQ=255)                         gggCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCAc                                                                                                                            <  1:496429/150‑1 (MQ=255)                               gccctggccctCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCCACATCGTGCAGTTTCAG‑GACCACCTGTGCCAGGCGCGCGACCGCCAGGCGCTGATCGATCACCTGCAc                                                                                                                      <  2:1043327/150‑1 (MQ=255)                                   tggccctCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATc                                                                                                                  <  1:703208/150‑1 (MQ=255)                                      ccctCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGAACTGTGCCAGGTGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGtt                                                                                                               <  1:788359/150‑1 (MQ=255)                                       cctCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGACGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTa                                                                                                              >  1:1116888/1‑150 (MQ=255)                                        ctCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTAt                                                                                                             >  2:147657/1‑150 (MQ=255)                                         tCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATg                                                                                                            >  2:700852/1‑150 (MQ=255)                                                   acaGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTAcg                                                                                                  >  2:1172947/1‑150 (MQ=255)                                                   acaGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTAcg                                                                                                  >  1:1000228/1‑150 (MQ=255)                                                         cggcgATACGGTGGTCAATGGCAGTGGCAGCAGCCAGCCGCAGGCAGCCCACCCCCACATCGTGCA‑TACGAGCTTCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCcat                                                                                            <  1:677421/150‑1 (MQ=255)                                                                        cAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTc                                                                             <  2:1000228/150‑1 (MQ=255)                                                                           tGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCt                                                                          >  1:667898/1‑150 (MQ=255)                                                                             ccAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCTTGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTAt                                                                        <  1:509476/150‑1 (MQ=255)                                                                             ccAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTAt                                                                        >  2:110928/1‑150 (MQ=255)                                                                               aGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATgg                                                                      <  2:1174466/150‑1 (MQ=255)                                                                                      gcaTCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTcgacg                                                               >  2:934757/1‑150 (MQ=255)                                                                                         tCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGacgacg                                                            <  2:1034814/150‑1 (MQ=255)                                                                                            aGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCTGGCGGTcgcccgcctgg                                                                                                                                          >  1:537485/1‑66 (MQ=37)                                                                                            aGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCTGGCGGTCGCCCGCCTgg                                                                                                                                          <  2:537485/69‑1 (MQ=37)                                                                                            aGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCCTCCACAACGCCTCGCTATGGCTCGACGACGGcc                                                         >  2:550096/1‑150 (MQ=255)                                                                                            aGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACCACGCCTCGCTATGGCTCGACGACGGcc                                                         <  2:625613/150‑1 (MQ=255)                                                                                            aGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGcc                                                         >  1:402922/1‑150 (MQ=255)                                                                                              ctgcaggcagCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTc                                                       >  1:949188/3‑150 (MQ=255)                                                                                               cGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCCAGAGTAGATACGATGGGTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCg                                                      >  1:1114062/1‑150 (MQ=255)                                                                                                      agCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCt                                               >  2:930100/1‑150 (MQ=255)                                                                                                       gCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTc                                              >  2:224176/1‑150 (MQ=255)                                                                                                        cccACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATg                                                                                                            >  2:217657/1‑87 (MQ=255)                                                                                                        cccACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATg                                                                                                            <  1:217657/87‑1 (MQ=255)                                                                                                         ccACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCga                                            <  2:1027003/150‑1 (MQ=255)                                                                                                            tcccAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGCGCTACCCCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGAc                                         <  2:460651/149‑1 (MQ=255)                                                                                                             cccAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGAcc                                        >  1:78654/1‑150 (MQ=255)                                                                                                              ccAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCa                                       >  1:460651/1‑150 (MQ=255)                                                                                                               cAACATCGTGCAGTTTCAG‑GCCGACCTGTTGCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCAt                                      >  1:1125507/1‑150 (MQ=255)                                                                                                               cAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCAt                                      >  1:461739/1‑150 (MQ=255)                                                                                                                      gTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGCCGGACTACGCGCCATTCACGACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGc                               <  1:147657/150‑1 (MQ=255)                                                                                                                      gTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGc                               >  1:225754/1‑150 (MQ=255)                                                                                                                      gTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGc                               <  1:14414/150‑1 (MQ=255)                                                                                                                                   gCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCGCGAGACCATGTTCCGCCTGCACGTCGaa                   <  1:550096/150‑1 (MQ=255)                                                                                                                                     cGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAgc                 <  1:700852/150‑1 (MQ=255)                                                                                                                                       aCCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAgcgc               <  1:930100/150‑1 (MQ=255)                                                                                                                                       aCCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAgcgc               <  2:402922/150‑1 (MQ=255)                                                                                                                                             gCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTacc         <  2:78654/150‑1 (MQ=255)                                                                                                                                              cctggcggtcgCCCGCCTGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAAACTCGAGACCATGTTCCGCCTGAACGTCGAAGCGCCCTacca        >  1:1033896/4‑150 (MQ=255)                                                                                                                                               cAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTaccac       <  2:225754/150‑1 (MQ=255)                                                                                                                                                    gggcgACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATACACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCaa  >  2:670194/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                    GTTACCGGAGCCAGCCAGCGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAG‑GCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAA  >  NC_002947/5255582‑5255877

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑