BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,127,252	A→G	6.4%	intergenic (+27/+34)	rpmE → / ← yiiX	50S ribosomal subunit protein L31/putative lipid binding hydrolase, DUF830 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,127,252	0	A	G	6.4%	702.0 / 25.7	220	intergenic (+27/+34)	rpmE/yiiX	50S ribosomal subunit protein L31/putative lipid binding hydrolase, DUF830 family protein
Reads supporting (aligned to +/- strand): ref base A (98/107); new base G (7/7); total (105/115)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-02

Only 100 of 224 total aligned reads displayed.

GGCTGGAAGAGCGACGCGGCCTTAAACCGAGGTTTTCCCATGAAAAAAGATATTCACCCGAAATACGAAGAAATTACTGCTAGCTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCA‑GCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTT‑GTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTA  >  NC_000913/4126974‑4127516

ggCTGGAAGAGCGCCGCGGCCTTAAACCGAGGTTTTCCCATGAAAAAAGATATTCACCCGAAATACGAAGAAATTACTGCCAGCGGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGCTCATGACCTGAACCTCGACGTGTGCA‑GCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCa                                                                                                                                                                                                                                                                           <  1:187391/279‑1 (MQ=255)

    ggAAGAGCGACGCGGCCTTAAACCGAGGTTTTCCCATGAAAAAAGATATTCACCCGAAATACGAAGAAATTACTGCTAGCTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCA‑GCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAACTAAGACCGGATTTCCGAAAAAAGCGCCGCa                                                                                                                                                                                                                                                                           >  2:673980/1‑275 (MQ=255)

                                        tGAAAAAAGATATTCACCCGAAATACGAAGAAATTACTGCTAGATGATCTTGCGGTAACGTAATGAAAATCCGCTACACCGTTGGTCATGACCTGAACCTCGACGTGTGCA‑GCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCAGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGAACGGCGCTTTTTTTTGGCCTGAATGTTGGTAGTAATa                                                                                                                                                                                                                                       >  2:559481/1‑275 (MQ=255)

                                                       cccagAAATACGAAGAAATTACTGCTAGCTTCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCA‑GCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCAt                                                                                                                                                                                                                    <  1:131602/275‑1 (MQ=255)

                                                          acagaTACGAAGAAATTACTGCTAGCTGGTCTTGCGGTAACGTAATGAAAATCCGCGCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCA‑GCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGAAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGcc                                                                                                                                                                                                                      <  2:409029/270‑1 (MQ=255)

                                                                                                aaCGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCA‑GCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTAc                                                                                                                                                                                                                                      <  2:377135/220‑1 (MQ=255)

                                                                                                            aTCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCA‑GCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGc                                                                                                                                                                    <  2:80755/274‑1 (MQ=255)

                                                                                                                gCTCCACCGGTGGTCATGAGCTGAGACTTG‑CGGGTGGACGCAATGGCCACCCTTTCTTCACTGTCAAAAAGCGTGATGTTGCTACCGGTGTCCGTGTTGACCGCTTTAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAGAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGGCGTATTACCTCGAAAACAGCAGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGac                                                                                                                                                               <  4:563080/275‑1 (MQ=255)

                                                                                                                cgcctgccGTTGGTCATGACCTGATCCTGGACGTGGGCA‑GCAAGTCCCACCCGTTCTTCCCTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGac                                                                                                                                                               <  4:205935/269‑1 (MQ=255)

                                                                                                                                   ccTGAACCTCGACGTGTGCA‑GCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCAc                                                                                                                                                                                                                                                                          >  5:778211/1‑149 (MQ=255)

                                                                                                                                       aaCCTCGACGTGTGCA‑GCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCgtggt                                                                                                                                                                                                   <  1:114459/216‑1 (MQ=255)

                                                                                                                                         ccTCGACGTGTGCA‑GCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCt                                                                                                                                                                                                                                                                    >  7:244223/1‑149 (MQ=255)

                                                                                                                                           tCGACGTGTGCA‑GCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCttt                                                                                                                                                                                                                                                                  >  5:1527471/1‑149 (MQ=255)

                                                                                                                                             ggcGTGGGCA‑GGAAGTTCCACCCGTTCTTCAATGGCAAAAAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGt                                                                                                                               <  1:27707/276‑1 (MQ=255)

                                                                                                                                             gACGTGTGCA‑GCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCttttt                                                                                                                                                                                                                                                                >  6:1560497/1‑149 (MQ=255)

                                                                                                                                                gtgtGCA‑GCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCAAAAAAAGCGCCGTGCGGCGCttttttcg                                                                                                                                                                                                                                                             <  7:761440/149‑3 (MQ=255)

                                                                                                                                                 tgtgCA‑GCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTtgt                                                                                                                                                                                                                                                            <  5:1471600/149‑1 (MQ=255)

                                                                                                                                                 tgtgCA‑GCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTtgt                                                                                                                                                                                                                                                            <  5:1185867/149‑1 (MQ=255)

                                                                                                                                                         caAGTGCCGCCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGaa                                                                                                                                                                                                                                                     <  5:1211048/149‑1 (MQ=255)

                                                                                                                                                              gCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTg                                                                                                                                                                                                                                                <  7:484065/149‑1 (MQ=255)

                                                                                                                                                                 aCCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTc                                                                                                                                                                                                                                             <  5:995330/149‑1 (MQ=255)

                                                                                                                                                                    cGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTa                                                                                                                                                                                                                                          <  8:157063/149‑1 (MQ=255)

                                                                                                                                                                     gTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTAt                                                                                                                                                                                                                                         <  6:1502301/149‑1 (MQ=255)

                                                                                                                                                                     gTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTAt                                                                                                                                                                                                                                         >  6:861903/1‑149 (MQ=255)

                                                                                                                                                                       tcttcACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAACAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCCAAAGCGGACTATTCTTGCCGAACCGTTCTTTTAGCTTTACCt                                                                                                          >  4:178187/1‑274 (MQ=255)

                                                                                                                                                                       tcttcACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGCGTTGTGGCGCATCGAAAACGGCTTGTGGCAAGCCCACCGTTTCCTCAAGCGGAAGATTCTTCCCGTAACGTTCTTTGATCttt                                                                                                              >  2:620634/1‑270 (MQ=255)

                                                                                                                                                                        cttcACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTg                                                                                                         <  2:525982/274‑1 (MQ=255)

                                                                                                                                                                         ttcACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGATCGTATTAcc                                                                                                                                                                                                                                     >  8:1502724/1‑149 (MQ=255)

                                                                                                                                                                           cACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCTGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCAc                                                                                                                                                                                                                                   <  6:366054/149‑1 (MQ=255)

                                                                                                                                                                            aCTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACg                                                                                                                                                                                                                                  >  6:1063686/1‑149 (MQ=255)

                                                                                                                                                                              tGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGt                                                                                                                                                                                                                                              >  8:932274/1‑135 (MQ=255)

                                                                                                                                                                              tGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAAAAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGt                                                                                                                                                                                                                                              <  7:1254188/135‑1 (MQ=255)

                                                                                                                                                                              tGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAAAAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGt                                                                                                                                                                                                                                              >  8:1254188/1‑135 (MQ=255)

                                                                                                                                                                                gCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACTGTCTTTGACCTTTGCCTTAACTT‑GTgg                                                                                               >  4:504145/1‑275 (MQ=255)

                                                                                                                                                                                    aCAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGc                                                                                                                                                                                                                          >  7:461198/1‑149 (MQ=255)

                                                                                                                                                                                      aGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTT‑GTGGgctgct                                                                                         <  2:205085/275‑1 (MQ=255)

                                                                                                                                                                                        cGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGa                                                                                                                                                                                                                                 >  4:710115/1‑138 (MQ=255)

                                                                                                                                                                                        cGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTg                                                                                                                    >  2:152563/1‑247 (MQ=255)

                                                                                                                                                                                          tGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCAt                                                                                                                                                                                                                    >  7:1474488/1‑149 (MQ=255)

                                                                                                                                                                                            aTGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTc                                                                                                                                                                                                                  >  7:1623327/1‑149 (MQ=255)

                                                                                                                                                                                               ttGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAAtttt                                                                                                                                                                                                                                                 <  6:1642890/115‑1 (MQ=255)

                                                                                                                                                                                               ttGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGc                                                                                                                                                                                                                       <  5:424432/141‑1 (MQ=255)

                                                                                                                                                                                               ttGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCttt                                                                                                                                                                                                               <  5:1128842/149‑1 (MQ=255)

                                                                                                                                                                                               ttGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCttt                                                                                                                                                                                                               <  8:881680/149‑1 (MQ=255)

                                                                                                                                                                                                tGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGTCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTg                                                                                                                                                                                                              <  7:1116140/149‑1 (MQ=255)

                                                                                                                                                                                                tGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCATATTACCACGAAAACAGCGGCCATTCTTTg                                                                                                                                                                                                              <  6:128943/149‑1 (MQ=255)

                                                                                                                                                                                                 gCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAAtttt                                                                                                                                                                                                                                                 >  6:24190/1‑113 (MQ=255)

                                                                                                                                                                                                 gCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTgg                                                                                                                                                                                                             >  8:1315425/1‑149 (MQ=255)

                                                                                                                                                                                                  cTACCGGTGGCCTTTTTGACCGCTTCAACATGCGTTTCAACATCCCGGGCAGCAAATAAGACCCGATTTCCCAACAAAGCGCCGCACGTCGCTTTCTTTCTACCCCAACTTTTGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTCGAGAGTTGTGGCGCATCGAACACGGATTGTGGCCAGACCACCGTTTCCTCATGCGGAACCACCGTGCCGTACCGTTCTTTACCCTTTGCCTGACCTT‑TTTCgctgcccgcggcgcact                                                                              >  4:54790/1‑262 (MQ=255)

                                                                                                                                                                                                       ggTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCg                                                                                                                                                                                                       >  8:914841/1‑149 (MQ=255)

                                                                                                                                                                                                          ggCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCgtgg                                                                                                                                                                                                    <  5:260572/149‑1 (MQ=255)

                                                                                                                                                                                                             cGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGa                                                                                                                                                                                                 >  6:556337/1‑149 (MQ=255)

                                                                                                                                                                                                                   gACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTtgtg                                                                                                                                                                                           <  5:1655898/149‑1 (MQ=255)

                                                                                                                                                                                                                    aCCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTgg                                                                                                                                                                                          >  5:840728/1‑149 (MQ=255)

                                                                                                                                                                                                                    aCCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTgg                                                                                                                                                                                          >  6:201617/1‑149 (MQ=255)

                                                                                                                                                                                                                    aCCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTgg                                                                                                                                                                                          >  5:1691619/1‑149 (MQ=255)

                                                                                                                                                                                                                    aCCGCCTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTCTTTTTGTGCCGGAATTTTGGTCGTATTACCACGACAACAGCGGCCGTTCTTTTGCGCGCGTGNTGAGTTGTCGCGCATCGAAAACGGCTTATCGCGAGACTACCGTGTCCTCAGCGGGAATCTGCTTGCGGTTACGTTCTTTTNGTGTTGCCTCAACTTCGGGCCCCGGNACGGTCG‑ACTCTTTCAATTTCtgctgccccgc                                                        >  3:173393/1‑273 (MQ=255)

                                                                                                                                                                                                                      cGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGgcg                                                                                                                                                                                        >  7:573398/1‑149 (MQ=255)

                                                                                                                                                                                                                         ttCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCAt                                                                                                                                                                                     <  8:1623327/149‑1 (MQ=255)

                                                                                                                                                                                                                         ttCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTT‑GTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTgcccgccc                                                       >  2:569818/1‑274 (MQ=255)

                                                                                                                                                                                                                         ttCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGGTCTTTGAGCtt                                                                                                               >  2:305868/1‑219 (MQ=255)

                                                                                                                                                                                                                            aacaaGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGa                                                                                                                                                                                  >  7:1272802/1‑149 (MQ=255)

                                                                                                                                                                                                                             aaaaacGTTTCGACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCCCGGCGCTTTTTTTTTGCCTGGATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGcc                                                                                                           <  4:305864/214‑1 (MQ=255)

                                                                                                                                                                                                                              caaGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGaaa                                                                                                                                                                                >  8:1026837/1‑149 (MQ=255)

                                                                                                                                                                                                                              caaGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTT‑GTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACgcgc                                                 <  2:178190/275‑1 (MQ=255)

                                                                                                                                                                                                                               aaGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTT‑GTGGgc                                                                                             <  4:377465/230‑1 (MQ=255)

                                                                                                                                                                                                                                aGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAAc                                                                                                                                                                              <  8:225176/149‑1 (MQ=255)

                                                                                                                                                                                                                                 gCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTc                                                                                                                                                                                                                  >  5:1600081/1‑112 (MQ=255)

                                                                                                                                                                                                                                  cGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACgg                                                                                                                                                                            <  8:1272802/149‑1 (MQ=255)

                                                                                                                                                                                                                                   gTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCttt                                                                                                                                                                                                               <  7:850509/113‑1 (MQ=255)

                                                                                                                                                                                                                                    tttCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCt                                                                                                                                                                          <  7:1026837/149‑1 (MQ=255)

                                                                                                                                                                                                                                    tttCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCt                                                                                                                                                                          >  7:614198/1‑149 (MQ=255)

                                                                                                                                                                                                                                      tCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTtg                                                                                                                                                                        <  8:461198/149‑1 (MQ=255)

                                                                                                                                                                                                                                      tCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTtg                                                                                                                                                                        >  8:624672/1‑149 (MQ=255)

                                                                                                                                                                                                                                       cAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTtgt                                                                                                                                                                       <  6:413129/149‑1 (MQ=255)

                                                                                                                                                                                                                                       cAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTtgt                                                                                                                                                                       <  6:778211/149‑1 (MQ=255)

                                                                                                                                                                                                                                       cAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTtgt                                                                                                                                                                       <  8:1474488/149‑1 (MQ=255)

                                                                                                                                                                                                                                         aCATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCt                                                                                                                                                                                                                 <  5:546401/105‑1 (MQ=255)

                                                                                                                                                                                                                                         aCATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTtgt                                                                                                                                                                       >  5:233556/1‑147 (MQ=255)

                                                                                                                                                                                                                                          cATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGc                                                                                                                                                                    <  7:665596/149‑1 (MQ=255)

                                                                                                                                                                                                                                            tCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTAt                                                                                                                                                                                                                                         >  8:1134844/1‑78 (MQ=255)

                                                                                                                                                                                                                                            tCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCga                                                                                                                                                                  >  6:505637/1‑149 (MQ=255)

                                                                                                                                                                                                                                            tCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCga                                                                                                                                                                  <  8:820984/149‑1 (MQ=255)

                                                                                                                                                                                                                                             cccGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTT‑GTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTCCTCGCCCCAGCGCATTCCCAGCGCGtt                                   >  4:55716/1‑274 (MQ=255)

                                                                                                                                                                                                                                                 ggCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGtc                                                                                                                                                                                                                                          >  4:458941/1‑71 (MQ=255)

                                                                                                                                                                                                                                                     gcaAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTAc                                                                                                                                                                                                                                      <  2:458945/72‑1 (MQ=255)

                                                                                                                                                                                                                                                      caAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCgtggt                                                                                                                                                                                                   <  1:65396/106‑1 (MQ=255)

                                                                                                                                                                                                                                                      caAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCgtggt                                                                                                                                                                                                   >  3:65396/1‑106 (MQ=255)

                                                                                                                                                                                                                                                       aaaTAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGttt                                                                                                                                                       <  5:1382345/149‑1 (MQ=255)

                                                                                                                                                                                                                                                         aTAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTcc                                                                                                                                                     <  7:19997/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                  gATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCgg                                                                                                                                             <  3:116394/148‑1 (MQ=255)

                                                                                                                                                                                                                                                                    tttCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAat                                                                                                                                          >  6:464087/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                       acaaaaaaaGCGCCGTGCGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCga                                                                                                                                                                  >  7:918284/4‑122 (MQ=255)

                                                                                                                                                                                                                                                                         gAAAAAAGCGCCGCACGGCGCTTTTTTCGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCgtg                                                                                                                                                                                                     <  6:533394/85‑1 (MQ=255)

                                                                                                                                                                                                                                                                          aaaaaaGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGaa                                                                                                                                           <  7:949327/142‑1 (MQ=255)

                                                                                                                                                                                                                                                                           aaaaaGCGCCGTGCGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTg                                                                                                                                   >  6:1452695/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                           aaaaaGCGCCGTGCGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTAGCGACCGTGGTGGGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTGTCCTCAAGCGGAATATTCTTg                                                                                                                                   >  6:724554/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                            aaaaGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTCGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGc                                                                                                                                  >  5:1408247/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                             aaaGCGCCGTGCGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGcc                                                                                                                                 >  8:460667/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                                       aCGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCt                                                                                                                       >  5:1587956/1‑149 (MQ=255)

GGCTGGAAGAGCGACGCGGCCTTAAACCGAGGTTTTCCCATGAAAAAAGATATTCACCCGAAATACGAAGAAATTACTGCTAGCTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCA‑GCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTT‑GTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTA  >  NC_000913/4126974‑4127516

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 10 ≤` `ATCG/ATCG` `< 22 ≤` `ATCG/ATCG` `< 35 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`