Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	65,381	T→C	T134A (ACT→GCT)	polB ←	DNA polymerase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	65,381	0	T	C	98.5%	203.0 / ‑5.0	66	T134A (ACT→GCT)	polB	DNA polymerase II
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (25/40);  total (26/40)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

TCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCAAGCGAGGAGGCGTCGCCATTCTCCGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCTATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGTGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑ACGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGC‑TGGCGGTG‑AAAA‑TCCTTTAACGCCAGCGG‑TGTCAGGCGAAAGCCTTGTTCACCCTGCAAAATATGCTGAGCGCG  >  NC_000913/65134‑65627                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                               cacccgacgctcctccgatctttGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAAGCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGcc                                                                                                                                                                                                                                                             <  2:225542/230‑1 (MQ=255)    tAGTTGGCAAACCAGGCGTGTATTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCa                                                                                                                                                                                                                                                            <  2:96670/249‑1 (MQ=255)           cacaCCAGCCGTTTAGTTTTTCAAGTAGCTGCGGGCGGCTGGACCCGTAGTCCAGTATAAAAGCGAGCGCGTAGCCGTCGCCATTCGCTGGCCCCAGCATATAAACGATGCGCTGCTCGCAGCCTTCCAGGCCGATGCAGTACTGCTCACCGTGGCGGGTTGTTTCAAGATCAATAGAAACCGGCTTTAGTGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCata                                                                                                                                                                                                                                                  <  2:269291/248‑1 (MQ=255)                      ttGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGGTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCCTTAAAGATAGCGcccgatcgggagagcacacgtcc                                                                                                                                                                                                                                        >  1:225542/1‑230 (MQ=255)                                 cagtagcttcgggcggCTGGCGACGTAGACCAGTTCAAAAGCGAGCGCGGAGGCGTCTCCATTCTCTGGCCCCAGCATATAAACGATGCTCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTTGACCCACACC‑ggtg                                                                                                                                                                                                                           <  2:470575/242‑1 (MQ=255)                                          cgggcggCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGGTTCAATATCAATAGAAAACCACTTGAGCGGCGGACGGTAGTCGGGATGCGGTTTCAGACGGGCCTTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCCCCAC‑GGTGAGGGGATa                                                                                                                                                                                                                   >  1:575412/1‑250 (MQ=255)                                                gcTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCt                                                                                                                                                                                                            <  1:368271/251‑1 (MQ=255)                                                   ggCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTcc                                                                                                                                                                                                          <  1:431363/250‑1 (MQ=255)                                                           attccatgtcaaaaTCCAGCGCGGAGGCGCCGCCATTCTCTGGCCCCCGCATATAAACGAGGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGGGGCGGGTGGTTTCAATATCTATAGAAACCCACTGGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACTGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATa                                                                                                                                                                                                  <  2:376303/239‑1 (MQ=255)                                                                          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)                                                                                              ttCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCATTCCAGGCCGCTGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAAAAAACCACTTGAGCGGCGGACGATCATCGGGAGGTGGGTTCAGGCGGGCATTAACGATAGCGCCA‑ATGTGCATATTACCCTCGACCCACCCC‑GGTGAGGTGGTAAACCACTCCATACGATAGCGTTGTGGCggtggcacacagagcacgttgaatgg                                                                                                                                                              >  1:593445/1‑227 (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:167115/250‑1 (MQ=255)                                                                                                                                         ccAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTg                                                                                                                  <  1:405540/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)                                                                                                                                                gATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑ggtgag                                                                                                                                                                                                                         <  2:147216/142‑1 (MQ=255)                                                                                                                                                gATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑ggtgag                                                                                                                                                                                                                         >  1:147216/1‑142 (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)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                               TCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCAAGCGAGGAGGCGTCGCCATTCTCCGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCTATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGTGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑ACGC‑ACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGC‑TGGCGGTG‑AAAA‑TCCTTTAACGCCAGCGG‑TGTCAGGCGAAAGCCTTGTTCACCCTGCAAAATATGCTGAGCGCG  >  NC_000913/65134‑65627

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