Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,926,292	C→T	intergenic (+96/‑16)	queF → / → ppnN	7‑cyano‑7‑deazaguanine reductase/nucleotide 5'‑monophosphate nucleosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,926,292	0	C	T	93.8%	84.8 / ‑3.1	32	intergenic (+96/‑16)	queF/ppnN	7‑cyano‑7‑deazaguanine reductase/nucleotide 5'‑monophosphate nucleosidase
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (17/13);  total (19/13)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAA‑GTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATT  >  NC_000913/2926047‑2926538                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        gAACGCATCTTTAATGCCCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATa                                                                                                                                                                                                                                                     >  1:355330/1‑250 (MQ=255)     gCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAg                                                                                                                                                                                                                                                <  1:202911/251‑1 (MQ=255)                      ttACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTacac                                                                                                                                                                                                                                >  1:295814/1‑250 (MQ=255)                         cGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACatat                                                                                                                                                                                                                            >  1:190779/1‑251 (MQ=255)                           cTTCTGCCCGCCAGAAAAATTGCGCGTTTACGCACGATATACCCGTCGTGGCGGTATGGAGATTAACCCGTGGCGCAGAAATAGCGATTTTGGCCCAACGACCACACGACTGGTTCGGAAATAAATTTTTTCTCAATTTTGCGTGCTGGATTTACGCAGAAGGATGTGAAAGGTCATCAGGCGGGGCTATTGTAATCAAAGGGAATGAAGATATTGGTACAATAAGGAATTTTTCTTGATTAcacaaata                                                                                                                                                                                                                           >  2:187573/1‑246 (MQ=255)                                        gAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTc                                                                                                                                                                                                             <  1:66380/251‑1 (MQ=255)                                              ttGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATgg                                                                                                                                                                                                        >  1:285653/1‑250 (MQ=255)                                                        aCGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCg                                                                                                                                                                                             >  1:351232/1‑251 (MQ=255)                                                                 ataCCCGTCGTGGCGGTCTGGACATTAACCCGTGGAGCAGTAAAAGCGCTGTTGTCCCCTCGACCAAAAAACTGGATCGGCCATAAATTTTTTCTCACTTTTGCGTGCTGGGTTCACGCAGAAGGTTTTGAAAGGGCTTCAGGCAGGGCTATTTTAATCAAAGGGAATGTCGATATTCGTCTCCTAAGGAGTTTTTCTTTATTACACATACTCGCCCGCTTTGCTCCATGGATTTGTTGTCggcactgga                                                                                                                                                                                     >  2:118816/1‑242 (MQ=255)                                                                            ggctgTCTGGACATTAACCCGTGGCGCAGTAATATCGATGTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATTTTGTCGCAGCTGGAA‑GTGGATATGc                                                                                                                                                                         <  1:251060/246‑1 (MQ=255)                                                                            ggCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATTAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAA‑GTGGATATGCt                                                                                                                                                                        >  1:355895/1‑251 (MQ=255)                                                                                                          aaTAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAA‑GTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCa                                                                                                                                          <  2:190779/251‑1 (MQ=255)                                                                                                             aGCGATTTTGGCCCATCGCCCACAAGACTGGTTCGGCAATAAATTGTTTCTCAATTTTGCTTGCTGGATGCACGCAGAAGGTTGTGAAAGTTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATCCTTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAA‑GTGGATATGCTCAAACGCACCGCCAGCAGCGACCTCTATCAACt                                                                                                                                       <  1:80634/251‑1 (MQ=255)                                                                                                                     tGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAA‑GTGGACATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGc                                                                                                                                >  1:90219/1‑250 (MQ=255)                                                                                                                              gACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAA‑GTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCAc                                                                                                                      >  1:58004/1‑251 (MQ=255)                                                                                                                                           gTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAA‑GTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAAc                                                                                                                                        >  1:222958/1‑220 (MQ=255)                                                                                                                                                                           gCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTTTGGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAA‑GTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCTGTAGTTTGACCGATAACAGCAAAGAAt                                                                         <  1:118816/251‑1 (MQ=255)                                                                                                                                                                              gtaTTCACGCTGAAGGTTGTGAAAGGTCATCAGGTAGGTTTATGGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTGGATTACACATATGATCCCGTTTGGCTCCATGGTTATGGGGCCGCTGCTGGAA‑GTGTATATGCTTAACCGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGGAGTTTGACCGATAACAGCAAAGAATTGc                                                                      <  2:244001/249‑1 (MQ=255)                                                                                                                                                                                    aCGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAC‑GTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGtctc                                                                >  2:368976/1‑251 (MQ=255)                                                                                                                                                                                     cGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAA‑GTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCg                                                               >  2:191548/1‑251 (MQ=255)                                                                                                                                                                                       aaGAAGGTGGTGAAAGGCCATCAGGCATGGCTATTTTAATAAAAGGGAATGACGTTATTCGTTTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATTTGTTGCCGCAGCTGGAA‑GTGGATATGCTTAACCGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGtt                                                             <  2:285653/250‑1 (MQ=255)                                                                                                                                                                                         gAAGGTTGTGAAAGGCCATCAGGCAGGGCTATGGTAATCAATGGGAAGGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACAGATTAGCCCTCTTGGCTCCATGGATATGTGGCCGCTG‑TGGAAAGTGGATATGCTTAAACGCCCCGCCAGCAGCGCCCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTTCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGtttt                                                           <  2:331679/251‑1 (MQ=255)                                                                                                                                                                                                        tcatcaGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAA‑GTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAACTTTCGATAtt                                            >  2:219484/1‑251 (MQ=255)                                                                                                                                                                                                         catcaGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATAGGTTGTCGCAGCTGGAA‑GTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTa                                           <  1:121737/251‑1 (MQ=255)                                                                                                                                                                                                                 aggGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAA‑GTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTg                                   >  2:396315/1‑251 (MQ=255)                                                                                                                                                                                                                               cAAAGGGAATGCCGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAA‑GTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGCACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACgcgg                     >  2:22433/1‑251 (MQ=255)                                                                                                                                                                                                                                aaaGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAA‑GTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGCTAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAAAgcggc                    >  2:10685/1‑251 (MQ=255)                                                                                                                                                                                                                                aaaGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAA‑GTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATAtt                                            <  2:26031/227‑1 (MQ=255)                                                                                                                                                                                                                                aaaGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAA‑GTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATAtt                                            >  1:26031/1‑227 (MQ=255)                                                                                                                                                                                                                                               atTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAA‑GTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTg     >  1:65536/1‑251 (MQ=255)                                                                                                                                                                                                                                                  cGTCTCATAAGGATTTTTTCTTGATTACCAATATTAGCCCGCTTGGCCCCAGGGATATGTTGTCGCAGCTGGAA‑GTGGATATGCTAAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGAtt  <  2:351232/251‑1 (MQ=255)                                                                                                                                                                                                                                                  cGTCTCATAAGGAGTTTTTTTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAA‑GTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCTATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGAt   <  2:355895/250‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        GAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAA‑GTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATT  >  NC_000913/2926047‑2926538

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