Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,627,641	C→T	F109F (TTC→TTT)	nadR →	DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,627,641	0	C	T	96.9%	96.7 / NA	32	F109F (TTC→TTT)	nadR	DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (14/17);  minor base G (0/1);  total (14/18)
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.98e-01

CCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTAT‑CCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACAC‑CTCGGAAGAAGCCGATGCGCCACAGTATATGGAAC  >  NC_000913/4627397‑4627879                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 ccAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGCGAGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGATTGAAGAc                                                                                                                                                                                                                                            >  1:679337/1‑251 (MQ=255)           tttAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTACCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGt                                                                                                                                                                                                                                  >  1:624925/1‑251 (MQ=255)                gCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCgcag                                                                                                                                                                                                                             <  1:31203/251‑1 (MQ=255)                     ttACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGcc                                                                                                                                                                                                                        >  1:142257/1‑251 (MQ=255)                            aTGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTg                                                                                                                                                                                                                 <  2:541912/251‑1 (MQ=255)                                 aaaaTCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGaa                                                                                                                                                                                                                                               >  1:624783/1‑216 (MQ=255)                                 aaaaTCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGaa                                                                                                                                                                                                                                               <  2:624783/216‑1 (MQ=255)                                             cccAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCAGTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCgctg                                                                                                                                                                                                <  1:280212/251‑4 (MQ=255)                                               cAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGt                                                                                                                                                                                              >  1:30919/1‑251 (MQ=255)                                                  cgtgcAAAAGCTGGAGGCGTTTCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCGTTGGTGTGGTAGGCGGTAAGTTCTACCCGCTGCATACCGGACATATCTGCCTTAGCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATGCCCTCGATCGCGCGTTGGTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTAt                                                                                                                                                                                           <  1:624189/248‑1 (MQ=255)                                                              ggAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACtttt                                                                                                                                                                                >  1:346832/1‑250 (MQ=255)                                                                       gCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAg                                                                                                                                                                       <  1:89061/250‑1 (MQ=255)                                                                        cACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCCTTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGCTACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGa                                                                                                                                                                      >  1:288179/1‑250 (MQ=255)                                                                                                       agaagaAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTTCCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCTCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAagag                                                                                                                                      <  1:433462/251‑1 (MQ=255)                                                                                                             aaaccattggtttggTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATg                                                                                                                                <  2:30919/237‑1 (MQ=255)                                                                                                                       tgtggTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTAT‑c                                                                                                                     >  1:60292/5‑251 (MQ=255)                                                                                                                            tATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGATTTGACGATACACGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATCTCAGAAAAATATTCGCATTCATGCTTTCAACGACGAGGGCAAGGAGCAGTAT‑CCGCAc                                                                                                                >  2:676904/1‑251 (MQ=255)                                                                                                                             aTTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTAT‑CCGCACg                                                                                                               <  1:88480/251‑1 (MQ=255)                                                                                                                              ttCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACt                                                                                                                                                                                   >  1:678857/1‑183 (MQ=255)                                                                                                                              tgcGGTAAGGTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAATACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACt                                                                                                                                                                                   <  2:678857/181‑1 (MQ=255)                                                                                                                                  gTAAGTTCATCCCGCTGCATATCTCACATATCTACCTTATTCAGCGCGCCGGTCGCCAGTGTGACGAACTGCAGATCCTTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAATCTTTTAAATATCAGAAAAATATTCGCTTGCATGCTTTCAACGAAGAGGGCATGGAGCCGTAT‑CCGCACGGCTgg                                                                                                          <  2:486964/251‑1 (MQ=255)                                                                                                                                             ccgcTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGCAGAGGGCATGGAGCCGTAT‑CCGCACGGCTGGGATGGGTgg                                                                                                 >  2:544492/4‑249 (MQ=255)                                                                                                                                                  gcatatctgACATATCTACCTTATCCAGCGCGCCTGTAGCCATGTGGACGAACTGCATAGCATTATGTGTTGTGACGATCCTCGCGATCGCGCGGTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCTCCTGGCCTGATCGTCTGCTCTGGTTATTACAATCTTTTAAATATCAGAAAAATATTCGCATGCATGCTTTCAACGAAGAGGTCATGGAGCC‑TATTCCGTACGGCTGGGATGTGTGGatcatcg                                                                                          <  2:679337/243‑7 (MQ=255)                                                                                                                                                    aTACCGGACATATCTACCTTAGCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTAT‑CCGCACGGCTGGGATGTGTGGAGCAACGGc                                                                                        <  1:146225/251‑1 (MQ=255)                                                                                                                                                             atatCTACCTTATCCGGCGTGCCGGTTGCCACGGTGACGAACTGCATATCATTATGGGTTGTGACGACACACGTGATTGTCCGTTTTTGGAAGACAGCGCTAGTTCGCAGCAGCCCACCGTGCCGGATCTTCTTCTCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATGCATGCTTTCAACGAAGAGGGCATGGAGCCGAAT‑CCGTACGGCTGTGATGTGTGGAGCAACGGCATCaaaaat                                                                               <  2:673933/251‑2 (MQ=255)                                                                                                                                                                   aCCTTATCCAGCGCGCCTGGTGCCAGGTTGACGAGCTGCATATCACTATGGGATTTGACGATACCCGCGATCGCGCGGTGTTTGAAGACATCGCTATGTCGCAGCAGCCCACCGTGCCGGAACGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAAGATTCTCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTAT‑CCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATg                                                                         <  2:288179/251‑1 (MQ=255)                                                                                                                                                                   aCCTTAGCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTAT‑CCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTa                                                                           <  2:60292/249‑1 (MQ=255)                                                                                                                                                                            aGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTAT‑CCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGaaa                                                                  >  1:672892/1‑249 (MQ=255)                                                                                                                                                                                               tGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGGTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTAT‑CCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGgacctg                                              <  1:629942/250‑1 (MQ=255)                                                                                                                                                                                                                  aTGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCTATATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGATCCGTAT‑CCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGTGATTCAGCCGGACCTGATCTACAC‑CTCGGaagaag                          <  1:51231/250‑1 (MQ=255)                                                                                                                                                                                                                              gATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGAATGGAGCCGTAT‑CCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACAC‑CTCGGCAGAAGCCTATGCGCcac              >  2:202391/1‑250 (MQ=255)                                                                                                                                                                                                                                          cgcgcaTTGTTTGACGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTCTTCCAAACTTTTACATATCAGCAACATATTCGCATTCATGCTTTCAACGAAGAGGGAATGGAGCCGTAT‑CCGAACGGCTGGGTTGTGTGGAACAACGGCATCACAAAGTTTATGGCTGCAAAAGGGATTCGGCCGGGCCTGAT‑AAAACACTCGGAAGAAGCCGATGACCCCCAGTATATggcac  >  2:426006/1‑247 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 CCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTAT‑CCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACAC‑CTCGGAAGAAGCCGATGCGCCACAGTATATGGAAC  >  NC_000913/4627397‑4627879

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