Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,926,772	C→A	G155G (GGC→GGA)	ppnN →	nucleotide 5'‑monophosphate nucleosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,926,772	0	C	A	97.1%	109.7 / ‑6.0	34	G155G (GGC→GGA)	ppnN	nucleotide 5'‑monophosphate nucleosidase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (14/19);  total (15/19)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

GCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAACGAATTGATCATC  >  NC_000913/2926526‑2927009                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              gCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTTCTGGGGAGGTc                                                                                                                                                                                                                                           <  1:256599/251‑1 (MQ=255)             taagcttccGGAAGAGGCTTGTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTGCCCAACCTGAATCGCGACATCGCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAAc                                                                                                                                                                                                                                <  1:201201/243‑1 (MQ=255)                            ggCTTTTGTCGATGGGCTAATTATTCGCGCTTTTCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGt                                                                                                                                                                                                                <  1:140026/251‑1 (MQ=255)                                tttGTCGATTGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCTGTCCTGCGTGACATGCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATg                                                                                                                                                                                                             <  2:4339/250‑1 (MQ=255)                                tttGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTGGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTAt                                                                                                                                                                                                              <  1:116202/249‑1 (MQ=255)                                   gTCGATGGGCGAATTATTCGCGCTGTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATGCTCTTCGTTTACGGGCAAATCCATAACACCTTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCAGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTATCGAAAACGAGTATTTGTATGccca                                                                                                                                                                                                         <  2:210620/251‑2 (MQ=255)                                               attattCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGaaa                                                                                                                                                                                             <  1:257485/251‑1 (MQ=255)                                                  attCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCTTCGGAAACCa                                                                                                                                                                                          <  2:194416/251‑1 (MQ=255)                                                  attCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAAcc                                                                                                                                                                                           >  2:249187/1‑250 (MQ=255)                                                                       ttGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGa                                                                                                                                                                      >  1:257295/1‑250 (MQ=255)                                                                                   cTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCgg                                                                                                                                                         >  2:110848/1‑251 (MQ=255)                                                                                            aTTCTCTTCGTTTACGGGCAAATCCATAACGCCGTTCGTTTTCCCAACCTGAAGCGCGACAACTCCGTGCACAGCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTTGGCCTGCGTGAGCTGGATAGCTGGAACGGCTGTGGTcc                                                                                                                                                <  2:142755/251‑1 (MQ=255)                                                                                                      tttACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATg                                                                                                                                      >  1:192746/1‑251 (MQ=255)                                                                                                           gggCAAATCCATAACACCGTTCGTTTTCCCAACCGGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGCATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAgc                                                                                                                                 >  1:102772/1‑251 (MQ=255)                                                                                                                           ccGTTCGTTTTCCCAACCTGAATCGCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCGGCTGGGGAGGTCACTCAATTAACGAAAGCGAGTTTTTGTATGCCCGTCGCGGCGGAAACCAGCGGGGGCTGCGGGAGCTTAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGtgct                                                                                                                 <  2:35403/251‑1 (MQ=255)                                                                                                                               tcgttTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGGCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTTGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCg                                                                                                              <  2:90167/250‑1 (MQ=255)                                                                                                                                              gAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCgcag                                                                                               <  1:163035/250‑1 (MQ=255)                                                                                                                                                   tcGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGAGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGCTGAAAGGTGCTGCGGTCGGACACGCGCAGAAGCGt                                                                                         >  1:122189/1‑251 (MQ=255)                                                                                                                                                     gACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTa                                                                                       <  1:194946/251‑1 (MQ=255)                                                                                                                                                           cccGTCCACATCACTAACCGGATCGTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCAAAATATGGTGGTCTGCTGGGGAGGTCGCTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCTCCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGa                                                                                 <  2:102772/250‑1 (MQ=255)                                                                                                                                                                      cACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGtttt                                                                       <  1:110848/250‑1 (MQ=255)                                                                                                                                                                                 tCTTTTCCATCTCGCGTAACGCTAGCGAGCTGCATGTGGATGAAGCGCCAAATATGCTGGTCTGCTGCAGAGGTCACTCAATTAACGCAAACGAGTATTTGTATGCCCGCCACGTCGGAAACCAGATGGGCCAGCGTGTGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGCTGCCATCGCCGATGAAAGGTGCTGCGGGCGGACAAGCGCAGCGGCGTGTCAAAGActggagtgttcttgggatgaa                                                            >  2:24135/1‑230 (MQ=255)                                                                                                                                                                                    tttCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGTGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGAGAGCTGAATATCTGCACCGGCTGTGGTCCGGGTGCGATGGAAGCGCCGATGAATGGTGCTGCGGGCGGACACGCGCAGCAGTGTTACAGAGACAGTGGTTTTATTGGTATGACaga                                                         <  2:272183/250‑1 (MQ=255)                                                                                                                                                                                    tttCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACa                                                                               >  1:28547/1‑228 (MQ=255)                                                                                                                                                                                    tttCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACa                                                                               <  2:28547/228‑1 (MQ=255)                                                                                                                                                                                        cATCTTGCGTAACGCTCGCGCGCTGCATGAGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGATTTATTGGTATGACAGATCCGt                                                    >  1:87898/1‑251 (MQ=255)                                                                                                                                                                                              gCGTAAGGCTCGCGCGCTGCATGTGGGAGAAGAGCCACATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGATGGGCCTGCGTGAGCTGAATAACTGCACCGGCTGAAGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGACGGACACGCGCAGACGCGTTACAGAGACAGTCGTTTTTTTGGTATGACAGAGCCGACGAtt                                               >  1:240790/1‑250 (MQ=255)                                                                                                                                                                                                    cGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGAGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGAGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATcgccg                                        >  1:136477/1‑251 (MQ=255)                                                                                                                                                                                                        cgcgcgCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTg                                     >  1:104734/1‑250 (MQ=255)                                                                                                                                                                                                             gcTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACGCAATTAACGAAAACGAGTAGTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGGCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGc                               <  1:42815/251‑1 (MQ=255)                                                                                                                                                                                                                gCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCt                             >  1:11255/1‑250 (MQ=255)                                                                                                                                                                                                                       ggTGAAGCGCCAAAGATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCt                     <  1:188594/251‑1 (MQ=255)                                                                                                                                                                                                                                       tggtggTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAACGAATTGatc     >  1:220316/1‑251 (MQ=255)                                                                                                                                                                                                                                           ggtCTGCTGGGGCGGTCACTCAATTAACGCAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACACAGACAGACGTTTTCTTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCACAGAATTGatcatc  >  2:186308/1‑250 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              GCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAACGAATTGATCATC  >  NC_000913/2926526‑2927009

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