Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,410,289	T→A	8.0%	I320N (ATT→AAT)	nadA →	quinolinate synthase iron‑sulfur cluster subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,410,289	0	T	A	8.0%	144.8 / 2.3	50	I320N (ATT→AAT)	nadA	quinolinate synthase iron‑sulfur cluster subunit
Reads supporting (aligned to +/- strand):  ref base T (18/28);  new base A (2/2);  total (20/30)
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.86e-01

TGCAGCAGTTGTGCCCGGACAAGGAATTCGTCGAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAATCGCATGCTGGACTTTACCCAG‑GCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGACGCCAGACCGCGTCATACTCTTCGCGGGCTCGCCCACTCCCACAG  >  NC_002947/1410142‑1410434                                                                                                                                                    |                                                                                                                                                     tGCAGCAGTTGTGCCCGGACAAGGAATTCGTCGAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAAtt                                                                                                                                                     <  5:7551/149‑1 (MQ=255)  gcagcaGTTGTGCCCGGACAAGGAATTCGTCGAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAAttt                                                                                                                                                    >  5:94234/1‑149 (MQ=255)    agcagTTGTGCCCGGACAAGGAATTCGTCGAAGCCCCCCCCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAAATTTc                                                                                                                                                  >  7:299582/1‑149 (MQ=255)    agcagTTGTGCCCGGACAAGGAATTCGTCGAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTc                                                                                                                                                  <  5:152296/149‑1 (MQ=255)    agcagTTGTGCCCGGACAAGGAATTCGTCGAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTc                                                                                                                                                  >  3:124269/1‑149 (MQ=255)             gCCCGGACAAGGAATTCGTCGAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCg                                                                                                                                         <  5:10640/149‑1 (MQ=255)               ccGGACAAGGAATTCGTCGAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAAATTTCGTTggtgccgt                                                                                                                                       <  8:299582/149‑8 (MQ=255)                 ggACAAGGAATTCGTCGAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGgcgc                                                                                                                                     <  6:221270/149‑1 (MQ=255)                       ggAATTCGTCGAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAAATTTCGTTggtgcc                                                                                                                                         <  8:256779/139‑6 (MQ=255)                       ggAATTCGTCGAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAAATTTCGTTggtgcc                                                                                                                                         >  7:256779/1‑134 (MQ=255)                         aaTTCGTCGAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGcc‑c                                                                                                                            <  3:285786/149‑1 (MQ=255)                                 gAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGAt                                                                                                                    >  4:108837/1‑149 (MQ=255)                                  aaGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATc                                                                                                                   <  4:13465/149‑1 (MQ=255)                                  aaGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATc                                                                                                                   <  2:168968/149‑1 (MQ=255)                                        ccACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGgcgc                                                                                                                                     <  8:271136/126‑1 (MQ=255)                                         cACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑Gc                                                                                                           >  3:279813/1‑149 (MQ=255)                                         cACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑Gc                                                                                                           <  1:49641/149‑1 (MQ=255)                                           ccgccgGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCtt                                                                                                         <  3:108837/149‑1 (MQ=255)                                              ccgGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAAt                                                                                                      <  3:176501/149‑1 (MQ=255)                                              ccgGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAAt                                                                                                      <  7:264794/149‑1 (MQ=255)                                              ccgGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAAt                                                                                                      <  3:308562/149‑1 (MQ=255)                                                 gCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAATCg                                                                                                    >  7:108102/1‑148 (MQ=255)                                                    aCGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAATCGCATg                                                                                                <  4:205394/149‑1 (MQ=255)                                                       gcgcGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGa                                                                                                                     >  1:231448/1‑126 (MQ=255)                                                       gcgcGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGa                                                                                                                     <  2:231448/126‑1 (MQ=255)                                                       gcgcGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAATCGCATGCTg                                                                                             <  4:124269/149‑1 (MQ=255)                                                         gcgGCATGCCGCAGCTGTGCGCAGTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGACTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAATCGCATGCTGGa                                                                                           <  6:250853/149‑1 (MQ=255)                                                          cgGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAATCGCATGCTGGAc                                                                                          <  4:279813/149‑1 (MQ=255)                                                          cgGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAATCGCATGCTGGAc                                                                                          <  4:89517/149‑1 (MQ=255)                                                             cATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAATCGCATGCTGGACttt                                                                                       <  6:94234/149‑1 (MQ=255)                                                                       cTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAATCGCATGCTGGACTTTACCCAg‑gcgg                                                                            >  2:217154/1‑149 (MQ=255)                                                                                         gATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAATCGCATGCTGGACTTTACCCAG‑GCGGCGCGGCTGAAGCTGTCt                                                           <  7:52092/148‑1 (MQ=255)                                                                                              cGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAATCGCATGCTGGACTTTACCCAG‑GCGGCGCGGCTGAAGCTGTCTGGTAAt                                                     >  4:270979/1‑149 (MQ=255)                                                                                               gATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAATCGCATGCTGGACTTTACCCAG‑GCGGCGCGGCTGAAGCTGTCTGGTAATg                                                    >  1:213171/1‑149 (MQ=255)                                                                                                    acacCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAATCGCATGCTGGACTTTACCCAG‑GCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGa                                               >  5:47121/1‑149 (MQ=255)                                                                                                      acCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAATCGCATGCTGGACTTTACCCAGTGCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGAc                                              <  1:217154/149‑1 (MQ=255)                                                                                                        ccTGGTGCGGGTGGTGGATTGCCTGCGTAACGGCAACAACTAAATTTTCTTTTACCCGTCGCTGGTGCC‑GAAGGCTATCAAGCCTATTTAATCGCATGCTGGACTTTACCCAG‑GCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGACGc                                            <  1:271678/149‑1 (MQ=255)                                                                                                        ccTGGAGCGGGTGCTGGATTGCCTGCATAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAATCGCATGCTGGACTTTACCCAG‑GCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGACGc                                            >  1:133899/1‑148 (MQ=255)                                                                                                            gAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAATCGCATGCTGGACTTTACCCAG‑GCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGACGCCAGAc                                       >  4:52964/1‑149 (MQ=255)                                                                                                                        gATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAATCGCATGCTGGACTTTACCCAG‑GCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGACGCCAGACCGCGTCATAct                            <  3:270979/148‑1 (MQ=255)                                                                                                                             ccTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAATCGCATGCTGGACTTTACCCAG‑GCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGACGCCAGACCGCGTCATACTCTTcgc                      <  6:47121/149‑1 (MQ=255)                                                                                                                                gCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAATCGCATGCTGGACTTTACCCAG‑GCGGCGCGGCTGAAGCTGTCTGGTaa                                                      >  3:1416/1‑114 (MQ=255)                                                                                                                                gCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAATCGCATGCTGGACTTTACCCAG‑GCGGCGCGGCTGAAGCTGTCTGGTaa                                                      <  4:1416/114‑1 (MQ=255)                                                                                                                                 cGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAATCGCATGCTGGACTTTACCCAG‑GCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGACGCCAGACCGCGTCATACTCTTCGCGGGc                  <  6:167796/149‑1 (MQ=255)                                                                                                                                  gTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCGCAAGGCGATCAAGCC‑GCTTAATCGCATGCTGGACTTTACCCAG‑GCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGACGCCAGACCGCGTCATACTCTTCGCGGGc                  >  1:65331/1‑149 (MQ=255)                                                                                                                                     aCGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAATCGCATGCTGGACTTTACCCAG‑GCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGACGCCAGACCGCGTCATACTCTTCGCGGGCTCGc              >  3:223718/1‑149 (MQ=255)                                                                                                                                      cGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAATCGCATGCTGGACTTTACCCAG‑GCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGACGCCAGACCGCGTCATACTCTTCGCGGGCTCGcc             >  8:150161/1‑149 (MQ=255)                                                                                                                                        gCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAATCGCATGCTGGACTTTACCCAG‑GCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGACGCCAGACCGCGTCATACTCTTCGCGGGCTCGCCCa           >  4:181394/1‑149 (MQ=255)                                                                                                                                             aaCGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAATCGCATGCTGGACTTTACCCAG‑GCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGACGCCAGACCGCGTCATACTCTTCGCGGGCTCGCCCACTcc       >  2:189455/1‑148 (MQ=255)                                                                                                                                               cGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAATCGCATGCTGGACTTTACCCAG‑GCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGACGCCAGACCGCGTCATACTCTTCGCGGGCTCGCCCACTCCca     >  8:215844/1‑148 (MQ=255)                                                                                                                                                gAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAATCGCATGCTGGACTTTACCCAG‑GCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGACGCCAGACCGCGTCATACTCTTCGCGGGCTCGCCCACTCCcaca   <  7:215844/149‑1 (MQ=255)                                                                                                                                                 aaaTTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAATCGCATGCTGGACTTTACCCAG‑GCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGACGCCAGACCGCGTCATACTCTTCGCGGGCTCGCCCACTCCCACAg  <  3:181394/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     TGCAGCAGTTGTGCCCGGACAAGGAATTCGTCGAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCC‑CAAGGCGATCAAGCC‑GCTTAATCGCATGCTGGACTTTACCCAG‑GCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGACGCCAGACCGCGTCATACTCTTCGCGGGCTCGCCCACTCCCACAG  >  NC_002947/1410142‑1410434

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