BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,209,826	T→G	14.8%	noncoding (54/76 nt)	gltV →	tRNA‑Glu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,209,826	0	T	G	14.8%	196.7 / 3.7	110	noncoding (54/76 nt)	gltV	tRNA‑Glu
Reads supporting (aligned to +/- strand): ref base T (31/55); new base G (18/2); total (53/57)
Fisher's exact test for biased strand distribution p-value = 1.07e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 114 total aligned reads displayed.

CTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTC‑GGGTTGTGAGGTTAAGCGACTAAGCGTACAC  >  NC_000913/4209588‑4210063

cTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGt                                                                                                                                                                                                                                                  >  1:866060/1‑239 (MQ=31)

    ctcTTCCGGAGGGCGTTTACCACTTTGTGATTCATGACTGGTGTGAAGGCGTAACAAGGTGACCTTAGGTGAACCTGCTGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑Ga                                                                                                                                                                                                                                             <  2:1064739/237‑1 (MQ=11)

     ccTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑Gaa                                                                                                                                                                                                                                            <  1:1031658/239‑1 (MQ=31)

     ccTTCGGGAGGACGCTTACCACGTTGTGATTCATGTCTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGAAACCTGCGGTTGGATCACCTCCTTTCCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑Gaa                                                                                                                                                                                                                                            <  2:23758/239‑1 (MQ=9)

         cGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GGAtac                                                                                                                                                                                                                                         >  1:646557/1‑236 (MQ=14)

         cGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATcc                                                                                                                                                                                                                                         >  1:1009851/1‑238 (MQ=31)

         cGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATcc                                                                                                                                                                                                                                         >  1:459891/1‑238 (MQ=31)

                    ttACCACGTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACg                                                                                                                                                                                                                              <  2:1107314/238‑1 (MQ=14)

                     tACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGcc                                                                                                                                                                                                                            <  1:885973/239‑1 (MQ=31)

                           tttGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACtt                                                                                                                                                                                                                        >  1:32697/1‑237 (MQ=31)

                             tgtgATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGGAGCCACttgttg                                                                                                                                                                                                                    >  1:1079337/1‑239 (MQ=14)

                                    cATGCCTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGGTCTGATACCCCAGGGGGCCCCCACTTCCGGGTTtgt                                                                                                                                                                                                              >  2:137882/1‑239 (MQ=11)

                                    cATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCCCGGGGGAGCCACTTGCTGGTTtgtg                                                                                                                                                                                                             >  2:344730/1‑239 (MQ=9)

                                      tGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGCTAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGA‑AACAGGAGCAC‑GAATCCCCCAGGGGGCGGCACTTGCTGGGTTgtgag                                                                                                                                                                                                           >  1:1123154/1‑239 (MQ=11)

                                              gTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGt                                                                                                                                                                                                    <  1:979714/238‑1 (MQ=31)

                                               tGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGG‑AACAGGGGCTC‑GCATCCCCTCAGGGGCCCCACTTCGTAGGTTgtgaggggaaatc                                                                                                                                                                                                   >  2:209285/1‑230 (MQ=11)

                                                gAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCAc                                                                                                                                                                                                 <  1:427766/239‑1 (MQ=31)

                                                gAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCAc                                                                                                                                                                                                 <  2:236257/239‑1 (MQ=31)

                                                gAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCAc                                                                                                                                                                                                 <  1:201225/239‑1 (MQ=31)

                                                gAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCAc                                                                                                                                                                                                 <  1:1083635/239‑1 (MQ=31)

                                                gAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCGTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGTGGG‑AACAGGGGTTC‑GAATCGCTTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGGCAc                                                                                                                                                                                                 <  2:688173/239‑1 (MQ=9)

                                                 aaGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCAc                                                                                                                                                                                                 <  2:667102/238‑1 (MQ=31)

                                                  aGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCAcc                                                                                                                                                                                                <  1:165295/238‑1 (MQ=31)

                                                      gTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATcccc                                                                                                                                                                                                                                       >  1:980266/1‑195 (MQ=31)

                                                      gTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCAc                                                                                                                                                                                                 <  2:576462/233‑1 (MQ=31)

                                                      gTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGC‑AAAAGGGGTTC‑GAATCCCCCTGGGGGAGCCCCTTTCTGGTTTGTGAGTGGAAGTTACCTGc                                                                                                                                                                                             >  1:1054486/1‑237 (MQ=11)

                                                       tAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAAGGGACGCCACTTGCTGGTTTGTGAGTGGAAGACACCTGCCtt                                                                                                                                                                                          >  2:954896/1‑239 (MQ=14)

                                                              gTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACCGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTTATAtc                                                                                                                                                                                    >  1:596896/1‑238 (MQ=14)

                                                               tAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATAtct                                                                                                                                                                                   <  2:459891/238‑1 (MQ=31)

                                                                 aCCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCaa                                                                                                                                                                                <  1:1154619/239‑1 (MQ=31)

                                                                    ggatgtgaagttGCGGTTGGATCACCTCCTTGCCTTAAAGACGTGTTCTTGGCAGTGCTCACGCAGATTGACTGATAGAAAGTGAAGAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTACGTCTATAGGCCCAGGACACCGCCCTTTCCCGGCGGT‑ATCAGGGATTC‑GAATCCCCTAGGGGACGCCACTTGCTTGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑Ac                                                                                                                                                                            <  2:1066031/228‑1 (MQ=11)

                                                                    gTAGGGTAACGTGCGGTTGGATCACCTCCTTACGTTAAAGCAGCGTTCTGTGCAGTGCGCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCGTGCGAAGCAGACTGATAGGTCCCCTTCGTCTAGAGGCCCGGGACACCGCCCGTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATACCTCAAA‑Ac                                                                                                                                                                            <  2:637734/239‑1 (MQ=11)

                                                                       ggggAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGC‑AACCGGGGGTC‑GAATCCCCCAGGGGACGCCACCTGCTGGTTTGTGAGGGAAAGTCCCACGCCTTAATATCCTaaa‑aat                                                                                                                                                                           >  2:770175/1‑235 (MQ=11)

                                                                             ccTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGGTC‑GCATCCCCCAGGGGAAGGCCCCTGCTGGTTTGGTGGGGAAAGTCACCCTGCTTAAATTCTCAAA‑ACTCCTCtttg                                                                                                                                                                   >  2:911313/1‑237 (MQ=11)

                                                                             ccTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAAACAGGGG‑TC‑GAATCCCCTAGGGGACGCCACCTGCTGGTTTGTGTGGGAAAGTCACCCGCCTTAATCACTCAAA‑ACTACTCTTCg                                                                                                                                                                   >  2:1105806/1‑239 (MQ=11)

                                                                             ccTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGGTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCg                                                                                                                                                                   >  1:601501/1‑239 (MQ=14)

                                                                               tGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCggg                                                                                                                                                                 <  1:215707/239‑1 (MQ=31)

                                                                                gctgTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGATGCAGACTGATACGTCCCCTTCGTCTAGGGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGt                                                                                                                                                                <  1:918535/236‑1 (MQ=14)

                                                                                gCGGTTGGATCACCTCCTTACCTTAAAGAAGCGCTCTTTGCAGTGCGCACACAGATTGTCTTATAGAAAGTGAAAAGCATGGCGTCTTGCGAATCAGACTGATGCGTGCCCTTCGTCTAGAGGCCCAGGACAGCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGt                                                                                                                                                                <  2:382184/239‑1 (MQ=11)

                                                                                  ggTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTg                                                                                                                                                               <  1:11577/238‑1 (MQ=31)

                                                                                        aTCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTg                                                                                                                                                               <  1:214772/232‑1 (MQ=31)

                                                                                         tCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTga                                                                                                                                                       <  1:954896/239‑1 (MQ=31)

                                                                                              tcctTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAAACCCCTAGGGGACGCCACTTGCTGGTGTGTGGGTGAAAGGCACCTGCCTTAATATCTCAAA‑ACTCATCTTTGGGGGAGGTTTggagtgt                                                                                                                                                  >  2:116390/1‑233 (MQ=11)

                                                                                               ccttaccttaAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCAc                                                                                                                                                                                                 <  1:643372/192‑1 (MQ=31)

                                                                                                 ttaccttaAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTg                                                                                                                                               <  1:1143198/239‑1 (MQ=31)

                                                                                                  taccttaAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGGTC‑GAATCCCCTAGGGGGCGCCCCTTCCTGGTTTGTGTGTGAAAGACACACGCCTTATTATCTCAAA‑ACTCATCTTCGGGGGATGttttgttgtcttct                                                                                                                                              >  2:214772/1‑228 (MQ=11)

                                                                                                   accttaAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGGTC‑GAATCCCCTAGGGGGCGCCACTTGCTGGTTTGTGAGTGAAAATCACCTGCCTTAATATCTCAAA‑ACTCATATTCGGGTGATGTTTGAGATATTTGct                                                                                                                                             >  1:416868/1‑239 (MQ=9)

                                                                                                       taAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGGTC‑GAATCCCCTAGGGGACGCCACTTGCTGGGTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGGGATGTTTTAGATATTTGCTCttt                                                                                                                                         >  1:75887/1‑239 (MQ=9)

                                                                                                         aagaagCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTaa                                                                                                                                       >  1:1066023/1‑239 (MQ=31)

                                                                                                          agaagCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGGTC‑GAAACCCCTAGGGGGCGCCACTTGCTGGGTTGTGAGGGAAAGTCACCCGCCTTTATATCTCAAA‑ACTTATCTTCTGGTGGTGTTTGGGATATTTGCTCTTTaaa                                                                                                                                      >  2:717536/1‑239 (MQ=11)

                                                                                                               cGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGGTC‑GAATCCCCTAGGGGGCGCCCCTTGCTGGTTTTTGAGGGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTTTTCTTTTAAAATCt                                                                                                                                 >  1:455401/1‑239 (MQ=11)

                                                                                                                           tGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTTCTGGGTTGTGAGAGAAAGTCACCTGCCTTAATATCTCAAA‑AATCATCTTCGGGTGTTGTTTTAGATATTTGATCTTTTAAAATCTGGATCAAGCTGa                                                                                                                     >  1:616645/1‑239 (MQ=11)

                                                                                                                           tGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGGTGTTTGAGATATTTGCTCTTTTAAAATTTGGATTAAACTGa                                                                                                                     >  1:744086/1‑239 (MQ=9)

                                                                                                                            gCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑CACAGGGGGTC‑GAATACCCTAGGGGCCACCCCTTGCTGGTTTCTGCGTGAAAGGTACCTGGCTTAATCTCCCACA‑AATCATCTTTGGGCGTAGTTTGTGAGATTCACTCTTTTAAAATcatgatacaaccca                                                                                                                     >  2:594348/1‑225 (MQ=11)

                                                                                                                            gCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑Ga                                                                                                                                                                                                                                             <  2:112752/119‑1 (MQ=32)

                                                                                                                            gCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑Ga                                                                                                                                                                                                                                             >  1:112752/1‑119 (MQ=32)

                                                                                                                            gCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGaa                                                                                                                    <  1:752565/239‑1 (MQ=31)

                                                                                                                            gCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGaa                                                                                                                    <  1:608844/239‑1 (MQ=31)

                                                                                                                             cTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCCCTTGCTGGTTTGTGAGTGAAAGGCACCTGCCTTAATATCTCAAA‑AATCATCTTCGGGGGATGTTTTAGATATTTGCTCTTTTAAAATCTGGATCAAGCTGaaa                                                                                                                   >  1:343577/1‑239 (MQ=9)

                                                                                                                             cTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGGTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGGCACCTGCCTTAATATCTCCAA‑ACTCATCTTCGGGTGATGTTTTTGATATTTTCTCTTTTAAAATCTTGTTCAAGCTGaaa                                                                                                                   >  1:325699/1‑239 (MQ=11)

                                                                                                                              tCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGaaaa                                                                                                                  >  1:1111742/1‑239 (MQ=31)

                                                                                                                                acacaGATTGTGTGATAGAACGTGACAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAAGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAAGGGTTC‑GCATCCCCAAGGGGCCGACACTTGCTTGGTTGTGAGTGAAAATCCAACCCCTTGTTCTCTCAAACACTCATCTT‑GGCAGATGTGTGTGGCATTTTATCTTTATAAATCTgcgtctcactggacatt                                                                                                                >  2:511933/1‑223 (MQ=2)

                                                                                                                                 cacaGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTg                                                                                                               <  1:515630/239‑1 (MQ=31)

                                                                                                                                 cacaGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTg                                                                                                               <  1:412125/239‑1 (MQ=31)

                                                                                                                                 cacaGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTg                                                                                                               <  1:344730/239‑1 (MQ=31)

                                                                                                                                 cacaGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTg                                                                                                               <  1:264438/239‑1 (MQ=31)

                                                                                                                                 cacaGATTGTCTGAGAGAAAGGGAAAAGCGAGCCTCCTTGCTAAGCAGACTGATCCGTCCACTTCGTCTAGAGGCCCAGGACACCGCCGTTGCACGGAGGT‑AACAGGGCGAC‑GAATCCCCTAGGGGACTCCGCTTGCTGGTTTGTGAGTGAAAGTCAACTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGATATAGTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTg                                                                                                               <  2:297410/239‑1 (MQ=11)

                                                                                                                                  acaGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTg                                                                                                                                               <  1:631613/206‑1 (MQ=31)

                                                                                                                                     gATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGACGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AAAAGGGGTTC‑GCAACCCCTAGGGGACGCCACTTTCTGGTTTGTGAGTGAAAGTCACCCGCCCTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTGTTCTTTTAAAAACTTGATCAAGCTGAAAATTGaaaa                                                                                                           >  2:764208/1‑238 (MQ=11)

                                                                                                                                       ttGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAaca                                                                                                          <  2:866060/238‑1 (MQ=31)

                                                                                                                                       ttGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAaca                                                                                                          <  2:183682/238‑1 (MQ=31)

                                                                                                                                          tCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGCTACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCCCGGCGGT‑TACCAGGGGTC‑TAAATCCCTACGGGACGCCACTTGCTGGTTTGTGCTTTAACGTCACCTCCCTTTATTTCTCCCA‑ACTCCTCCTTGGGGTATGTTTACGCTATTTTCTCTTTTAACATCTTTATCaaacatcgaattttaacactga                                                                                                      >  2:396984/1‑219 (MQ=2)

                                                                                                                                          tCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACCAGGGGTC‑GAATCCCCTAGGGGACGCCACTTGCGGGTTTTTGAGGGAAAATCACCTGCCTTATTATCTCCAA‑ACTCATCTTCGGGGGATGTTTTGGATATTTTCTCTTTTAAAATCTGGATAAAACTGAAAATTTAAACACTGa                                                                                                      >  2:944899/1‑239 (MQ=11)

                                                                                                                                           cTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGa                                                                                                      <  1:137137/238‑1 (MQ=31)

                                                                                                                                              aTAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTT‑GAATCCCCTAGGGGACGCCACCTGCTGGTTTGTGAGGGAAAGGCACCTGCCTTAATATTTCAAA‑AATCATCTTTGGGGGATGTTTGGGATATTTGCTCTTTAAAAATCTGGGTCAAGCTTGAAATTGAAAAAATGaacaa                                                                                                  >  1:164860/1‑239 (MQ=11)

                                                                                                                                              aTAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGC‑AACAAGGGCGC‑GAAACCCCTATGGGGAGCCCCTTGCTGGGTTGTGAGTGAAAGGCACCCGCCTTAAATTTTCCAA‑AATCAGCTTCTGGGGTTGGTTGGGATATTTTCTCTTTAAAAATCTTGATCAAACTGAAACTTTGAACACTGaaaaa                                                                                                  >  2:603651/1‑239 (MQ=11)

                                                                                                                                               tAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGGGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGaacaac                                                                                                 <  1:694246/239‑1 (MQ=14)

                                                                                                                                                                                                 cGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GGATCCCCTCGGGGGCGCCACTTGCTGGTTTGTGAGGGAAAGGAACCTGCCTTTATATCTCAAA‑ACTCATCTTCGGGTGGTGTTTTAGATATTTGCTCTTTAAAAATCTGGGTCAAGCTGAAAATTTAAAAACTGAACAACGAAAGTTGTTCGTGAGTCTCCCCAATTTTCGCAACCCGATGATGAATCGt                                               >  2:574832/1‑239 (MQ=11)

                                                                                                                                                                                                 cGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGGGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTTAAACACTGAACAACGGAAGTTGTTCGTGAGTCTCTCAAATTTTTGCCAAACGATGATGAATCGt                                               >  1:1069572/1‑239 (MQ=2)

                                                                                                                                                                                                         ggCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGaaa                                         <  1:907983/237‑1 (MQ=31)

                                                                                                                                                                                                         ggCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAAGCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGaaa                                         <  2:17118/237‑1 (MQ=14)

                                                                                                                                                                                                            ccAGGACGCCGCCCTTTCACGGCGGT‑AACAGCGGGAC‑GAATCGCCTAGGCGGCGGCGCTCGCTGGTTCGTTAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATAGTTGCTCGTTAAAAATCTGGAGCAAGCTGAAAATGGAAACACTGAACAACGAAAGTTGTACGTGAGTCTCTCAAAGTGTCGCAACACGATGATGAATCGTAAGAAACATCt                                    <  2:416868/239‑1 (MQ=11)

                                                                                                                                                                                                            ccAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTCGGGGACGCCACCTGGTGGTTTGTGAGTGAAAGTGAACTGCCTTAATATCTCAAA‑ACTCATCTTTGGCTGATGTTTGAGCTATTTTCTCTTTAAAAATCTGGAATTAGCTGAAAATTGAAACACTGTACAACGAAAAATGTTCGTGAGACTCTCAAATTTTTGCAACACGCTGATGAATTGTAATAAACATCt                                    >  1:897673/1‑239 (MQ=11)

                                                                                                                                                                                                             cAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCtt                                   <  1:586433/239‑1 (MQ=31)

                                                                                                                                                                                                             cAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAATCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAATCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCt                                    >  1:100124/1‑238 (MQ=14)

                                                                                                                                                                                                             cAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGGTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGGGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCt                                    >  1:52933/1‑238 (MQ=14)

                                                                                                                                                                                                                 accccGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTC‑g                                <  1:60368/234‑1 (MQ=14)

                                                                                                                                                                                                                     cGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTC‑GGGTtgt                          <  1:603651/239‑1 (MQ=31)

                                                                                                                                                                                                                      gCCCTTTCACGGCGGT‑AACAGGGGGTC‑GAATCCCCTAGGGGGCGCCACCTGCTGGGTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGGGGTGTTTTGGATATTTGCTCTTTTAAAATCTGGATCAAGCTGAAAATTGAAAAACTGAACAACGGAAGTTGTTCGTGAGTCTTCCACATTTTTGTAACACAATGATGAATCGTAAGAAAAATCTTC‑GGGAtgtg                         >  2:197985/1‑239 (MQ=11)

                                                                                                                                                                                                                          tttCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCCCTTGCTGGTTTTTGAGGGAAAATCACCTGCCTTTATATTTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTTAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAATTGTTTGTGAGTCTTTCAAAATTTTCCAAAACGAATATTAAACGTAAGAAACATCTTCGGGGTTGTgggg                      >  2:609253/1‑239 (MQ=11)

                                                                                                                                                                                                                          tttCACGGCGGT‑AACAGGGGGTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGGCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGGGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAATTTTTCGTGAGTCTCTCAAATTTTCGCAACCCGATGATGAATCGTAAAAAACATCTTC‑CGGATGTGAGGt                     >  1:308195/1‑239 (MQ=11)

                                                                                                                                                                                                                           ttaaCGGCTGT‑AACAGGGGTTC‑GAATCCCCTAGGAGACGCCTCTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTC‑GGGTTGTGAGGtt                    <  2:522179/236‑1 (MQ=9)

                                                                                                                                                                                                                                  cggT‑AACAGGGGGTC‑GAATCCCCTAGGGGACGCCCCTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTTCGATATTTGCTCTTTAAAAATCTGGATCAAACTGAAAATTGGAACACTGAACAACGAAAGGTGTTCGGGAGTTTCTCAAATTTTCTCCACACGATGATGAATCGTAAGAAACATCTTT‑GGGTTGTGAGGTTAAACg               >  2:399312/1‑237 (MQ=11)

                                                                                                                                                                                                                                     t‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTC‑GGGTTGTGAGGTTAAGCGACTa           <  1:574832/238‑1 (MQ=31)

                                                                                                                                                                                                                                       aaCAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTC‑GGGTTGTGAGGTTAAGCGACTaa          <  1:209285/238‑1 (MQ=31)

                                                                                                                                                                                                                                            ggtttC‑GAATCCCCTAAGGGATGCCACTTGGTGGTTGGTGAGTGAAAGGAACGTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTGCTATTTAAAAATCTGTATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAAATTTCGCAACACGATGATGAATCGTAAGAAGCATCGTC‑GGGTTGTGAGGTTAAGCGACTAAGCGTa     <  2:343577/235‑1 (MQ=11)

                                                                                                                                                                                                                                             ggTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTC‑GGGTTGTGAGGTTAAGCGACTAAGCGTa     <  2:596896/237‑1 (MQ=31)

                                                                                                                                                                                                                                              gTTC‑GAATCCCCTAGGGGATGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTC‑GGGTTGTGAGGTTAAGCGACTAAGCGTaca   <  2:455401/238‑1 (MQ=14)

                                                                                                                                                                                                                                              gTTC‑GAATCCCCTAGGGGACGGCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTC‑GGGTTGTGAGGTTAAGCGACTAAGCGTaca   <  2:325699/238‑1 (MQ=14)

                                                                                                                                                                                                                                              gTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTC‑GGGTTGTGAGGTTAAGCGACTAAGCGTaca   <  1:420907/238‑1 (MQ=31)

CTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGT‑AACAGGGGTTC‑GAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAA‑ACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTC‑GGGTTGTGAGGTTAAGCGACTAAGCGTACAC  >  NC_000913/4209588‑4210063

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