Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	903,026	T→G	10.5%	N158T (AAC→ACC)	nadC ←	quinolinate phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	903,026	0	T	G	10.5%	100.1 / 4.3	38	N158T (AAC→ACC)	nadC	quinolinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (17/17);  new base G (2/2);  total (19/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 = 9.06e-01

TTCGTCCAGGCTTTCGACCTCGATTTCCACCGGCTTGCCCGGCGCAATGCGGTGTGCTGCGGCCACCGCCTCGGCTACGCCGCCGCTGGCGGCGATGTGGTTTTCCTTGATCAGGAAGGCGTCGTACAGGCCTATGCGGTGGTTGTCGCACCCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGCGGGTGTCCAGCAGGCGCACCTGGGTGCCGGCTACCAGGTCGGCGAGAAAGCGCGCACGGGTGGCGACCCCTGACAGCAGCTG  >  NC_002947/902884‑903171                                                                                                                                               |                                                                                                                                                  ttCGTCCAGGCTTTCGACCTCGATTTCCACCGGCTTGCCCGGCGCAATGCGGTGTGCTGCGGCCACCGCCTCGCCTACGCCGCCGCTGGCGGCGATGTGGTTTTCCTTGATCAGGAAGGCGTCGTACAGGCCTGTGCGGTGGTTGTCg                                                                                                                                              >  3:184315/1‑148 (MQ=255)    gTCCAGGCTGTCGACCTCGATTTCCACCGGCTTGCCCGGCGCAATGCGGTGTGCTGCGGCCACCGCCTCGGCTACGCCGCCGCTGGCGGCGATGTGGTTTTCCTTGATCAGGAAGGCGTCGTACAGGCCTATGCGGTGGTTGTCGcac                                                                                                                                           <  6:307058/148‑1 (MQ=255)         ggCTTTCGACCTCGATTTCCACCGGCTTGCCCGGCGCAATGCGGTGTGCTGCGGCCACCGCCTCGGCTACGCCGCCGCTGGCGGCGATGTGGTTTTCCTTGATCAGGAAGGCGTCGTACAGGCCTATGCGGTGGTTGTCGCACCCACCg                                                                                                                                     <  4:246073/149‑1 (MQ=255)            tttCGACCTCGATTTCCACCGGCTTGCCCGGCGCAATGCGGTGTGCTGCGGCCACCGCCTCGGCTACGCCGCCGCTGGCGGCGATGTGGTTTTCCTTGATCAGGAAGGCGTCGTACAGGCCTATGCGGTGGTTGTCGCACCCACCg                                                                                                                                     >  3:139266/1‑146 (MQ=255)              tCGACCTCGATTTCCACCGGCTTGCCCGGCGCAATGCGGTGTGCTGCGGCCACCGCCTCGGCTACGCCGCCGCTGGCGGCGATGTGGTTTTCCTTGATCAGGAAGGCGTCGTACAGGCCTATGCGGTGGTTGTCGCACCCACCGCa                                                                                                                                   <  4:139266/146‑1 (MQ=255)               cGACCTCGATTTCCACCGGCTTGCCCGGCGCAATGCGGTGTGCTGCGGCCACCGCCTCGGCTACGCCGCCGCTGGCGGCGATGTGGTTTTCCTTGATCAGGAAGGCGTCGTACAGGCCTATGCGGTGGTTGTCGCACCCACCGCAGGt                                                                                                                                <  3:305867/148‑1 (MQ=255)                    tCGATTTCCACCGGCTTGCCCGGCGCAATGCGGTGTGCTGCGGCCACCGCCTCGGCTACGCCGCCGCTGGCGGCGATGTGGTTTTCCTTGATCAGGAAGGCGTCGTACAGGCCTATGCGGTGGGTGCGACAACCACCGCata                                                                                                                                 >  4:294831/1‑140 (MQ=255)                    tCGATTTCCACCGGCTTGCCCGGCGCAATGCGGTGTGCTGCGGCCACCGCCTCGGCTACGCCGCCGCTGGCGGCGATGTGGTTTTCCTTGATCAGGAAGGCGTCGTACAGGCCTATGCGGTGGGTGCGACAACCACCGCata                                                                                                                                 <  3:294831/142‑3 (MQ=255)                                       ccGGCGCAATGCGGTGTGCTGCGGCCACCGCCTCGGCTACGCCGCCGCTGGCGGCGATGTGGTTTTCCTTGATCAGGAAGGCGTCGTACAGGCCTATGCGGTGGTTGTCGCACCCACCGCAGGTCACTGCGTATTTTTGCGCCAGgcgc                                                                                                       <  4:58667/149‑1 (MQ=255)                                          gcgcAATGCGGTGTGCTGCGGCCACCGCCTCGGCTACGCCGCCGCTGGCGGCGATGTGGTTTTCCTTGATCAGGAAGGCGTCGTACAGGCCTATGCGGTGGTTGTCGCAGCCAGCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAgg                                                                                                    <  2:76474/149‑1 (MQ=255)                                                   ggTGTGCTGCGGCCACCGCCTCGGCTACGCCGCCGCTGGCGGCGATGTGGTTTTCCTTGATCAGGAAGGCGTCGTACAGGCCTATGCGGTGGtt                                                                                                                                                  <  3:46025/94‑1 (MQ=255)                                                   ggTGTGCTGCGGCCACCGCCTCGGCTACGCCGCCGCTGGCGGCGATGTGGTTTTCCTTGATCAGGAAGGCGTCGTACAGGCCTATGCGGTGGtt                                                                                                                                                  >  4:46025/1‑94 (MQ=255)                                                         ctgcGGCCACCGCCTCGGCTACGCCGCCGCTGGCGGCGATCTGGTTTTCCTTGATCAGGAAGGCGTCGTACAGGCCTATGCGGTGGTTGTCGCACCCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGtttt                                                                                      <  5:314131/148‑1 (MQ=255)                                                          tgcGGCCACCGCCTCGGCTACGCCGCCGCTGGCGGCGATGTGGTTTTCCTTGATCAGGAAGGCGTCGTACAGGCCTATGCGGTGGTTGTCGCACCCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGc                                                                                    >  1:156631/1‑149 (MQ=255)                                                            cGGCCACCGCCTCGGCTACGCCGCCGCTGGCGGCGATGTGGTTTTCCTTGATCAGGAAGGCGTCGTACAGGCCTATGCGGTGGTTGTCGCACCCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGCgg                                                                                  <  6:23478/149‑1 (MQ=255)                                                            cGGCCACCGCCTCAGCTACGCCGCCGCTGGCGGCGATGTGGTTTTCCTTGATCAGGAAGGCGTCGTACAGGCCTATGCGGTGGTTGTCGCACCCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGCgg                                                                                  <  8:32326/149‑1 (MQ=255)                                                              gCCACCGCCTCGGCTACGCCGCCGCTGGCGGCGATGTGGTTTTCCTTGATCAGGAAGGCGTCGTACAGGCCTATGCGGTGGTTGTCGCACCCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGCGGgt                                                                                <  6:166952/149‑1 (MQ=255)                                                                    gCCTCGGCTACGCCGCCGCTGGCGGCGATGTGGTTTTCCTTGATCAGGAAGGCGTCGTACAGGCCTATGCGGTGGTTGTCGCACCCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGCGGGTGTCcag                                                                          >  8:193568/1‑149 (MQ=255)                                                                       tCGGCTACGCCGCCGCTGGCGGCGATGTGGTTTTCCTTGATCAGGAAGGCGTCGTACAGGCCTATGCGGTGGTTGTCGCACCCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGCGGGTGTCcagcag                                                                       >  1:81689/1‑149 (MQ=255)                                                                           cTACGCCGCCGCTGGCGGCGATGTGGTTTTCCTTGATCAGGAAGGCGTCGTACAGGCCTATGCGGTGGTTGTCGCACCCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGCGGGTGTCCAGCAGgcgc                                                                   <  2:81689/149‑1 (MQ=255)                                                                            tACGCCGCCGCTGGCGGCGATGTGGTTTTCCTTGATCAGGAAGGCGTCGTACAGGCCTATGCGGTGGTTGTCGCACCCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGCGGGTGTCCAGCAGGCGCa                                                                  >  3:249993/1‑149 (MQ=255)                                                                               gccgccgcTGGCGGCGATGTGGTTTTCCTTGATCAGGAAGGCGTCGTACAGGCCTATGCGGTGGTTGTCGCACCCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGc                                                                                    >  7:310200/1‑128 (MQ=255)                                                                               gccgccgcTGGCGGCGATGTGGTTTTCCTTGATCAGGAAGGCGTCGTACAGGCCTATGCGGTGGTTGTCGCACCCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGc                                                                                    <  8:310200/128‑1 (MQ=255)                                                                               gccgccgcTGGCGGCGATGTGGTTTTCCTTGATCAGGAAGGCGTCGTACAGGCCTATGCGGTGGTTGTCGCACCCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGCGGGTGTCCAGCAGGCGCACCt                                                               >  1:125685/1‑149 (MQ=255)                                                                                 cgccgcTGGCGGCGATGTGGTTTTCCTTGATCAGGAAGGCGTCGTACAGGCCTATGCGGTGGTTGTCGCACCCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGCGGGTGTCCAGCAGGCGCACCTgg                                                             <  7:193568/149‑1 (MQ=255)                                                                                                       ttCCTTGATCAGGAAGGCGTCGTACAGGCCTATGCGGTGGTTGTCGCACCCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGCGGGTGTCCAGCAGGCGCACCTGGGTGCCGGCTACCAGGTCGGc                                         >  4:769/1‑147 (MQ=255)                                                                                                          cTTGATCAGGAAGGCGTAGTACAGGCCTATGCGGTGGTTGTCGCACCCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGCGGGTGTCCAGCAGGCGCACCTGGGTGACGGCTACCAGGTCGGCgaga                                     >  3:242284/1‑148 (MQ=255)                                                                                                           ttGATCAGGAAGGCGTCGTACAGGCCTATGCGGTGGTTGTCGCACCCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGCGGGTGTCCAGCAGGCGCACCTGGGTGCCGGCTACCAGGTCGGCgaga                                     >  4:114608/1‑147 (MQ=255)                                                                                                               tCAGGAAGGCGTCGTACAGGCCTATGCGGTGGTTGTCGCACCCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGCGGGTGTCCAGCAGGCGCACCTGGGTGCCGGCTACCAGGTCGGCGAGAAAgcgc                               <  4:249993/149‑1 (MQ=255)                                                                                                               tCAGGAAGGAGTCGTACAGGCCTAGGCGGTGGTGGTCGCACCCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGCGGGTGTCCAGCAGGCGCACCTGGGTGCCGGCTACCAGGTCGGCGAGAAAgcgc                               <  4:42799/149‑1 (MQ=255)                                                                                                                cAGGAAGGCGTCGTACAGGCCTATGCGGTGGTTGTCGCACCCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGAGGGTGTCCAGCAGGCGCACCTGGGTGCCGGCTACCAGGTCGGCGAGAAAgcgcg                              >  1:55917/1‑149 (MQ=255)                                                                                                                   gAAGGCGTCGTACAGGCCTATGCGGTGGTTGTCGCACCCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGCGGGTGTCCAGCAGGCGCACCTGGGTGCCGGCTACCAGGTCGGCGAGAAAGCGCGCAc                           >  6:144475/1‑149 (MQ=255)                                                                                                                         gtcgtACAGGCCTATGCGGTGGTTGTCGCACCCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGCGGGTGTCCAGCAGGCGCACCTGGGTGCCGGCTACCAGGTCGGCGAGaaa                                   <  2:113654/135‑1 (MQ=255)                                                                                                                         gtcgtACAGGCCTATGCGGTGGTTGTCGCACCCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGCGGGTGTCCAGCAGGCGCACCTGGGTGCCGGCTACCAGGTCGGCGAGaaa                                   >  1:113654/1‑135 (MQ=255)                                                                                                                             tACAGGCCTATGCGGTGGTTGTCGCACCCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGCGGGTGTCCAGCAGGCGCACCTGGGTGCCGGCTACCAGGTCGGCGAGAAAGCGCGCACGGGTGGCGAc                 >  3:113096/1‑149 (MQ=255)                                                                                                                                aGGCCTATGCGGTGGTTGTCGCACCCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGAGGGTGTCCAGCAGGCGCACCTGTGTGACGGCTACCAGGTCGGCGAGAAAGCGCGCACGGGTGGCGAcacc              >  7:140267/1‑149 (MQ=255)                                                                                                                                   ccTATGCGGTGGTTGTCGCACCCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGCGGGTGTCCAGCAGGCGCACCTGGGTGCCGGCTACCAGGTCGGCGAGAAAGCGCGCACGGGTGGCGACCCCTGa           >  1:43306/1‑149 (MQ=255)                                                                                                                                       tGCGGTGGTTGTCGCACCCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGCGGGTGTCCAGCAGGCGCACCTGGGTGCCGGCTACCAGGTCGGCGAGAAAGCGCGCACGGGTGGCGACCCCTGAca         <  1:26005/147‑1 (MQ=255)                                                                                                                                        gCGGTGGGTGCGACAACCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGCGGGTGTCCAGCAGGCGCACCTGGGTGCCGGCTa                                                   >  3:114591/1‑104 (MQ=255)                                                                                                                                        gCGGTGGGTGCGACAACCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGCGGGTGTCCAGCAGGCGCACCTGGGTGCCGGCTa                                                   <  4:114591/104‑1 (MQ=255)                                                                                                                                            tggtTGTCGCACCCACCGGAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGCGGGTGTCCAGCAGGCGCACCTGGGTGACGGCTACCAGGTCGGCGAGAAAGCGCGCACGGGTGGCGACCCCTGACAGCAGCTg  <  4:242284/149‑1 (MQ=255)                                                                                                                                            tggtTGTCGCACCCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGCGGGTGTCCAGCAGGCGCACCTGGGTGCCGGCTACCAGGTCGGCGAGAAAGCGCGCACGGGTGGCGACCCCTGACAGCAGCTg  <  3:114608/149‑1 (MQ=255)                                                                                                                                              gtTGTCGCACCCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGCGGGTGTCCAGCAGGCGCACCTGGGTGCCGGCTACCAGGTCGGCGAGAAAGCGCGCACGGGTGGCGACCCCTGACAGCAGCTg  >  1:217670/1‑147 (MQ=255)                                                                                                                                               |                                                                                                                                                  TTCGTCCAGGCTTTCGACCTCGATTTCCACCGGCTTGCCCGGCGCAATGCGGTGTGCTGCGGCCACCGCCTCGGCTACGCCGCCGCTGGCGGCGATGTGGTTTTCCTTGATCAGGAAGGCGTCGTACAGGCCTATGCGGTGGTTGTCGCACCCACCGCAGGTCACTGCGTATTTTTGCGCCAGGCGCAGGCCCGGCAGGGTTTTGCGGGTGTCCAGCAGGCGCACCTGGGTGCCGGCTACCAGGTCGGCGAGAAAGCGCGCACGGGTGGCGACCCCTGACAGCAGCTG  >  NC_002947/902884‑903171

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