BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,065,603	T→A	15.3%	I55F (ATC→TTC)	yihO ←	putative sulphoquinovose importer

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,065,603	0	T	A	15.3%	196.4 / 3.7	96	I55F (ATC→TTC)	yihO	putative sulphoquinovose importer
Reads supporting (aligned to +/- strand): ref base T (57/20); new base A (6/11); total (65/31)
Fisher's exact test for biased strand distribution p-value = 3.75e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTATTACGCCGTGCGTTTACGGGCTTCT  >  NC_000913/4065367‑4065832

catcaGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑a                                                                                                                                                                                                                                      >  2:979634/1‑238 (MQ=255)

catcaGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑a                                                                                                                                                                                                                                      >  2:978606/1‑238 (MQ=255)

catcaGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTAAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑a                                                                                                                                                                                                                                      >  2:616001/1‑238 (MQ=255)

   caGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATc                                                                                                                                                                                                                                   >  2:518374/1‑238 (MQ=255)

   caGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGGTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCCTATAAAATAAAGGGACGAAATTTGCCTTTAGGCACAATATTTT‑T‑ACGCGAGTCGAGTAAAAAA‑CAGGTG‑AGCATACCGGTAAACG‑CGGTAAAAAAATGCGACACCAGAAAGa‑a‑aacc                                                                                                                                                                                                                                   >  1:361018/1‑236 (MQ=255)

        tGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCa                                                                                                                                                                                                                             >  2:148813/1‑239 (MQ=255)

        tGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAAATTCGCCACCAGCAAGA‑T‑AATCCCGCCa                                                                                                                                                                                                                             >  2:710119/1‑239 (MQ=255)

        tGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAACATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGGG‑AGCATATCGGTAAACG‑CGGGAAAAAACTTGTCCACACGAAAGA‑A‑AATCCCGCCa                                                                                                                                                                                                                             >  1:225152/1‑239 (MQ=255)

        tGTATGAAAGTCCAAACATCATAAACAGCGCCGTGGAAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACAAAAATTGCCTTTAGGCCCAATATTTT‑T‑TACCAAGTCGAGTAAAAAT‑CAGGCG‑AGCATACCGGTAAAAT‑CGGCAAAAAACTTCGCCACCAGAAAGA‑A‑AATCCCGCCa                                                                                                                                                                                                                             >  1:626595/1‑239 (MQ=255)

          tATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCAACCAGAAAGA‑T‑AATCCCGCCATa                                                                                                                                                                                                                           >  2:1018001/1‑239 (MQ=255)

           aTGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGGGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCGGAAAGA‑T‑AAACCCGCCATAg                                                                                                                                                                                                                          >  2:941898/1‑239 (MQ=255)

                aGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑CTAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAAATTCGCAACCAGAAAGA‑T‑AATCCCGCCTTAGTagcc                                                                                                                                                                                                                     >  2:112242/1‑236 (MQ=255)

                  tCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCa                                                                                                                                                                                                                    >  2:452165/1‑238 (MQ=255)

                  tCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAg                                                                                                                                                                                                                   >  2:55382/1‑239 (MQ=255)

                  tCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑ACGGTG‑AGCATACCGTTAAACC‑CGGTAAAAAACTTCACCACCAAAAAGA‑A‑AAACCCCCCaaaaaaaaaag                                                                                                                                                                                                                   >  1:837019/1‑229 (MQ=255)

                    cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGATTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAAAG‑CGGGAAAAAACTTCGCAACCAGAAAGA‑A‑AATCCCGCCAAAGaaaaaaagg                                                                                                                                                                                                                 >  1:846173/1‑233 (MQ=255)

                    cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAGAT‑CCGGTG‑AGCATATCGGAAAACG‑CGGTAAAAAAATTCGCCACCAGAAAGA‑A‑AATCCCGCCATAGGAAACAGGc                                                                                                                                                                                                                 >  2:598395/1‑239 (MQ=255)

                    cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGc                                                                                                                                                                                                                 >  2:922083/1‑239 (MQ=255)

                    cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGAAAACG‑CGGAAAAAAACTTTGCCACCAGAAGGA‑T‑AATCCCGCCATAGTAAGCAGGc                                                                                                                                                                                                                 >  2:721744/1‑239 (MQ=255)

                    cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTA‑AGCATATCGGTAAACG‑CGGTAAAAACCTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGc                                                                                                                                                                                                                 >  2:9148/1‑239 (MQ=255)

                    cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGGG‑AGAATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCCGCAAGA‑T‑AATCCCGCCATCGAAAGCAgg                                                                                                                                                                                                                  >  2:319770/1‑238 (MQ=255)

                    cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCAGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCAAATCGGCAAACG‑CGGGAAAAAACTACGCCACCCAAAAGA‑G‑AATccaacaaaaaaaagaagga                                                                                                                                                                                                                 >  1:1086048/1‑222 (MQ=255)

                            atAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCTCTGCATATAAAATAAAGGGACGAAATTTGACTTTAGGCCCCATATTTT‑T‑CCGAGAGACGAGTAAAAAT‑CCGTTG‑TGCCAATAGGTAACCC‑TGGTAAAAACCTTCGCAACCAGACAAA‑A‑AACCCCCCCAAATAAAACAGGCATACCaaa                                                                                                                                                                                                         >  1:1204148/1‑239 (MQ=255)

                            atAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑GGGTAAAAAAATTCGCCACACGAAGGA‑T‑AATCCCGCCATAGTAACCAGGCATCCCt                                                                                                                                                                                                           >  2:408583/1‑237 (MQ=255)

                            atAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTaa                                                                                                                                                                                                         >  2:667635/1‑239 (MQ=255)

                              aaaCAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGAAAAAAACTTTGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCAAAtc                                                                                                                                                                                                       >  2:291777/1‑238 (MQ=255)

                                  ggccgggtcactaggtttgtTTTCACCGACAAAAAAAAGGGGGTCGAA‑AAAAACGGGGGCTTTCAAATATCCCCCGCGGAAA‑CCGCTTAATAAAAAAAAAAGGGAAAAATGTTGCCTTTAGCCCAAATATTTC‑T‑CCGCGATGCGATTAAAACA‑CCGGTG‑AGCATTCCGGTAAACG‑AGGTAAAAACCTCCGCCACAAGAAAAA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATc                                                                                                                                                                                                   <  1:452165/221‑1 (MQ=255)

                                   gcgcGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCCGAAAGA‑T‑AATCCCGCCATATTAAGCAGGCATCACTAATTCATCg                                                                                                                                                                                                  >  2:865917/1‑239 (MQ=255)

                                   gcgcGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGAAACCAGGCATCCCTAATTcctg                                                                                                                                                                                                   >  2:311397/1‑235 (MQ=255)

                                    cgcgGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCgg                                                                                                                                                                                                 >  2:261733/1‑239 (MQ=255)

                                    cgcgGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCgg                                                                                                                                                                                                 >  2:811831/1‑239 (MQ=255)

                                     gcgGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAAtt                                                                                                                                                                                                       >  2:904744/1‑232 (MQ=255)

                                       ggTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGtat                                                                                                                                                                                              >  2:585503/1‑239 (MQ=255)

                                         tggtcagttgttgtttAACCCGAAAAAAAGAAGGGGTGGAGAAAAAACTGGAGCAGTACAATGAACGCTGAGGGAC‑CCGTTGAATAAAAAAAAAAGGGACAAATTTCGCTTTTAGGCCTAATTTTTT‑C‑CCGCGAGTCGAGAAAAAAC‑CGGGAG‑ATCATCTCGGAAAACG‑CGGAAAAAATCTCCGCCACCAGAAAAA‑A‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGtat                                                                                                                                                                                              <  1:746146/225‑1 (MQ=255)

                                                 gtcgtTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCTCGAGTAGGGTAAAAAT‑CAGGTG‑CGCATATCGGTAAACG‑GGGTAAAAACCTTCGCAACCAGAAAGA‑G‑AATCCCGCCATATGAAGCAGCAAACCTAAATTAATAGGTAAAAAATTTgag                                                                                                                                                                                    >  1:261854/1‑239 (MQ=255)

                                                   cgtTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagga                                                                                                                                                                                  >  2:277425/1‑239 (MQ=255)

                                                             ggCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAACGAGGCACCCCTAATTCATAGGGAAAAAAATTtagaaaaaaaaaaga                                                                                                                                                                        >  1:713553/1‑227 (MQ=255)

                                                                      aaGGTGGGGGAG‑AAAAACGGGAGTGTTGCAATTAACCCCCCCGGAA‑CCGCTGAATATAAAATAAAGGGACAAAATTTGCCTTTAGCCCCAATATTTT‑C‑CCGCGAGTCGAGTAAAAAC‑CCGGTG‑AGCATATCGGAAAACG‑CGTTAAAAACCTTCGCCACCAGAAAGA‑T‑AACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCa                                                                                                                                                               <  2:1086048/239‑1 (MQ=255)

                                                                       tggtggtGGAT‑ATAAAGTGGGGTTTTTCATTTCACCCCGCCGGAC‑CCTGTGAAAAAAAGAAAAAGGGACAAAATTTGCTTTTGGCCCAAATATTTT‑T‑CCGCTGGTGGAGAAAAAAC‑CGGGTG‑AGCATATCGATAAACG‑CGTAAAAAACCCCCGCCACAAAAAAAA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCag                                                                                                                                                              <  1:667635/238‑1 (MQ=255)

                                                                         gtggtgGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAATG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGAATCCCTAATTCATCGGCAAAAAAGTTGAGGAGAAACAAAGTACCGATGCaaag                                                                                                                                                            >  1:888658/1‑239 (MQ=255)

                                                                                     aaaCTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGATTCGAGTAAAAAT‑CCGGTG‑AGCTTATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGAATAAAATTTGAGTAGATACAGAGAACCGATACacagactacaattca                                                                                                                                                 >  1:760861/1‑228 (MQ=255)

                                                                                           gcggttgacAATTACCGCCGGCGGAC‑CCGTTG‑AAAAAAAATAAAAGGGCCAAATTTGGCCTTAGGGCCAATTTTTTCT‑CGGCGAGCTGTGAAAAAAC‑CCGGTG‑ACAATTGCGGTAAACG‑CGGAAAAAACCTGCACCACAGGAAAAA‑T‑AACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCg                                                                                                                                            <  1:922083/230‑1 (MQ=255)

                                                                                                ttGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTc                                                                                                                                       >  2:521077/1‑238 (MQ=255)

                                                                                                ttGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTAATCGGTATAAAATTTGAGGAGTTACAAAGAACCGCTGCAGAGCATTAAATTCGACCCGAcacc                                                                                                                                       >  2:289208/1‑235 (MQ=255)

                                                                                                ttGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCACAAAGA‑T‑AATCCCGCCATAGTAAGCAGGAATCCCTAATTCATCGGGATACAATTTGAGGAGATAAAGAGTAACGATACAGAGCATTAAATTCGAAcagcaggc                                                                                                                                       >  1:555281/1‑231 (MQ=255)

                                                                                                 tcctttttaccccgccgGAC‑CCTCTGAAATAAAAAAAAAAGGAAAAAATCTGCCTTGAGGCACAATTTGTT‑T‑CCGCGAGTAGAGAAAAAAC‑CCGGAG‑GGCATACCGGTAAACG‑AGGTAAAAACCCGCCCCACAAAAAAGA‑A‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGc                                                                                                                                     <  1:710119/228‑1 (MQ=255)

                                                                                                   caaatttacccgcccGAACCCGCTGAAAAAAAAAAAAAGTGACGACATGTGCCTTTAGGCCAAATTATTG‑C‑CCGCGAGTACAATAAAAAG‑CCGGTG‑AGCTTACCGGTAAACG‑CGGTAAAAAACTTCGCCACAAGAAAAA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGcc                                                                                                                                    <  1:55382/230‑1 (MQ=255)

                                                                                                     aTTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCcca                                                                                                                                  >  2:834494/1‑238 (MQ=255)

                                                                                                        aaCGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGGG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCCAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGACTTAAATTCGAACCGACGGCGCccataa                                                                                                                               >  2:650847/1‑237 (MQ=255)

                                                                                                                              tataAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTc                                                                                                               >  2:411483/1‑233 (MQ=255)

                                                                                                                                 aaaaaaagGGGAAAAATTTTGCTTTTAGGCCAAATATTTT‑C‑CCGCGAGTGGAGTAAAACT‑CCGGTG‑AGCATATGGGAAAACG‑CGTAAAAAACCTTCGCCCCAAAAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGAtt                                                                                                      <  2:837019/239‑1 (MQ=255)

                                                                                                                                        aGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTGAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTACTTCATCGGCATAAAATTTGAGGAGATACAGGGTACCGCCGCAAACCATTAAATTCGAACCAGCGCCCCCCCACCAAAAAGCGATTTTTTccccacaaataaggtta                                                                                               >  1:487036/1‑225 (MQ=255)

                                                                                                                                        aGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGACACGACGCCGCCCATACCTAAGGGGATTTTTTccccaaggttacgtttt                                                                                               >  1:1019316/1‑225 (MQ=255)

                                                                                                                                        aGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGtttt                                                                                               >  2:470634/1‑239 (MQ=255)

                                                                                                                                        aGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGCACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGtttt                                                                                               >  2:131707/1‑239 (MQ=255)

                                                                                                                                        aGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGAATAAAATTTGAGTAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGACGATTCAGgttt                                                                                               >  2:722927/1‑239 (MQ=255)

                                                                                                                                          ggccGAAATTTCCCTTGAGGCCCAATCTTTT‑C‑CCGGGAATCAAATAAAAAC‑CGGGTG‑AGCATACGGGAAACCG‑CGGTAAAAACCTTCACCACCAAAAAGA‑A‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGtttt                                                                                               <  2:831254/234‑1 (MQ=255)

                                                                                                                                             cGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTCTAAAATTTGAGGAGATACAGAGTACCGATCCAGAGCATTAAATTCGCACCGACGTCGCCCATACCATAGGCGATTTTTTCCCCCAGCTTCAGTTTTAATGt                                                                                          >  2:513184/1‑239 (MQ=255)

                                                                                                                                             cGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGt                                                                                          >  2:363511/1‑239 (MQ=255)

                                                                                                                                             cGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACAGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCCCCGATTCAGTTTTAATGt                                                                                          >  2:1262948/1‑239 (MQ=255)

                                                                                                                                             aaaaaaTTGTCTTTGGGGCCAATTTTTT‑TCCCGGGAGGCGAAAAAAACC‑CCGGGG‑AGCATAGCGATAAACG‑CGGAAAAAAACTTCGCCACAAAAAAGA‑A‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAAt                                                                                            <  2:853689/236‑1 (MQ=255)

                                                                                                                                                     gCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAACAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATCAAATTTATGGAGATAACAAGTACCGATGCAGAGCATTAAATTGGAGCAGAGGTAGCCCAAACCACAGAGGATTTTTTTCCGCGGAATCATTTTTTATGTTAGTggca                                                                                  >  1:43970/1‑237 (MQ=255)

                                                                                                                                                     gCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAACGTCACGGGAt                                                                                  >  2:524240/1‑239 (MQ=255)

                                                                                                                                                          tAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCAACGCCGCCCATACCTAAGGCGGTTTTTTCCCCCAAATTAAATTTTAATGTCAGTGGATTAAgg                                                                             >  1:519614/1‑239 (MQ=255)

                                                                                                                                                           gtgcctattttttt‑t‑CCGAGCGTCGATAAAACAC‑CCGGTG‑GCAATAT‑GGAAAAAG‑CAGTAAAAAAATTTCCCCCAAAAGAAA‑TAAATCCCGCCATAGTAAGCAGTCATCCCTAATTCATCGGTATAAAATGTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGt                                                                            <  1:979634/229‑1 (MQ=255)

                                                                                                                                                           gggCCCAATTTTTT‑C‑CCGGGAGGCGGAAAAAAAC‑CCGGTG‑AGCATAGGGGTAAACG‑GGGAAAAAATCTGCCCCACAAAAAAGA‑A‑AACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGt                                                                            <  2:846364/238‑1 (MQ=255)

                                                                                                                                                           agcccaaaTTTTTT‑C‑CCGCGATTCGAAAAAAACC‑CGGGTA‑AGAAAATCGGTAAACG‑CGGTAAAAATCTCCGCCACCAGAAAAA‑T‑ACTCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGt                                                                            <  2:626595/233‑1 (MQ=255)

                                                                                                                                                            ggCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGCCAGTGGATTATGGTc                                                                           >  2:531357/1‑239 (MQ=255)

                                                                                                                                                            gaccaaaTTTTTT‑T‑CCGCGGGGCGAGTAAAAAT‑CGGGTA‑AGCATAGGGATAAGCG‑GGGAAAAAAATGGCAACACAAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGt                                                                            <  2:589438/233‑1 (MQ=255)

                                                                                                                                                            cagcacatttttcgt‑CCGAGAGTTTAATAAAGAG‑CGGGGG‑AGTATAGCGGTAAACG‑CGAAAAAAACCTTCCCCCACAAAAAGA‑A‑AATCCCGCCATAGTAAGCAGGCATTCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGt                                                                            <  1:148813/225‑1 (MQ=255)

                                                                                                                                                               cacatttttt‑c‑CCGGGAGTAGAATAAAAAC‑CCGGGG‑AGCATCGCGG‑AAACGTAGATAAAAAACTTCCCAACCAAAAAAA‑A‑AACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGa                                                                        <  2:1248463/233‑1 (MQ=255)

                                                                                                                                                                  catatttttgcccgaaggacagtAAAAATCCCGGTC‑AACATCCGGGTAACCG‑GGGAAAAAACCTTCCCCCACAAAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGa                                                                        <  2:555281/219‑1 (MQ=255)

                                                                                                                                                                  atatTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGAATCCCTAATTCATCGGAATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTAGCCCATACAAAAGGCGATTTTTTCCCCCAAATTAAATTTTAATGTCAGTGGATTAAGGtaaaaaaa                                                                     >  1:622078/1‑232 (MQ=255)

                                                                                                                                                                       gt‑t‑CCGCGAGACGGATAAAACT‑CCGGTG‑AGAATTTCGGAAAACG‑GGGAAAAAAACTCCCCCCCCAAAAAGA‑T‑AACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCa                                                                <  2:1019316/238‑1 (MQ=255)

                                                                                                                                                                        t‑t‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGCCGTTGCCCATACCATAGGGGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGAGGATTTGAGTCAGACATAATCa                                                                >  2:581406/1‑238 (MQ=255)

                                                                                                                                                                         ttcccGGGAGGCGATTAAAAAC‑CCGGTG‑AAAATTGCGGAAAACG‑CGTAAAAAACCTTCGCCACCAAAAAGA‑A‑AACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCa                                                                <  2:888658/236‑1 (MQ=255)

                                                                                                                                                                               cgagtcgagtAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCtt                                                          >  2:808970/1‑239 (MQ=255)

                                                                                                                                                                               cgagtcgagtAAAAAT‑CCGGAG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGAAGAGCATTAAATTCGAACCGACGTCACCCATACCATAGGCGATTTTTTCCACCCGGTTCCGTTTTACATTAAGGGGATTAATGTCAGACATCATCACTCCtt                                                          >  1:783654/1‑239 (MQ=255)

                                                                                                                                                                                      agtAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACg                                                    >  2:817149/1‑238 (MQ=255)

                                                                                                                                                                                      aatAAAAAC‑CCGGTG‑AGTATAGCGGAAAACG‑GGAAAAAAACCTTCCCCACAAAAAATA‑T‑AACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACg                                                    <  2:760861/236‑1 (MQ=255)

                                                                                                                                                                                       attaaaaT‑CCGGGGTAGAATGTCGGAAACCG‑CGGAAAAAACCTTCCCCACAAAAAAGA‑T‑AACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACg                                                    <  2:846173/235‑1 (MQ=255)

                                                                                                                                                                                                        gCATAGCGGTAAGCG‑CGTAAAAAACCTTCCCCAACAAAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGcgtcgt                                    <  2:951204/238‑1 (MQ=255)

                                                                                                                                                                                                         cATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGGGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGttt                                  >  2:717117/1‑239 (MQ=255)

                                                                                                                                                                                                           tctgggaaaaCG‑GGGTAAAAACCTTCCCCACCAGAAAAA‑A‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTttat                                <  2:198133/232‑1 (MQ=255)

                                                                                                                                                                                                           tatCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATACTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCACGGGATTATGGTACGACATAATCCCTCCTTATAACGCCCCGCAGGGCGTCGttttt                                >  1:109180/1‑237 (MQ=255)

                                                                                                                                                                                                            gttgggtacgtggggaaaaaaaCCTCCCCA‑CAGAAAAATT‑ACTCCCGCCATAGTAAGCAGGCATCCCTAATGCATCGGTATAAAATTTGAGGAGATACAGATTACCGAGGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATCCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTttat                                <  1:470634/223‑1 (MQ=255)

                                                                                                                                                                                                                 tAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTttattattat                          >  2:111583/1‑239 (MQ=255)

                                                                                                                                                                                                                   aaCG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGCTTTTTTCCCGCAGACTCAGTTTTAATTTCAGTGGCTTTGGGTCAGACATAACCCCTCCTTCTAACCCCCCCCAGGAGCTCGTttttttttatt                         >  1:705927/1‑238 (MQ=255)

                                                                                                                                                                                                                      cgggtaaaaaaaCTTCCCCACCAAAAAAA‑T‑AACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTATTACg                      <  2:225152/233‑1 (MQ=255)

                                                                                                                                                                                                                           tAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGGTTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCCCAGGGCGTCGTTTATTTTTATTAcccccgg                 >  1:913996/1‑236 (MQ=255)

                                                                                                                                                                                                                                 ccTTCGCAACAAAAAAGA‑T‑AACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTATTACGCCGTGCGTTTa           <  2:622078/238‑1 (MQ=255)

                                                                                                                                                                                                                                 catccgcagctaGAAAAA‑T‑ACCCCCGCCGTAGGAAGGATACATCCCTAATGCATAGGTATAAAATTTGAGGAGATCCAGAGTACCGATGCAGCGAATTAAATGCGAACAGACGTCGCCCATACCATAGGCGACTTTTTCCCGCAGATTCAGTTTTAATGGCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGTAGGGCGTCGTTTATTATTATTACGCCGTGCGTTTa           <  1:112242/228‑1 (MQ=255)

                                                                                                                                                                                                                                   tccccccccaAAAATA‑A‑AACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTATTACGCCGTGCGTTTACg         <  2:261854/232‑1 (MQ=255)

                                                                                                                                                                                                                                   ctCGCCACAAAAAAAA‑A‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTATTACGCCGTGCGTTTACg         <  2:361018/238‑1 (MQ=255)

                                                                                                                                                                                                                                     cgccccaaaaaaaA‑T‑AACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTATTACGCCGTGCGTTTACg         <  2:783654/228‑1 (MQ=255)

                                                                                                                                                                                                                                     cGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTATTACGCCGTGCGTTTACgg        >  2:979239/1‑238 (MQ=255)

                                                                                                                                                                                                                                          caaGAAAGA‑T‑AACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTATTACGCCGTGCGTTTACGGGCTTCt  <  2:1204148/237‑1 (MQ=255)

CATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCG‑ATAAACTGGAGCGTTGCAATTAACGCCGCCGGAA‑CCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTT‑T‑CCGCGAGTCGAGTAAAAAT‑CCGGTG‑AGCATATCGGTAAACG‑CGGTAAAAAACTTCGCCACCAGAAAGA‑T‑AATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTATTACGCCGTGCGTTTACGGGCTTCT  >  NC_000913/4065367‑4065832

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