BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,062,716	A→C	14.2%	T623T (ACA→ACC)	scpA →	methylmalonyl‑CoA mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,062,716	0	A	C	14.2%	266.2 / 3.0	113	T623T (ACA→ACC)	scpA	methylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand): ref base A (61/28); new base C (3/20); total (64/49)
Fisher's exact test for biased strand distribution p-value = 1.93e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 116 total aligned reads displayed.

CCAAGCCAGTGTGTTACCGGCGTGATTGCGCAAAGCTATCATCAGTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGC  >  NC_000913/3062480‑3062954

ccAAGCCAGTGTGTTACCGGCGTGATTGCGCAAAGCTATCATCAGTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGAACGCGGCGCGAAA‑GAGATCGCCAGCGCCTATTCCGATCTCGGTTTCTACGTAGATTTAAA‑CCCGATGTTCTCTaca                                                                                                                                                                                                                                                >  2:1028346/1‑237 (MQ=255)

  catccgtttttttgccgggttGATTGCGAAAAACTTGCATCAGGTGGAGAAATCCGCCTCCGAGGTCGATGCCATTGTTGCGGAAACGGGGCAGTTCCTTGACGACAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACTCGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTg                                                                                                                                                                                                                                            <  1:523010/219‑1 (MQ=255)

         tgtgtTACCGGCGTGATTGCGCAAAGCTATCATCAGTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGGTTCGACGTAGGTTTTAG‑CCCGAGGTTCTCTAAACCTGAAaaga                                                                                                                                                                                                                                      >  2:420565/1‑238 (MQ=255)

          gtgtTACCGGCGTGATTGCGCAAAGCTATCATCAGTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGGCACGATCGCGGCGCGAAA‑GTGATCGCCAACGCCTATTCCGATCTCGGTTTCGCCGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAagagct                                                                                                                                                                                                                                     >  2:929473/1‑236 (MQ=255)

                ccGGCGTGATTGCGCAAAGCTACCATCAGTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGAAAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATcgcccgc                                                                                                                                                                                                                              <  2:891929/239‑1 (MQ=255)

                    cGTGATTGCGCAAAGCTATCATCAGTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GGGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTTAAGAGATCGCCCGCCTg                                                                                                                                                                                                                           >  2:369507/1‑238 (MQ=255)

                     gTGATTGCGCAAAGCTATCATCAGTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCCGCGCGAAA‑GTGATCCCCAGCGCCTATTCCGATCTCGGTTTCGGCGGAAATTTAAA‑CCCGCAGTGCTCTACACCTCAagcggttcttttaaaaca                                                                                                                                                                                                                         >  1:218849/1‑223 (MQ=255)

                        aTTGCGCAAAGCTATCATCAGTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGCAA‑GTGATCCCCAGCGCCTATTCCGATCTCGGTTTCGACGGAGATTTAAG‑CCCGCAGTTCTCTCCAcccatgccctcatccacatctgacg                                                                                                                                                                                                                       >  1:371425/1‑215 (MQ=255)

                                 agcttcccttaGTCTGAAAAACTGGCCTCCGAGTTCGATGCCATTGTTTGCCAAACGGGGGAGTTCCTTGCCTACAATGGT‑CGTCG‑CCC‑GCGCAGTCCGATCGCTAAGATGGGCCCGGATGGAAACGATTGCGGCGCGAAA‑GTGATCCCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGACGTAGAAAacga                                                                                                                                                                                                             <  1:1308785/229‑1 (MQ=255)

                                 aGCTATCATCAGTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACCCCTGAAGAGATCGCCCGCCTGGCCGTAGAAAacg                                                                                                                                                                                                              >  2:1308785/1‑238 (MQ=255)

                                 aGCTATCATCAGTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAacg                                                                                                                                                                                                              >  2:462162/1‑238 (MQ=255)

                                 aGCTATACTCAATCTGAGAAAGCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTGGCCGAAAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAacga                                                                                                                                                                                                             <  2:1316235/239‑1 (MQ=255)

                                             tCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGt                                                                                                                                                                                                 >  2:14543/1‑239 (MQ=255)

                                             tCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATAGCCCGCCTGGCCGTAGAAAACGACGTTCActaagt                                                                                                                                                                                                 >  2:973524/1‑234 (MQ=255)

                                             tCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GGGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGCAGAGATCGCCCGCCGGGCCGTAGCAAACGCGGTTCACGGAgg                                                                                                                                                                                                 >  2:1180189/1‑238 (MQ=255)

                                             tCTGAGAAATCGGCCTACGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCCCGAAA‑GTGTTCGCCAGCGCCTATTCCGATCCCGGTTTCTCCGTACTTTTAGC‑CCCAGTGTCCCCTACCCCTAAAGCAATAGCCCCCCTGGACGGAGAAAAAGACGATCCCggatg                                                                                                                                                                                                 >  1:1140013/1‑235 (MQ=255)

                                              ctgggaatcctgccctcGAGTTTGAGGCCATTTGTGCGGAACGGGAGCAGGTCCTGGCCGAAAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGAGCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTa                                                                                                                                                                                                   <  1:171874/220‑1 (MQ=255)

                                              cTGGGAAAGCGGCCCCCGAGTTGGCTCCCCTTGTTGCGCAAACGGACGATTTCCTTGCCGACAAGGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGGTGGGCCAGGAGGGACCCGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTa                                                                                                                                                                                                   <  1:1044682/236‑1 (MQ=255)

                                                gagaactccgtccccggttgtgATGCCGTTGTTGGACCAACGGGGCCGTTCCCTGCCGACAATGGT‑CGGCC‑CCC‑GCGCATTCCGAGCGGTAAGATGGGCCAGGATGGCAGCGACCGCGGCGTGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACCCCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTggg                                                                                                                                                                                              <  1:230253/218‑1 (MQ=255)

                                                gagaAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGCCGTTCACGTAGTggg                                                                                                                                                                                              >  2:327775/1‑239 (MQ=255)

                                                gagaAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTggg                                                                                                                                                                                              <  2:595133/239‑1 (MQ=255)

                                                gagaAATCCGCCTCCGGGGTCGAGTCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACCCCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTggg                                                                                                                                                                                              <  2:1030824/239‑1 (MQ=255)

                                                  gaAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAACCGACGTTAACGTAGTGGgc                                                                                                                                                                                             >  2:778153/1‑238 (MQ=255)

                                                                            ttgttgCGCAAACGGAGCAGTTCCTTGCCGACAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATaaa                                                                                                                                                                  >  2:261755/1‑239 (MQ=255)

                                                                                 tgcccaaCGGGGCGTTTCGTTGCCGACAATGGTCCGACG‑CCC‑G‑GCATTCTGATCGCGAAAGTGGGCCAGGCGGGACGCGATCGCGGCCCGAAAGGTGA‑CGCCAGCGCCCATCCCGCTCCCGGTTCCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCt                                                                                                                                                             <  1:420565/234‑1 (MQ=255)

                                                                                 ggccAACCGGGGCGGTT‑CTTTCCGGCAATGGG‑CGGCCCCCC‑GCGCATTGTGATCGCTAAGATGGGCCAGGATGGACCGGAGCGGGCCTCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCGGATGTTCTCTCCACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCt                                                                                                                                                             <  1:1120904/236‑1 (MQ=255)

                                                                                 gcgcAACCGGAGCAGTTCCTTGCCGCAAATGGG‑CGTGG‑CCC‑GCGCATTCTGAGCGCTAAGATGGGCCAGGATGGACCCGTTCGCGGCGCGAAA‑GGGATCGCCAGCCCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACCCCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCt                                                                                                                                                             <  1:369507/239‑1 (MQ=255)

                                                                                 gcgcAAACGGAGGAGTTCCTGGCCGACAGTGGT‑CCTCG‑CCC‑GAGCTTTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCt                                                                                                                                                             <  2:563906/239‑1 (MQ=255)

                                                                                 gcgcAAACGGAGCAGGTCCGTGCCGACAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCt                                                                                                                                                             <  2:1130020/239‑1 (MQ=255)

                                                                                    cAAACGGAGCAGTTCCTTGCCGACAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCGGa                                                                                                                                                           >  2:1226333/1‑238 (MQ=255)

                                                                                       aCGGAACAGTTCCTTGCCGACAATCGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GGGATCGCCAGCGCCTATTCCGATCCCGGTTTCGACGTGGATTTAAG‑CCCGCTGTCCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAAAGACGATCACCGAGTGGCCGGATCCCCCCGTGCCGCCGGCCATACAACGcccacccc                                                                                                                                                       >  1:362640/1‑232 (MQ=255)

                                                                                                     tGCCGACAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTCTTCCGATCTCGGTTTCGACGTAGATTTAAG‑CACGATGTTCTCTCCAACTGGAGAGATCGCCCGCCTGGCCGTAGGAAACACCGTTTACGCTATGGGGGAACCACCACTGGCCGCGCGGCAACAAACGCTGGCCACGGAACCggggcgag                                                                                                                                         >  1:983672/1‑233 (MQ=255)

                                                                                                     tGCCGACAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACCCCTGAAGAGATCGCCCGCCTCGCCGTAGTAAAAGCAGTCTCCGAAGAGGGCGCATCCTCACTGGCTGCCGGCCATAAAACGCTGACCCCGGAACTTGTCgga                                                                                                                                          >  1:1266601/1‑236 (MQ=255)

                                                                                                     tGCCGACAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTg                                                                                                                                                                                                                                            >  1:417351/1‑140 (MQ=255)

                                                                                                     tGCCGACAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTg                                                                                                                                                                                                                                            <  2:417351/140‑1 (MQ=255)

                                                                                                     tGCCGACAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGaa                                                                                                                                          >  2:1252408/1‑238 (MQ=255)

                                                                                                     tGCCGACAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGGCATAAAACGCTGATCCCGGAACTGGTCGaa                                                                                                                                          >  2:976574/1‑238 (MQ=255)

                                                                                                         gaaaagggT‑CGTCG‑CCCTGCG‑ATTTTGATCGCTAAGGTGGGCCAGGATTGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCt                                                                                                                                     <  2:426180/233‑1 (MQ=255)

                                                                                                          aCAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCCGATCCCGGAACTGGTCGAAGCGCTg                                                                                                                                    >  2:1230103/1‑239 (MQ=255)

                                                                                                          aCAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGCCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTg                                                                                                                                    >  2:765067/1‑239 (MQ=255)

                                                                                                                       ccccc‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATggg                                                                                                                          <  2:436807/236‑1 (MQ=255)

                                                                                                                             gcgcATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGGcgcg                                                                                                                    >  2:1004790/1‑239 (MQ=255)

                                                                                                                                  ttCTGATGGCAAAGTGGGGCCCGGCGGGACCCGATCCGGGCGGGAAA‑GTTGTCCCCATGCCCTTTTCCGTTCCCGGTTTCGACGTAGATTTAAG‑CCCGATGGTCCCTACCCCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCAGCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGat                                                                                                               <  1:323392/239‑1 (MQ=255)

                                                                                                                                       aTCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAC‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCTTAGAAACCGCCGTTTCCGTAGTGGCACCACCCTACCCGGCTGCCTGTCATAAAACGCTGACCCCGGAACCGGCCTCACGCCCGAAAAAAAGGGGAAGGGAAGagaccgg                                                                                                          >  1:643701/1‑233 (MQ=255)

                                                                                                                                           cTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATAACTggttg                                                                                                      >  2:667334/1‑235 (MQ=255)

                                                                                                                                                   ggCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGCTCCCGGGAATGGTCGCAGCGCTGACAAAATGGGGACCCGAGAATATCTGCGGGGGCGCggg                                                                                              >  1:54283/1‑239 (MQ=255)

                                                                                                                                                   ggCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGaa                                                                                                                                                   >  2:447416/1‑186 (MQ=255)

                                                                                                                                                    gCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATG‑TCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCt                                                                                                                                                                                      <  1:124720/149‑1 (MQ=255)

                                                                                                                                                    gCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATG‑TCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCt                                                                                                                                                                                      >  2:124720/1‑149 (MQ=255)

                                                                                                                                                        ggATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGc                                                                                          >  2:861727/1‑238 (MQ=255)

                                                                                                                                                         gATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCGGCCTGGCCGTAGCAAACGCCGTTCACGTAGTGGGCGACTCCCAACTGACTCCGGGGCCTAAAACGCTGACCCCGGAACTGGTGGAACCCCCGAAAAAATGGGGCCCCGAAGAAAT‑TGCGAGGCTGCggggggggg                                                                                        >  1:50276/1‑235 (MQ=255)

                                                                                                                                                              acacGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCCCGTAGCGGGCCCATCCTCACTGGCTGCCGGTCATAAACCGCTGTTCCCGGAACTGGCCGAACCGCGGAAAAAATGGGGACCCGAAAATATCAGAGGGGGCGCggggggggagaatc                                                                                   >  1:678248/1‑231 (MQ=255)

                                                                                                                                                              acacGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCg                                                                                                        >  2:858342/1‑218 (MQ=255)

                                                                                                                                                              acacGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTGAtt                                                                                    >  2:467673/1‑238 (MQ=255)

                                                                                                                                                              acacGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCAtt                                                                                    >  2:782786/1‑238 (MQ=255)

                                                                                                                                                              acacGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGGCAtt                                                                                    >  2:555447/1‑238 (MQ=255)

                                                                                                                                                              acacGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCCCGGGTGGCGTCAtt                                                                                    >  2:437500/1‑238 (MQ=255)

                                                                                                                                                              acacGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGGCATAAAACGCTGATCCCGGAACTGGTCGGACGGCTGAAAAAATGGGCACGGGAAGATACCGTCGTGGACGCGGGGGGCGTCAtt                                                                                    >  1:1267689/1‑238 (MQ=255)

                                                                                                                                                                  gATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGACGAAACCTCCGTGGTCGCGGGGGCGGTCATTccgcc                                                                               >  2:1289908/1‑239 (MQ=255)

                                                                                                                                                                   gtgggggcgCACAA‑GGGATCGCCACGGGCTTTTCGGATCCCGGTGTTGACGGAGTCTTAAC‑GCGGATGTTTCCTCACCCTGAAGAGAGCGCCCGCCGGGCCGTAGAAAACGACGTTAACGTAGTGGGCGCATCCTCACGGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTccgccg                                                                              <  1:1028346/234‑1 (MQ=255)

                                                                                                                                                                   gtCGCGGGGCGAAA‑AGGGTCGCCAGCGCCTATCCCGAACTCGGTTTCGACGGAGATTTAAG‑CCCAATGTTGCCTACCCCGGAAGAGATCGCGCGCCTGGGCGTAGAAAGCGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTccgccg                                                                              <  1:1289908/238‑1 (MQ=255)

                                                                                                                                                                    gcggggcgCGAAA‑GTGATCGCCAGCGCCTATTCCGATCGCGGTTTCGACGTAGATTTAAG‑CCCGATGTCCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTccgccg                                                                              <  2:1267689/234‑1 (MQ=255)

                                                                                                                                                                               aaagttctctcccgcctcttttCCTATTTCGGTTTCGGCGTAGATTTAAG‑CCCGATGTTCTCTACCCCTGAAGGGATCCCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGc                                                                   <  1:555447/219‑1 (MQ=255)

                                                                                                                                                                                a‑GTGATCGCCAGCCCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATTTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCt                                                                 <  2:1140013/239‑1 (MQ=255)

                                                                                                                                                                                   tgaccgccgggccttATCCCGTTCGCGGTTTCGGCGTAGTTTTAAG‑CCGGATGTTCTTTACACCTGACGAGACCGCCCGCCTGGCCGTAGAAAACGACGTTCCCGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCt                                                                 <  1:778153/223‑1 (MQ=255)

                                                                                                                                                                                          cAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGTGGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAag                                                        >  2:1064731/1‑239 (MQ=255)

                                                                                                                                                                                             cgcCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGCAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGGTTACGCCTTCCTGCAagacc                                                     >  2:158396/1‑237 (MQ=255)

                                                                                                                                                                                             cgcCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAagag                                                      >  2:1191187/1‑238 (MQ=255)

                                                                                                                                                                                                    tCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGCCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGt                                              >  2:818172/1‑239 (MQ=255)

                                                                                                                                                                                                        aTCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTCGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGGCTTGAAAAAATGGGGACCCCAAGATATCTGCGTGGTCGCGGGGGGCGTCATACAGCCCCCGGATTACGCCTTCC‑GCAAGAGCgcggccgggggg                                         >  1:27751/1‑232 (MQ=255)

                                                                                                                                                                                                         tctcGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCCGGATTACGCCTTCCTGCAAGAGCGCggcggggcgg                                         >  2:1017427/1‑239 (MQ=255)

                                                                                                                                                                                                                    gcgttgtttccaG‑CCCGGTTTTCCCCACCCGTGAAGAGCCTCCCCGCCGGGCCGTAGAAAACGCCGTTTCCGTGGTGGGGGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATTTCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGt                              <  1:1230103/228‑1 (MQ=255)

                                                                                                                                                                                                                    gcGTCGATTTAAG‑CCCGATGTTCCTGACCCCTGCAGCGATCGCCCGCCTGGACTGAGAAAACGGCGTTCACGGAGGGGGCGCATCCTCACTGGCTGCCGGTCCTAAAACGCTGTTCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGCCGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGt                              <  1:1017909/238‑1 (MQ=255)

                                                                                                                                                                                                                     cgtgattgtAAG‑CCCGATGTTCTCTACACCTGAAGAGACCGCCCGCCTGTCCGTAGAAAACGACGTTCACGTAGTGGGCGCAGCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGt                              <  1:462162/230‑1 (MQ=255)

                                                                                                                                                                                                                     cggaaatgtAAG‑CCCGTTTGTCTCTCCCCATGAAGAGAGCGCCCGCCTGGCCGTAGAAACCGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGt                              <  2:935520/230‑1 (MQ=255)

                                                                                                                                                                                                                     cGTAGATTTAAG‑CCCGATGTTCTTTACCCCTGAGGAGATCGCCCGCCTGGCCGTAGAAAACGACGTGCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGGCGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGt                              <  2:689696/238‑1 (MQ=255)

                                                                                                                                                                                                                     cGTAGATTTAAG‑CCCGATGTTCTCTACACTTGAAGAGATCGCCCGACTGGCCGTAGCAAACGACGTTCACGTAGTGGGCGCCTCCTCGCTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACTCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGt                              <  2:176908/238‑1 (MQ=255)

                                                                                                                                                                                                                      gtttagtgacaagccgagtttCACTACTCCGGAGGAGTCCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTCGGGGGCGCATCCTCACTGGCTGACGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGGGAAGATATCTGCGTGGTCGCGGGTGGCGTCAGTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGt                              <  1:902953/219‑1 (MQ=255)

                                                                                                                                                                                                                       gtgATTTTAGCCCCGATCTTCTCTCCACCCGAGGAGATCGCCCGCCTGGCCGGTAAAAACGACGTTCAGGTAGTGGGCGCATCCCCACGGGGTGCCGGTCAAAAAACGCTGAACCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAAATATCTGCGTGGGCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGGAAGAGCGCGGGGTGGCGGCGATTTATGGTcc                            <  1:14543/237‑1 (MQ=255)

                                                                                                                                                                                                                         gttttAAG‑CCCGATGTTCGCTCCCCCTGAAGGGACCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAgg                         <  2:362640/237‑1 (MQ=255)

                                                                                                                                                                                                                         gATTTAAG‑CCCGAGTTTCCCATCACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGAATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAgg                         <  2:638208/239‑1 (MQ=255)

                                                                                                                                                                                                                           tttAAG‑CCCGATTTTCCCTACCCCTGAAGAGATCGCCCGGCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTa                       <  2:1022619/239‑1 (MQ=255)

                                                                                                                                                                                                                           tttAAG‑CCCGATGTTCACTACACCTGAAGAGATCGGCCACCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAgg                         <  2:1126647/237‑1 (MQ=255)

                                                                                                                                                                                                                           tttAAC‑CCCGATGTTCCCAACCCCTGAAGAGACCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTa                       <  2:1068593/239‑1 (MQ=255)

                                                                                                                                                                                                                            tttacgccctttttcctccACCCCGGAAGAGACCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCAGCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTa                       <  1:782786/224‑1 (MQ=255)

                                                                                                                                                                                                                              aaG‑CCGGATGTTCTGTACACCTGAAGAGGTCGCCCGGCTGGCCGTAGAAAACGACGGTCACGTAGTGGGCGCAACCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATTTCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTacac                    <  2:389695/239‑1 (MQ=255)

                                                                                                                                                                                                                               aacccctgtGTTCTCTACCTCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTGGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAAGGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTacac                    <  2:773559/231‑1 (MQ=255)

                                                                                                                                                                                                                               aG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGGGCGGGGGGgcggcgggtttggggcaaggaaaacc                   >  1:1013496/1‑219 (MQ=255)

                                                                                                                                                                                                                                   ccaaTGTTCCCTACCCCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCt                  <  2:1265187/234‑1 (MQ=255)

                                                                                                                                                                                                                                     gATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAACAATGGGGCAGCGCAGAATTCTCCGCGGTGCCGGGGGGCGTCATTCCCCCGCGGGGTTACCCTTTCCGGCATGAGCGCGCGGTGGGAGCGCTTTagtgttcggcggacatcagac              >  1:430882/1‑219 (MQ=255)

                                                                                                                                                                                                                                     gATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGGGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTACGGTCCAGGTCAACCTAt                >  2:374903/1‑237 (MQ=255)

                                                                                                                                                                                                                                     gATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGACGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCCGCGTGGTAGCGGGTGGCGTCATTCCGCCGCGGGGTTCAGCCTTCCTGCGAGAGGCCGGCGGGGGGGCGGTTTATGGCTCAGGGAccccatgc               >  1:824453/1‑234 (MQ=255)

                                                                                                                                                                                                                                      aTGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGGGGCGTCATTCGGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGCGGTGGCGGGTGTTTATGGCCCAGGAACACCTATgt              >  2:1326160/1‑237 (MQ=255)

                                                                                                                                                                                                                                       tGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGTGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTTCGTGGTCGCGGGTGGCGGCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGGGATTTTTGGTCCAGGAAAAccaaagccc            >  1:1062213/1‑232 (MQ=255)

                                                                                                                                                                                                                                       tGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATgtt             >  2:827463/1‑236 (MQ=255)

                                                                                                                                                                                                                                         ttCTCTACCCCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGa          <  2:54283/239‑1 (MQ=255)

                                                                                                                                                                                                                                           ctctACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGCTTTATGGTCCAGGt                        >  2:115070/1‑223 (MQ=255)

                                                                                                                                                                                                                                           ctctACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACa        >  2:769823/1‑239 (MQ=255)

                                                                                                                                                                                                                                           ccctccacCTGAAGGGATCCCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACa        <  2:678248/234‑1 (MQ=255)

                                                                                                                                                                                                                                              tACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACAccctgccccttataca     >  2:902953/1‑225 (MQ=255)

CCAAGCCAGTGTGTTACCGGCGTGATTGCGCAAAGCTATCATCAGTCTGAGAAATCGGCCTCCGAGTTCGATGCCATTGTTGCGCAAACGGAGCAGTTCCTTGCCGACAATGGT‑CGTCG‑CCC‑GCGCATTCTGATCGCTAAGATGGGCCAGGATGGACACGATCGCGGCGCGAAA‑GTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAG‑CCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGC  >  NC_000913/3062480‑3062954

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