BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,173,363	Δ2 bp	intergenic (‑490/+918)	gatD ← / ← gatB	galactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,173,363	0	C	.	82.8%	370.8 / 80.8	116	intergenic (‑490/+919)	gatD/gatB	galactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand): ref base C (0/20); new base . (60/36); total (60/56)
Fisher's exact test for biased strand distribution p-value = 5.59e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	NC_000913	2,173,364	0	C	.	82.8%	377.0 / 74.6	116	intergenic (‑491/+918)	gatD/gatB	galactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand): ref base C (0/20); new base . (60/36); total (60/56)
Fisher's exact test for biased strand distribution p-value = 5.59e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 125 total aligned reads displayed.

GCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGAC‑ACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCG  >  NC_000913/2173123‑2173613

||

  aTTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGAC‑ACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTccc                                                                                                                                                                                                                                               >  1:749463/1‑250 (MQ=255)

                  gtgtcaaaaCCAGTAGTTTGCGTCGCGACCCACAGGGGGATGCTCATAATGA‑GACCACCCGAGATTAGGGTGCGGAACAGATTTCCACGAGGCGCGGCGACCGCCATCGCCACGAAGAAGCCGATGGGGGCAAGATC‑GCCAAACGGCAGCACCGGATTCCCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATc                                                                                                                                                                                                                              <  2:10616/245‑1 (MQ=255)

                    gtgCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGAC‑ACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCa                                                                                                                                                                                                                             >  2:631396/1‑250 (MQ=255)

                     ggCAAACCAATAGTTGGGGTCGCGATCAACAGGGGGAGTCTAATAATGAGGCC‑ACCCGAGATAAAGGTGCGGAACAGATTCCCACGATGCCGGGCGACCGCCATCGCCCCGAAGAGGCCGATGGGGGCAAGATC‑GCAAAACGGCAGCACCTGATTCCCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTGGCGGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGcctgttcttg                                                                                                                                                                                                                           <  2:102247/250‑9 (MQ=255)

                       cAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGAC‑ACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATAAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGc                                                                                                                                                                                                                         <  1:739702/251‑1 (MQ=255)

                        aaaCCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGAC‑ACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGCTACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                                                                                                                                                                                                                        >  2:540098/1‑251 (MQ=255)

                          aCCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGAC‑ACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCa                                                                                                                                                                                                                       >  1:799411/1‑250 (MQ=255)

                           ccAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGAC‑ACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAt                                                                                                                                                                                                                     <  1:1158416/251‑1 (MQ=255)

                               gaGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGAC‑ACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg                                                                                                                                                                                                                 >  2:219564/2‑251 (MQ=255)

                                     gCGTCGCGATCCACAGGGTGATGCTCATAATGATGAC‑CCCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTc                                                                                                                                                                                                            <  1:527359/250‑1 (MQ=255)

                                      cgtggcgATCCACAGGGTGAGCCACATAATGAGGAC‑ACCCGAGATTAAGGGGCGGAACAGTTTTCCACGATGCAGGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCt                                                                                                                                                                                                          <  2:1099788/247‑1 (MQ=255)

                                      cGTCGCGATACACAGGGTGATGCTCATAATGATGAC‑ACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GCGAGTGGGATAAAAATCAGGCTTGCCGATACAACC‑‑GACGTATGTCCCAGCAGCACCGCCGGAACAAGGCCAATCAGGAACTCCt                                                                                                                                                                                                          >  2:601671/1‑251 (MQ=255)

                                       gTCGCGATCCACAGGGAGATGCTCATAATGATGAC‑ACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTg                                                                                                                                                                                                         >  1:535226/1‑251 (MQ=255)

                                         ggcgATCCACAGGGTGATGCTCAAAATGATGAC‑ACCCGAGATTAAGGTGCGGAACAGATTCCCACGATGCACGGCGCCCGCCATCGCCACGAAGAAGCCGATGGGGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GAAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAc                                                                                                                                                                                                       <  2:749463/250‑1 (MQ=255)

                                            gATCCACAGGGTGATGCTCATAATGATGAC‑ACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGAGCAAGGCCAATCAGGAACTCCTGAccgc                                                                                                                                                                                                    >  2:278839/1‑251 (MQ=255)

                                             aTCCACAGGGTGATGCTCATAATGATGCC‑ACCCGCGATTAAGGTGCGCAACAGATTTCCACGATGCCCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTTGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGCGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCACTCAGGAACTCCTGAccgcc                                                                                                                                                                                                   >  1:300996/1‑251 (MQ=255)

                                             aTCCAAAGGGGGATGCTCAAAATGAGGAC‑CCCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACT‑GCAATTAAAATG‑GTGAGGGGGATAAAAACCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgcc                                                                                                                                                                                                   <  1:820812/251‑1 (MQ=255)

                                               ccACAGGGTGATGCTCATAATGATGAC‑ACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgccg                                                                                                                                                                                                  >  1:431512/1‑250 (MQ=255)

                                                        gATGCTCATAATGATGAC‑ACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGAAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCAGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTAGCCt                                                                                                                                                                                        >  1:69549/1‑251 (MQ=255)

                                                              cATAATGATGAC‑ACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCCGATTCCCCGGCAAACAAACA‑GCAATTAAACTG‑GTGAGTGGGATAAAAATCAGGCTTCCCGATACCACC‑‑GCCGTATGTCCCAGCCGCAACGCCGGATCAAGGCCAATCAGCAACTCCTGACACCCGCCCTTCGCCTGTAAAc                                                                                                                                                                                  >  2:1159509/1‑251 (MQ=255)

                                                                  atgatgAC‑ACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACt                                                                                                                                                                              >  2:734245/1‑251 (MQ=255)

                                                                   atgataccacacGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGGCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAAATAAAATG‑GTGCGGGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGAACAAGGCCACTCAGCAACTCCTGACCGCCGGACTTCGCCTGTAAACGACt                                                                                                                                                                              >  2:84108/12‑251 (MQ=255)

                                                                     atgaC‑ACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACg                                                                                                                                                                           >  2:274955/1‑251 (MQ=255)

                                                                       gaC‑ACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACg                                                                                                                                                                           >  1:17424/1‑249 (MQ=255)

                                                                            cccaagaTTAAGTCGCGGAACAGATTCCCCCGATGCGCGCCGCCGGCCATCGCCAGAAAGAAGCCGATGGGGGCAAGATC‑GCCAAACGGCATCACTTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGAATAAAAACCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACTAGCCTg                                                                                                                                                                     <  2:343556/247‑1 (MQ=255)

                                                                               gagaTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCGGACCGCGGAACTTCCCCTGAAAACGACTACGAGCCTGCtt                                                                                                                                                                  >  2:574512/1‑251 (MQ=255)

                                                                                   ttAAGGTGCGGAACGGATTTCCCCGATGCACGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACTTGATTACCCGGCACACAAACA‑GAAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGc                                                                                                                                                               <  1:37605/250‑1 (MQ=255)

                                                                                    cctgtccgggAACAGATGCCCACGATGCCGGGCGACCGCCATCGCCACGAAGAAGCCGAGGGTGGCAAGATC‑CCCAAACGGCAGCACCTGTTTACCCGCCACACAAACA‑CCAATAAAAAGG‑GGGAGGGGGAAAAAAACCAGGCTCGCCTATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCGGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCg                                                                                                                                                              <  2:490378/242‑1 (MQ=255)

                                                                                     aaGGTGCGGAACAGATTTCCACGATGCACGGCGCCCGCCATCGCCACGAAGAAGCCGAGGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACAAAAACA‑GCAATTAAAAGG‑GTGAGTGGGATAAAAATCAGGCTTGCCAATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                                                                                                                                            <  1:960449/251‑1 (MQ=255)

                                                                                       ggTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGCACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATgg                                                                                                                                                          >  2:806434/1‑251 (MQ=255)

                                                                                        gTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATggg                                                                                                                                                         >  1:896435/1‑251 (MQ=255)

                                                                                         ggCGGAACAGATTCCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGc                                                                                                                                                        <  1:835573/250‑1 (MQ=255)

                                                                                          gCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCg                                                                                                                                                       >  2:93385/1‑251 (MQ=255)

                                                                                            ggAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCg                                                                                                                                                              >  1:490378/1‑242 (MQ=255)

                                                                                             aaacaaatccccgcgacgcgcGGCGCCCGCCAGCGCCCCGAAGAACCCGATGGTGGCAAGCTC‑GCCAAAGGGCGGCCCCTGATTCCCGGGCACACAAACA‑GCAATAAAAAGG‑GTGAGGGGAAAAAAAATCAGGCTTGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTa                                                                                                                                                    <  2:457658/231‑1 (MQ=255)

                                                                                               aaataTTTCCACGTTCCCCGGCGACCGCAATCGCCCCGAAGAAGCCGATGGTGGCAA‑ATCCGCAAAACGGCAGCACCTGATTACCCGGCACACAAAAA‑GCAATTAAAATG‑GTGAGGGGGAAAAAAACCAGGCTTGCCGAT‑CC‑CCCCGCCTTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                                                                                                                                                  <  2:732397/247‑1 (MQ=255)

                                                                                               aCAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaa                                                                                                                                                   >  2:147419/1‑250 (MQ=255)

                                                                                                  gATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCGACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAcc                                                                                                                                                >  1:471351/1‑250 (MQ=255)

                                                                                                  gATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAcc                                                                                                                                                >  1:422521/1‑250 (MQ=255)

                                                                                                       ccACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatccat                                                                                                                                          <  1:666386/251‑1 (MQ=255)

                                                                                                        cccgaggCACGGCGACCGCCATCGCCACGAAGAAGCCGAGGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatg                                                                                                                                         <  1:540098/245‑1 (MQ=255)

                                                                                                             tGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATggg                                                                                                                                    >  2:248703/1‑251 (MQ=255)

                                                                                                             ggCACGGCGACCGCCATCGCCACGAAGAAGCGGATGGTGGAAAGATC‑GCCAAACGGCAGCACCTGATTCCCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGAATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATggg                                                                                                                                    <  1:534352/250‑1 (MQ=255)

                                                                                                                  ggcacccGCCATGGCCACGAAGAAGCAGATGGGGGAAAGACC‑GCAAACCGGCAGCCCCTGATTCCCCGGCACACAAACA‑GCAATAAAAAGG‑GTGAGTGGGATAAAAATCAGGCTTCCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTaa                                                                                                                               <  2:684886/246‑1 (MQ=255)

                                                                                                                     gcccGCCATCGCCCCGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCAACTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGAAACGCCGGATCAAGGCCAATCAGGAACTCCCGACCGCCGAACTTCGCCTGTAAACGACTCCGAGTCTGCTTAGCGATGGGCGGTAAACCATCAATGATGGGTTTAATcc                                                                                                                            >  2:602898/3‑250 (MQ=255)

                                                                                                                       ccGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAcc                                                                                                                          <  1:274955/251‑1 (MQ=255)

                                                                                                                          ccATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGt                                                                                                                       <  1:390852/251‑1 (MQ=255)

                                                                                                                            cgcGCCCCGAAGAAGCCGAGGGTGGCAAGATC‑GCAAAACGGCAGCTCCTGATTACCCGGCACACAAAAA‑GCAATTAAATGG‑GTAAGTGGGATAAAAATCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTgg                                                                                                                     <  1:513880/249‑1 (MQ=255)

                                                                                                                             ccgccagaagaaacccGATGGTGGAAAGATC‑GCAAACCGGCAGCACATGATAACCGGGCACACAAACA‑GCAATTAAAATG‑GTGAGGGGGATAAAAACCGGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGc                                                                                                                    <  2:384941/237‑1 (MQ=255)

                                                                                                                                ccACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCAtt                                                                                                                 <  1:1081215/251‑1 (MQ=255)

                                                                                                                                   cGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCGGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagc                                                                                                              <  1:489600/251‑1 (MQ=255)

                                                                                                                                       gaagCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACACACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagcagc                                                                                                           >  2:173385/1‑250 (MQ=255)

                                                                                                                                       gaagCAGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCCGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGACGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagcagc                                                                                                           >  2:404308/1‑250 (MQ=255)

                                                                                                                                           cgtgtttgggtaataTC‑GCCAAACGGAGGCGCGTGTTTACCGGGCACACAAACA‑GCAATTAAA‑TGTGTTGTTGAAAAAAAAAGCAG‑CTTGCCCGTACC‑CCCCGCTGTATGTCCCAGCAGCAACGCCGGATAAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGTCTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCAc                                                                                                      <  2:673719/237‑1 (MQ=255)

                                                                                                                                            cGATGGGGCCAAGATC‑GCCAAACGATAGCTCCTAATTACCCGCCACAAAAACA‑GCAATTAAATTG‑GTCAGTGGAAAAAAAATAAGCCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGAAACTCCTAACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACt                                                                                                     <  1:1137315/251‑1 (MQ=255)

                                                                                                                                               tggtggCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgcc                                                                                                  >  1:501532/1‑251 (MQ=255)

                                                                                                                                                  gggaaaGATC‑GCAAAACGGCAGCCCCTGATTCCCCGGCACACAAACA‑GAAATTAAAATG‑GTGAGTGGAAAAAAAATCAGGCTGGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgccgcc                                                                                               <  2:353260/247‑1 (MQ=255)

                                                                                                                                                  gggaaaGAGC‑GCCAACCGGCAGCCCCTGATTACCCGTCACACAAACA‑GCAATAAAAATG‑TGGTGTGGGATAAAAATAAGGCTTCCCGA‑ACC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgccgcc                                                                                               <  2:586814/247‑1 (MQ=255)

                                                                                                                                                           c‑GCCAAACGGCAGCACCTGATTCCCCGGCACAAAAACA‑GCAATAAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCGGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTAcc                                                                                      <  1:1150358/251‑1 (MQ=255)

                                                                                                                                                             gCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCAGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAAATCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCg                                                                                     >  2:1000304/1‑251 (MQ=255)

                                                                                                                                                                aaaCGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGcc                                                                                   >  1:677959/1‑250 (MQ=255)

                                                                                                                                                                          aCCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATac                                                                        >  1:477972/1‑251 (MQ=255)

                                                                                                                                                                          aCCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATac                                                                        >  1:734271/1‑251 (MQ=255)

                                                                                                                                                                            cTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATacac                                                                      >  1:1081629/1‑251 (MQ=255)

                                                                                                                                                                                 tccccGGCACACAAACA‑GCAATTAAAAGG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCttt                                                                  <  1:631396/248‑1 (MQ=255)

                                                                                                                                                                                    ccGGCACA‑AAACAAGCAATAAAAAGG‑GTGGGGGGGAAAAAAGCCGGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACa                                                              <  1:102566/251‑1 (MQ=255)

                                                                                                                                                                                       gccccaaaaacacaaataaaaaGG‑GTGAGGGGAAAAAAAATCAGGCTTGCGG‑‑ACCCCCCCGCCGTATGTCCCAGCCGCAACGCCGGATAAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACTCCTTTGACATc                                                            <  2:896435/233‑1 (MQ=255)

                                                                                                                                                                                        cacacaAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGt                                                          >  1:179783/1‑251 (MQ=255)

                                                                                                                                                                                           acaAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAAc                                                       >  1:47324/1‑251 (MQ=255)

                                                                                                                                                                                            aaaaCA‑GCAATAAAAAGG‑TTGAGGGGAAAAAAAATCAGGCTGGCCGA‑ACC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAAcc                                                      <  1:147419/250‑1 (MQ=255)

                                                                                                                                                                                             aaccA‑GCAATTAAAATG‑GTGAGTGGGATAAAAACCAGGTTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAAccc                                                     <  1:806434/248‑1 (MQ=255)

                                                                                                                                                                                              aaCA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGCCTGGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCcg                                                    <  1:248703/251‑1 (MQ=255)

                                                                                                                                                                                                cA‑GCAATTAAAATG‑GGGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGGTAAACCATCCATGAGGGGGTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACATTTGACATCGTAACCcgcg                                                  >  2:471466/1‑251 (MQ=255)

                                                                                                                                                                                                     aaTTAAAATG‑GTGAGTGGGAAAAAAATCAGGCTTCCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGa                                              <  1:33606/251‑1 (MQ=255)

                                                                                                                                                                                                      ataaaaaGG‑GTGAGTGGGAAAAAAATCAGGCTTGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGAt                                             <  1:284271/248‑1 (MQ=255)

                                                                                                                                                                                                        aaaaaGG‑GTGAGTGGAATAAAAATCAGCCTTCCCGATAC‑‑CTCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCTTAACCCGCGAGGATGc                                           <  1:561338/250‑1 (MQ=255)

                                                                                                                                                                                                        aaaaaGG‑GTGAGGGGGATAAAAATCAGGCTTGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATg                                            <  1:173385/249‑1 (MQ=255)

                                                                                                                                                                                                         aaaaTG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGcc                                          >  1:638381/1‑251 (MQ=255)

                                                                                                                                                                                                           aaTG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCg                                         >  2:522308/1‑250 (MQ=255)

                                                                                                                                                                                                            aTG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGAt                                       >  1:1184354/1‑251 (MQ=255)

                                                                                                                                                                                                              g‑gTGAGTGGGATAAAAACCAGGCTTGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATaa                                     <  1:742412/251‑1 (MQ=255)

                                                                                                                                                                                                                 tGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATc                                   >  1:630497/1‑251 (MQ=255)

                                                                                                                                                                                                                  gAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATc                                   >  1:41996/1‑250 (MQ=255)

                                                                                                                                                                                                                      gggATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGcc                                          >  2:325337/1‑239 (MQ=255)

                                                                                                                                                                                                                      gggAAAAAAATCGGGCTTGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGcc                                          <  1:325337/239‑1 (MQ=255)

                                                                                                                                                                                                                       ggATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAg                                                                                 >  2:884034/1‑199 (MQ=255)

                                                                                                                                                                                                                        gATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGcccc                             >  1:1030358/1‑250 (MQ=255)

                                                                                                                                                                                                                        gATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGcccc                             >  2:459268/1‑250 (MQ=255)

                                                                                                                                                                                                                         aTAAAAAGCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTTACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCAt                           <  1:219564/251‑1 (MQ=255)

                                                                                                                                                                                                                            aaaaTCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACACGCGAGGATGCCGATAATCAGCCCCATcac                        >  1:349652/1‑251 (MQ=255)

                                                                                                                                                                                                                              aaaaacaggttgcccgacccccccGCCGTATGTCCCAGCAGCAACGCGGGATCAAGCCCAATCAGGAACTCCTGCCCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATcac                        <  2:734271/231‑1 (MQ=255)

                                                                                                                                                                                                                               aaaccggcttgcccgatcccaccgcctTATGTCCCAGCAGCAACGCCGGATCAAGGCCAACCAGGACCTCCTGCCCGCCGACCTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATcaca                       <  2:799411/230‑1 (MQ=255)

                                                                                                                                                                                                                                      ttGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACg              >  2:811077/1‑251 (MQ=255)

                                                                                                                                                                                                                                         tcccgaaccccccGCCGTAGTTCCCGGCAGCACCGCGGGACAAAGGCCACTCAGGAACCCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCAc               <  2:422521/240‑1 (MQ=255)

                                                                                                                                                                                                                                             gtcacccccGCCGTAGGTCCCAGCAGCACCGCCGGATAAAGGCCAATCAGAAACCCCTGCCCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGAc         <  2:638381/246‑1 (MQ=255)

                                                                                                                                                                                                                                             gataccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTgcc         >  2:938261/8‑249 (MQ=255)

                                                                                                                                                                                                                                             gataccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGAc         >  1:684269/8‑251 (MQ=255)

                                                                                                                                                                                                                                                ccgcccGCCGGTTGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGACTTCCTGACCCCCGAACTTCGCCTGTAAACGACTGCGAGCCTGCTTGGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACAGCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAgg      <  2:192296/248‑1 (MQ=255)

                                                                                                                                                                                                                                                   cccGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTc   <  1:93385/251‑1 (MQ=255)

                                                                                                                                                                                                                                                    ccgccgTATGTCCCAGCAGAAACGCCGGATCAAGGCCAATCAGGACCCCCGGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCg  <  2:179783/251‑1 (MQ=255)

||

GCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGAC‑ACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATC‑GCCAAACGGCAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATG‑GTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCG  >  NC_000913/2173123‑2173613

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