Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,818,362	G→A	intergenic (‑88/+111)	argY ← / ← argV	tRNA‑Arg/tRNA‑Arg

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,818,362	0	G	A	94.5%	186.1 / ‑5.3	78	intergenic (‑88/+111)	argY/argV	tRNA‑Arg/tRNA‑Arg
Reads supporting (aligned to +/- strand):  ref base G (0/4);  new base A (36/33);  total (36/42)
Fisher's exact test for biased strand distribution p-value = 1.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCCACCTCTTCCATACCTTCATCGCTGAAAGTACGTAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCAAATGGCGGTGAGGCGGGGATTCGAACCCCGGATGCAGCTTTTGACCGCATAC  >  NC_000913/2818115‑2818602                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  aTCCAGCTGAGCTGCGTATGCATCGGGAAAGTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCt                                                                                                                                                                                                                                                <  1:22333/251‑1 (MQ=255)   ccagatgatgttctgatccatggggAAATTTTTTTTACTGCTGATATTGATGACGCTAATAACGCCAAATCAAGTAAGATATGGTGCACCCGGGAGGATGCGAACCTCCGACCGCGCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTTCTTTTCAGCAGAGTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTg                                                                                                                                                                                                                                              <  2:43890/229‑1 (MQ=18)     aGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCg                                                                                                                                                                                                                                            >  1:178504/1‑251 (MQ=255)         ttgttACGGATGCCTCGGGAAATTTCTCTATCTGCCGATACTGAGACCACTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGAGGCACCGGGAAACTTTCTTTACTGCAGATATTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCgg                                                                                                                                                                                                                                           <  2:557979/244‑1 (MQ=18)         gAGCTACGGATGCATCGGGAAATTTATTTAACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTACTTGACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGc                                                                                                                                                                                                                                          <  1:407047/249‑1 (MQ=255)           tctgcGGATGCATCGGGAAATTTAGTTGACTGATGAGATTGATACCGCGACTAAAGCCATATCAAGTAAGAGTTGGGGCAGCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCTTCGCCGAGTACACTATCCAGCGGAG‑CTACGGAGGCATCGGGAAACGTACTTTACAGCAGATTTTTGGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGAGTACTCACCGTGCGGc                                                                                                                                                                                                                                          <  2:452827/243‑1 (MQ=18)           gCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGCTCGTAGCCGAGTACTCTATCCAGCTGAC‑CTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGc                                                                                                                                                                                                                                          >  1:452827/1‑247 (MQ=255)            cTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTCCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCCGCTCGcc                                                                                                                                                                                                                                     >  1:372487/1‑251 (MQ=255)             tACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGACTCGccc                                                                                                                                                                                                                                    >  1:557978/1‑251 (MQ=255)                ggacccatcgtgaaaggtaatttACTGGTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCAAGTACTCTACCCAGCTGAG‑CTACGGATGCATCGGGAAACATACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTc                                                                                                                                                                                                                                 <  2:247888/231‑1 (MQ=21)                ggATGCATCGGGAAATTTTTTTTACGGCTGATCTGGATTCCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTc                                                                                                                                                                                                                                 <  1:507945/251‑1 (MQ=255)                ggATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTc                                                                                                                                                                                                                                 >  1:181438/1‑251 (MQ=255)                 gATGCATCGGGAAATTTCTTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCCCCTAAAGCCGTATCAAGTAAGAGATGATGCATCCGGGAGGTTTACTAACCTTGCGGCTCCCCCTTcc                                                                                                                                                                                                                                >  2:300227/1‑250 (MQ=25)                  aTGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGCTTACTCACCTTGCGGCTCGCCCTTCgg                                                                                                                                                                                                                               >  1:70765/1‑251 (MQ=255)                    gCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCACGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGc                                                                                                                                                                                                                             <  1:296956/251‑1 (MQ=255)                    ccattgggAAATGTATTTTACTGCTGATATTTAAACCGCTACTAAAGCCATATCAAGTAAGATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCggg                                                                                                                                                                                                                              <  1:497343/245‑1 (MQ=21)                        cGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGtt                                                                                                                                                                                                                         >  1:515499/1‑251 (MQ=255)                          ggtaaTTTATTTTGCGTCTGTTATTGTTTCCGCGACTAAAGCAATATTAATTATGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTACTGTACTGCAGATTTTTTGATACCGCTACTATAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCTGGCCGTTG‑c                                                                                                                                                                                                                      <  2:665536/248‑1 (MQ=17)                                    tttACTGCTTATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAAc                                                                                                                                                                                                             <  1:18289/250‑1 (MQ=255)                                    tttACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTACTTTACTACAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGGGCATCCGGGAGGATTACTCACCTTGCGGCACGCCCTTcatgcggtggctaaagaaa                                                                                                                                                                                                               >  1:167622/1‑231 (MQ=21)                                       actgtgtatatTGAGACCTCTACTAAATCCAGATCAAGGAAGAGATGGTGCAGACGGGGGGAATCGAACCTCCGACCGGTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑GTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGATATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGTCGTTG‑CTAAAGTAACGTTa                                                                                                                                                                                                         <  1:604545/244‑1 (MQ=17)                                             gATATTGATGCCGCTACATAAGCCATAGTAAATAAAAGATGGGGCATCCGGGAGTATTCGAACCTCCGACCGCTCCGTTCGTAGCCGAGTACTCTATCCAGCTTAG‑CTACGGATGCCTCGGGGAAATTATTTTACATCAGATATTTTGATTCCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCGTCGGGCCGTTG‑CTAAAGCAACGTTAtcctc                                                                                                                                                                                                    <  1:289568/250‑1 (MQ=14)                                                         gCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGCAACTTACTGTACGGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATGCGGGAGGATTACACACCTTGCGGCTCGCCCCTCGGGCCTTTGGCTAA‑GGAACGTTTACCTCCCTgggcccctcg                                                                                                                                                                                       >  2:296013/1‑243 (MQ=18)                                                         gCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATAGGGAAACTTACTTTAATGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAACTAAGAGTTGGTGCATCCGGGAGGCTTACTCACCTTGCGGCTCGCCCTTATGACCGTTT‑CTAAAACAACATTATCCTCCCTGGTCCCTGct                                                                                                                                                                                       >  1:123463/1‑250 (MQ=18)                                                           tactaAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCggg                                                                                                                                                                                     >  1:447712/1‑250 (MQ=255)                                                               aaaGCCATATCAAGTAAGATATGGTGCATCCGGGAGGATTCGAACCCCCGCCCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCCCCCTGGTGCCTGCGGtta                                                                                                                                                                                   <  1:486033/249‑1 (MQ=255)                                                               aaaGCCAGATCAAGTCAGAGATGGGGCATCCGGGGGGATGCGAACCGCCGACCGCTCGGTTCGGAGCCGAGTCCGCTACCCAGCTGAG‑CTACGGATGCAGCGGGAAACTTTCTTTACTGCAGATTTTTTGATACCGCTCCTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGttat                                                                                                                                                                                  <  1:507683/250‑1 (MQ=18)                                                                 agctatatcttgtaagcgattgggCATCCGGTAGGTTACTAACCTCCGACCGATCGGTTCGTTGCCGAGTACTCTATCCAGCTGAT‑CTACGGATGCATCTGGAAACTTACTTTTCTGCAGATTTTTTGAGACCGCTTCTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTAc                                                                                                                                                                               <  2:410779/228‑1 (MQ=17)                                                                    aaggttgcggaggataatggggaTCCGGGAGGGTTCGAACCTCCGACCAGTCGGTGCGTAGCCGAGTACACTATCCAGCTTAG‑CTACGGATTCATCGGGAAACGTACTTGACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTGG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGc                                                                                                                                                                            <  2:147411/229‑1 (MQ=17)                                                                        gcAAGTAAGAGATGTGGCATCCGGGAGAAGTCGAACCTCCGCCCGCTCGGTTCGTAGCCGAGGTCTCTACCCAGCTGAG‑CTACGGATGCATCGGGAAACTTTCTTTTCTGCAGTTTTGTTGAAACCCCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCCCCTTGCGGCCCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCCAcc                                                                                                                                                                        <  2:705315/250‑1 (MQ=18)                                                                          aaGTAAGAGATGGTTCATTCGGGAGTATTCGAACCTCCGACCCCGCGGTTCGTAGCCGAGTACTCTATCCATCTGTG‑CTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTGGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTGGCGGCTTGCCCTTCGGGCCGGTG‑CTAAAGCAACGTTATCCTCTCTGGTTCCTGCGGTTATTACCTCCACctc                                                                                                                                                                      <  2:372487/251‑1 (MQ=21)                                                                          aaGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCCACctc                                                                                                                                                                      <  1:78740/251‑1 (MQ=255)                                                                              aaGAGATGGTGCATACGGGAGGATTCGAAACTCTGACCGCTCGGTTCGTAGCCGAGTACTCGATCCAGCTGAG‑CTACGGATGCAACGGGAAACTGACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCCACctctgc                                                                                                                                                                   <  1:654807/250‑3 (MQ=18)                                                                                gagaTGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTCTCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTACAGCAACGTTATCCTCCCTTGTGCCTGCGGTTATTACCGCCCCCTCTTCCAt                                                                                                                                                                >  1:266360/1‑251 (MQ=255)                                                                                 agaTGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGACCTTCGGGCCGTTG‑CTAAAGCAACGTTAACCTCCCTGGTGCCTGCGGTTTTTACCGCCCCCTCTTCCAt                                                                                                                                                                >  2:416232/1‑250 (MQ=255)                                                                                   aTGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTACTTTACTGCCGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCCTCGGGCCGGTG‑GTAAAGCAACGATATCCTACCTGGTGGCTGCGGGTATTACCACCACCTCTTCCATAcc                                                                                                                                                             >  1:34150/1‑251 (MQ=18)                                                                                          aTCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCCACCTCTTCCATACCTTCATCg                                                                                                                                                      <  1:416232/251‑1 (MQ=255)                                                                                                gaggaTTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTCAAACAACGTTATTCTCCCTGGTGCCTTCGGGTATTACCGCCACCTCTTCCATACCTTCACCGCttaa                                                                                                                                                 >  1:231445/1‑247 (MQ=18)                                                                                                 aggaTTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATTCTCCCTGGTGCCTGCGGTTATTACCGCCACCTCTTCCCTACCTTCCTCGCTGaa                                                                                                                                                 >  1:440564/1‑249 (MQ=255)                                                                                                    aTTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAGACTTACTTTACAGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGGCGATT‑ATAAAGCAACGTTTTCCTCCCTGGGGCCTTCGGTTATTACCCCCACCTCTTACATCCCTTCCTCGATGAAAGTAc                                                                                                                                            >  2:158069/1‑251 (MQ=18)                                                                                                         aaCCTCCGACCGCCCGGTTCGTAGCCGAGTACTCTATCCGGCTGAG‑CTCCGGATGCATCGGGAAACTTACTTTACTGCAGATTGTTTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCCACCTCTTCCATACCGTCATCGCTGAAAGTACGTAGa                                                                                                                                       <  1:240880/251‑1 (MQ=18)                                                                                                                    gCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTACTTTATTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCCACCTCTTCCATACCTTCATCGCTGAAAGTACGTAGAATATGGTGCAt                                                                                                                            <  1:696896/251‑1 (MQ=255)                                                                                                                     cTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGCGCCTCCGGTTATTACCGCCACCTCTTCCATACCTTCATCGCTGTAAGTACGGAGAATATGGTGCATc                                                                                                                           >  2:202307/1‑251 (MQ=255)                                                                                                                       cggtgcgttgtcGAGTACTCCTTCCATCTGAGTCT‑CGGACGCATCTGGATACTTACTTTATTGCATATTTTTTGATCCCGCTTCTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTGGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCTACCTCTTCCATACCTGCAGCGCTGAAAGGACGTAGAATATGGTGCAt                                                                                                                            <  2:624903/237‑1 (MQ=18)                                                                                                                       cGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCCACCTCTTCCATACCTTCATCGCTGAAAGTACGTAGAATATGGTGCATc                                                                                                                           <  1:300227/249‑1 (MQ=255)                                                                                                                         gTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGTGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCa                                                                                                                                                                                                               <  1:201081/163‑1 (MQ=255)                                                                                                                                       aCTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCGCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCCACCTCTTCCATACCTTCATCGCTGAAAGTACGTAGAATATGGTGCATCCGGGAGGATTCGAACCTc                                                                                                         <  2:447712/251‑1 (MQ=255)                                                                                                                                                  cTGAG‑CTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCCGGTGCCTGCGGTTATTACCGCCACCTCTTCCATACCTTCATCGCTGAAAGGACGTAGAATATGGTGCATCCGGGAGGATTCGAACCTCCCACCGCTCgg                                                                                              >  2:228946/1‑251 (MQ=255)                                                                                                                                                  cTGAG‑CTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGATCGCCCTTCGGGCCGTAG‑CTAAAGCAAAGTTATCCTCCCTGGTGCCTTCGGTTATTACCGCCAACTCTTCCATACCTTCATTGCTGAACGGACGTAGAATATGGTACAACCGGGAGGATTCGAACCTCCgccccctagg                                                                                              >  1:708438/1‑242 (MQ=18)                                                                                                                                                        cTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCCACCTCTTCCATACCTTCATCGCTGAAAGTACGTCGAATATGGTGCATCCGGGAGGATTCGAACCTCCGCCCGCTCGGTTCGt                                                                                         >  2:77885/1‑251 (MQ=255)                                                                                                                                                          aCGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCCACCTCTTCCATACCTTCATCGCTGAAAGCACGTAGAATCTGGTGCATCCGGGAGGATTCGCACCTCCGACCGCTCGGTTCGt                                                                                         >  1:370930/1‑249 (MQ=255)                                                                                                                                                              aTGCATCGGGACACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAATCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCCACCTCTTCCCTACCTTCATCGCTGAAAGTACGTAGAATATGGTGCATCCGGGAGGATTCGAACCTCCCACCGCTCGGTTCGTAGCCgc                                                                                   >  2:110581/1‑250 (MQ=255)                                                                                                                                                               tGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCCACCTCTTCCATACCTTCATCGCTGAAAGTACGTAGAATATGGTGCATCCGggag                                                                                                                     <  1:425781/216‑1 (MQ=255)                                                                                                                                                               tGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCACTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCCACCTCTTCCCTACCTTCATCGCTGAAAGGACGTAGCATATGGTGCATCCGggag                                                                                                                     >  2:425781/1‑216 (MQ=255)                                                                                                                                                                    cGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTTTCCTCCCCGGTGCCTGCGGTTATTAACGCCACCACTTCCATACCTTCATCGCTGAAAGTACGTAGAATTTGGTGCCTCCGGGATGATTCGAACCTCCGACCGCTCCGTTCTTAGCCCCGTActc                                                                             >  1:440414/1‑251 (MQ=21)                                                                                                                                                                         aaCTTACTTTACTGCAGATTTTTTGATACCTCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCTCCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCCACCTCTTCCATACCTTCATCGCTGAAAGTACGTAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATcc                                                                        <  1:553352/251‑1 (MQ=255)                                                                                                                                                                           cttacttTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCCACCTCTTCCCTACCTTCATCGCTGAAAGTACGTAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGAACGCTCGGTTCGCAGCCGAGTACTCTATcc                                                                        >  2:374082/1‑249 (MQ=255)                                                                                                                                                                                          aTTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACATTGCGGCTCGCACTTCGGGCCGTTG‑CTAAAACAAAGTTATCCTCCCTGATGCCTGAGGTTATTACCGCCAACTCTTCCATACCTTCATCGCTGAAAGTACGTAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATACAGCTGAGCTACGGATGc                                                       >  1:664661/1‑251 (MQ=17)                                                                                                                                                                                             ttttGATACCGCTACTAAAGCCTTATCAAGTAAGAGATGGTGCACCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCTTTT‑CTAAAGCAACGTTATCCCCCCTGGTGCCTGCGGTTATTACCGCCACCTCTTCCATACCTTCATCGCTGAAAGTACGTAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCaa                                                     <  1:546516/250‑1 (MQ=255)                                                                                                                                                                                               ttGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAAAGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCCACCTCTTCCATACCTTCATCGCTGAAAGTACGTAGACTATGGTTCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGAACTCTATCCAGCTGAGCTACGGCTGCCAATg                                                  >  1:520011/1‑251 (MQ=17)                                                                                                                                                                                                    aCCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCCACCTCTTCCATACCTTCCTCGCTGAAAGTACGTAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGa                                                                  >  1:151055/1‑230 (MQ=255)                                                                                                                                                                                                    aCCGCTACTAAAGCCGTATCAAGTAAGAGAGGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCGTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCCACCTCTTCCATACCTTCATCGCTGAAAGTACGTAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGa                                                                  <  2:151055/230‑1 (MQ=255)                                                                                                                                                                                                     ccGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCCACCTCTTCCATACCTTCATCGCTGAAAGTACGTAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCAAATGGCGGTg                                            <  1:158069/251‑1 (MQ=255)                                                                                                                                                                                                        ctactaAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCCACCTCTTCCATACCTTCATCGCTGAAAGTACGGAGGATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGACGGAAATGGCGGTGAgg                                         >  1:565951/1‑251 (MQ=255)                                                                                                                                                                                                                gCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCGCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCCACCTCTTCCATACCTTCATCGCTGAAAGTACGTAGAATATGGTGCATCCGGGAGGATTCGAGCCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCAAATGGCGGTGAGGCGGGGAtt                                 <  2:440414/251‑1 (MQ=255)                                                                                                                                                                                                                   gTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCCACCTCTTCCATACCTTCATCGCTGAAAGTACGTAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCAAATGGCGGTGAGGCggggggttg                               >  1:576332/1‑245 (MQ=255)                                                                                                                                                                                                                   gTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCCACCTCTTCCATACCTTCATCGCTGAAAGTACGTAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCAAATGGCGGTGAGGCGGGGATTCGa                              <  1:296013/251‑1 (MQ=255)                                                                                                                                                                                                                        aaGTAAGAGATGGTGCATCCGGGAGGATTACTCTCCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCCACCTCTTCCATACCTTCATCGCTGAAAGTACGTAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCAAATGGCGGTGAGGCGGGGATTCGAAcccc                         <  2:576332/251‑1 (MQ=255)                                                                                                                                                                                                                         aGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCCACCTCTTCCATACCTTCATCGCTGAAAGTACGTAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCAAATGGCGGTGAGGCGGGGATTc                                <  1:130394/243‑1 (MQ=255)                                                                                                                                                                                                                         aGTAAGAGATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCCACCTCTTCCATACCTTCATCGCTGAAAGGACGTAGAATATGGGGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCAAATGGCggggggggggggggtctagttcgg                        >  2:130394/1‑229 (MQ=17)                                                                                                                                                                                                                         aGTAAGAGAGGGTGCGTCCGGGATTATTACTCCCTTTACGGCTCGCCCTTCGGGTCGTTG‑ATAAAGCAACGTTATCCTCCCTGGTGACTGCGGTTATTACCGCCACTTCTTCCATACCTTCATCGCGGAAAGTACGTAGAATATGGTGCATCCGGGAGTATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCAAATGGCGGTGAGGCGGGGATTCTAACCCCg                        <  1:283669/251‑1 (MQ=18)                                                                                                                                                                                                                             agagATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGGCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCCGGGGCCTGCGGGTATTACCGCCACCTCTTCCATACCTTCATCGCTGAAAGAACGTAGAATATGGTGCATCCGGGAGGAGTCGAAACTCCGCCCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTCCGGCTGCCAATGGCGGGGAGGCggggggtttcaccacccgtg                    >  2:281935/1‑235 (MQ=17)                                                                                                                                                                                                                             agagATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTTGTGCCTGCGGTTATTACCGCCACCTCTTCCATACCTTCATCGCTGAAAGTACGTAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCAAATGGCGGTGAGGCGGGGATTCGAACCCCGGATg                    <  1:561101/251‑1 (MQ=255)                                                                                                                                                                                                                             agagATGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCCACCTCTTCCATACCTTCATCGCTGAAAGTACGTAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCCAATGGCGGTGGGGCgggggggtgagacccgga                      >  1:128728/1‑235 (MQ=17)                                                                                                                                                                                                                               agaTGGTGCATCCGGGAGGATTACTCACCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCCACCTCTTCCATACCTTCATCGCTGAAAGTACGTAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCAAATGGCGGTGAGGCGGGGGTTTGAACCCCGGGTGcc                  >  1:595752/1‑250 (MQ=255)                                                                                                                                                                                                                                    ggtcagccGGGAGGTTTACTCACCTTGCTGCTCGCTCTTCGGGCATTGG‑CTAAAGCAACGTTATACTCCCTGGTGCTTGCTGTTATTTCCGCCACCTTTTCCATATCTTCAACGCTGAAAGTACGTAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCAAATGGCGGTGAGGCGGGGATTCGAACCCCGGATGCAGCttt             <  2:572419/245‑1 (MQ=18)                                                                                                                                                                                                                                      gCATCCGGGATGATTACTCCCCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCACCGTTACCCTCCCTGGTGCTTGCGGTTATTACCGCTACCTCTTCCATACCTTCATCGCTGAAAGTACGTAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCAAATGGCGGTGAGGCGGGGATTCGAACCCCGGATGCAGCTTTTg           <  1:673197/251‑1 (MQ=255)                                                                                                                                                                                                                                               aggaTTACTCACCGTGCGGCTCGGCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCCACCTCTTCCATACCTTCATCGCTGAAAGTACGTAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCAAATGGCGGTGAGGCGGGGATTCGAACCCCGGATGCAGCTTTTGACCGCATAc  <  1:703317/251‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  ATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAG‑CTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTG‑CTAAAGCAACGTTATCCTCCCTGGTGCCTGCGGTTATTACCGCCACCTCTTCCATACCTTCATCGCTGAAAGTACGTAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCAAATGGCGGTGAGGCGGGGATTCGAACCCCGGATGCAGCTTTTGACCGCATAC  >  NC_000913/2818115‑2818602

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