BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,521,124	A→C	12.3%	noncoding (76/76 nt)	valX →	tRNA‑Val

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,521,124	0	A	C	12.3%	217.4 / 10.1	88	noncoding (76/76 nt)	valX	tRNA‑Val
Reads supporting (aligned to +/- strand): ref base A (22/54); new base C (7/5); total (29/59)
Fisher's exact test for biased strand distribution p-value = 5.48e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACTCTCCCTATAATGCGACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTAT  >  NC_000913/2520896‑2521362

gAACTCTCCCTATAATGCGACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGcc                                                                                                                                                                                                                                      <  2:1736229/239‑1 (MQ=255)

 aaCTCTCCCTATAATGCGACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGGTTAACAAGAAAATGAGAAGTGGGTGATTTGCTCAACCGGGATAGCACCTCCCTTACAAGGAGGGGGGCGGG‑GGTTAGA‑TCCCC‑TCATTCCCCCCCCCTTTTTTGCCa                                                                                                                                                                                                                                     >  1:1590544/1‑239 (MQ=25)

   ttgtcgctagattgtgACTTCACTCAGCGTGTGGCATTAGCCCGGCTGGGAGACCACCTTCCTTACAAGGAGGGGGTCGGCGGTTCGATCTCGTCATCACCCCCCAACTCCTTTATTTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCCCTTTCTCGCCAGc                                                                                                                                                                                                                                   <  1:967832/224‑1 (MQ=21)

                gCGACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTTTGTAGGCCCCGCCGGGTAGCAAGAAATTTAGAAGTGGGTGGTTAGCGCAGCTGGTAGAGCACCTACCCCACAAGGAGGGGGGCGGC‑GGGTAGACGCCCG‑TTTTCACCCCCAACttttcgacacataatactcgcttg                                                                                                                                                                                                                       >  1:1072939/1‑218 (MQ=21)

                gCGACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGACCCCGCCGTGGAGCAAGAAATTGAGAATTGGGGGATTAGCGCAGCCGGGAGAGCACCTCCCCTCCAAGGAAGGGGGCGGC‑GGGTCGG‑ACCCG‑TCTTCACCCCCCCCCTTCTTCCCCTCTAAAATTCTTggt                                                                                                                                                                                                                      >  1:1317211/1‑237 (MQ=21)

                 cGACTCCGCACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCTGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTa                                                                                                                                                                                                                      <  2:596092/238‑1 (MQ=255)

                  gACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTATCCGCCGTGTAGAAAGAAAATGAGCAGTTGGTGATTAGCCCAGCTTGGCGAGCACCCCCCTTATAAGGAGGGGGGCGGCGGGTTCG‑‑TCCCGCTCTTAACCC‑CCACTTTCTCTCCCGCTAGCTTTTTTCTaat                                                                                                                                                                                                                    >  1:596092/1‑238 (MQ=21)

                       cacacaGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCTACTTTCTTTAGTCCCCGCCGTGTAGCAAGAAAATTTTAATTAGGTGGTTACCCCTTCTGGGAGTGCACCCCCCCCCCAAGGGAAGGGGCGTC‑GGTATGG‑CCCCC‑TCTTCACCCCCCCCCTTCTCCCCCGCtccacttcgtgggacccttg                                                                                                                                                                                                               >  1:405153/1‑220 (MQ=9)

                       cacacaGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑GCCTCACCCACCACCTTCTCGCCCGCTAAAATTCTTTTaaaaaaaa                                                                                                                                                                                                               >  2:1135294/1‑236 (MQ=255)

                        ccacaGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAAtgtg                                                                                                                                                                                                              <  1:1564829/238‑1 (MQ=255)

                        actcaGCGGGGGTGAGTAGCTCAGCTGAGATGGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAAtgtg                                                                                                                                                                                                              <  1:301689/236‑1 (MQ=255)

                          acaGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGaa                                                                                                                                                                                                            <  2:1590544/239‑1 (MQ=255)

                               gggggTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTa                                                                                                                                                                                                         <  2:1009833/237‑1 (MQ=255)

                                   gggtttAGCTCAGCTGGTTGAGCACCTACCTTACAAGGAGGGGGTCGCGGGTTCGATCCCGTCATCACCCACCAACTACCTTATGTAGACTCCGCCGTGTAGCGAGAAATTGAGAAGTGGGTTATTAGCTCAGCGGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGaa                                                                                                                                                                                                    <  1:569962/234‑1 (MQ=37)

                                   gTGATTAGCTCAGCTGGGAGACGATCTTCCTTACAAGGAGGGGGTCGGCGGTTCGATTCCGTCATCACCCACCAACAACGTTATGGAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGaa                                                                                                                                                                                                    <  1:133546/238‑1 (MQ=37)

                                       ttAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTgg                                                                                                                                                                                                <  2:1557278/238‑1 (MQ=255)

                                            tCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGAt                                                                                                                                                                                           <  2:1074926/238‑1 (MQ=255)

                                                cTGGGAGAGGAGCTCCCTTCCAAGGAGGGGGTCGGCGGTTCGAACTCGTCATCACCCACCAACTACTGTATGGAGCCGCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGc                                                                                                                                                                                       <  2:576602/238‑1 (MQ=37)

                                                cTGGGAGAGCACCTGCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGc                                                                                                                                                                                       <  1:366736/238‑1 (MQ=255)

                                                cTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGGAGTCTCCGCCGTGTAGCAAGGAATTGAGAAGGGGGTGAGTAGCTCAGCTGGGAGAGCACCTCACTTACCACGGGGGGGTCGGC‑GGGTTGG‑TTCCC‑TCTTCAACCCCCCCTCTTTTGCCCGCTTAAATTTTTTTAAAAATGGGAAATAACCCAGTggggggtttgc                                                                                                                                                                                       >  1:1557278/1‑230 (MQ=19)

                                                cTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCAACCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGGGATTAGCTCAGCTGGGAGAGCACCTCCCTTAAAAGGAGGGGGTCGGC‑GGTTCGA‑TGCCG‑TCATCACCCACCCCCTTCTCGCCAGCTAAATTTTTTGGAAAAATTGGGAGGAACGGAGTggggtaatagc                                                                                                                                                                                       >  2:1268179/1‑230 (MQ=25)

                                                 tGGGTGAGCACATCCGTGAAAAGAAGGGGGTCGGCGGATCGATCTCTTCATCACCCACCAGCTGCGTTATGTAGTCTCCGCCGTGTAGCAAGAGATTGAGAAGTGGGTGAGTAGCTCAGCGGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCt                                                                                                                                                                                      <  1:734683/238‑1 (MQ=37)

                                                 gtggAGAGAACCGCCTTTACAAGGAGGGGGTCGGCGGTTCGATCGCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTTAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCt                                                                                                                                                                                      <  2:941105/236‑1 (MQ=255)

                                                          aCCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGG‑TCCCC‑TCATCACCCACCCCTTTCTCCCCAGCTAAATTTTTTGTAAAAATGGTAAGTACCGAAAAGGGGGAATAGCGCCGCTGGgag                                                                                                                                                                            >  2:1646815/1‑239 (MQ=37)

                                                            cTCCCTGACTAGGAGGGGGTCGGCGGTTCGATCCCGTTATAACAGACCAACTAGTTTATGTAGTTCCCGCCGTATAGCGTCAACTGGAGCGGCGGGTGTGTAGAAGAGCTAGGAGAGCCCATCCTCTTCATGGAGTGGGTGGCG‑GGTGCGC‑TCTCT‑CCATAGCACACCCCTTGCTCGATAGTTAAATTTCTCGAAAAAGCGTGAAGTCCCGAAGTGGGTGATTAGCTCAGCTGGGagag                                                                                                                                                                          <  1:26866/239‑1 (MQ=255)

                                                               gcTTACGAGGAGGGGGTCGGCGGTTCGTGCCCGGCCTCACCCACCAACTACTTTATGGAGTGTCCGCCGTGTAGCAAGAAATTGAGAAGGGGGTGATTAGCGCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCAc                                                                                                                                                                       <  1:1135294/238‑1 (MQ=38)

                                                               ccTTAGAAGGATGGGGTCGGCGGTTCGATCCCGTCAGCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCAc                                                                                                                                                                       <  1:758713/239‑1 (MQ=255)

                                                               ccTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑TGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCAc                                                                                                                                                                       <  2:867163/239‑1 (MQ=255)

                                                                 ttctaaGGAGGGGGTCGACCGTTCGCTCCCGTGATCACCCACCAACTTCTTTATGTAGTCTCCGCCGTGTAGCAATAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTCCCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCt                                                                                                                                                                     <  1:891231/235‑1 (MQ=37)

                                                                 ttACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTgg                                                                                                                                                                                                <  2:279680/212‑1 (MQ=255)

                                                                 ttACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTGTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCt                                                                                                                                                                     <  2:787516/239‑1 (MQ=255)

                                                                    cAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTccc                                                                                                                                                                  <  2:527584/239‑1 (MQ=255)

                                                                         aggGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTAc                                                                                                                                                              <  2:1146358/238‑1 (MQ=255)

                                                                         aggGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTAc                                                                                                                                                              <  2:1710021/238‑1 (MQ=255)

                                                                                                  atcaCCCACCAACTGCTTTATGGAGTCTCCCCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCACCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTAGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGAt                                                                                                                                       <  2:290508/236‑1 (MQ=255)

                                                                                                     acccaccAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGc                                                                                                                                                                                  >  2:1343228/1‑190 (MQ=255)

                                                                                                     acccaccAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGc                                                                                                                                                                                  <  1:1343228/190‑1 (MQ=255)

                                                                                                            aaCTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGGGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGGCGGCGGTTTGATTCCGTCATcacc                                                                                                                           >  2:980921/1‑238 (MQ=255)

                                                                                                              ctactTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCCCTTTCTCGCCCGCTCAAATTCTTGTAAAAATGTGAAGTACCGAAATGGGGGGTTAGCTCCGCTTGGCGAGCACCCCCCCTACAAGAAAGGGGTCGGCGGTTTCATACCCTGATcacccc                                                                                                                         >  1:695590/1‑237 (MQ=25)

                                                                                                               tactTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCt                                                                                                                                                                                      >  1:1552406/1‑176 (MQ=255)

                                                                                                               tactTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCt                                                                                                                                                                                      <  2:1552406/176‑1 (MQ=255)

                                                                                                                 ctTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTgg                                                                                                                                                                               <  2:1453388/181‑1 (MQ=255)

                                                                                                                        gaggttctgccgTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCAAGACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCgg                                                                                                                <  1:1040455/229‑1 (MQ=255)

                                                                                                                           tctcCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCCCCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGGGATTAGGTTCGCTGGGAGAGCACCCCCCTTACAAAAAAGGGGTCCGCGGATCGGCTCCCTCATCACCCCCCACCTccgcgcgc                                                                                                           >  1:922153/1‑232 (MQ=25)

                                                                                                                           tctcCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTgg                                                                                                                                                                               <  2:1645645/171‑1 (MQ=255)

                                                                                                                           tctcCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTgg                                                                                                                                                                               >  1:1645645/1‑171 (MQ=255)

                                                                                                                           tctcCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGGTTGGTCCCGGCATCACCCACCCCTTTGGGGCgg                                                                                                           >  2:1040455/1‑238 (MQ=255)

                                                                                                                           tctcCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGCAAAAATGTGGAGTACCGAATAGGGTGGTTAGCTCCGGTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGGTTCCGTCGTCCACccccacattccgcgcg                                                                                                           >  1:761008/1‑229 (MQ=25)

                                                                                                                            ctcCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGAGGGTGATTAGCTCAGCTGTGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATACCGTTATCCCCCCCCAATTCCGGTCttt                                                                                                          >  2:154307/1‑239 (MQ=37)

                                                                                                                               cgccgTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGT‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCAc                                                                                                                                                                       >  2:273650/1‑175 (MQ=255)

                                                                                                                               cgccgTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGT‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCAc                                                                                                                                                                       <  1:273650/175‑1 (MQ=255)

                                                                                                                                gccgTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGAC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGAGTCGGCGGGTCGATTCCGTTATCACCCACCACTTCGGGGCGTTAGc                                                                                                       >  2:1519371/1‑238 (MQ=255)

                                                                                                                                    tgtAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TACCG‑TCATCACCCACCAATTTCTCGCCAGATAAATTTTTTGTAAAAAAGTGAAGTACCGAAGTGGGGGGTTAGGTCAGTTGGGAGAGCACCTCCCTTAACAGGAGGGGGTCGGCGGTTCGGTCCCCTGATcacccacccccttccggggttttctccgt                                                                                                  >  1:1623700/1‑218 (MQ=25)

                                                                                                                                          aaGAAATGGAAACGTGGGGGATTAGGTCTGCAGAGGGAGCACCTCCCTTATACAGAGGTGGTCTGC‑GGTTCGA‑TACCG‑TCATCACCCACCACTTTCACGCCAGCTAAACTCCTTGTAAAAATGAGAAGTGCCGGAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTa                                                                                             <  1:1406218/238‑1 (MQ=21)

                                                                                                                                              aaTTGAGAAGTCGGGGATGAGCGCCGCTGGGAGATCGCGTCCCTTCTAAGGAGGGGGTCGGG‑GGTTTGA‑TCCCG‑TCATCACCCACCCCTTTCTCGACAGGTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATGAGCTCAGCTGGGAGAGCACCCCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTagag                                                                                          <  1:76368/237‑1 (MQ=25)

                                                                                                                                                ttGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCAACCCCCTTACAAGGAGAGGGGCGGCGGTTGGATCCCTTCATTAACCACCACCTCCGGAGGTTAgacccgctcggagaggagt                                                                                      >  2:293202/1‑221 (MQ=25)

                                                                                                                                                 ttagaAGTGGGTGATTAGCTCAGCTGGGAGAGCGCCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTTCCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGtt                                                                                     <  2:460470/237‑1 (MQ=255)

                                                                                                                                                          ggTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGGGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCAGACCCCCCCACCTCCGGACGTTTGCACACGTGGGAGAGCAGGTTACTTTTaa                                                                            >  1:919674/1‑239 (MQ=25)

                                                                                                                                                          ggTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGGGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCCGCATCACCCCCCCCTTCGGGGCGTTAGCTCCGTTGGGAGGGCAGTTGACctttta                                                                            >  2:1735606/1‑237 (MQ=25)

                                                                                                                                                                      aGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGC‑TCCCG‑TCATCACCCCCCACTTTCTCGCCAGCTAAATTTTTTGTAAAAATGTGAAGTACCGAAATGGGGGCTTAGGTCAGCTGCGGGAGCACCCCCCTTACAAGGAGGGGGGCGGCGGGTTGGGCCCGTCATCACCCACCACTTCCGGTCGTTGGGTCCGTTGGGAGGGCGGTTGGCTTTTATTccctttgtgg                                                                 >  1:1509560/1‑229 (MQ=25)

                                                                                                                                                                       gttGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCAACACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTATCGAAGTGGGTGATTAGCTCAGCTGTGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACtt                                                                                <  2:823791/220‑1 (MQ=255)

                                                                                                                                                                         tGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCa                                                               <  2:1623700/237‑1 (MQ=255)

                                                                                                                                                                                 cACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCCGCTAAATTTCTTGTAAAAAGGTGAAGTACCGCAGTGGGGGATTAGCTCCGCTGGGAGAGCAACTCCCCTACAAGGAAGGGGGCGGCGGGTTGGGTCCGTTAACACCCCCCCCTTCACGGCGTTTGCTCCGGTGGGAGAGCCGTTTGCTTTTTATCAATTTGGCCCCGGTTTGAAt                                                      >  1:1205066/1‑238 (MQ=19)

                                                                                                                                                                                    cTCCCTTGGAAGGAGGGGGTCGGC‑GGTTCAA‑TCTCG‑CCCTCGCCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTg                                                  <  1:1691931/239‑1 (MQ=255)

                                                                                                                                                                                       ccTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTaga                                                                                           <  2:845683/195‑1 (MQ=255)

                                                                                                                                                                                       ccTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTGCTTGTAAAAATGTGAAGTTCCGAAGTGGGTGATTAGCTCAGCTGATAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTg                                                  <  2:82319/236‑1 (MQ=255)

                                                                                                                                                                                              aggaggTGGTGGGC‑GGTCCGA‑TGCCG‑TCATCACCAACAACTTCATCGCCAGCTAAATTTCTTGGAATAATGTGAAGTACCGAAGTGGGTGATGAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACccac                                        <  1:154307/239‑1 (MQ=18)

                                                                                                                                                                                              aggaggGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACccac                                        <  2:1491674/239‑1 (MQ=255)

                                                                                                                                                                                              aggaggGGGGCGGG‑GGTTCGC‑AATCG‑TCATCGCCCACGCCTTTCTCGCCAGCTAAATTTCTTGTAAAAATGAGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACccac                                        <  1:1438801/239‑1 (MQ=18)

                                                                                                                                                                                                 aggGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACccacca                                      <  2:695590/238‑1 (MQ=255)

                                                                                                                                                                                                 aggGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACccacca                                      <  2:557031/238‑1 (MQ=255)

                                                                                                                                                                                                 aggGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTATCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACccacca                                      <  2:1300622/238‑1 (MQ=255)

                                                                                                                                                                                                  gggggTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCGGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGGTTGATTCCGTCAGCACCCCCCACTTCTGGTTGTTAGAGCAGTTGGTAGAGGCGGTGTCTTTTTATCAATTTGTCCCCGGGTCTAATGCTGCCCGGccccaccca                                    >  1:1198224/1‑233 (MQ=21)

                                                                                                                                                                                                    gggTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCa                                                               <  2:488491/210‑1 (MQ=255)

                                                                                                                                                                                                             ggTTCGA‑TCCCG‑TCATTGCCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAg                           <  2:922153/238‑1 (MQ=255)

                                                                                                                                                                                                             ggTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAg                           <  2:658306/238‑1 (MQ=255)

                                                                                                                                                                                                             ggTTCGA‑TCCCG‑CCATCACCCACCCCTTTCTCGACAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAGGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAg                           <  2:1556822/238‑1 (MQ=255)

                                                                                                                                                                                                              ggttatagcccG‑TCATAACCCACCACTTTCTCGCCAGGTAAATTTCGTGTAAAAATGTGAAGTAGCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAg                           <  1:980921/230‑1 (MQ=255)

                                                                                                                                                                                                                        cG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACGGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAgg                <  2:1027570/239‑1 (MQ=255)

                                                                                                                                                                                                                        cG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAgg                <  2:1205066/239‑1 (MQ=255)

                                                                                                                                                                                                                        cG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAgg                <  2:405153/239‑1 (MQ=255)

                                                                                                                                                                                                                            catcgCCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcg              <  2:1072939/238‑1 (MQ=255)

                                                                                                                                                                                                                            catcaCCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcg              <  2:1502813/238‑1 (MQ=255)

                                                                                                                                                                                                                            catcaCCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcg              <  2:1331339/238‑1 (MQ=255)

                                                                                                                                                                                                                            caacacccaccACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcg              <  2:106399/235‑1 (MQ=255)

                                                                                                                                                                                                                                acccaccACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCCCCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACtt                                                                                >  1:1321854/1‑168 (MQ=255)

                                                                                                                                                                                                                                acccaccACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACtt                                                                                <  2:1321854/168‑1 (MQ=255)

                                                                                                                                                                                                                                acccaccACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtttt         >  2:1638767/1‑239 (MQ=255)

                                                                                                                                                                                                                                   caccacTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTa      <  2:1116702/239‑1 (MQ=255)

                                                                                                                                                                                                                                   caccacTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTa      <  2:1317211/239‑1 (MQ=255)

                                                                                                                                                                                                                                    accacTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTat     >  1:1129353/1‑238 (MQ=255)

                                                                                                                                                                                                                                      cacTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCCCTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACCTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTa   >  2:1206572/1‑239 (MQ=255)

                                                                                                                                                                                                                                       acTTGCTCGCCTGCTAAATTTCGTGTAAAAGCGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGAGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTAt  <  1:1453105/239‑1 (MQ=255)

GAACTCTCCCTATAATGCGACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGC‑GGTTCGA‑TCCCG‑TCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTAT  >  NC_000913/2520896‑2521362

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