Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	566,651	0	T	C	54.7%	19.1 / 86.5	65	intergenic (‑673/‑33)	intD/renD	putative integrase/protein RenD
Reads supporting (aligned to +/- strand):  ref base T (16/13);  new base C (17/18);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.57e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAACT‑GACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTG  >  NC_000913/566412‑566893                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   gcggagTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCC‑ATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGt                                                                                                                                                                                                                                          <  1:355193/247‑1 (MQ=255)  cggagTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCC‑ATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTc                                                                                                                                                                                                                                        <  1:220646/249‑1 (MQ=255)                 ggAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAGTGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCC‑ATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGc                                                                                                                                                                                                                          >  2:173790/1‑250 (MQ=255)                 ggAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGAGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCC‑ATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGc                                                                                                                                                                                                                          >  1:334060/1‑250 (MQ=255)                        gAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAACT‑GACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTAt                                                                                                                                                                                                                  >  1:174494/1‑251 (MQ=255)                            gcgcACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCC‑ATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACg                                                                                                                                                                                                              >  1:12193/1‑251 (MQ=255)                             cgcACTACTGGCTCGTTACGAACCTGTATCAGAACATGCGGGCCAATTCACTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCC‑ATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGTCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGCCTGAAAGGAGCAAGTGTATGACgg                                                                                                                                                                                                             <  1:173427/251‑1 (MQ=255)                                actACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCC‑ATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCa                                                                                                                                                                                                          <  2:170102/249‑1 (MQ=255)                                     tggctggttACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCC‑ATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAag                                                                                                                                                                                                       >  1:370670/5‑249 (MQ=255)                                      ggctggttACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCC‑ATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTGCGGGCTGAAAGGAGCAAGTGTATGACGGGCAAAGAgg                                                                                                                                                                                                    <  1:188769/248‑1 (MQ=255)                                         tggttACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCC‑ATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCa                                                                                                                                                                                                  >  1:401620/1‑250 (MQ=255)                                                 aaCTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAACT‑GACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCAt                                                                                                                                                                                         >  2:332242/1‑251 (MQ=255)                                                  acctGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCC‑ATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCAtt                                                                                                                                                                                        <  2:173078/248‑1 (MQ=255)                                                  acctGTAGCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCC‑ATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCGCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGGATGACGTGCAAAGAGGCAATTATTCAtt                                                                                                                                                                                        <  2:22241/248‑1 (MQ=255)                                                    ctGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCC‑ATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTAc                                                                                                                                                                                      >  1:150847/2‑251 (MQ=255)                                                     tGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAACT‑GACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTAcc                                                                                                                                                                                     >  2:264371/1‑251 (MQ=255)                                                         ccAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAACT‑GACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTggg                                                                                                                                                                                 >  2:171725/1‑251 (MQ=255)                                                           aGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCC‑ATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTgggg                                                                                                                                                                                >  2:314249/1‑250 (MQ=255)                                                           aGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCC‑ATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGa                                                                                                                                                                               <  1:242741/251‑1 (MQ=255)                                                            gAACATGCGGGCAAATGCGCTGGCTGACGCGGAATTACGGCTCAAGGCTGCCGATGAACT‑GACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATCCATTACCTGGGGAc                                                                                                                                                                              <  2:35379/251‑1 (MQ=255)                                                              aCATGCTGGCCACTGCGCTTACTGACGCGGAGTTACGTCGTAAGGCCGCAGATGAGCT‑TGTCAATATGACTGCGAGAATTAACCGTGGTGGGGGGATGCCTGAACCAGAAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCAAAGGATCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAACGCGGCAATTATTCGTTACCTGGGGACGc                                                                                                                                                                            >  2:183264/1‑251 (MQ=255)                                                               cATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAACT‑GACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGGTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAATATCGCAGAAATTAAAGTTAAGTTCGGACTGAAAGGATCAAGTGTATGACGTGCAAAGAGGCAATTATTCATTACCTGGGGACGCa                                                                                                                                                                           <  1:185484/251‑1 (MQ=255)                                                                    gggCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCC‑ATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATaat                                                                                                                                                                       >  2:346247/1‑249 (MQ=255)                                                                       ccAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAACT‑GACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCt                                                                                                                                                                   >  2:352500/1‑251 (MQ=255)                                                                         aaTGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCC‑ATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGCCGCATAATAGCTTc                                                                                                                                                                 >  2:215132/1‑251 (MQ=255)                                                                             cgcTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAACT‑GACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCtgtg                                                                                                                                                             >  1:131589/1‑251 (MQ=255)                                                                              gctggctgACGCGGAATTACGGCGCAAGGCTGCCGATGAACT‑GACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGc                                                                                                                                                            <  1:299413/251‑1 (MQ=255)                                                                                               tACGGCGCAAGGCTGCCGATGAACT‑GACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCcgcg                                                                                                                                            <  1:171725/250‑1 (MQ=255)                                                                                                        aGGCTGCCGATGAACT‑GACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGgcg                                                                                                                                  <  1:90686/251‑1 (MQ=255)                                                                                                         ggCTGCCGATGAACT‑GACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGgcgc                                                                                                                                 <  2:260555/251‑1 (MQ=255)                                                                                                                   gAACT‑GACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATc                                                                                                                       >  2:31445/1‑251 (MQ=255)                                                                                                                   gAACT‑GACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATAAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAAt                                                                                                                        <  2:286550/250‑1 (MQ=255)                                                                                                                               tATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTAccc                                                                                                            >  1:78382/1‑251 (MQ=255)                                                                                                                                aTGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatc                                                                                                             >  1:39409‑M1/1‑237 (MQ=255)                                                                                                                                aTGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTTCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCa                                                                                                           <  2:160813/251‑1 (MQ=255)                                                                                                                                aTGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCa                                                                                                           <  2:160797/251‑1 (MQ=255)                                                                                                                                  gACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGc                                                                                                         >  2:98448/1‑251 (MQ=255)                                                                                                                                                aCCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatgg                                                                                           <  2:384937‑M1/251‑31 (MQ=255)                                                                                                                                                   gtggtgAGACGATACCTGAACCAGGAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGTCAATTATTCATTACCTTGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCACTAGTGGACAAGCGGc                                                                                        >  2:110198/1‑251 (MQ=255)                                                                                                                                                      gtgAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacggg                                                                                     <  1:29547‑M1/251‑37 (MQ=255)                                                                                                                                                      gtgAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGAGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAcaagcataaatcaggccgcggctaaaatggcacggg                                                                                     >  1:67678‑M1/1‑215 (MQ=255)                                                                                                                                                         aggcgatCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGTCGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggc                                                                                    >  2:290207‑M1/4‑212 (MQ=255)                                                                                                                                                               tACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTaa                                                                             >  1:307244/1‑250 (MQ=255)                                                                                                                                                                        cAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCtcagtc                                                                   >  1:168069/1‑251 (MQ=255)                                                                                                                                                                          gTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCtcagtcag                                                                 <  1:400559/251‑1 (MQ=255)                                                                                                                                                                            aaaaCAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCaa                                                                                                                               <  1:42114/187‑1 (MQ=255)                                                                                                                                                                            aaaaCAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCaa                                                                                                                               >  2:42114/1‑187 (MQ=255)                                                                                                                                                                                  acTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACTCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCTACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggt                                                         <  1:410441‑M1/251‑65 (MQ=255)                                                                                                                                                                                    ttCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTc                                                        >  2:246678/1‑250 (MQ=255)                                                                                                                                                                                        tgtcaTGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTg                                                                                                                                                 <  1:69671/153‑1 (MQ=255)                                                                                                                                                                                        tgtcaTGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTg                                                                                                                                                 >  2:69671/5‑157 (MQ=255)                                                                                                                                                                                           catGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCGTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggcgaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctggcgaac                                                <  2:150847‑M1/250‑74 (MQ=255)                                                                                                                                                                                              gggCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaa                                              >  2:220757/1‑250 (MQ=255)                                                                                                                                                                                                gCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctggcgaacggtgt                                           <  2:401620‑M1/251‑79 (MQ=255)                                                                                                                                                                                                             aaaTCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctggcgaacggtgtattaccggtttgc                              >  1:100622‑M1/1‑160 (MQ=255)                                                                                                                                                                                                              aaTCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGt                              >  1:156522/1‑250 (MQ=255)                                                                                                                                                                                                              aaTCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTa                             <  1:264371/251‑1 (MQ=255)                                                                                                                                                                                                                     cacaGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctggcgaacggtgtattacctgtttgctaccaggg                      <  1:118675‑M1/251‑100 (MQ=255)                                                                                                                                                                                                                         ggCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGt                  <  1:98448/251‑1 (MQ=255)                                                                                                                                                                                                                               ggCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctggcgaacggtgtattaccggtttgctaccagggaagaacgg              <  1:290207‑M1/249‑108 (MQ=255)                                                                                                                                                                                                                                    aGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctggcgaacggtgtattaccggtttgctaccagggaagaacgggaagga        >  1:410713‑M1/1‑137 (MQ=255)                                                                                                                                                                                                                                     gATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTc       >  2:301082/1‑250 (MQ=255)                                                                                                                                                                                                                                     gATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctggcgaacggtgtattaccggtttgctaccagggaagaacgggaaggaaa      <  2:39409‑M1/251‑116 (MQ=255)                                                                                                                                                                                                                                      aTCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGc     <  1:31445/251‑1 (MQ=255)                                                                                                                                                                                                                                         gCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATGCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCGTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAATCTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTg  <  1:46156/251‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   GCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAACT‑GACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTG  >  NC_000913/566412‑566893

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