Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	783,000	C→A	7.5%	Q306K (CAG→AAG)	nadA →	quinolinate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	783,000	0	C	A	7.5%	113.7 / 2.1	55	Q306K (CAG→AAG)	nadA	quinolinate synthase
Reads supporting (aligned to +/- strand):  ref base C (28/20);  new base A (2/2);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.43e-01

ATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGAATAAT  >  NC_000913/782859‑783136                                                                                                                                              |                                                                                                                                         aTTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCa                                                                                                                                         >  2:174513/1‑143 (MQ=255)     tgGCAACCGATCGGGGTATTTTCTACAGAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGcc                                                                                                                                     <  1:711273/143‑1 (MQ=255)     tgGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCGTTCAGGcc                                                                                                                                     >  1:340626/1‑143 (MQ=255)          aCCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCGGGCCATCGc                                                                                                                                >  1:81580/1‑143 (MQ=255)            cGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATCGTCCGTGGATGGCCATGCATGGCCTTCAGGCCCTCGccg                                                                                                                              >  2:485926/1‑141 (MQ=255)              aTCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCagag                                                                                                                            <  1:833869/143‑1 (MQ=255)               tCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCTGAgg                                                                                                                           >  1:856159/1‑143 (MQ=255)                    gTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCAGCGCGCATTGGCCGTGGATGGCCATGCATGGCCTTCAGGCCATCGCAGACGCATTa                                                                                                                      >  2:124735/1‑143 (MQ=255)                    gTACTTCCTACAAAAGGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTa                                                                                                                      <  1:4806/143‑1 (MQ=255)                     tATTGTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAGCCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAg                                                                                                                     <  2:726736/143‑1 (MQ=255)                     gaTTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACGAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAg                                                                                                                     <  2:244049/142‑1 (MQ=255)                        tttCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAAc                                                                                                                  <  2:712263/143‑1 (MQ=255)                             aCAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGTACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGCCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACaggaa                                                                                                             >  1:666144/1‑143 (MQ=255)                                aaaTGCAGCAGGCGATGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGAAGCAGCGCGCATTGTCAGTGGATGGCCATGAATGGCATTCAGTCGATCGCAGAGGCATTAGAACaggaagga                                                                                                          >  2:325850/1‑143 (MQ=255)                                      agcagGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAAt                                                                                                    <  2:631303/143‑1 (MQ=255)                                      agcagGCGGTGCCAGAGAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAAt                                                                                                    <  2:414040/143‑1 (MQ=255)                                            cGGTGTGAGATATAGCGTTACTGGAAGCGCCGAACGCAGGTGGGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCGGTGGATGGCCATGAATGGCCGTGAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAg                                                                                              <  2:813827/143‑1 (MQ=255)                                                         agagTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGCTCATGTtgatg                                                                                 >  2:947212/1‑143 (MQ=255)                                                                cTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAggat                                                                          >  1:868784/1‑141 (MQ=255)                                                                cTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCAGCGGCTGCGCGCGGTGGCCGTGGATGGCCATGAATGGCCTTCAGGCAATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGGTCATGTTGACGGAAGGCt                                                                          >  1:229490/1‑143 (MQ=255)                                                                   gAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCCTTAGAACAGGAAGGAAGCAATCACGAGGGCCATGTTGATGAACGGCTGCg                                                                       >  2:56966/1‑143 (MQ=255)                                                                    aaGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGCGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGACAGGCTGCga                                                                      >  1:176736/1‑143 (MQ=255)                                                                        gccAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCgagaga                                                                  <  2:17273/142‑1 (MQ=255)                                                                        aCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCgagaga                                                                  >  1:437673/1‑143 (MQ=255)                                                                        aCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCgagaga                                                                  <  1:550528/143‑1 (MQ=255)                                                                        aCCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCgagaga                                                                  <  1:750262/143‑1 (MQ=255)                                                                             ccGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTGCGTGGATGGCCATGAATGGCGTTCAGGCCATCGCAGAGGAATTAGAACAGGAAGGAAGCAAGCACGAGGTTCAGGTTGAGGAAAGGCGGCGAGGGAGGgcg                                                             >  1:61336/1‑143 (MQ=255)                                                                             ccGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGCTGAAAGGCTGAGCGCGAGGgcg                                                             >  2:655259/1‑143 (MQ=255)                                                                                        ggTGCAACCTGCCGCCGCTGCGCGCATTGGCCGTGGATGGCCATGAGTGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCt                                                  <  2:937570/143‑1 (MQ=255)                                                                                          tGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGCACGGCTGCGAGAGAGGCCGCTGGcgccgcccg                                                >  2:420507/1‑140 (MQ=255)                                                                                          tGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGCGCGAGGGCGCTGGTGCCGCTCa                                                >  1:771572/1‑143 (MQ=255)                                                                                          tGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGGGCATGTTGATGAAGGGCTGCGAGAGAGGGCGCTGGTGCCGCGCa                                                >  2:527626/1‑143 (MQ=255)                                                                                           gCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCaa                                               <  1:174513/143‑1 (MQ=255)                                                                                             aaCCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCaagc                                             >  1:580340/1‑141 (MQ=255)                                                                                               ccTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGGATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGt                                           >  2:802075/1‑143 (MQ=255)                                                                                               ccTGCCGCAGCCGCGCGGGTGGGCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCAGGAGGTGCAGGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGt                                           <  2:420660/143‑1 (MQ=255)                                                                                                   ccGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAAGACGCTGATTCCGCTCAATCATATga                                       >  1:795657/1‑142 (MQ=255)                                                                                                    cGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGATGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCt                                      <  1:279024/143‑1 (MQ=255)                                                                                                     gcggttGCGCGCATTGTCCGGTGGTGGCCATGGATGGCCTTCAGGCCATCGCAGAGGCATTCGAACAGGAAGTAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGACAGGGCGCTGGTGCCGCTCAATCGTTTGCTg                                     <  2:820075/138‑1 (MQ=255)                                                                                                     gCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTg                                     <  1:539852/143‑1 (MQ=255)                                                                                                       ggCTGCGCGCATTGGCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGa                                   <  2:275849/142‑1 (MQ=255)                                                                                                             gcgcCTTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGc                             <  2:751868/143‑1 (MQ=255)                                                                                                             gcgcATTGGCGGTGGATGGCCCGGAATGCCCTCGCGGTTGTGCCCGATCCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGc                             <  1:451089/143‑1 (MQ=255)                                                                                                                cATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCTTATGCTGGCTTTTGCgcc                          >  2:453901/1‑141 (MQ=255)                                                                                                                     tCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTacac                     >  1:296224/1‑143 (MQ=255)                                                                                                                            aTGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAGGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGGATGCTGGAATTTGCGGCTACAATACCTgg              >  1:783794/1‑143 (MQ=255)                                                                                                                            aTGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACCTgg              >  2:736364/1‑143 (MQ=255)                                                                                                                                tctgcgaTGGCCTGAAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATaa          >  1:168766/7‑143 (MQ=255)                                                                                                                                ccATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTGCACGACGTGGATaa          >  1:884983/1‑143 (MQ=255)                                                                                                                                ccATGAATGGCCTTCAGGCCATCGCAGAGGCAATAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTGTGCTGGATTTTGCGGCTACACTACGTGGATaa          >  2:322345/1‑143 (MQ=255)                                                                                                                                  aTGAATGGCCTTGAGGCCAACGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACg        <  1:131065/143‑1 (MQ=255)                                                                                                                                     gaTGGCCTGAAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGaat     <  2:168766/142‑1 (MQ=255)                                                                                                                                      aTGGCCTGAAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCgagaga                                                                  <  2:180727/81‑1 (MQ=255)                                                                                                                                      aTGGCCTGAAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCgagaga                                                                  >  1:180727/1‑81 (MQ=255)                                                                                                                                        ggCCTTCCGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGATGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGaataat  <  1:186890/143‑1 (MQ=255)                                                                                                                                              |                                                                                                                                         ATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGTGCAACCTGCCGCAGCTGCGCGCATTGTCCGTGGATGGCCATGAATGGCCTTCAGGCCATCGCAGAGGCATTAGAACAGGAAGGAAGCAATCACGAGGTTCATGTTGATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGAATAAT  >  NC_000913/782859‑783136

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