Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	65,454	A→G	R109R (CGT→CGC)	polB ←	DNA polymerase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	65,454	0	A	G	98.5%	196.0 / NA	66	R109R (CGT→CGC)	polB	DNA polymerase II
Reads supporting (aligned to +/- strand):  ref base A (1/0);  major base G (22/42);  minor base T (0/1);  total (23/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

CAAGCGAGGAGGCGTCGCCATTCTCCGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCTATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGTGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑ACGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGC‑TGGCGGTG‑AAAA‑TCCTTTAACGCCAGCGG‑TGTCAGGCGAAAGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGG‑AATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTGCAA  >  NC_000913/65207‑65687                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                cgaGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGGTTTAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGGGAGGTGGTAAACCGCTCCATCAGATAGCGGTCTGGCGG‑gcgc                                                                                                                                                                                                                                              >  1:119869/3‑251 (MQ=255)   ggcgcggaggCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGc‑ac                                                                                                                                                                                                                                           <  1:316419/246‑1 (MQ=255)   aGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTTAATATCAAAAAAAACCCACTTGAGCGGCGGACGGGAGTCGGGGAGCGGGTTCTGAAGGGCGTTAACGATAACGCCC‑TTTTGGATATCCCCCCCGCCCCCCACC‑CGGGGGGGGAAAAAACGCTCCCTCAGATAACGTTCTGGCGG‑agcgcc                                                                                                                                                                                                                                            >  2:675583/1‑246 (MQ=255)     ggcggatgtgTCGCCATTCTCTGGCCCAAGCAGATAAACGATGCTTGGCCCGCAGCCGACCAGGCCGAGGAAGTGCATCTCACCGGGGCGGGTGGTTTCAATATCAATAGAAACCCGCGTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGGGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACAt                                                                                                                                                                                                                                         <  2:429507/242‑1 (MQ=255)      gcggaggCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATc                                                                                                                                                                                                                                        <  1:99061/249‑1 (MQ=255)            gCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTc                                                                                                                                                                                                                                  <  1:483219/251‑1 (MQ=255)              gTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCCATATCTATTGAAACCCACTTGAGCGGCGGACGCTAGTCGGGATGCGGTTTCAGTCCGGCCTTCACGATAGCGCCC‑TTGTGCGTAACACCCTCTACCCACCCC‑GGCGGGGTGATAAACCGCTCCCTCAGATAGCGTTCTCGGGG‑GCGC‑CCATCGGCCTCgg                                                                                                                                                                                                                                >  1:9722/1‑250 (MQ=255)               gcgtcattcgcgggcccaatcATATTAACGAGGCGCTGCCAGCGGACTACCAGGCCGATGCTGTACAGCTCACCGTGGCGGGTGGTGTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGTGATGCGGTTTCAGACGGGCATGAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGCGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGt                                                                                                                                                                                                                                <  2:634338/230‑1 (MQ=255)                cGCCCTTCTCTGGCCCCAGCATATAAACGTAGCGATGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATACCACCCGCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGt                                                                                                                                                                                                                                <  2:255463/249‑1 (MQ=255)                   cATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACg                                                                                                                                                                                                                           <  1:640415/251‑1 (MQ=255)                         ctggCCCCAGCATATAAACTAGGCGTTTCCCGCATGCTTCCAGGCCCATGCAGTACAGCTCACCGTGGCGGGTGGTTTTAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCAGTTAT‑CATATCACCCTCGATCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTa                                                                                                                                                                                                                        <  2:702365/246‑1 (MQ=255)                              cccAGCATATAAATGTTGCGCGGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCCGCGGGGCGGGTGGTTTCAATATCAAGAGAATCCCACTTGAGCGGCGGACGATAGTCTGGATGCGGTTTCAGACGGGCATTGACGATAGCGCCA‑TTGTGCATAGCACCCTCGACCCACACC‑GGTGGGGTGATAGACCGCTGCATCAGATAGCGCTCTGGCGG‑TGGC‑ACATCGGCCTCGTAGACTGTAACGCCAcc                                                                                                                                                                                                                <  2:637742/251‑1 (MQ=255)                                  gCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCa                                                                                                                                                                                                            <  1:278024/251‑1 (MQ=255)                                   cATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCAc                                                                                                                                                                                                           <  1:509259/251‑1 (MQ=255)                                    atatAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAAACCACTTGAGCGGCGGACGATAGTCGGGATGCGGGTTCAGAAGGGCCTTAACGATAACGCCA‑TTGTGCATATCACCCTCGACCC‑CACCCGGGGGGGGGATAAACCGCTCCCTCAGGAAGCGTTCTGGCGG‑GGGC‑ACATAGGCCCCCTTAACGGTAACGCCACCtttccc                                                                                                                                                                                                          >  1:701050/1‑247 (MQ=255)                                      ataAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCa                                                                                                                                                                                                        <  1:458267/251‑1 (MQ=255)                                               gcgcTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCtttt                                                                                                                                                                                               <  1:485439/251‑1 (MQ=255)                                                    gCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTGACCGTGGCGGGTGGTTTCAATATCAAGAGAAACCCACTTGAGCGGCGGACGATAGGCGGGATGCGGTTTCAGACGTGCATTAACGATAGCGCCA‑TGGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCGCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCGTCACGCAGGCGCTTTTCGt                                                                                                                                                                                            <  2:351166/249‑1 (MQ=255)                                                        gcggcattCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGGGTCAATATCTATAGAAACCCACTTGCGCGGCGGGCGATAGTCGGGAGGCGGGTTAAGACGGGCCATAACGAGAGTGCCC‑TTGTGAATGTCACCCTCGGACCACAAC‑GGGGGGGTGGTAAACCCCCACATCCGGAAACGTTCTTGCGGGAC‑C‑ACACGGGCCCCCTCGACGGTAACGCCAACTTTACGCACGCGCttttgagtaattt                                                                                                                                                                                      >  2:206250/7‑242 (MQ=255)                                                                ccAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTg                                                                                                                                                                               <  1:405540/250‑1 (MQ=255)                                                                ccAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTg                                                                                                                                                                               <  1:167115/250‑1 (MQ=255)                                                                  gggCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCg                                                                                                                                                                            <  2:705545/250‑1 (MQ=255)                                                                     ccGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATg                                                                                                                                                                         <  1:675583/251‑1 (MQ=255)                                                                                 aGCTCACCGTGGCGGGTGGTTTCAATAACAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTa                                                                                                                                                             <  1:258972/251‑1 (MQ=255)                                                                                     cACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTaaa                                                                                                                                                           <  1:263617/249‑1 (MQ=255)                                                                                      cccGTGGCGGGTGGTTTCAATAGCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGc                                                                                                                                                        <  1:212348/250‑1 (MQ=255)                                                                                      cccGTGGCGGGTGGATTCAATATCAATCGAAACCCACTTGAGCGGCGGACGATAGTCGTGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCGCGACCCACACC‑GGTGAGGGGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGGAGACGGTAACGCCACCTTCACGCAGGCGCTGTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGc                                                                                                                                                        <  2:479084/250‑1 (MQ=255)                                                                                        cGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCa                                                                                                                                                      <  1:392874/251‑1 (MQ=255)                                                                                                   gtTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCTGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCTCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGC‑Tg                                                                                                                                          <  2:36639/251‑1 (MQ=255)                                                                                                     ttCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCCGC‑TGGc                                                                                                                                        >  1:642351/1‑251 (MQ=255)                                                                                                       cAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACCTCGGCCTCGTAGACGGTAACGCACCCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGCTGGGCGCGCCAGTAACAGCCACACACCGGcctgcgg                                                                                                                                       >  2:154289/1‑245 (MQ=255)                                                                                                        aaTATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGC‑TGGCGGt                                                                                                                                     <  1:465335/251‑1 (MQ=255)                                                                                                                   aaaCCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑AAATCGGCCTCGGCGAAGGTAACCCCACCTTCACGCCGGCGCTTTTTGGTATTTATCCATTTGCGATGGGCGAGACAATAAGGACCATACAACGGCGTGG‑GGTG‑aaaaaccctt                                                                                                                         >  2:89367/1‑245 (MQ=255)                                                                                                                    aaCCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCTGCCTCGTAGACGGTAACGCCCCCTTCACGCCGGCGCTTTTCGTAAATCATCAATTGGCGATGGGCGCGCCACGAAAGACCATACACCCgc‑gggcgggggaaaacactat                                                                                                                         >  2:166601/1‑238 (MQ=255)                                                                                                                    aaCCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGC‑TGGCGGTG‑AAAA‑TCCTTTa                                                                                                                       <  1:333534/251‑1 (MQ=255)                                                                                                                        cACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCCCGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGC‑TGGCGGTG‑GAAA‑ACCTTTAACGc                                                                                                                   >  1:138616/1‑251 (MQ=255)                                                                                                                            tGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGC‑TGGCGGTG‑AAAA‑TCCTTTAAAGCCAGc                                                                                                               >  1:75618/1‑251 (MQ=255)                                                                                                                            ggAGCGGCGGACGATAGTCGGGATGCGGTCTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCGCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGC‑TGGCGGTG‑AAAA‑TCCTTTAACGCAAGc                                                                                                               <  2:163984/250‑1 (MQ=255)                                                                                                                              aGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGCCGGTAACGCCACCTTCACGCAGGCGCCTTTCGTCATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCCTTCAACGGC‑TGGCGGGG‑GAAA‑TACTTTatacccaacgg                                                                                                             >  1:8225/1‑241 (MQ=255)                                                                                                                              aGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCAAACC‑GGAGAGGTGAAAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGCCACA‑CGGCCTCGTAGACGGAAACGCCACCTTCACGCAGGCGCTTTTTGTAATTTATCAAATGGGGATGAGGGGGACCATAAAGGCCCTACACCGGC‑GGTCGGGG‑aaaagactttttaaacagc                                                                                                                >  2:194059/1‑234 (MQ=255)                                                                                                                               gcggcgGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGC‑TGGCGGTG‑AAAA‑TCCTTTAACGCCAGCGg‑t                                                                                                           <  1:438738/251‑1 (MQ=255)                                                                                                                                cggcggACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGC‑TGGCGGTG‑GAAA‑TCCTTTAACGCCAGCgggg                                                                                                           >  1:114550/1‑249 (MQ=255)                                                                                                                                  gcggACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGC‑TGGCGGTG‑AAAA‑TCCTTTTACGCCAGCGG‑GGTc                                                                                                        >  1:379925/1‑251 (MQ=255)                                                                                                                                    ggACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGC‑TGGCGGGG‑AAAA‑TCCTTTAACGCCAGCGG‑TGTccg                                                                                                      >  1:196205/1‑249 (MQ=255)                                                                                                                                     gACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGCTAGCGCCA‑TTGTGCATATCACCCTCGACCCGCACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCGCGTCGCCGGTAACGCCACCTTCACGCAGGCCATTTTCGTAATTCATCAATTGGCGAAGGGCGCGACCGTAAAGGCAAAACACCGGC‑TGGGGGGG‑GAAA‑ACCTTTAACGCCACCggcgtac                                                                                                        >  2:211556/1‑244 (MQ=255)                                                                                                                                       cGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGTCAGTAAAGGCCATACACCGGC‑TGGCGGTG‑AAAA‑TCCTTTAACGCCAGCGG‑TGTCAGGCg                                                                                                   <  1:60542/251‑1 (MQ=255)                                                                                                                                            gTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGC‑TGGCGGGG‑AAAAACCCTTTAACGcacccgggtcagggggacgc                                                                                                >  2:139121/1‑232 (MQ=255)                                                                                                                                               gggATGCGGTTTCAGACGGGTATTACCGATAGCGCCA‑TAGGGCAGATCTCCCACGACCCACACC‑GGGGAGGTGATAAACCGCTCCATCAGATATCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTATCGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCGATTGGCTATGGGCGCGACAGTAAAGGCCATACGCCGGC‑GGGCGGTG‑AAAA‑TCCTTTAACGCCAGCGG‑TGTCAGGCGaccccctg                                                                                           <  2:150555/251‑8 (MQ=255)                                                                                                                                               gggAGGCGGTTTCAGACGGGCATTATAGATAGCACCA‑TTGTGCATACCACCCCCGACCCCCACC‑GGTGAGGTGATAAACCGCTCCATCAGTGAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCTCCTGCACGCAGGCGCTTTTCGTAATTCAGCAATGGGCGATGGGCGCGACAGTAAAGGCCATACACCGGC‑TGGCTGTG‑AAAA‑TCCTTTAACGCCAGCGG‑TGTCAGGCGAACGCCtt                                                                                           <  2:32553/251‑1 (MQ=255)                                                                                                                                                  aTGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGC‑TGGCGGGG‑AAAA‑TCCTTTAACGCCCGCGG‑TGTCAAGCGAACGCCTTGtt                                                                                        >  1:701041/1‑251 (MQ=255)                                                                                                                                                           cAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATAAACCGGC‑TGGCGGTG‑AAAA‑ACCTTTAAAGCCAGCGG‑TGTCAGGctgacgccattgtcaccccgac                                                                               >  1:481031/1‑230 (MQ=255)                                                                                                                                                           cAGACGGGCAGTAACGATAGCGCCA‑TGTTGCATATCACCCGCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGC‑TGGCGGTG‑AAAA‑TCCTTTAACGCCAGCGG‑TGTCAGGCGAACGCCTTGTTCACCCTGCa                                                                               <  2:123459/251‑1 (MQ=255)                                                                                                                                                                      tAACGCTAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGCGATACACCGATCCATCAGATCGCGATCTGGCGG‑GCGC‑ACATCGGACTCGTAGACGGTTACGCCACCTTCAAGCAGGCGCTTTTCGTAATTCCACCATTGGCCAGGGGCGCGAACGGAAAGGCCATACAACGGC‑CGGCGGGGTAAAC‑TCCTTAAACGCCAGGGG‑TG‑GAGGCGAAAGCCttatttcacacgccaaaaagcac                                                                      >  2:668198/1‑228 (MQ=255)                                                                                                                                                                           aTAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATAAACCCGC‑TGGCGGTG‑GAAA‑TCCTTTAAAGCAAAAGGATG‑CAAGCGAACGCCTTGTTAACCCTGCAAAataatgtgagctgg                                                               >  1:551647/1‑240 (MQ=255)                                                                                                                                                                                cgccttggcacatcacctccTCGACCCACCCC‑GGTTAAGTGATAACCCGCTCCATCAGATAGCGTTCTGGCGG‑GCTC‑ACATCGGCCTCTTAGACGGTCACGCCTCCTTCACGCACGCTCATTTCGTAATTCATCAATTGGCTATGGGGGCGACAGTAAAGGCCATCCACCGGC‑TGGCGGTG‑AAAA‑TCCTTTAACGCCAGCGG‑TGGCAGGCGAACTCCTTGTTCACCCTGTAAAATATGCTGTGCGCggtg                                                            <  2:498561/232‑3 (MQ=255)                                                                                                                                                                                         gCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGC‑TGGCGGTG‑AAAA‑TCCTTTAACGCCAGCGG‑TGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATcg                                                  <  1:531674/251‑1 (MQ=255)                                                                                                                                                                                           atatCCCCCGCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCTCCTTCACGCAGGCGCTGTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGC‑TGGCGGTG‑AAAA‑TCCTTTAACGCCAGCGG‑TGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATcggc                                                <  2:575412/251‑1 (MQ=255)                                                                                                                                                                                           atagcACCCCCGTCCCACACC‑GGTGGGGTGATAAACCGTGCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGC‑TGGCGGTG‑AAAA‑TCCTTTAACGCCAGCGG‑TGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCGGAGCGCGGTGAACCTGATcggc                                                <  2:8225/247‑1 (MQ=255)                                                                                                                                                                                                         cctaacC‑GTTGAGGTGATAAACCTCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGC‑TGGCGGTG‑AAAA‑TCCTTTAACGCCAGCGG‑TGTCAGGCTAACGCCGTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGG‑AATAAACggca                                 <  2:577322/247‑4 (MQ=255)                                                                                                                                                                                                                  gtgaggtgaTAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGC‑TGGCGGTG‑AAAA‑TCCTTTAACGCCAGCGG‑TGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGGAAAAAAC‑CCACGGGctct                         >  1:307017/1‑251 (MQ=255)                                                                                                                                                                                                                    gaggtgaTAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGC‑TGGCGGTG‑AAAA‑TCCTTTAACGCCAGCGG‑TGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGG‑AATAAACGCCACGGACTCtt                        <  1:64270/250‑1 (MQ=255)                                                                                                                                                                                                                       gtgaTAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGAAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGC‑TGGCGGTG‑AAAA‑TCCTTTAACCCCAGCGG‑TGTCAGGCGAACGACTTGTTCACCCTGCAAAATATGCTTAGGGCGGGGAACCTGACCGGCGGG‑AATAAAcaacccaggcgtttgtcg                    >  1:403798/1‑234 (MQ=255)                                                                                                                                                                                                                             aaCCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGC‑TGGCGGTG‑AAAA‑TCCTTTAACGCCAGCGG‑TGTCAGGCGAACGCCtt                                                                                           <  2:108316/175‑1 (MQ=255)                                                                                                                                                                                                                             aaCCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGC‑TGGCGGTG‑AAAA‑TCCTTTAACGCCAGCGG‑GGTCAGGCGAACGCCtt                                                                                           >  1:108316/1‑175 (MQ=255)                                                                                                                                                                                                                                     cATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCTCCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGTCATACACCGGC‑TGGCGGTG‑AAAA‑TCCTTTAACGCCAGCGG‑TGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGG‑AATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCt      <  1:80995/251‑1 (MQ=255)                                                                                                                                                                                                                                     cATCAGATAGCGTCCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCGTCACGCAGGCGCTTTTCGTACTTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGC‑TGGCGGTG‑AAAA‑TCCTTTAACGCCAGCGG‑TGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGG‑AATAAACGCCACGGACTCTTGCTGTGCAAGCGTAACCt      <  2:522529/251‑1 (MQ=255)                                                                                                                                                                                                                                      aTCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGC‑TGGCGGTG‑AAAA‑TCCTTTAACGCCAGCGG‑TGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGG‑AATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTg     <  1:154289/251‑1 (MQ=255)                                                                                                                                                                                                                                         aGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGC‑TGGCGGTG‑AAAA‑TCCTTTAACGCCAGCGG‑TGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGG‑AATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTGCaa  <  2:701041/251‑1 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                CAAGCGAGGAGGCGTCGCCATTCTCCGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCTATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGTGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑ACGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGC‑TGGCGGTG‑AAAA‑TCCTTTAACGCCAGCGG‑TGTCAGGCGAAAGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGG‑AATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTGCAA  >  NC_000913/65207‑65687

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