BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,621,987	C→G	10.3%	A73G (GCG→GGG)	yjjJ →	putative protein kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,621,987	0	C	G	10.3%	371.5 / 3.7	162	A73G (GCG→GGG)	yjjJ	putative protein kinase
Reads supporting (aligned to +/- strand): ref base C (71/69); new base G (13/9); total (84/78)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 164 total aligned reads displayed.

TATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTG‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGG‑GCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGA‑AGATGTGCTCTACGCCCTGACCGTATTTAACGGT  >  NC_000913/4621749‑4622222

tATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTG‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCTTCTGGCGGGTGGACGATA‑CCGGAAAGgc                                                                                                                                                                                                                                                >  2:849736/1‑239 (MQ=255)

   tCACGAATCTATTCATTAATGAGCGAGCTGACTGCTCTTTTACTGCAAGGGCCGCGTTCTG‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGCCGATA‑CCGGAAAGACGCa                                                                                                                                                                                                                                             >  2:1169883/1‑239 (MQ=255)

     aCGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTG‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTTAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGGGGACGATA‑CCGGAAAGGCGCATa                                                                                                                                                                                                                                           >  2:440085/1‑239 (MQ=255)

       gatccttttCATTAATGAGCAGGTGGGCGGATTTTTTGCTGCAAGGGCCGCGTTCTG‑CCCCGGAATTGCGCCGGGGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATaaa                                                                                                                                                                                                                                         <  1:1243732/232‑1 (MQ=255)

            tATTCATTAATGAGCGAGCTGACTGATCTTTTACTGGAAGGGCCGCGTTCTG‑TCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTTACGCCTTGTTGCCAGAGCAGATCGGGTGATTCGCTTTGGTCAAGCACGGGCAAGGCGCGCTGCACTGTTGCGTCCTTATCGCGGAATTGGGCGTATACCCGTATGGCGGGGGGGCCGTA‑CCGGGAGGGCGCAAAAAtatcc                                                                                                                                                                                                                                    >  1:1260426/1‑235 (MQ=255)

            tATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTG‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGTGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATCTGCACTGCTGCGTCCTTATCGCGGAATTGCGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGTAAAGCGCATAAAttag                                                                                                                                                                                                                                     >  2:1151547/1‑236 (MQ=255)

            tATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTG‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGGGGACGATA‑CCGGAAAGGCGCATAAATTCg                                                                                                                                                                                                                                     >  2:895528/1‑238 (MQ=255)

                 aTTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTG‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGGTATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGGCGAAA‑CCGGAAGGGGGCAAAAAATCGCCGaa                                                                                                                                                                                                                                >  1:61235/1‑237 (MQ=255)

                 aTTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTG‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGGGGACGATA‑CCGGAAAGGGGCAAAAATTCGCCGAc                                                                                                                                                                                                                                >  2:631600/1‑238 (MQ=255)

                   tAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTG‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGTCGGATGGACGTGA‑CCGGAAAGGCGCTTAAATTTGCCGAC‑ac                                                                                                                                                                                                                             >  2:836999/1‑237 (MQ=255)

                   tAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTG‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGGGGACGATA‑CCGGAAAGGGGCATAAATTCGCCGActt                                                                                                                                                                                                                              >  2:36065/1‑236 (MQ=255)

                   tAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTG‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCCCGGAAATGAGCGTATTCCCGCCTGGCGGGTGGACGTAA‑CCGGAAGGTGGCATAAAttttccgtcctc                                                                                                                                                                                                                             >  1:775623/1‑229 (MQ=255)

                      tGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTG‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGCAAGGCGCATAAATTCGCCGGC‑ATCCgg                                                                                                                                                                                                                         >  2:1066147/1‑239 (MQ=255)

                        cgcgagcTGACTGATCTTTTACTGCAAGGGCCGCGTTCTG‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGACGCGTCCTTATCGCGGAATTGAGCGATTTCCCGTCTGGCGGGTGGGTGCAA‑CCCGACAGGCACATCAATTGGCCCAC‑ACCCGGtg                                                                                                                                                                                                                       >  1:944850/2‑238 (MQ=255)

                        agcgagcTGACTGATCTTTTACTGCAAGGGCCGCGTTCTG‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGCAAGGCGCATAAATTCGCCGAC‑ATCCGGtt                                                                                                                                                                                                                       >  2:738422/1‑239 (MQ=255)

                        agcgagcTGACTGATCTTTTACTGCAAGGGCCGCGTTCTG‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGGATTCCCCTCTGGCGGCGGGA‑GCTACCCGGCAAGGGGCCTAAATTTGCCgcacaccgggt                                                                                                                                                                                                                        >  1:282426/1‑229 (MQ=255)

                        agcgagcTGACTGATCTTTTACTGCAAGGGCCGCGTTCTG‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGAAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCCTAAATTCGCCGAC‑ATCCGGtt                                                                                                                                                                                                                       >  2:304774/1‑239 (MQ=255)

                            agcTGACTGATCTTTTACTGCAAGGGCCGCGTTCTG‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGAATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTACCCGCC‑ATCCGGTtttgt                                                                                                                                                                                                                   >  2:1030134/1‑238 (MQ=255)

                             gcTGACTGATCTTTTACTGCAAGGGCCGCGTTCTG‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTtgt                                                                                                                                                                                                                     >  2:925018/1‑236 (MQ=255)

                                gactgaTCTTTTACTGCAAGGGCCGCGTTCGG‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGAGCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGc                                                                                                                                                                                                               <  2:64217/239‑1 (MQ=255)

                                 actgaTCTTTTACTGCAAGGGCCGCGTTCTG‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTGGTTAAAGCACGGGCAACGCGATATGCCCTGCTGCGTCCTTATCGCGGAATTGAGCGTCTGCCCGACTGGCGGTGGGCAGCTA‑CCGGTAAGGGGCCGAAATTCGCCggccaccgggttggggggcc                                                                                                                                                                                                               >  1:945821/1‑220 (MQ=255)

                                   tgaTCTTTTACTGCAAGGGCCGCGTTCTGCCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCg                                                                                                                                                                                                             <  2:439498/239‑1 (MQ=255)

                                   tgaTCTTTTACTGCAAGGGCCGCGTTCTG‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCg                                                                                                                                                                                                             <  2:564455/238‑1 (MQ=255)

                                       cTTTTACTGCAAGGGCCGCGTTCTG‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGCGGACGATA‑CCGGAAAGGGGCAAAAATTTCCCGAC‑ACCCGGTAGTGCGGGCCCCCggg                                                                                                                                                                                                        >  1:64217/1‑239 (MQ=255)

                                          ttccgggaaGGGCCGCGTTCGG‑CCCCGGAATTGCGCCAGGGGTTGGCAATCAGGCAGGGGACGGTCTCACGCCTTGTTGCCAGAGAAGAACGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAg                                                                                                                                                                                                     <  1:1122601/232‑1 (MQ=255)

                                          ttACTGCAAGGGCATCGTTCTG‑CCCCGGAATTGTGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTTTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTTAAGCACGGGCAACGAGCTATGCACTGCTGCGTCCTTATCGAGGAATTGAGTGTATTCCGTGGTGGTGGGGGGACGCGA‑CCGCTAAAGGGCACAAATGCGCCGGC‑AACCGGTGGGGCTGGGCACAgggggg                                                                                                                                                                                                     >  1:691560/1‑236 (MQ=255)

                                           taccgaaaGGGCCGCGTTCTG‑CCCCGGAATTGCGCCAGCCTCTGGCAATCAGGCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGt                                                                                                                                                                                                    <  2:422904/233‑1 (MQ=255)

                                           tACTGCAAGGGCCGCGTTCTG‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGCAGATCGGGTGATTCGCTTTGTTAACGCACGGGCAACGCGCTATGCACTGCTGAGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGTCGGGAGGACGATA‑CCGGAAAGGCGAATAAATTCGCCGAC‑ACACGGTTGGGCCGGCCGCAGGGagag                                                                                                                                                                                                    >  1:489080/1‑236 (MQ=255)

                                           tACTGCAAGGGCCGCGTTCTG‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTTTCGCGGAATTGAGCGAATTCCCGTCTGCCGGGGGGACGGTA‑CCGGAAAGGCGCCTAAATTTACCGCA‑ATCCGGTTGGGCTGGCCGCAGGGAAgg                                                                                                                                                                                                    >  1:295571/1‑238 (MQ=255)

                                           tACTGCAAGGGCCGCGTTCTG‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCTCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGt                                                                                                                                                                                                    >  2:66894/1‑239 (MQ=255)

                                               gCAAGGGCCGCGTTCTG‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGCGCGTATTCCCGTCTGGCGGGGGGACGACA‑CCGGAAAGGCGCATAAATTCGCCAGC‑ACTCGGTTGTGCTGGAGGCGGGGAAGTTggt                                                                                                                                                                                                >  1:1154082/1‑237 (MQ=255)

                                                cAAGGGCCGCGTTCTG‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTg                                                                                                                                                                                                  >  2:1122601/1‑236 (MQ=255)

                                                  aGGGCCGCGTTCTG‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTAGGCGTATTCCCGTCTGGCGGGGGGTCGATA‑GCGGGAAGGCGGATAAATTCTCCGGC‑CTCCGGGTGTTGTGGCCGCCGGGAAGTTggattg                                                                                                                                                                                             >  1:1139158/1‑234 (MQ=255)

                                                      ccGCGTTCTG‑CCCCGGAATTGCGCCAGCGTTTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGATTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTaa                                                                                                                                                                                          <  2:259746/238‑1 (MQ=255)

                                                       cgcgTTCTG‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAACGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGGaaaa                                                                                                                                                                                        >  2:368996/1‑237 (MQ=255)

                                                               g‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCCTAAATTCGCCCGC‑AACCGGTTGTGCTGGCCCCAGGGAAGTTGTCTGGTAACAGGGGCCGa                                                                                                                                                                                >  2:53910/1‑239 (MQ=255)

                                                               g‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGa                                                                                                                                                                                <  2:1251050/239‑1 (MQ=255)

                                                                     ggAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGcg                                                                                                                                                                           <  2:61235/239‑1 (MQ=255)

                                                                      gAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAGGTTGCCTGGTAACAGGGGCAGTTGGcg                                                                                                                                                                           >  2:1231979/1‑238 (MQ=255)

                                                                             gcCGGCGTCTGGCATTCAGTCAGGGGGCGTTCCCCCCGCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATAGGCACTGGTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAg                                                                                                                                                                   <  1:1086347/239‑1 (MQ=255)

                                                                             gcCAGGGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACa                                                                                                                                                                    <  2:665259/238‑1 (MQ=255)

                                                                              ccAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAg                                                                                                                                                                   <  2:282426/238‑1 (MQ=255)

                                                                                 gCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCTACGCGATATGCACTGCTGCGTCCTTATAGCGGAACTGCGCGCATTCCCGTCTGGCGGGGGGTCGATA‑CCGGAAGGGGCCATCAATTCGCCGCC‑ATCCGGTTGTGCTTGCCTCCGGGAAGTTTTCGGGTAACAGGTCCAGTTGGcgcaggacaaaggg                                                                                                                                                               >  1:1078919/1‑227 (MQ=255)

                                                                                 gCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGGGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGCGGt                                                                                                                                                               >  2:266379/1‑239 (MQ=255)

                                                                                     cTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTg                                                                                                                                                            <  2:1059867/238‑1 (MQ=255)

                                                                                          aaTCAGGCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGt                                                                                                                                                       <  2:1078919/238‑1 (MQ=255)

                                                                                           aTCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGGC‑ATCCGGTTGGTCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCGGTTGGCGCCGACAAGTGGTTTGATGGttt                                                                                                                                                     >  2:598560/1‑239 (MQ=255)

                                                                                               ttcagtcctcGTTCTCCCGCCTTGTTCCCCGGTAAGACCGGGGTCTTCGCTTTGGGAAAGCCCGGGCAACGCGATATGCCCTGCTGCGTCCTTACGGCGGAATTGAGCGTTTTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGccc                                                                                                                                                 <  1:895528/230‑1 (MQ=255)

                                                                                               gtcagGCGACGTTCTCACGCCTTGGTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGccc                                                                                                                                                 <  2:148872/239‑1 (MQ=255)

                                                                                                 cagGCGACGTTCTCCAGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGGAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTAGCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATC‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGccc                                                                                                                                                 <  2:405674/237‑1 (MQ=255)

                                                                                                          tttttccgcttttttgggatagaagaTCGGTTGATTCTTTTTGGTAAAGCACGGGCAACGCCATATGCACTGCTGCGTCCTTATGGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGAc                                                                                                                                      <  1:517228/219‑1 (MQ=255)

                                                                                                           tcccACGCCTTGTTGCCAGGGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑ACGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGAc                                                                                                                                      <  2:775623/235‑1 (MQ=255)

                                                                                                                     tgttgCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGAAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTACGa                                                                                                                            >  2:1014369/1‑238 (MQ=255)

                                                                                                                       ttgCCAGAGAAGAGCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCg                                                                                                                         <  2:295571/239‑1 (MQ=255)

                                                                                                                         tccgaagatggttgggttgttCGCTTTGGTAACGGCCGGGCACCGCGATATGCTCTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCa                                                                                                                       <  1:1030898/220‑1 (MQ=255)

                                                                                                                         tccaggaaatgtggggTGGTTTGCTTTGGTAACGCCCGGGCAACGCCATTTCCACTGCTGGGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCa                                                                                                                       <  1:452904/226‑1 (MQ=255)

                                                                                                                                     cGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCCGCAGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGAACTCCGACCGCAGGGCTTTTTAgg                                                                                                           >  2:197735/1‑239 (MQ=255)

                                                                                                                                       gtttgttCGCTTTTGTAAAGTCCGGGCAAGTGCATTTGCACTGGTGCGTTCTTATCGCGGTATTGGGCGTATTCCCGTGTGGCGGGTGGACGATA‑CCGGAAAGTGTCATAAATTCGCCGAC‑AACCCGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTTGTATTTGACCGATCTCCTACCGCAGGGCTTTTTAGG‑gc                                                                                                        <  1:197735/234‑1 (MQ=255)

                                                                                                                                        gTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGAATTTAACCGATCTCCGACCCCCGGGCCTTTTAGG‑gc                                                                                                        >  2:1015420/1‑238 (MQ=255)

                                                                                                                                          gATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGCC‑ATCCGGTTGTGCTGGCAGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGCTTGATGGTTTGCCCTGGGATTTGACCGATCTCCGACCGCAGGGCTTTTTAggcatg                                                                                                       >  2:1111099/1‑234 (MQ=255)

                                                                                                                                          gATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGG‑acggg                                                                                                     >  2:632052/1‑237 (MQ=255)

                                                                                                                                          gATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTAACCGATCTCCGACCGCAGGGCTTTTTAGG‑gcgcg                                                                                                     >  2:246921/1‑239 (MQ=255)

                                                                                                                                                 tttGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGCC‑ATCCGGTGGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGAGCAGATGGCGAACACCAGGGGGTTGGTGGGTTGGCCTGGGAGTTTACCCTCGTCCGACCGCAGGGGTTTTTTgg‑gggggggggggg                                                                                              >  1:711106/1‑228 (MQ=255)

                                                                                                                                                  ttGGAAAAGCACGGGCAACGGGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGGCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCTCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGG‑GCGCGCGTGGGGc                                                                                             <  2:229326/239‑1 (MQ=255)

                                                                                                                                                            ccGGGCAGCGCGATTGGCCCTGCTGCGTCCTTACCGGGGAATTGAGCTTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGG‑GCGCGCGTGGGGCAGGAAGTTAg                                                                                   <  1:246921/238‑1 (MQ=255)

                                                                                                                                                                          atGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGG‑GCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGaa                                                                      <  2:1283802/238‑1 (MQ=255)

                                                                                                                                                                              aCTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGGATTTGACCGGTTCCCGACCCACGGGCTTTTTCGG‑GCCCGGGTGGGGCGGGAAGGTAGCcgacaaacggaatttg                                                                  >  1:525553/1‑224 (MQ=255)

                                                                                                                                                                               ctggtgGGTCCGTAACGCGGAATCTGGCGTGTTCCCGGCTGGGGGGTGGCCGATA‑CCGAAAAGGCGACTAGATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTGTTCTGGTAACAGGCGCGGATGGCGACGAACAGTGGTTTGAGGGGTTGCACTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGG‑GCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATctgac                                                                <  1:368996/239‑1 (MQ=255)

                                                                                                                                                                                tggtgcGTCATTTTGGGGGAATTGAGCGTTTCCCCGTGTGGCGGGTGGACGATA‑CCGGGAAGGCGCATAAATGCGCCGAC‑ATCCGGTTTGGCTGGGCGCAGGGAAGTTGTCTGGGAACAGGCGGAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGTGCTTTTTAGG‑GCGCGCGTGGGGAAGGAAGTTAGCCGCGCAACCGAATctgac                                                                <  1:1169883/235‑1 (MQ=255)

                                                                                                                                                                                tgctgcGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGCAAGGCGCATAAATTCGACGAC‑ATACGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACCGGCGCAGGTGGCGACGAACAGTGGGTGGATGGTTTGCCCTGGTATTTGACCGATCTCCGCCCGCCGGGACTTTTTGG‑GCGCGCGGGGGGCAGGAAGTTGCACGCGAAACTGaaatgattt                                                               >  1:546370/1‑232 (MQ=255)

                                                                                                                                                                                tgctgcGTCCTTATCGCGGAATTGAGCGTATTCCCGTATGGCGGGTGGACGATA‑CCGGCAAGGCGCATAAATTCGCCCAC‑ATCCGGTTGTGCTGGCCGCAGGGGAGTTGTCTGGAAACAGGCTCAGAGATGTCCGAACAGGGGGTTGATGGGTTGACCTGGTTTTAGACCGATCTCCGACCCCAGGGCTTTTTAGGCGGCCGCGT‑GGGCTGGAAGtgtgcggggaaagcgaaacgagca                                                               >  1:993853/1‑216 (MQ=255)

                                                                                                                                                                                tgcggggttcttgtgggggAATTGAGCGTTTTCCCGGTGGGGGGGGGGACGATA‑CCGAAAAGGCTCATAATTTCGCCGAC‑AGCCGGTTGTGCTGGCCGCGGGTAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGG‑GCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATctgac                                                                <  1:560752/221‑1 (MQ=255)

                                                                                                                                                                                  ctgcGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGG‑GCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATctgactga                                                             >  2:317958/1‑239 (MQ=255)

                                                                                                                                                                                   gtcGTCCTTATCGGGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGG‑GCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATctgactga                                                             <  2:1204983/236‑1 (MQ=255)

                                                                                                                                                                                       tctttttcccgGAATTGGGGGTATGCCCCTTTGGCGGGGTGGCGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGGGCTGGCCGCAGTGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGG‑GCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGata                                                        <  1:393139/230‑1 (MQ=255)

                                                                                                                                                                                       tCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGG‑GCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATATGACTGATGata                                                        >  2:69/1‑239 (MQ=255)

                                                                                                                                                                                       tCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGAATTTGACCGATCTCCGACCGCAGGGCTTTTTAGG‑GCGCGCGTGGGGCAGGAAGTTAGCCGCGCCACTGAATCTGACTTGTGata                                                        >  2:1209848/1‑239 (MQ=255)

                                                                                                                                                                                       tCCTTAGCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGG‑GCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGata                                                        <  2:288154/239‑1 (MQ=255)

                                                                                                                                                                                            aTCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGGACAGTGGGTTGAGGGGTTGCCCTGGGATTTGACCGATGTCCTCCCGCCGGGCTTTTTATG‑TGGCGCGTGGGGGCGACAGTTACCCGCGCACCTGATTCTGACGTCAGATATACgg                                                   >  1:1204983/1‑238 (MQ=255)

                                                                                                                                                                                            aTCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACATGGGTTTGATGGTTTGCCCTGGGATTTGACCGATCTCCGACCGCAGGGCTTTTTAGG‑GCCGCCGTGGGGCAGGAAGTTAGCAGCGCACCTGAATCTTACtgatgtattccc                                                    >  1:297291/1‑231 (MQ=255)

                                                                                                                                                                                            aTCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGG‑GCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACg                                                    >  2:206144/1‑238 (MQ=255)

                                                                                                                                                                                              cgcgGAATTGGGCGTATTCCCGGCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGG‑GCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGtct                                                 <  2:944850/239‑1 (MQ=255)

                                                                                                                                                                                                         gCGTATTCCCGTCGGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGG‑GCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATctgactga                                                             <  2:957990/216‑1 (MQ=255)

                                                                                                                                                                                                             tttccggggtGGCGGGGGGCCGATC‑CGGTAAAGGCGCATAAAGTGGGCGGC‑CGCCGGTTGTGCCGGCCGGAGGGAAGTTGTCGGGTAACAGGTGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGG‑GCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGaaga‑ag                                  <  1:849736/229‑1 (MQ=255)

                                                                                                                                                                                                             aTTCCCGGCTGGCGGGTGGGCGAAC‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGG‑GCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGaaga‑ag                                  <  2:825928/238‑1 (MQ=255)

                                                                                                                                                                                                                cctgttgggggggTGGACGATA‑CCGGAAAGGCGCATAAATGCGCCGAC‑AGCCGGGTGGTCTGGCGGGAGGGAGGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGGTCTCCGACCGCAGGGCTTTTTAGG‑GCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGA‑AGAtgt                              <  1:26241/229‑1 (MQ=255)

                                                                                                                                                                                                                    tCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGGCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGG‑GCGCGCGTGGGGCAGGCAGTTAGCCGCGCAACTGAATCTGACTGATCAAATACGTCTCTGGCCGgcaga‑aaaatggct                           >  2:1163686/1‑230 (MQ=255)

                                                                                                                                                                                                                      tGGGGGGTGGACGATC‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGG‑GCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGA‑AGATGTGCTCTa                        <  2:546370/239‑1 (MQ=255)

                                                                                                                                                                                                                              tttcGATG‑CCGGAAAGTGGCATTATGTCCCCGGC‑ATCTGGTTGTGGTGGCCGGAGGGAAGTTGTCTGGGAACAGGCGCAGATGTCGAAGAACAGTTGTTGGATGTTTCGCCCTGGTATTGTACCGATCTCCGACCGCAGGGCTTTTTAGG‑GCGCGCGTGGGGCAGGAAGTTAGCCGCGGAACTGAATCTGACTGATGCTATGCGTCTCTGGCAGGAAGA‑AGATGTGCTCTACGCCCTGa                <  1:32271/236‑1 (MQ=255)

                                                                                                                                                                                                                               ggcttggcctggAAGGGCGTATAAATGCGCCCAC‑TTGCGGTTGTGCTGGCCGCAGGGAGTTTGTTTGGTAACAGGCGCATATGGCCACGAACAGTGGTTGGATGGTTTGCCCTGGTATTTGTCCGATCGCCGTCCGCAGGGCTTTTTAGG‑GCGGGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGA‑AGATGTGCTCTACGCCCTGa                <  1:819734/229‑1 (MQ=255)

                                                                                                                                                                                                                               gACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGG‑GCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGA‑AGATGTGCTCTACGCCCTGa                <  2:705556/238‑1 (MQ=255)

                                                                                                                                                                                                                               gACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGGGCTGGCCGGAGGGAAGTTGGCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGG‑GCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGA‑AGATGTGCTCTACGCCCTGa                <  2:1253659/238‑1 (MQ=255)

                                                                                                                                                                                                                                gataagacgggAAAGGGACATTAATTCGCCGGC‑ATCCGGTTGGGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGGAGATGGCGACGAACAGTGGGTTGATGGTTTGCCGTGGTATTTGACGGATCTCCGACCGCAGGGCTTTTGAGG‑GCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGA‑AGATGTGCTCTGCGCCCTGa                <  1:181068/229‑1 (MQ=255)

                                                                                                                                                                                                                                  gATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCCCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTGGACCGATCTCCGACCGCAGGGCTTTTTAGG‑GCGCTCGTGGGGCAGGACGTTAGCCGCCCAACCGAATCTGACTGATCATTTCCGTCGTTGGCAGGAAGACAAGTGGGC‑CTACACCCTGACCGt            >  1:690625/1‑239 (MQ=255)

                                                                                                                                                                                                                                         ggAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGG‑GCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGA‑AGATGTGCTCTACGCCCTGACCGTATTTaa      <  2:167449/239‑1 (MQ=255)

                                                                                                                                                                                                                                           aaaGGGGCATAAATTCGCCGGC‑ATCCGGTGGGGTTGGCCGCAGGGAAGTTGTCTGGGAACAGGCGCAGATGGCGACGAACAGTGGTTTGAGGGTTGGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGG‑GCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGA‑AGATGTGCTCTACGCCCTGACCGTATTTAAc     <  1:625307/238‑1 (MQ=255)

                                                                                                                                                                                                                                           aaaGGCGCATAAACTCGCGGAC‑ACCCGGTTGTGTGGGGCGCAGGGAAGTTGTCTGGTAACAGGCGGAGAGGGGTGCCAACAGGGGTTTGGTGGTTTGCCCTTGTATTTGACCGATCCCCGACCGCAGGGCTTTTTAGG‑GCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGA‑AGATGTGCTCTACGCCCTGACCGTATTTAAc     <  1:893195/238‑1 (MQ=255)

                                                                                                                                                                                                                                           aaaGGCGAAGAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGG‑GCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGA‑AGATGTGCTCTACGCCCTGACCGTATTTAAc     <  2:690625/238‑1 (MQ=255)

                                                                                                                                                                                                                                             aGGCGCATAATGTCG‑CGACAATCCGGTTGTGGGGGCCGGAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGCCGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGG‑GCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGA‑AGATGTGCTCTACGCCCTGACCGTATTTAACGGt  <  1:1000910/239‑1 (MQ=255)

                                                                                                                                                                                                                                             aGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGG‑GCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGA‑AGATGTGCTCTACGCCCTGACCGTATTTAACGGt  <  2:525553/239‑1 (MQ=255)

TATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTG‑CCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATA‑CCGGAAAGGCGCATAAATTCGCCGAC‑ATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGG‑GCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGA‑AGATGTGCTCTACGCCCTGACCGTATTTAACGGT  >  NC_000913/4621749‑4622222

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