Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	113,695	G→A	W84* (TGG→TGA)	guaC →	GMP reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	113,695	0	G	A	96.0%	78.0 / NA	26	W84* (TGG→TGA)	guaC	GMP reductase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (13/11);  total (13/13)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTATTGAAGAAGATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACGCTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTG‑CATAAACACTATTCTGTCGAAGAGTGGCAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTGACTGGTGAAATGTGT  >  NC_000913/113448‑113938                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    ggaTTGAAGAAGATCTGAAGTTAGTTTTTAAAGACGTGCTCATCCGCCCTAAACGCTCCACGCTTAAAAGCCGTGCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATGGCCGCAAACATGGACACCGTAGGCACAGTTTCTATGGCCTCGGCGTTGGCTGCTTGTGATATTTTGAGTGC‑GTGGAATAAACACTGTTCTGCTTAAGAGTGGCaa                                                                                                                                                                                                                                                  <  2:304550/249‑1 (MQ=255)                 taaGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACGCTCCCCTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCATGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTTGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGTCTGCTGTG‑CATAAACCCTATTCTGTCGAAGAGTGACAAGCGTTTATcaaca                                                                                                                                                                                                                                     <  2:281090/247‑1 (MQ=255)                 gAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACGCTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTG‑CATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATcaaca                                                                                                                                                                                                                                     <  2:281098/248‑1 (MQ=255)                   aGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACGCTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTG‑CATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAAttct                                                                                                                                                                                                                                <  1:36189/251‑1 (MQ=255)                                                   aaaCGCTCCACTCTTAAAAGCCGGTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGGGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTG‑CATAAACACTATTCGGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGtt                                                                                                                                                                                                  <  2:297687/249‑1 (MQ=255)                                                                   aaaGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTG‑CATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGAt                                                                                                                                                                                <  1:69227/251‑1 (MQ=255)                                                                       ccGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACGCCGTATGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTG‑CATAAACACTATGCTGTCGAAGACTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTTAAACATGTGATGGTTTCTACCGGTCCGTCTGATGctg                                                                                                                                                                            <  2:149417/251‑3 (MQ=255)                                                                                          gAACGGCAATTCACCATCAAACAATCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGGCTGCTGTG‑CATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACa                                                                                                                                                         <  2:50968/251‑1 (MQ=255)                                                                                                 aaTTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTG‑CATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTAGCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATtctc                                                                                                                                                  <  1:139495/251‑1 (MQ=255)                                                                                                                 ttCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTG‑CATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCAt                                                                                                                                  >  1:153768/1‑251 (MQ=255)                                                                                                                     ggTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTG‑CATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTaa                                                                                                                               >  1:306443/1‑250 (MQ=255)                                                                                                                                               tATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTG‑CATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATg                                                                                                    >  2:290795/1‑251 (MQ=255)                                                                                                                                                                  acCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTG‑CATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCgg                                                                                 >  2:232980/1‑250 (MQ=255)                                                                                                                                                                  acCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTG‑CATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGt                                                                                 >  2:205026/1‑251 (MQ=255)                                                                                                                                                                        gggccaTTTTCTATGGCCTATGCGCTGGCTTCTTTTGATATTTTGACTGCTGTG‑CATAAACACTATTCTGTCGAAGAGTGCCAAGCGTTTACCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGTATTTCGAAAAAACTAAACAGATTCTAGACCTGAGCACCGCATTAAACTTCGTTTGTATTCACGTGGCGAATGGTTATACCGAATACTTCGTGCAGtt                                                                           <  2:165974/247‑1 (MQ=255)                                                                                                                                                                            caTTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTG‑CATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGtt                                                                       <  1:215624/251‑1 (MQ=255)                                                                                                                                                                                                        ttttGATATTTTGACTGCTGTG‑CATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCATTTCGTTGCGAAAGCGCGTGAAGAGTGGCCGAcc                                            >  1:213126/1‑250 (MQ=255)                                                                                                                                                                                                                                aTAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGt                     >  2:4318/1‑250 (MQ=255)                                                                                                                                                                                                                                   aaCACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGGACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACGTCGTTTGTATTGTCGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATGTGTGCTGGTAAGGTAGTg                 <  2:196144/251‑1 (MQ=255)                                                                                                                                                                                                                                    acacTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTg                 >  1:8691/1‑250 (MQ=255)                                                                                                                                                                                                                                      acTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTGACt              >  2:117289/1‑251 (MQ=255)                                                                                                                                                                                                                                          ttCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTGACTGGTg          >  1:210869/1‑251 (MQ=255)                                                                                                                                                                                                                                                cGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTAATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTGACTGGTGAAAtgt    >  2:157438/1‑251 (MQ=255)                                                                                                                                                                                                                                                  aaGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTGACTGGTGCAAtgtgt  >  2:1708/1‑251 (MQ=255)                                                                                                                                                                                                                                                     agTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGt                          <  2:277638/224‑1 (MQ=255)                                                                                                                                                                                                                                                     agTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGt                          >  1:277638/1‑224 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    GTATTGAAGAAGATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACGCTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTG‑CATAAACACTATTCTGTCGAAGAGTGGCAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTGACTGGTGAAATGTGT  >  NC_000913/113448‑113938

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