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	96.0%	73.0 / NA	25	G155G (GGC→GGA)	ppnN	nucleotide 5'‑monophosphate nucleosidase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (15/9);  total (16/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACTGATTAATCCCCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAACGAATT  >  NC_000913/2926530‑2927002                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                       gAACTGATTAATCCCCCGGAAGAGGCTTTTGCCGAAGGGCGAATTATTCACACTTTGAAGGCCAAATTGTTAGCGGTCTTGAGTGCCATTCTCTTCGTTTACGGGCAAATACATAACACAGTTAGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCAACTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATCTGGTAGTCTGCTGGGGAGGTAACTc                                                                                                                                                                                                                                >  1:127262/1‑251 (MQ=255)         tAATCCCCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAAc                                                                                                                                                                                                                         <  1:207694/250‑1 (MQ=255)             cccccGGCAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGaaa                                                                                                                                                                                                                     >  2:110035/1‑250 (MQ=255)                 cGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAg                                                                                                                                                                                                                <  2:299669/251‑1 (MQ=255)                                                                            gTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCAc                                                                                                                                                     >  1:170111/1‑251 (MQ=255)                                                                                                  tttACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATg                                                                                                                               >  1:250515/1‑251 (MQ=255)                                                                                                  tttACGGGCAAATCCATAACACCGTTCGTTTTCCCAACATAAGTATAGACAACCCCGTCCACATCACTACCCTGGTCTCTTCCATCTTCCGTAAAGCTAGCGCGCTGCATGTGGGTGAAGCGCCCAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTACGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATg                                                                                                                               >  1:304176/1‑251 (MQ=255)                                                                                                                        cgttcgttTTCCCAACCTGAATAGCGACATCTCCGTCAACATCACTAACCTGGTCTTTTCCATCTTGCGTACCGCTCGCGCGCTGCATGTGGGTGAAGCGCCCAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGAATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCTGCTGAAGTCCGGGAGCGATGGAATCGCCGCTGACAGGagctg                                                                                                         >  2:125633/1‑251 (MQ=255)                                                                                                                          ttcgttTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCg                                                                                                       >  2:220507/1‑251 (MQ=255)                                                                                                                                   ccAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGacac                                                                                              <  1:80362/251‑1 (MQ=255)                                                                                                                                             tctcGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCagcagc                                                                                    >  1:157095/1‑251 (MQ=255)                                                                                                                                              ctcGACAACTCAGTCCACATCACTAACCTGGTTTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGAGGTCTGCTGGGGAGGTCACTCAATTAACGCAGACGAGTATTTGTATGCCCGTAGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGACCGGGAGCGGTGTAAGCGCAGATGAAAGGTGCTGCGCACGGACACACGCAGAAGCg                                                                                   >  2:83099/1‑251 (MQ=255)                                                                                                                                                cGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTGGCGTACCGCTCGCGCGCTGCATGTGGGGGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGTATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCagcagc                                                                                    <  2:158601/248‑1 (MQ=255)                                                                                                                                                 gACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTa                                                                                <  2:59486/251‑1 (MQ=255)                                                                                                                                                        ccGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGAc                                                                         <  1:95919/251‑1 (MQ=255)                                                                                                                                                         cGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGAc                                                                         >  2:86321/1‑250 (MQ=255)                                                                                                                                                                                         tGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGACAAATATGGTGGTCTGCTGGGGAGATCACTCAATTAACGAAAACGAGTATTTGGATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATATGCACCGGCTGTGGACCGGGAGCGATGGAAGAGACGCTGAAGGGAGCTGCGGACGGACACGCACAGAAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAACGCCGCCGAtt                                        >  2:115029/1‑251 (MQ=255)                                                                                                                                                                                           cGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTAt                                      >  1:152510/1‑251 (MQ=255)                                                                                                                                                                                              aaCGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCAc                                                                                                                                                     >  2:287876/1‑137 (MQ=255)                                                                                                                                                                                              aaCGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCAc                                                                                                                                                     <  1:287876/137‑1 (MQ=255)                                                                                                                                                                                               aCGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATAATcgcc                                  <  1:86321/251‑1 (MQ=255)                                                                                                                                                                                                           tGCATGTGGGTGAAGCGCAAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGGTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGCTTATCGCCGCTGAACCGCCt                      >  2:216730/1‑251 (MQ=255)                                                                                                                                                                                                              aTGTGGGCGAAGCGCCAAATATGATGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGAGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAAc                   >  1:89157/1‑251 (MQ=255)                                                                                                                                                                                                                              aaaTATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAAACGCCTAACCCGCTGGTCAACGAAt   >  1:23097/1‑251 (MQ=255)                                                                                                                                                                                                                               actatgtgtgtatgcgggggAGGTCACTCAATTACCGAAAACGAGTATGTGTATGCCCGTCGCGTCGAAAACCAGCGGGGCCTGCGTGAGCTGAATATCTGCACCGGCGGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCGGCGGTCGGACACGCGCAGCAGCGTTACAAAGACATACGTTTTATTGTTAGGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAACGAAtt  <  2:261497/235‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                       GAACTGATTAATCCCCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAACGAATT  >  NC_000913/2926530‑2927002

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