BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,662,349	A→T	6.1%	intergenic (‑220/+232)	iscR ← / ← trmJ	isc operon transcriptional repressor, suf operon transcriptional activator, oxidative stress‑ and iron starvation‑inducible, autorepressor/tRNA mC32,mU32 2'‑O‑methyltransferase, SAM‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,662,349	0	A	T	6.1%	245.7 / 4.1	97	intergenic (‑220/+232)	iscR/trmJ	isc operon transcriptional repressor, suf operon transcriptional activator, oxidative stress‑ and iron starvation‑inducible, autorepressor/tRNA mC32,mU32 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand): ref base A (39/45); new base T (2/4); total (44/53)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.98e-01

Only 100 of 108 total aligned reads displayed.

CTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGTTTTTTATTGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGT‑TCACTGTAGGCCTGATAAGACG‑CATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCG  >  NC_000913/2662111‑2662587

cTTTAGATGTTCGTCTCATGTCTTATTTCATCTTAAACTCGCCCCGGCCCGGGGTTTTTTATTTTAAAAGGGGGGGTTTTGTATAGCAGGGTCAAGTCTGACCTTCCCGAGTAAA‑TTGGTCAACTATTTACGTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCTGGGCGTTTGTCCCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTa                                                                                                                                                                                                                                                  <  2:221729/239‑1 (MQ=255)

cTTTAGATGTCATTCTCATGTCTTCATGCACCGCAACCTCGCCCCTGCACGGGGTTTTTTATTGTAAAAGTGGGGGTATTGCATTGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTa                                                                                                                                                                                                                                                  <  1:1106694/239‑1 (MQ=255)

cTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGTTTTTTATTGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTa                                                                                                                                                                                                                                                  <  1:176360/239‑1 (MQ=255)

 tttAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGTTTTTTATTGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAAGGCCATATCCG‑‑CCCTGGGGTTCTGTa                                                                                                                                                                                                                                                  >  1:791377/1‑238 (MQ=255)

 tttAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGTTTTTTATTGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTac                                                                                                                                                                                                                                                 >  1:221729/1‑239 (MQ=255)

     gttgtcagttctgtGTCTTAATTCGCCTCAACCCCGCCCCGTCCCGGGGTTTGTTATTGTAAAAGTGGGGGTTTTGCATAGCAGGGTTAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGGCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCTTATCCG‑‑ACCTACGGTTCTGTTCACTg                                                                                                                                                                                                                                             <  2:744658/226‑1 (MQ=255)

         tCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGTTTTTTATTGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTaaccctgt                                                                                                                                                                                                                                           >  1:207476/1‑230 (MQ=255)

          cAGTCTCATGGCTTACTTGCCCTCAAACTCGCCCCTGCCCGGGGTTTGTTAATGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGc                                                                                                                                                                                                                                        <  1:355773/239‑1 (MQ=255)

           aGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGTTTTTTATTGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTTTAgg                                                                                                                                                                                                                                         >  1:686868/1‑237 (MQ=255)

           aGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGTTTTTTATTGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGGATTTAACCTTCAATGCCATTTTTTTATAGTTGCGTTTGTACGTGTTGTCGGATGCGGCGTGAGCGCCCTATCCG‑‑CCCTTGGGTTCTGTACACCGTatg                                                                                                                                                                                                                                         >  1:658664/1‑235 (MQ=255)

           aGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGTTTTTTATTGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGGCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGGATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTTTACCTGTTGTCGGAGCCGGGCTGAACGCCCAACCCG‑‑CCCCGGGGGTCGTTAACCTGTCGGcc                                                                                                                                                                                                                                       >  2:794572/1‑239 (MQ=255)

              ctcATGTCTTACTTCCCCTTAAACCCGCCCCTGCCCGGGGTTTTTTATTGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGcc                                                                                                                                                                                                                                       <  1:642371/236‑1 (MQ=255)

                  tGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGTTTTTTATTGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTAAACTGTAGGCCTGATAAg                                                                                                                                                                                                                                >  1:150488/1‑239 (MQ=255)

                   aagttACTTCCCCTCAAACTCGCCCCTGCCCGGGGTTTTTTATTGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAg                                                                                                                                                                                                                                <  1:1030980/235‑1 (MQ=255)

                       tACTTCACCTCAAACTCGCCCCTGCCCGGGGTTTTTTATTGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGTTAAGACGCa                                                                                                                                                                                                                           >  1:553576/1‑239 (MQ=255)

                        aCTTCACCTCAAACTCGCCCCTGCCCGGGGTTTTTTATTGTATAAGTGGGGGTATTGCATAGCAGGGCCAAGTCTGACATTCCCGAGTCAA‑TTGGTCAACTATTTACTTGACTGATTTCGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTAACTGTTGTCGGATGCGGCGTGAACGCGATATCCG‑‑TCCTGGGCTTTTGTACCCTGTAGGCCGGACAACACGCAt                                                                                                                                                                                                                          >  1:709473/1‑239 (MQ=255)

                             aCCTCAAACTCGCCCCTGCCCGGGGTTTTTTATTGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTTGGATGCGGCGTGAACGCCATATTCG‑‑CCCTGGGGTTTTGTACACTTTAGGCCTGATAAAACGCATTAcg                                                                                                                                                                                                                      >  2:388566/1‑238 (MQ=255)

                             aCCTCAAACTCGCCCCTGCCCGGGGTTTTTTATTGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATCCGGCGTGAACGCCCTATCCG‑‑CCCTGGGGTTCTTTACACTGTAGGCCCTATAAAACGCCTTAcc                                                                                                                                                                                                                      >  2:62030/1‑237 (MQ=255)

                              ccTCCAACCCGCCCCTGCCCTGGGTTTTTTTTTGTAAAAGGGGGGGTATTGCTAAGCAGGGTAAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTTTTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTAcgcg                                                                                                                                                                                                                    <  2:559639/239‑1 (MQ=39)

                              ccTCAAACTCGCCCCTGCCCGGGGTTTTTTATTGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTa                                                                                                                                                                                                                                                  >  2:530296/1‑209 (MQ=255)

                              ccTCAAACTCGCCCCTGCCCGGGGTTTTTTATTGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTa                                                                                                                                                                                                                                                  <  1:530296/209‑1 (MQ=255)

                              ccTCAAACTCGCCCCTGCCCGGGGTTTTTTATTGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTAcgcg                                                                                                                                                                                                                    <  1:387729/239‑1 (MQ=255)

                                 cAAACTCGCCCCTGCCCGGGGTTTTTTATTGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAAGCCCATATCCC‑‑CCCTGGGGTTTTGTCCACTTTAGGCCTGGTAAGACGCAtttccccg                                                                                                                                                                                                                   >  2:1030980/1‑231 (MQ=255)

                                    aCTCGCCCCTGCCCGGGGTTTTTTATTGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTa                                                                                                                                                                                                                                                  >  2:528834/1‑203 (MQ=255)

                                    aCTCGCCCCTGCCCGGGGTTTTTTATTGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTa                                                                                                                                                                                                                                                  <  1:528834/203‑1 (MQ=255)

                                           ccTGCCCTGGTTTTTTTTTTGTAAAAGTGGGGTTTTTTCATTGCAGGGTCAAGGCTGACATTCCCGAGTAAA‑TTTGGCAAATATTTGCTTGGCTGATTTAGTCGGGGATTTAACCTTCAATGCCATTTTTTTATCGGTGCGTTTGTTCCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTTGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCaa                                                                                                                                                                                                       <  2:552243/239‑1 (MQ=25)

                                                    ggTTTTTTATTGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACCCCATATCCG‑‑CCCTGGGGTCTTTTACACTGTAGGCCTATTAAGCCGCATTACCCGCCCCATCAGGCAACgccctgcg                                                                                                                                                                                              >  2:42814/1‑232 (MQ=255)

                                                          ttAGTGTAAAAGGGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATgcg                                                                                                                                                                                        <  1:70612/239‑1 (MQ=255)

                                                               ttAAAAGTTGGGGTTTTGCAATGGTGGGTCAAGGTTGGCATTCCCGAGTAAA‑TTGGTCAAATTTTTTTTTGACTTATTTAGTCGGGTATTTAACCCTCAGTGCCATTTTTTTATCGTGGCGCTTGTACCTTTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTg                                                                                                                                                                                   <  2:881040/238‑1 (MQ=25)

                                                                tAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGCCTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCCTATCCG‑‑CCCTGGGGGTCTGTTCACTTTAGGCCCGATAAGACGCATTTCCCGTCGCATAAGGGAACgccgtgccgttggagcgggg                                                                                                                                                                                  >  2:422078/1‑220 (MQ=255)

                                                                tAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAAACGCATTACCCGCCGCCTCAGGCAACGGCGGTGGGATggcacgcgg                                                                                                                                                                                  >  2:347151/1‑231 (MQ=255)

                                                                   aaGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATgcgacgcgtaat                                                                                                                                                                               <  1:347151/239‑7 (MQ=255)

                                                                   aaGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGGTCGGTACCCTGTAGGCCTGCAAAGACGCATTACGCGTCGCATAAGGAACAGGCTGTCGGGTggcgggggaac                                                                                                                                                                                >  2:388561/1‑228 (MQ=255)

                                                                   aaGGGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACg                                                                                                                                                                               <  1:422078/239‑1 (MQ=255)

                                                                     ttGGGGGTTTTTCTTATCCGTGTTAAGTGTGACTTTCCCGTGTAAA‑TTGGTCAACGTTTTTCTTTACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTAGCGGGGCGTTTGTACCTGTTGTCGGATGCGGCGTTAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGcc                                                                                                                                                                             <  2:502930/238‑1 (MQ=21)

                                                                          ggTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGAGTCGGCGTGAACGCCtgtt                                                                                                                                                                         >  1:282233/1‑235 (MQ=255)

                                                                             acggtattgcgtggttaatgtgtaCATTCCCGGTTAAA‑TTGGTCAATTATTTACTTGGCTGATTTTGTCGGGTATTTACCCTTCAGTGCCATTTTTTTATCGGTGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                                                      <  2:474575/215‑1 (MQ=21)

                                                                               tGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTCCACTGTAGGCCGGGTAAGACGCATTACCCGCCGATGCAGGCAACGGCGGTCGGATgcgggcgaaaccaccttatcaac                                                                                                                                                                    >  2:38330/1‑219 (MQ=255)

                                                                                 tattgtagggttaatgttGACATTCCCGATTAAA‑TTGGTCAACTATTTTCTTGACTGATTTAGGTGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCGGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTa                                                                                                                                                                 <  2:847587/222‑1 (MQ=25)

                                                                                  aTAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACCGGTCGCATCAGGCAACGGCTGTCTGATGAGGGCTGAAAGCCTTATCCGAcccaa                                                                                                                                                                >  1:196981/1‑236 (MQ=255)

                                                                                   tAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCCGATGGGGCGTGAAACGCTTATCCGcccaccg                                                                                                                                                               >  1:646877/1‑235 (MQ=255)

                                                                                    aGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCTTGAAAGCCTTTTCCGAccaccgg                                                                                                                                                              >  1:847587/1‑234 (MQ=255)

                                                                                     gCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATTCGGCGTGAAACCCTTATCCGACCTACGGt                                                                                                                                                             >  1:211407/1‑239 (MQ=255)

                                                                                      cAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTCGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGTTGCGGCGTGAACACCTTATCCGCCCCACCGtt                                                                                                                                                            >  1:459007/1‑239 (MQ=255)

                                                                                       aGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGtt                                                                                                                                                            <  1:698446/238‑1 (MQ=255)

                                                                                       aGGGTAAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGtt                                                                                                                                                            <  1:424072/238‑1 (MQ=255)

                                                                                        gggTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCTTATCCG‑‑ACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                                                      >  1:11905/1‑227 (MQ=255)

                                                                                        gggTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCTTATCCG‑‑ACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                                                      <  2:11905/227‑1 (MQ=255)

                                                                                         ggTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGCACGCCCTATCCG‑‑CCCTGTGGTTCTGTACACTGTAGGCCTGGTAAGACGCATTCAGCGTCGCATCCTGCAACCGCCGTCGGTGCCGGCGGCAACCCCTTTTcccccccccagcctct                                                                                                                                                         >  2:617193/1‑225 (MQ=255)

                                                                                          gTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACCCGTCGCATCAGGCAACGGCTTTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCttt                                                                                                                                                        >  1:944109/1‑237 (MQ=255)

                                                                                           tCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACCCTTTAGGCCTGATAAGACGCATTACCCGCCGCATCTGGCAACGGCTTTCGGATCGGGGGGGAACGCCCTTTCCCAcccaagggctcggg                                                                                                                                                       >  2:527492/1‑227 (MQ=255)

                                                                                           tCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGAc                                                                                                                                                                                                                              >  1:356968/1‑168 (MQ=255)

                                                                                           tCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGAc                                                                                                                                                                                                                              <  2:356968/168‑1 (MQ=255)

                                                                                           tCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGTTGCGGCGTGAACGCCTTAACCGGCCCACGGGTcccgt                                                                                                                                                       >  1:474575/1‑235 (MQ=255)

                                                                                            caattttggcATTTCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTTGTCGGGTATTTAACCTTCAGGGCCATTTTTTTATCGTGGCGTTTGTACCTTTTGTGGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGt‑t                                                                                                                                                      <  2:635775/229‑1 (MQ=25)

                                                                                            aaagttttgcATTCCCAGGTAAATTTGGTC‑ACTTTTTTCTTTACTGATTTAGGCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGGGGCCTTTGTACCTGGTGGCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGt‑t                                                                                                                                                      <  2:944109/229‑1 (MQ=18)

                                                                                            aaagtgagccattcctagatAAA‑TTGGTCAACTTTTTTCTTGGCTGTTTTTGTGGGGTTTTGAACCTTGAATGTCATTTTTTTATCGTGGTGTTGTTCCCTGTTGTCGGTTGCGGCGGGAACGCCATACCCG‑‑GCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGt‑t                                                                                                                                                      <  2:791377/219‑1 (MQ=17)

                                                                                              aGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGGTGCGGCGTGACACCCTTATccccccaccgggttctttta                                                                                                                                                     >  2:633595/1‑220 (MQ=37)

                                                                                                tCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGAAACGGCCGTCGGATGCGGGGGGAAACCCCTATCCCGCCTACGGGTTTGt‑ttc                                                                                                                                                    >  2:248654/1‑234 (MQ=37)

                                                                                                  tGACATTCCAGAGTAAA‑GTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGT‑TCACTGTa                                                                                                                                               <  1:42814/239‑1 (MQ=255)

                                                                                                  ctgcgttCCCGAGTAAA‑TGGGTCAGTATTTTACTGGAGTGTTTTATGCGGTTATTTAACCTTCAGTGCCATTTTTTTTTCGTGGGGTTTGTACCTGTGGTCGGATGCGGCGTGAACGCCGTATCCGCCCCCTGGGGTTCTGTACACTGGAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGT‑TCACTg                                                                                                                                                 <  2:22890/234‑1 (MQ=18)

                                                                                                       attaccctgttaatgtgTCAACTATGTCCTGGACCGTTTTAGGTGGGTTTTTACCCTTTAGTGCCATTTTTTTATCGTGGCGTTGTTTCCGGTTGTCGGATGCGGCGTGAACGCCATAGCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGT‑TCACTGTAGGc                                                                                                                                            <  2:150488/222‑1 (MQ=17)

                                                                                                         cccGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAAAGCCTTATCCGACCTACGGGTCTGT‑TCACTGTAGGCCTgt                                                                                                                                        >  1:191705/1‑238 (MQ=255)

                                                                                                               tAAA‑TTGGTCAACTATTTTCTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGT‑TCACTGTAGGCCTGATAAGa                                                                                                                                   <  1:388561/238‑1 (MQ=255)

                                                                                                                          atcTATTTGCTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGT‑TCACTGTAGGCCTGATAAGACG‑CATTAcgcg                                                                                                                       <  1:990510/237‑1 (MQ=255)

                                                                                                                           tcgatttttttGACTTGTTTAGTCGGTGATTTAACCTTCATGGTCATTTTTGTAGCTGGGCGTTTGTTCCTGTTGTCGTATTGGGTGGTGACGCCATATGCG‑‑CCTTGGTGTGCGGTACCTGGTAGGCCTGTTAAGACGAATTGCGCGACGCACCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGT‑TCACTGTAGGCCTGATAAGACG‑CATTACGCGt                                                                                                                      <  2:709473/230‑1 (MQ=9)

                                                                                                                            cTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGt‑t                                                                                                                                                      <  1:981519/207‑1 (MQ=255)

                                                                                                                                  aCTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACgg                                                                                                                                                                                                    >  2:953380/1‑156 (MQ=255)

                                                                                                                                  aCTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACgg                                                                                                                                                                                                    <  1:953380/156‑1 (MQ=255)

                                                                                                                                  aCTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATCCGGCGGGAACGCCCTTTCCGCCCTCCGGTTCTTT‑TCCCTTTAGGCCCGGTAAAAGC‑CATTAccccccccctca                                                                                                               >  2:873291/1‑228 (MQ=21)

                                                                                                                                    ttGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTCCCGTTCTGT‑TCACTGTAAGCCGAGTAAAGAG‑CATTACGCGACGCATCAgg                                                                                                             >  1:170893/1‑239 (MQ=21)

                                                                                                                                       aCTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCCACGGCCGTCGGAGGCGGCGCGAACACCTTATCCTACCTCCGGTTATGT‑GCCCTTTAGGCCCGATAAAAccccttcccccgccccacaggcaa                                                                                                           >  2:852676/1‑216 (MQ=21)

                                                                                                                                       aCTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGTTCTGT‑TCCCTGTAGGCCTGATAAGACG‑CCTTACGCGTCGCCTCAGGCa                                                                                                           >  1:652451/1‑238 (MQ=255)

                                                                                                                                       aCTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGT‑TCACCGTAGGCCCGATAAAACC‑CCTTACCCGGCCCATCAGGCac                                                                                                          >  1:381215/1‑238 (MQ=21)

                                                                                                                                       aCTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCCTTTCCGACCAACGGTTCTGT‑TGCCGTTAAGCCTGATAAGACG‑GATTAcggctgccatccggga                                                                                                           >  1:375650/1‑224 (MQ=21)

                                                                                                                                                     gTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCACCGGCTGTCGGATGCGGCGTGAACCCCCTTTCCGCCCTTCGGTTCTGTCTCCATTTA‑GCCTGCGAA‑ACGACATTAACCGTCGCCTCATGCAACGGCTGTCggcagc                                                                                            >  2:150943/1‑235 (MQ=18)

                                                                                                                                                               ttCAGTGCCATTTTTTTATCGTGGCGTTTGTCCCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGT‑TCACTGTAGGCCTGATAAGACG‑CATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACg                                                                                  <  1:248654/239‑1 (MQ=255)

                                                                                                                                                               ttCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                                                      >  1:424547/1‑157 (MQ=255)

                                                                                                                                                               ttCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGGAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                                                      <  2:424547/157‑1 (MQ=255)

                                                                                                                                                               ttCAGGGCCATTTTTTTATCGTGGCGTTTGTTCCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGT‑TCACTGTAGGCCTGATAAGACG‑CATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACg                                                                                  <  1:633595/239‑1 (MQ=255)

                                                                                                                                                                  cggtccATTTTTTGATGGTGGCGTTTGGTCCTGTTGGCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGT‑TCACTGTAGGCCTGATAAGACG‑CATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCt                                                                               <  1:150943/235‑1 (MQ=255)

                                                                                                                                                                  aGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGT‑TCACTGTAGGCCTGGTAAGACC‑CATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGGGTGAACGCCt                                                                               >  1:157630/1‑239 (MQ=255)

                                                                                                                                                                       cATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGGTCTGT‑TCACTGTAGGCCTGATAAGACG‑CATTACGCGTTGCATCAGGCAACGGCTGTCGGATGCGGCGTGGACGCCCTATc                                                                           >  1:776875/1‑238 (MQ=255)

                                                                                                                                                                        aTTTTTTTATCGGGGCGTTTTTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGT‑TCACTGTAGGCCTGATAAGACG‑CATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCg                                                                         <  1:527492/239‑1 (MQ=255)

                                                                                                                                                                        aTTTTTTTAGGGTGGTGTTTGTCCCTGTTGTCGGGTGGGGGGTGAACGCCTTTGCCG‑‑CCCTGGGGTTTTGTGCGCTGTAGGCCTGATAAGACGCATTACGCGTCGCATCCGGCAACTGCTGTCGGATGCGGCGTGAACGCCTTATCCGCCCTGCGGTTCTGT‑TCACTGTAGGCCTGATAAGACG‑CATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCg                                                                         <  2:652451/239‑1 (MQ=17)

                                                                                                                                                                        aTTTTTTTAACGTGGCGTTTGTTCCTGTTGTCGGATGCGGCGTGAACGCCATACCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGT‑TCACTGTAGGCCTGATAAGACG‑CATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCg                                                                         <  1:876110/239‑1 (MQ=255)

                                                                                                                                                                                  cGTGGTGTTTGTACCTGTTGGCGGATGCGGTGCGAACCCCATATCCG‑‑CCCTTGGGTTCTTTAACTTGTAGGCCTTATAAGCAGCATTACTCGTCGCCGCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGT‑TCACTGTAGGCCTGATAAGACG‑CATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGtt                                                               <  2:776875/239‑1 (MQ=18)

                                                                                                                                                                                   gatgggtaaaagagacaggtCGGATGCGGCGTGAACCCCTTATCCG‑‑CCCTGGGGTTCTGTACACTTTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGT‑TCACTGTAGGCCTGATAAGACG‑CATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGtt                                                               <  1:245011/220‑1 (MQ=17)

                                                                                                                                                                                       gtttggttgcgtttttCGGGTGAGTCGTGACCGCCTTTCCCG‑‑CGCTGGGGTTCTGTGCTATGTAGGCCTGAGGAGACGCTGTACGCTTTGCCGCTGGCAACGGTTGGCGGATGCGGCGTGAACGCCTTAGCCGACCTAGGGTTCTGT‑TCACTGTAGGCCTGATAAGACG‑CATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGt                                                           <  2:282233/223‑1 (MQ=14)

                                                                                                                                                                                        gtttggcacttttgtCGGATGCGGGGTGAACGCCATTTCCG‑‑CCCTGGGGTTCTGTGCACTGTAGGCCTGATAAGCAGCATTGCGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGT‑TCACTGTAGGCCTGATAAGACG‑CATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGtt                                                          <  1:312593/227‑1 (MQ=17)

                                                                                                                                                                                        gTTTGTGCCTGTTGTCGGATGCGGTTGGAACGCCATTGCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGT‑TCACTGTAGGCCTGATAAGACG‑CATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGtt                                                          <  1:70228/238‑1 (MQ=255)

                                                                                                                                                                                         ttttttccTGTTTTCGGGTGCGTGTTGAACGGCTTTTCCC‑‑CCCTGGGGTGCTGTACTGTGTAGCCCTGATAGCCCGCTTTCGCCGGCGCATCAGGTAACGGGTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGT‑TCCCTGTAGGCCTGATAAGACG‑CATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCa                                                        <  2:686868/233‑1 (MQ=17)

                                                                                                                                                                                             gaCCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGT‑TCACTGTAGGCCTGATAAGACG‑CATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGt                                                    <  1:874703/238‑1 (MQ=255)

                                                                                                                                                                                                   ttgtCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGT‑TCACTGTAGGCCTGGACAGACG‑CCTTACGCGTCGCATCAGGCAACGGCTGGCGGGTGCGGCGTGAACGCCCTATCCGACCCAccggtagtttgaccgtaggccg                                              >  2:610003/1‑218 (MQ=14)

                                                                                                                                                                                                                  tGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGT‑TCACTGTAGGCCTGATAAGACG‑CATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTa                                <  1:610003/238‑1 (MQ=255)

                                                                                                                                                                                                                     aCGCCTTATCCC‑‑CCCTGGTGTTCTGTACACTGTAGGCCTGATAAGGCGAATTTCGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGT‑TCACTGTAGGCCTGATAAGACG‑CATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTAcgcg                            <  1:38330/239‑1 (MQ=255)

                                                                                                                                                                                                                                  tccttGGGGTGCTGTGCAATTTAGCCCTGATAAGGCTCTTTTCGTGTCGCCTCAGGCACCGGCTGTCGGATGCGGCGTGAACGCCTTATCCGGCCTACTGTTCTGT‑TCACTGTAGGCCTGATAAGACG‑CATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGc                 <  2:191705/235‑1 (MQ=14)

                                                                                                                                                                                                                                     cTGGGGTTCTGGTCACTGTAGGCCTGATAAGACGCATTTCGCGTCGCAGCAGGCAACGGCTGTCGGAAGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGT‑TCACTGTAGGCCTGATAAGACG‑CATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCaa               <  1:617193/238‑1 (MQ=255)

                                                                                                                                                                                                                                        gggTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGT‑TCa                                                                                                                                                    >  2:125648/1‑103 (MQ=33)

                                                                                                                                                                                                                                                 acacTGTAGGCCTTATAAGCCGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGT‑TCACTGTAGGCCTGATAAGACG‑CATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCg  <  1:471602/239‑1 (MQ=255)

CTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGTTTTTTATTGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAA‑TTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCG‑‑CCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGT‑TCACTGTAGGCCTGATAAGACG‑CATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCG  >  NC_000913/2662111‑2662587

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