Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,828,654	C→T	G361R (GGG→AGG)	astD ←	aldehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,828,654	0	C	T	90.0%	52.0 / ‑5.7	20	G361R (GGG→AGG)	astD	aldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (13/5);  total (13/7)
Fisher's exact test for biased strand distribution p-value = 1.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGTTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAG‑TTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAAT  >  NC_000913/1828417‑1828886                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                          cgGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGTTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCTCGGCGTCAGCAAc                                                                                                                                                                                                                              >  1:290476/1‑251 (MQ=255)                                            ccccGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCTCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGgtcgtc                                                                                                                                                                                   <  2:293002/251‑1 (MQ=255)                                             cccgcccgCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCTCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTcc                                                                                                                                                                                  <  2:74081/251‑1 (MQ=255)                                                                                     gCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCTCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAG‑TTGCTGCCATGCAGTAACCACctgc                                                                                                                                         >  1:28016/1‑251 (MQ=255)                                                                                               aaCCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCTCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAG‑TTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCt                                                                                                                               >  2:92439/1‑251 (MQ=255)                                                                                                 ccAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCTCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAG‑TTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGt                                                                                                                             <  1:179899/251‑1 (MQ=255)                                                                                                                        cGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCTCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAG‑TTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATaa                                                                                                       >  2:87424/1‑250 (MQ=255)                                                                                                                                 ttCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCTCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAG‑TTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCgg                                                                                             >  1:42105/1‑251 (MQ=255)                                                                                                                                             atagCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGGCCGAACACATATTCATCTGGTACGCCAGAAAAGCCTGTCATTTCAATGATCCTCGGCGTCAGCAACGATGGCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGATTCCAG‑TTGCTGACATGAAGTAACCACCTGCTGTGCGGCCTGTTAAGAACTCAGACCGCCAATAAACGGCTGCGGTTCGTCATccc                                                                                  >  1:243104/4‑250 (MQ=255)                                                                                                                                             aTCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCTCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAG‑TTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATcccc                                                                                 >  1:243098/1‑250 (MQ=255)                                                                                                                                               cGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCTCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAG‑TTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTTATCCAAGt                                                                               >  2:247322/1‑251 (MQ=255)                                                                                                                                                       cGAAAGCATCATCAAGCCAGACGCACAATAACGGTCAGAACACCTATTCATCAGGTCCGCAAGCCACGACTGTCATTTCAATAATCCTCGGCGTCAGAAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCAATCACTTCCAG‑TTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCAAATAAACGGCTGAGGTTCGTCATCCCAGTTGCCCGGc                                                                       >  2:295812/1‑251 (MQ=255)                                                                                                                                                                      gCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGAACGCCAGCAACGCCTGTCATTTAAATGATCCTAGGCGTCAGCAACGATGTACCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAG‑TTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCgctggctg                                                        >  1:282242/1‑251 (MQ=255)                                                                                                                                                                                    aaCGGTCCGACCCCCTCTTCATCGGGCACGCCACCAACCCCTGTCATTTCAATGACCCTCGGCGTCAGCAACGATGTCCCTGCTGTTAAAAGGCGCGGCGCAAGAGGGGGCCGCCCGCCAATCGCT‑CCAGGTTGCTCCCAGGCAGTACCCCCCTGCGGTGCGGCCTGTTCAAAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGAc                                                      <  2:289526/239‑1 (MQ=255)                                                                                                                                                                                    aaCGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCTCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAG‑TTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGAc                                                      >  1:289526/1‑239 (MQ=255)                                                                                                                                                                                        gTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCTCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAG‑TTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGc                                       >  2:47561/1‑250 (MQ=255)                                                                                                                                                                                                acCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCTCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAG‑TTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGc                              >  2:71920/1‑251 (MQ=255)                                                                                                                                                                                                      tcagctggaaccccgcaaccGCCTGTCATTTCAAGGACCCCGGGCGTAAGCAACGTTGCCCCTGCTTGTAAAAAGCGGGGCGAAGGCAGGGGCCGCCCGCCCATCGCTCCCAG‑TTCCGGCCATGCAGTACCCCCCTGCTGTGGGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATc                           <  2:164628/229‑1 (MQ=255)                                                                                                                                                                                                        accTGGTCCGCCAGCAACGCCGGTCATTTCAATGATCCTCGGCGTCAGCAACGATGTCCCTGTTTGAAATAAGCGCGGCGCAAGCGGGGGTCGCCCGCCCATCGCTTCCAG‑TTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCt                      <  1:71920/249‑1 (MQ=255)                                                                                                                                                                                                                              gttaTTTAAATGATCCCCGGCGTCAGCAACTATGTCCCTGCTTGAAAAAAGCGCGGCGAAAGGAGGGGTCTCCCGCCCACCCCTCCCAG‑TTGCGGCAATGCGGTACCCCCCTGCTCGGCGGCCTGTTAAAAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAAt  <  1:92439/246‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                          CGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGTTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAG‑TTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCAAT  >  NC_000913/1828417‑1828886

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