Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	984,504	C→T	intergenic (+207/+15)	gloC → / ← aspC	hydroxyacylglutathione hydrolase GloC/aspartate aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	984,504	0	C	T	92.6%	73.3 / ‑0.7	27	intergenic (+207/+15)	gloC/aspC	hydroxyacylglutathione hydrolase GloC/aspartate aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (16/9);  total (18/9)
Fisher's exact test for biased strand distribution p-value = 5.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.13e-01

GCATAATCCCTTCCTGCAAGACGAAATGCCCGTCTGGTAAGGCACATAAAAAAGCCCGCTTTTAATGCTGGCCTGGATTTCTGGCAAAGTGCGCTTTGTTTATGCCGGATGCGGCACGAGCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAAATATTGCAGGGGACGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTCTCAGCCCGCTTTTCAGCGGGCTTCATTGTTTTTAATGCTTACAGCACTGCCACAATCGCTTCGCACAGCGGAGCCATGTTATCTGGTGTCATCCCGGCCACATTTACGCGACCAGAAGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAAC  >  NC_000913/984258‑984738                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                           gCATAATCCCTGCCTGCAAGACTAAATGCCCGTTTGGTAAGGCACATAAAAAAGCCCGCTTTTAATGCTGGCCTGGATTTCTGGCAAAGTGCGCTTTGTTTAGGCCGGATGCGGCACGAGCGCCGTATCCGGCCTACAAAATCGTGCAAATTCAAAATATTGCAGGGGACGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTCTCAGCCCGCTTTTCAGCGGGCTTTATTg                                                                                                                                                                                                                                        <  2:102996/251‑1 (MQ=255)                      cGAAATGCCCGTCTGGTAAGGCACATAAAAAAGCCCGCTTTTAATGCTGGCCTGGATTTCTGGCAAAGTGCGCTTTGTTTATGCCGGATGCGGCACGAGCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAAATATTGCAGGGGACGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTCTCAGCCCGCTTTTCAGCGGGCTTTATTGTTTTTAATGCTTACAGCACTg                                                                                                                                                                                                                   >  2:148262/1‑251 (MQ=255)                                    ggTAAGGCACATAAAAAAGCCCGCTTTTAATGCTGGCCTGGATTTCTGGCAAAGTGCGCTTTGTTTATGCCGGATGCGGCACGAGCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAAATATTGCAGGGGACGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTCTCAGCCCGCTTTTCAGCGGGCTTTATTGTTTTTAATGCTTACAGCACTGCCACAATCGCTTCg                                                                                                                                                                                                     <  1:42137/251‑1 (MQ=255)                                       aaGGCACATAAAAAAGCCCGCTTTTAATGCTGGCCTGGATTTCTGGCAAAGTGCGCTTTGTTTATGCCGGATGCGGCACGAGCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAAATATTGCAGGGGACGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTCTCAGCCCGCTTTTCAGCGGGCTTTATTGTTTTTAATGCTTACAGCACTGCAACAATCGCTTCGcac                                                                                                                                                                                                  >  1:212162/1‑251 (MQ=255)                                            acaTAAAAAAGCCCGCTTTTAATGCTGGCCTGGATTTCTGGCAAAGTGCGCTTTGTTTATGCCGGATGCGGCACGAGCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAAATATTGCAGGGGACGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTCTCAGCCCGCTTTTCAGCGGGCTTTTTTGTTTTTAATGCTTACAGCACTGCCACAATCGCTTCGCACAGCgg                                                                                                                                                                                             >  2:140981/1‑251 (MQ=255)                                                       cccGCTTTTAATGCTGGCCTGGATTTCTGGCAAAGTGCGCTTTGTTTATGCCGGATGCGGCACGAGCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAAATATTGCAGGGGACGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTCTCAGCCCGCTTTTCAGCGGGCTTTATTGTTTTTAATGCTTACAGCACTGCCACAATCGCTTCGCACAGCGGAGCCATGTTAt                                                                                                                                                                                  >  2:143482/1‑251 (MQ=255)                                                            gtttAATGCTGGCCTGGATTTCTGGCAAAGGGCGCTTTGTTTATGCCGGATGCGGCACGAGCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAAATATTGCAGGGGACGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTCTCAGCCCGCTTTTCAGCGGGCTTTATTGTTTTTAATGCTTACAGCGCTGCCACAATCGCTTCGCACAGCGGAGCCATGTTATCTgg                                                                                                                                                                              <  2:47589/249‑1 (MQ=255)                                                                           ggATTTCTGGCAAAGTGCGCTTTGTTTATGCCGGATGCGGCACGAGCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAAATATTGCAGGGGACGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTCTCAGCCCGCTTTTCAGCGGGCTTTATTGTTTTTAATGCTTACAGCACTGCCACAATCGCTTCGCACAGCGGAGCCATGTTATCTGGTGTCATCCCGGCcaca                                                                                                                                                              <  1:242695/251‑1 (MQ=255)                                                                             aTTTCTGGCAAAGTGCGCTTTGTTTATGCCGGATGCGGCACGAGCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAAATATTGCAGGGGACGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTCTCAGCCCGCTTTTCAGCGGGCTTCATTGTTTTTAATGCTTACAGCACTGCCACAATCGCTTCGCACAGCGGAGCCATGTTATCTGGTGTCATCCCGGCCACAtt                                                                                                                                                            >  1:160314/1‑251 (MQ=255)                                                                                     cAAAGTGCGCTTTGTTTATGCCGGATGCGGCACGAGCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAAATATTGCAGGGGACGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTCTCAGCCCGCTTTTCAGCGGGCTTTATTGTTTTTAATGCTTACAGCACTGCCACAATCGCTTCGCACAGCGGAGCCATGTTATCTGGTGTCATCCCGGCCACATTTACGCGAc                                                                                                                                                    >  1:323352/1‑251 (MQ=255)                                                                                      aaaGTGCGCTTTGTTTATGCCGGATGCGGCACGAGCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAAATATTGCAGGGGACGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTCTCAGCCCGCTTTTCAGCGGGCTTTATTGTTTTTAATGCTTACAGCACTGCCACAATCGCTTCGCACAGCGGAGCCATGTTATCTGGTGTCATCCCGGCCACATTTACGCGAcc                                                                                                                                                   >  1:307931/1‑251 (MQ=255)                                                                                             gcTTTGTTTATGCCGGATGCGGCACGAGCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAAATATTGCAGGGGACGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTCTCAGCCCGCTTTTCAGCGGGCTTTATTGTTTTTAATGCTTACAGCACTGCCACAATCGCTTCGCACAGCGGAGCCATGTTATCTGGTGTCATCCCGGCCACATTTACGCGACCAGAAGc                                                                                                                                             >  1:145994/1‑250 (MQ=255)                                                                                               tttGTTTATGCCGGATGCGGCACGAGCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAAATATTGCAGGGGACGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTCTCAGCCCGCTTTTCAGCGGGCTTTATTGTTTTTAATGCTTACAGCACTGCCACAATCGCTTCGCACAGCGGAGCCATGTTATCTGGTGTCATCCCGGCCACATTTACGCGACCAGAAGCAAc                                                                                                                                          >  1:406503/1‑251 (MQ=255)                                                                                                           ggATGCGACACGAGTGCCTGACCCGGCCTACAAAATCGTGCAAATGAAAAATATTGAAGGGGACGCGTTAGCCGGTTAAGCGGAGCGCATCAGGCATTGTGGCGTTTGTCATCATTCGCAGCCCGCTGTTCAGCGGGCTTTATTGTTTTTAATTCTTACAGCTCTGCTACAACCGCTTCGTACAGCGGGGCGATGTTATCTGGTTTCATCCCGGCCACATTTACGCGACCAGAAGCAACCGCATATACGcc                                                                                                                              <  2:40105/251‑1 (MQ=255)                                                                                                            gATGCGGCACGAGCGCCTTATCCTGCCTACAAAATCGTGCAAATTCAAAATATTGCAGGGGACGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTCTCAGCCCTCTTTTCAGCGGGCTTTATTGTTTTTAATGCTTACAGCACTGCCACAATCGCTTCGCACAGCGGAGCCATGTTATCTGGTGTCATCCCGGCCACATTTACGCGACCAGAAGCAACCGCATATACGCCa                                                                                                                             <  1:226543/251‑1 (MQ=255)                                                                                                                    aCGAGCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAAATATTGCAGGGGACGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTCTCAGCCCGCTTTTCAGCGGGCTTTATTGTTTTTAATGCTTACAGCACTGCCACAATCGCTTCGCACAGCGGAGCCATGTTATCTGGTGTCATCCCGGCCACATTTACGCGACCAGAAGCAACCGCATATACGCCAAACTCTTc                                                                                                                     <  2:190643/251‑1 (MQ=255)                                                                                                                       aGCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAAATATTGCAGGGGACGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTCTCAGCCCGCTTTTCAGCGGGCTTTATTGTTTTTAATGCTTACAGCACTGCCACAATCGCTTCGCACAGCGGAGCCATGTTATCTGGTGTCATCCCGGCCACATTTACGCGACCAGAAGCAACCGCATATACGCCAAACTCTTcgcc                                                                                                                  >  2:355640/1‑250 (MQ=255)                                                                                                                                                         cAAAATATTGCAGGGGACGCGTATGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTCTCAGCCCGCTTTTCAGCGGGCTTTATTGTTTTTAATGCTTACAGCACTGCCACAATCGCTTCGCACAGCGGAGCCATGTTATCTGGTGTCATCCCGGCCACATTTACGCGACCAGAAGCAACCGCATATACGCCAAACTCTTCGCGCAGGCGCAGCACTTGTTCTTTTGTCAGGCCACt                                                                                 <  1:342993/250‑1 (MQ=255)                                                                                                                                                              tatTGCAGGGGACGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTCTCAGCACGCTTTTCAGCGGGCTTTATTGTTTTTAATGCTTACAGCACTGCCACAATCGCTTCGCACAGCGGAGCCATGTTATCTGGTGTCATCCCGGCCACATTTACGCGACCAGAAGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAggg                                                                           >  1:51553/1‑250 (MQ=255)                                                                                                                                                                   cAGGGGACGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTCTCAGCCCGCTTTTCAGCGGGCTTTATTGTTTTTAATGCTTACAGCACTGCCACAATCGCTTCGCACAGCGGAGCCATGTTATCTGGTGTCATCCCGGCCACATTTACGCGACCAGAAGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCCGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAAc                                                                       >  2:153904/1‑250 (MQ=255)                                                                                                                                                                     ggggACGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTCTCAGCCCGCTTTTCAGCGGGCTTTATTGTTTTTAATGCTTACAGCACTGCCACAATCGCTTCGCACAGCGGAGCCATGTTATCTGGTGTCATCCCGGCCACATTTACGCGACCAGAAGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATg                                                                    >  1:138376/1‑251 (MQ=255)                                                                                                                                                                                    tGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTCTCAGCCCGCTTTTCAGCGGGCTTCATTGTTTTTAATGCTTACAGCACTGCCACAATCGCTTCGCACAGCGGAGCCATGTTATCTGGTGTCATCCCGGCCACATTTACGCGACCAGAAGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTtgatg                                                     >  1:24801/1‑251 (MQ=255)                                                                                                                                                                                      aTAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTCTCAGCCCGCTTTTCAGCGGGCTTTATTGTTTTTAATGCTTACAGCACTGCCACAATCGCTTCGCACAGCGGAGCCATGTTATCTGGTGTCATCCCGGCCACATTTACGCGACCAGAAGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTgatgat                                                   >  2:44042/1‑251 (MQ=255)                                                                                                                                                                                          gCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTCTCAGCCCGCTTTTCAGCGGGCTTTATTGTTTTTAATGCTTACAGCACTGCCACAATCGCTTCGCACAGCGGAGCCATGTTATCTGGTGTCATCCCGGCCACATTTACGCGACCAGAAGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATaa                                                 >  1:358556/1‑249 (MQ=255)                                                                                                                                                                                                     cAGGCAATGTTGCGTTTGTCATCAGTCTCAGCCCGCTTTTCAGCGGGCTTTATTGTTTTTAATGCTTACAGCACTGCCACAATCGCTTCGCACAGCGGAGCCATGTTATCTGGTGTCATCCCGGCCACATTTACGCGACCAGAAGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACAt                                                                     >  1:313618/1‑218 (MQ=255)                                                                                                                                                                                                                                   ccccGCTTTCCAGCGGGTTTTATTGTTTTTAATGCTTACAGCACTGCCACAAGCGCTTCGCACAGCGGAGCCATGTTATCTGGTGGCATCCCGGCCACATTTACGCGACCAGAAGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGAc      <  2:323352/250‑1 (MQ=255)                                                                                                                                                                                                                                        cTTTTCAGCGGGCTTTATTGTTTTTAATGCTTACAGCACTGCCACAATCGCTTCGCACAGCGGAGCCATGTTATCTGGTGTCATCCCGGCCACATTTACGCGACCAGAAGCAACCGCATATACGCCCAACTCTTCGCGCAGCCGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGACAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGaac  >  2:336407/1‑250 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                           GCATAATCCCTTCCTGCAAGACGAAATGCCCGTCTGGTAAGGCACATAAAAAAGCCCGCTTTTAATGCTGGCCTGGATTTCTGGCAAAGTGCGCTTTGTTTATGCCGGATGCGGCACGAGCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAAATATTGCAGGGGACGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTCTCAGCCCGCTTTTCAGCGGGCTTCATTGTTTTTAATGCTTACAGCACTGCCACAATCGCTTCGCACAGCGGAGCCATGTTATCTGGTGTCATCCCGGCCACATTTACGCGACCAGAAGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAAC  >  NC_000913/984258‑984738

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