BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,579,943	G→A	G103D (GGC→GAC)	tktB →	transketolase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,579,943	0	G	A	97.2%	284.5 / ‑5.6	110	G103D (GGC→GAC)	tktB	transketolase 2
Reads supporting (aligned to +/- strand): ref base G (1/2); major base A (61/42); minor base C (2/1); total (64/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 110 total aligned reads displayed.

ACTCTGGTCATCCCGGCGCGCCGATGGGCATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTGATTGG‑TTTTTACGATCACAACGG  >  NC_000913/2579700‑2580189

aCTCTGGTCATCCCGGCGCGCCGATGGGCATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACg                                                                                                                                                                                                                                                     <  1:210847/251‑1 (MQ=255)

  tctGGTCATCCCGGCGCGCCGATGGGCATGGCTGATATTGCCGCAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCAAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTCCATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGACGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACt‑atcatcc                                                                                                                                                                                                                                                   >  2:340289/1‑246 (MQ=255)

     ggTCATCCCGGCGCGCCGATGGGCATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAgg                                                                                                                                                                                                                                                >  1:942555/1‑251 (MQ=255)

     ggTCAGCCCGGCGCGCCGATGGGCATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCTGAGATTGACT‑ATACGCCAg                                                                                                                                                                                                                                                 <  1:371702/250‑1 (MQ=255)

          ccccGGCTCGCCGATGGGCATGGCTGATATTGCCGAAGTGCTGTGTAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATTATCGTG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATTCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTg                                                                                                                                                                                                                                           <  1:418422/250‑1 (MQ=255)

           cccGGCGCGCCGATGGGCATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTAACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCAACCCGGAGATTGACT‑ATACGCCAGGCGTTg                                                                                                                                                                                                                                           >  1:202964/1‑250 (MQ=255)

                cgcgcCGATGGGCATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAacca                                                                                                                                                                                                                                     <  1:932126/251‑1 (MQ=255)

                   gcCGATGGGCATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAaccaccc                                                                                                                                                                                                                                  >  2:706951/1‑250 (MQ=255)

                   gcCGATGGGCATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAaccacca                                                                                                                                                                                                                                  <  1:332283/251‑1 (MQ=255)

                        tGGGCATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGc                                                                                                                                                                                                                             >  1:425235/1‑251 (MQ=255)

                           gcagggcggATATTGCCGAATTGCTGTGAAACGAGTTTCTTAAACATACCCCTTCCGACCCAACCTGGTATGATCGCG‑ACCGCGTT‑ATTCTTTCCAACGGGCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCg                                                                                                                                                                                                                          <  2:716632/243‑1 (MQ=255)

                             aTGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCt                                                                                                                                                                                                                        >  1:698170/1‑251 (MQ=255)

                                 cTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAAAATAACACTACCGACCCAACCTGGTGTGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGCCACGCGTCGATGCTGCTCTACAGTTTGCTAAATCTGAAAGGTTACGACCTGCCGCTGGCAGCACTGAAGAACTTCCGTCAGTTGCATTCGAAAACACCAGGCCCCCCGGAGCTTGACT‑ATACGCCAGGAGTTGCAACCACCCACGGCCCCCTTgg                                                                                                                                                                                                                     >  2:954842/1‑250 (MQ=255)

                                    atatTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACa                                                                                                                                                                                                                  >  1:650185/1‑250 (MQ=255)

                                              aGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCg                                                                                                                                                                                                       >  2:641367/1‑250 (MQ=255)

                                              aGTGCTGGGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCg                                                                                                                                                                                                       >  2:890133/1‑250 (MQ=255)

                                               gTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCAGCTTGGACA‑CGGTTTGGCGa                                                                                                                                                                                                      >  1:184300/1‑250 (MQ=255)

                                                 gCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACg                                                                                                                                                                                                   >  1:545914/1‑251 (MQ=255)

                                                  ctgaggAACGCTTTTCTTAAACATAACCATACCGACCCAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAATCCCAGGCCACCAGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACAGGCCCGCTTGGACA‑AGGTTTGGCGAACGc                                                                                                                                                                                                  >  1:160584/5‑251 (MQ=255)

                                                   tgtgGAACGATTTTCTTAAACATAACCCTACCGACCAAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCATGGACA‑AGGTTTGGCGAACGcc                                                                                                                                                                                                 >  1:806958/1‑251 (MQ=255)

                                                    gtgGAACGATTTTCTTAAACATAACCCTACCGTCCCAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGcc                                                                                                                                                                                                 <  1:280180/250‑1 (MQ=255)

                                                      ggAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGt                                                                                                                                                                                               <  1:284085/250‑1 (MQ=255)

                                                      ggAACGATTTTCTTAAACATAACCCCACCGTCCCAACCTGTAATGATCGCG‑ACCGCTTT‑ATTCTTCCCACCTGGCACGCGTGGATGCTGCTCTATAGTTTGCTACATCTGGCCGGTTACGTCCTGCCGCTTGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTc                                                                                                                                                                                              <  2:583499/251‑1 (MQ=255)

                                                        aaCGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCgg                                                                                                                                                                                            >  1:812019/1‑251 (MQ=255)

                                                              tttCTTAAGCATAACCCTACAGACCCAACCTGGTATGATAGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTACTCTACAGTTTGCTCCATCTGACCGGATCGGCCCTGCCGCTGGAAGAACTGGAGAACTTCCGTCAGTTGCATTCGAAAACCCCCGACCACCCGGAGATTGCCT‑ATACGCCAGGCGTTGAAACCACCCCCTAACCTATTGGGCA‑AGTATTGGCGAACGCCGGCGGGcgggg                                                                                                                                                                                      >  2:551785/1‑247 (MQ=255)

                                                                  ttAAACATAACCCTACCGACCCAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGAt                                                                                                                                                                                   >  1:17022/1‑250 (MQ=255)

                                                                   tAAACATAACCCTACCGACCCACCCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATa                                                                                                                                                                                  <  1:715006/250‑1 (MQ=255)

                                                                    aaaCATAACCCTACCGACCCAACCTGGTATGATCGCG‑CCCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGc                                                                                                                                                                                >  1:803082/1‑251 (MQ=255)

                                                                    aaaCAAAACCCTCCCGACCCAACCTGGTATGATCCCG‑ACCGCTTT‑ATCCTTCCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGc                                                                                                                                                                                <  2:976301/251‑1 (MQ=255)

                                                                      aaaTAACCCTACCCACCCACACTGTTATGATCGCG‑ACCGTTTT‑ATTCGTTCCAACGGTCCCGCGTCGATCCTGCTCTACAGTTGGCTACATCGGCCCGGTTCCTCCCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAACCCCCAGGCCACCCGGATATTGACT‑ATACGCCAGGCGTGGAAACCACCATCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGGCGGTCTGACGATAGCag                                                                                                                                                                              <  2:141683/249‑1 (MQ=255)

                                                                       cATAACCCTACCGACCCAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGATGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑TGGTTTGGCGAACGCCGTCGAGCTGGCGATAGc                                                                                                                                                                                >  1:577218/1‑248 (MQ=255)

                                                                        atacaccttacGACCCAACCTGGTCTGA‑CCCGTACCGCTTT‑AGTCTGCCCAACGGTCACGCGTCGATGATTCACGACAGTTTTCTACATCTTACCGGTTACGACCTGCCGCTGGAAGATCTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGTCT‑ATACGCCAGGCGTTGAAACCGCCACCGGCCCGCTTGGACA‑AGGGTTGGCGAACGCCGTCGGGCTGGCGATAGCagag                                                                                                                                                                            <  2:952936/241‑1 (MQ=255)

                                                                        aTAACCCTACCGACCCAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCag                                                                                                                                                                              >  1:346202/1‑249 (MQ=255)

                                                                          aaCCCTACCGACCCAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCg                                                                                                                                                                          >  1:569158/1‑251 (MQ=255)

                                                                          aaCCCTACCGACCCAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCg                                                                                                                                                                          <  1:811250/251‑1 (MQ=255)

                                                                             ccTACCGACCCAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTac                                                                                                                                                                       <  1:890133/251‑1 (MQ=255)

                                                                             ccTACCGACCCAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCAGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGCTTGGCGAACGCCGTCGGGCTGGCGGTAGCaggacgtac                                                                                                                                                                       >  2:324058/1‑244 (MQ=255)

                                                                              cTACCGACCCAACATGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCACCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCAGTCGGGCTGGCGATAGCAGAGCGTac                                                                                                                                                                       >  1:779418/1‑250 (MQ=255)

                                                                                  cGACCCAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGACACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGGAGACCTTCCGTCAGTTGCATTCGAAAACCCCAGACCACCCGGAGATTGACT‑ATACGCCAGTCGTTGAGACCACCACCGGCCCGCTTGGACA‑AGGCTTGGCGAACACCGTCGGGCTGACGATAGCAGAGCGTACACTgg                                                                                                                                                                  >  1:540238/1‑251 (MQ=255)

                                                                                    aCCCAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCagagggaaaactgggg                                                                                                                                                                >  2:600223/1‑239 (MQ=255)

                                                                                      ccatccTGGTATGACCGCG‑‑CCGCTTTTATTCTTTCCAACGGTCGCGCGTCGATGCTGCTCTACAGTTTGCTACCTCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTggcggc                                                                                                                                                              <  2:350134/247‑1 (MQ=255)

                                                                                      ccAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACAGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGCGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTggcggg                                                                                                                                                              >  1:373877/1‑250 (MQ=255)

                                                                                        aaCCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGgc‑gc                                                                                                                                                           >  1:501908/1‑251 (MQ=255)

                                                                                           cTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGt                                                                                                                                                        >  1:343907/1‑251 (MQ=255)

                                                                                               aaTGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAACCCCCAGGCCAACCGGAGATTGACT‑ATACGCCAGACGTTGAAACCACCACCGGCCCGCTTGGAAA‑AGGATTGGCGAACGACGGCGGGTTGGCGATAGCAGAGGCTCAACTGGCGGC‑GCAATTTa                                                                                                                                                     >  2:256623/2‑250 (MQ=255)

                                                                                                aTGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAAc                                                                                                                                                   <  1:776306/251‑1 (MQ=255)

                                                                                                 tGATCGCT‑AGCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATTCTGATCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGACGAAATGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGCTTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGACGGGCTGGCGATAGCAACACGTCCCCTGGCGGC‑GCACTTTAAc                                                                                                                                                   >  2:110386/1‑250 (MQ=255)

                                                                                                   aTCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAAcc                                                                                                                                                  >  1:79046/1‑249 (MQ=255)

                                                                                                     ccccgcccGCTTT‑ATCTTTTCCACCGGTCACGCGTCGAGGCTGCTTTAAAGTTTGCTACATCTGACCGTTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAAccagc                                                                                                                                               <  2:698170/245‑1 (MQ=255)

                                                                                                        g‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAAccagcc                                                                                                                                              >  2:195521/1‑248 (MQ=255)

                                                                                                               ttt‑ATTCTTTCCAACGTCCCCGCGTCAATGCGTCTCTACAGTGTGCTACATCTGCCCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATg                                                                                                                                     <  2:677050/250‑1 (MQ=255)

                                                                                                                     tCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGAtc                                                                                                                                <  2:167458/250‑1 (MQ=255)

                                                                                                                     tCTTTCCAACGGGCACGCTTCGATGCTGCTCTACATTTTGCTACATCTGCCCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTAGACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGAtc                                                                                                                                <  2:866241/250‑1 (MQ=255)

                                                                                                                         tCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGC‑CCATGAGCtcgtcg                                                                                                                            >  2:777848/1‑250 (MQ=255)

                                                                                                                         tCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGAGAGCAGAGCGTACCCTGGCGGC‑CCAGttacaacaaccagaacctgagatcct                                                                                                                               >  2:267343/1‑224 (MQ=255)

                                                                                                                          ccaacgtgcacgctgcGATGCTGCTCTACAGTTTGCTACAGCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCAGGAGATCGTCGAt                                                                                                                          <  1:168411/236‑1 (MQ=255)

                                                                                                                          ccAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCCGTTTAAACAGCCAGA‑CCAtagaatagtcga                                                                                                                           >  2:860584/1‑239 (MQ=255)

                                                                                                                           cAACGGTCACGCTACGATGCTGCTCTCCAGTTTGCTACATCTGACCGGTTACGACATGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCACAGGCCACCAGGAGATAGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGATTGG‑CACAGGTTTAGCGAACGCCGTCGGGCAGGCGATAGCAGAGGGGACACAGGCGGCCGAAGTT‑AAccagccaagccctgggcacgttgaac                                                                                                                          >  2:632478/1‑232 (MQ=255)

                                                                                                                           cAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTTCATCTGACCGGTTACGGCCTTCCGCTTGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTTACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGATCGTCGATc                                                                                                                         <  1:221444/251‑1 (MQ=255)

                                                                                                                           cAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGATCGTCGATc                                                                                                                         >  1:452072/1‑251 (MQ=255)

                                                                                                                             acgttCACGGGTTGATGCTGCTCGTCATGTTGCTCCATCTGACCGGTTACGACCTTCCGTTGGAAGAACTGAAGAACTGCCGTCAGTTGCATTCGAAACCCCCAGGCCACCCGGAGATTTACT‑ATACGCCAGGCGTGGAATCCACCATCGGCCCTCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCTGC‑GCAGTTTACCCATCCAGA‑CCATGAGATCGTCGATCat                                                                                                                       <  2:840921/247‑2 (MQ=255)

                                                                                                                               ggTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGCTTAACCAGCCAGA‑CCATGAGAtactagatcacgt                                                                                                                     >  1:640812/1‑239 (MQ=255)

                                                                                                                                  cACGCGTCGATGCTGCTCTCCAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGCCA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGACGGC‑GCAGTTTAACCAGCC‑GATCCATGAGATCGTCGATCACTTccc                                                                                                                  >  2:400718/1‑249 (MQ=255)

                                                                                                                                  cACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACAAGCCAGA‑CCATGAGATCGTCGATCACTTCAc                                                                                                                  >  1:499921/1‑251 (MQ=255)

                                                                                                                                  cACGCGTCGATGCTGCTCTACAGTTTGCCACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGATCGTCGATCACTTCAc                                                                                                                  >  1:746750/1‑251 (MQ=255)

                                                                                                                                      cgTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGATCGTCGATCACTTCACCTAt                                                                                                              >  1:762517/1‑251 (MQ=255)

                                                                                                                                      cgTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGATCGTCGATCACTTCACCTAt                                                                                                              >  1:734625/1‑251 (MQ=255)

                                                                                                                                      cgTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTAGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGATCGTCGATCACTTCACCTAt                                                                                                              <  1:154073/251‑1 (MQ=255)

                                                                                                                                       gTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGATCGGCGCTCCCTTCACCtttg                                                                                                             >  2:936745/1‑248 (MQ=255)

                                                                                                                                           aTGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGTAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGATCGTCGATCACTTGACCTATGTGtt                                                                                                         >  2:262130/1‑251 (MQ=255)

                                                                                                                                               tgctCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCAAGA‑CCATTATATCGTCGATCACTTCACCTATGGGTTTATg                                                                                                     >  1:163751/1‑251 (MQ=255)

                                                                                                                                                gctCTACAGGTTGCTATATCTGACCGGTTACGCCCTGCCGCTGGAAGAACTGAAGAACTTACGTCAGTTGCATTCGAAATCCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGATCGTCGATCACTTCACCTATGTGTTTATg                                                                                                     <  1:849695/250‑1 (MQ=255)

                                                                                                                                                      cAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACg                                                                                              >  1:227066/1‑251 (MQ=255)

                                                                                                                                                       aagtggtaacatcggcccGGTTTCGTCCTGCCGCTGGAAGACCTGAAGAACTTCCGTCATTGGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACgg                                                                                             <  2:112577/235‑1 (MQ=255)

                                                                                                                                                            gCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTTAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GAAGTTTAACCAGCCAGA‑CCATGAGATCGTCGATCACTTTACCTATGTGTTTATTGGCGACGGCTGcc                                                                                        >  2:939801/1‑251 (MQ=255)

                                                                                                                                                                   ctgcccGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATg                                                                                   <  2:985600/245‑1 (MQ=255)

                                                                                                                                                                      cctggTTAAGACCTGCCGCTGGAATAACTGAAGAACTTCCGTCAGGTGCATTGGAAAACCCCAGGACACCCGGGTATTGACT‑AGACGCCAGGCGGTGAAACCACCACCGGCCCGCTTGGACA‑AGGCTTGGAGAACGCCGTCGGGCTGGCGAGAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCTAGT‑CCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAgg                                                                              <  2:361654/248‑1 (MQ=255)

                                                                                                                                                                               gACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCa                                                                     <  1:641367/251‑1 (MQ=255)

                                                                                                                                                                                       gCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCt                                                             <  1:340289/251‑1 (MQ=255)

                                                                                                                                                                                        cTGGCAGAACTGAAGAACTTCCGTCAGTTGCATTCGAACACCCCAGGCCACCCGGAGATTGCCT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGCTTAACCAGCCAGA‑CCATGAGAACGTCGATCCCTTCACCTATGTGTTTATGGACGACGGCTGCCTGATGGACGGAATTTCACACGCAGtatg                                                            >  2:889414/1‑248 (MQ=255)

                                                                                                                                                                                        cTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTg                                                            >  1:643743/1‑251 (MQ=255)

                                                                                                                                                                                                    aagaactTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTggcagg                                                >  1:470756/1‑251 (MQ=255)

                                                                                                                                                                                                          ttCCGTAAGTTGCATTCGAAAACCCCAGGCCACCCGGAGCTTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGGTTAAACCGCCAGA‑CCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGTTGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGAAGGGACGCt                                          >  2:626366/1‑251 (MQ=255)

                                                                                                                                                                                                                cctgcgccttccacaaCCCCAGGCCACCCGGAGATTGACT‑ATAGGCCAGGCGTGGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCATATTAACCAGCCAGA‑CCATGAGATCGTTGATCACTTCACCGATGTGTTTATGGGCGTCGGCTGCCTGAGGGAAGGTATTTCCCACTAAGTCTGGTCGCTGGCAGTCACGCTGGGACt                                    <  2:577218/237‑1 (MQ=255)

                                                                                                                                                                                                                     gCATTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGc                                <  1:600223/250‑1 (MQ=255)

                                                                                                                                                                                                                       aTTCGAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAg                             <  1:860584/251‑1 (MQ=255)

                                                                                                                                                                                                                         gcGAAATCCCCAGTCCACCCGGAGATTGACT‑ATTCGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGGGTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGATCGTCGATCTCTTCACCTATGTGTTTATGGGCGGCGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGTGCAAGc                            <  2:551101/249‑1 (MQ=255)

                                                                                                                                                                                                                           gAAAACCCCAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTgg                         >  1:980813/1‑250 (MQ=255)

                                                                                                                                                                                                                             aaaCCCCAGGCCAGCAGGAAATTGACT‑ATACGCCCGGCGTTGAAACCACCACCGGCCCGTCTGGTCA‑AGGTGTGGCGAACGCCGACGGGCTGGCGAGGGCAGTGCGTACCCTGGCGGC‑GCAGTTTCACCAGCCAGA‑CCAGGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGTGTGCCTGATGGAAGGTATTTCCCACGAATTCTGTTAGCTTGCAGGCACGCTGCGCCTGTGCATGCTGAtt                       <  2:362382/251‑1 (MQ=255)

                                                                                                                                                                                                                                ccccAGGCCACCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGCTAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGCCTGGGCAAGCTGCTTggg                    >  1:926293/1‑250 (MQ=255)

                                                                                                                                                                                                                                 cccAGGCCACCCGTAGAT‑GGATCATACCCCAGGCGTTGAACCCACCTCCGGCACGTTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGATCGTCGAGCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTGATTGG‑tt                  <  2:769924/251‑1 (MQ=255)

                                                                                                                                                                                                                                         aCCCGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTGATTGG‑TTTTTACGAt          >  1:71881/1‑251 (MQ=255)

                                                                                                                                                                                                                                            cGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCCCCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTTATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTGATTGG‑TTTTTACGATcac       <  1:87912/251‑1 (MQ=255)

                                                                                                                                                                                                                                            cGGAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCCGTTTAACCAGCCAGA‑CCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTGATTGG‑ATTTTACGATccc       >  2:57353/1‑249 (MQ=255)

                                                                                                                                                                                                                                             ggAGATTGGCT‑ATACGCCAGGCGTTGAAACCACCACCGTCCCGCTTGGACA‑ATGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGttt                                                                                                                                                      <  2:929584/109‑1 (MQ=255)

                                                                                                                                                                                                                                             ggAGATTGGCT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGCTTGGCGAACGCCGTCGGGCTGGCGATAGCCAAGCGTACACTGGCGGC‑GCAGttt                                                                                                                                                      >  1:929584/1‑109 (MQ=255)

                                                                                                                                                                                                                                             ggAGATTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTTTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTGATTGGGTTTTTagttcac        >  2:789259/1‑244 (MQ=255)

                                                                                                                                                                                                                                               agaTTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTGATTGG‑TTTTTACGATCACAAc    <  1:602921/251‑1 (MQ=255)

                                                                                                                                                                                                                                                gaTTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTGATTGG‑TTGTTACGATCACAACg   <  2:942555/251‑1 (MQ=255)

                                                                                                                                                                                                                                                 aTTGACT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTTGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGAATTTCCCACGAAGTCTGTTAGCTGGCAGGCACGCTGGGACTGGGCAAGCTGCTTGG‑TTTTTACTATCCCAACgg  >  2:766698/1‑251 (MQ=255)

ACTCTGGTCATCCCGGCGCGCCGATGGGCATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCG‑ACCGCTTT‑ATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCT‑ATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACA‑AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGC‑GCAGTTTAACCAGCCAGA‑CCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTGATTGG‑TTTTTACGATCACAACGG  >  NC_000913/2579700‑2580189

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