Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	65,316	T→A	I155I (ATA→ATT)	polB ←	DNA polymerase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	65,316	0	T	A	95.2%	189.1 / ‑3.6	64	I155I (ATA→ATT)	polB	DNA polymerase II
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base A (12/48);  total (15/49)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

GGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCAAGCGAGGAGGCGTCGCCATTCTCCGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCTATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGTGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑ACGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGT  >  NC_000913/65070‑65560                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        ggCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAg                                                                                                                                                                                                                                                     <  1:78215/251‑1 (MQ=255)      gtgtttGCAGCAATCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAAcc                                                                                                                                                                                                                                                <  1:52521/248‑1 (MQ=255)      gTTTGTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAAc                                                                                                                                                                                                                                                 <  1:710118/250‑1 (MQ=255)         tttGCATCATTCGCAGATCGAACTGCCACGCGGTCCAACCGATGATCACAGCAGGAGCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCATTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCa                                                                                                                                                                                                                                              <  1:664030/250‑1 (MQ=255)            gcagcaTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACtt                                                                                                                                                                                                                                           <  1:2371/250‑1 (MQ=255)             cagcaTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCCCACCGTGGCGGGGGGGTTCAATATCCACAGAAAACCACTTg                                                                                                                                                                                                                                          >  2:533494/1‑250 (MQ=255)              agcttatgCAGATCTAACTGCCCCACGGTACAACGGTGGACCACATCAAGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGGGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGAGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAg                                                                                                                                                                                                                                        <  2:228785/244‑1 (MQ=255)                caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCCTATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGGTTAAATATCAATACAAACACACTTGTgcg                                                                                                                                                                                                                                      >  1:695179/1‑251 (MQ=255)                   tCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAgcggc                                                                                                                                                                                                                                    <  1:486210/250‑1 (MQ=255)                     gCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGGTTCAAAATCAAAAAAAACCCACTTGAGcggcggg                                                                                                                                                                                                                                 >  1:579188/1‑250 (MQ=255)                               cTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAAACCACTTGAGCGGCGGACGATAGTCgg                                                                                                                                                                                                                       >  1:488484/1‑251 (MQ=255)                                 gccccacGTTCCAACCGAGGAGCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTGCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCCCTTGAGCGGCGGACGATAGTCgg                                                                                                                                                                                                                       <  2:137699/246‑1 (MQ=255)                                    ccacGGTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTGCCAGGAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCGCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGc                                                                                                                                                                                                                  <  2:309764/251‑1 (MQ=255)                                           ccAACAGATGATCACATCAGGATCGTAGTTGGCCAACCGGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGTCGACGTATTCCAGATCAAACGCGAGCGCGGAGGCGTCGCAATTCTCTGGCCCAAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGAGAGTCGGGATGCGGTTTCa                                                                                                                                                                                                           <  2:113001/251‑1 (MQ=255)                                                     atCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCAtt                                                                                                                                                                                                 <  1:263346/251‑1 (MQ=255)                                                     atCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGGTCAAAATCGAGCGCGGAGGCGGCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCAtt                                                                                                                                                                                                 <  2:710125/251‑1 (MQ=255)                                                         caTCAGGATCTGAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACg                                                                                                                                                                                             <  2:545036/251‑1 (MQ=255)                                                            cAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATa                                                                                                                                                                                          <  1:648453/251‑1 (MQ=255)                                                            cAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATa                                                                                                                                                                                          <  1:167830/251‑1 (MQ=255)                                                               gATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAgcg                                                                                                                                                                                       <  1:47021/251‑3 (MQ=255)                                                                 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‑GCGcac                                                                                                       <  1:316419/246‑1 (MQ=255)                                                                                                                                            aGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTTAATATCAAAAAAAACCCACTTGAGCGGCGGACGGGAGTCGGGGAGCGGGTTCTGAAGGGCGTTAACGATAACGCCC‑TTTTGGATATCCCCCCCGCCCCCCACC‑CGGGGGGGGAAAAAACGCTCCCTCAGATAACGTTCTGGCGG‑agcgcc                                                                                                       >  2:675583/1‑246 (MQ=255)                                                                                                                                              ggcggatgtgTCGCCATTCTCTGGCCCAAGCAGATAAACGATGCTTGGCCCGCAGCCGACCAGGCCGAGGAAGTGCATCTCACCGGGGCGGGTGGTTTCAATATCAATAGAAACCCGCGTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGGGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGCACAt                                                                                                     <  2:429507/242‑1 (MQ=255)                                                                                                                                               gcggaggCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGCACATc                                                                                                    <  1:99061/249‑1 (MQ=255)                                                                                                                                                     gCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGCACATCGGCCTc                                                                                              <  1:483219/251‑1 (MQ=255)                                                                                                                                                       gTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCCATATCTATTGAAACCCACTTGAGCGGCGGACGCTAGTCGGGATGCGGTTTCAGTCCGGCCTTCACGATAGCGCCC‑TTGTGCGTAACACCCTCTACCCACCCC‑GGCGGGGTGATAAACCGCTCCCTCAGATAGCGTTCTCGGGG‑GCGCCCATCGGCCTCgg                                                                                            >  1:9722/1‑250 (MQ=255)                                                                                                                                                        gcgtcattcgcgggcccaatcATATTAACGAGGCGCTGCCAGCGGACTACCAGGCCGATGCTGTACAGCTCACCGTGGCGGGTGGTGTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGTGATGCGGTTTCAGACGGGCATGAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGCGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGCACATCGGCCTCGt                                                                                            <  2:634338/230‑1 (MQ=255)                                                                                                                                                         cGCCCTTCTCTGGCCCCAGCATATAAACGTAGCGATGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATACCACCCGCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGCACATCGGCCTCGt                                                                                            <  2:255463/249‑1 (MQ=255)                                                                                                                                                            cATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGCACATCGGCCTCGTAGACg                                                                                       <  1:640415/251‑1 (MQ=255)                                                                                                                                                              ttCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCATTCCAGGCCGCTGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAAAAAACCACTTGAGCGGCGGACGATCATCGGGAGGTGGGTTCAGGCGGGCATTAACGATAGCGCCA‑ATGTGCATATTACCCTCGACCCACCCC‑GGTGAGGTGGTAAACCACTCCATACGATAGCGTTGTGGCggtggcacacagagcacgttgaatgg                                                                                      >  1:593445/1‑227 (MQ=255)                                                                                                                                                                  ctggCCCCAGCATATAAACTAGGCGTTTCCCGCATGCTTCCAGGCCCATGCAGTACAGCTCACCGTGGCGGGTGGTTTTAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCAGTTAT‑CATATCACCCTCGATCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGCACATCGGCCTCGTAGACGGTa                                                                                    <  2:702365/246‑1 (MQ=255)                                                                                                                                                                       cccAGCATATAAATGTTGCGCGGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCCGCGGGGCGGGTGGTTTCAATATCAAGAGAATCCCACTTGAGCGGCGGACGATAGTCTGGATGCGGTTTCAGACGGGCATTGACGATAGCGCCA‑TTGTGCATAGCACCCTCGACCCACACC‑GGTGGGGTGATAGACCGCTGCATCAGATAGCGCTCTGGCGG‑TGGCACATCGGCCTCGTAGACTGTAACGCCAcc                                                                            <  2:637742/251‑1 (MQ=255)                                                                                                                                                                           gCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCa                                                                        <  1:278024/251‑1 (MQ=255)                                                                                                                                                                            cATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCAc                                                                       <  1:509259/251‑1 (MQ=255)                                                                                                                                                                             atatAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAAACCACTTGAGCGGCGGACGATAGTCGGGATGCGGGTTCAGAAGGGCCTTAACGATAACGCCA‑TTGTGCATATCACCCTCGACCC‑CACCCGGGGGGGGGATAAACCGCTCCCTCAGGAAGCGTTCTGGCGG‑GGGCACATAGGCCCCCTTAACGGTAACGCCACCtttccc                                                                      >  1:701050/1‑247 (MQ=255)                                                                                                                                                                               ataAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCa                                                                    <  1:458267/251‑1 (MQ=255)                                                                                                                                                                                        gcgcTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCtttt                                                           <  1:485439/251‑1 (MQ=255)                                                                                                                                                                                             gCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTGACCGTGGCGGGTGGTTTCAATATCAAGAGAAACCCACTTGAGCGGCGGACGATAGGCGGGATGCGGTTTCAGACGTGCATTAACGATAGCGCCA‑TGGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCGCCATCAGATAGCGTTCTGGCGG‑GCGCACATCGGCCTCGTAGACGGTAACGCCACCGTCACGCAGGCGCTTTTCGt                                                        <  2:351166/249‑1 (MQ=255)                                                                                                                                                                                                 gcggcattCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGGGTCAATATCTATAGAAACCCACTTGCGCGGCGGGCGATAGTCGGGAGGCGGGTTAAGACGGGCCATAACGAGAGTGCCC‑TTGTGAATGTCACCCTCGGACCACAAC‑GGGGGGGTGGTAAACCCCCACATCCGGAAACGTTCTTGCGGGAC‑CACACGGGCCCCCTCGACGGTAACGCCAACTTTACGCACGCGCttttgagtaattt                                                  >  2:206250/7‑242 (MQ=255)                                                                                                                                                                                                         ccAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTg                                           <  1:167115/250‑1 (MQ=255)                                                                                                                                                                                                         ccAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTg                                           <  1:405540/250‑1 (MQ=255)                                                                                                                                                                                                           gggCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCg                                        <  2:705545/250‑1 (MQ=255)                                                                                                                                                                                                              ccGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATg                                     <  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‑GCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTa                         <  1:258972/251‑1 (MQ=255)                                                                                                                                                                                                                              cACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTaaa                       <  1:263617/249‑1 (MQ=255)                                                                                                                                                                                                                               cccGTGGCGGGTGGTTTCAATAGCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGc                    <  1:212348/250‑1 (MQ=255)                                                                                                                                                                                                                               cccGTGGCGGGTGGATTCAATATCAATCGAAACCCACTTGAGCGGCGGACGATAGTCGTGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCGCGACCCACACC‑GGTGAGGGGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGCACATCGGCCTCGGAGACGGTAACGCCACCTTCACGCAGGCGCTGTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGc                    <  2:479084/250‑1 (MQ=255)                                                                                                                                                                                                                                 cGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCa                  <  1:392874/251‑1 (MQ=255)                                                                                                                                                                                                                                            gtTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCTGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCTCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTg       <  2:36639/251‑1 (MQ=255)                                                                                                                                                                                                                                              ttCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCCGCTGGc     >  1:642351/1‑251 (MQ=255)                                                                                                                                                                                                                                                cAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGCACCTCGGCCTCGTAGACGGTAACGCACCCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGCTGGGCGCGCCAGTAACAGCCACACACCGGcctgcgg   >  2:154289/1‑245 (MQ=255)                                                                                                                                                                                                                                                 aaTATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑GCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGt  <  1:465335/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        GGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCAAGCGAGGAGGCGTCGCCATTCTCCGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCTATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGTGCCA‑TTGTGCATATCACCCTCGACCCACACC‑GGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGG‑ACGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGT  >  NC_000913/65070‑65560

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