BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,425,968	T→G	5.2%	intergenic (‑52/+38)	hisQ ← / ← hisJ	histidine ABC transporter permease/histidine ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,425,968	0	T	G	5.2%	694.6 / 25.3	199	intergenic (‑52/+38)	hisQ/hisJ	histidine ABC transporter permease/histidine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand): ref base T (102/85); new base G (7/4); total (109/90)
Fisher's exact test for biased strand distribution p-value = 7.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.35e-01

Only 100 of 202 total aligned reads displayed.

CAATCTGCAAACCGTAGAAAATCAGCAGCATCAACACTAAGTCCGGCACGCCACGAATCAGCGTGGTGTACCCTTCGAAAATCAGCCCCGAAAGCCGATTTTGCGAGAGCTTACCGCCAGCGCCAATTAAACCGATGATTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCA  >  NC_000913/2425694‑2426240

cAATCTGCAAACCGTAGAAAATCAGCAGCATCAACACTAAGTCCGGCACGCCACGAATCAGCGTGGTGTACCCTTCGAAAATCAGCCCCGAAAGCCGATTTTGCGAGAGCTTACCGCCAGCGCCAATTAAACCGATGATTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGt                                                                                                                                                                                                                                                                                  <  2:152306/275‑1 (MQ=255)

                                          ccGGCACGCCACGAATCAGCGTGGTGTACCCTTCGAAAATCAGCCCCGAAAGCCGATTTTGCGAGAGCTTACCGCCAGCGCCAATTAAACCGATGATTACAGCGAGCACTACAGCGCTGATAGCCAGCCCCAGCGTGACGAGCGCACCCTGTAAACTAACACCTGAAAACCCATCCAACATGCTGCCTGTCCTGTCGTGTGTGGCTAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGGAAGGCGCGGGCTTTCTTTCACCGATCGCCATTAGc                                                                                                                                                                                                                                        >  4:219041/1‑275 (MQ=255)

                                                                                     ccccGAAAGCCGATTTTGCGAGAGCTTACCGCCAGCGCCAATTAAACCGATGATTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTCGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTCGCCACCATAAACATCAAAATCGAAGTACCTTTTCGCTACTTtct                                                                                                                                                                                              >  1:482/1‑274 (MQ=255)

                                                                                                                               tAAACCGATGATTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGtt                                                                                                                                                                                                                                                                                 <  5:165977/149‑1 (MQ=255)

                                                                                                                                aaaCCGATGATTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTa                                                                                                                                                                                                                                                                                <  7:738742/149‑1 (MQ=255)

                                                                                                                                aaaCCGATGATTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTa                                                                                                                                                                                                                                                                                <  7:612944/149‑1 (MQ=255)

                                                                                                                                    cGATGATTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCg                                                                                                                                                                                                                                                                            >  8:109207/1‑149 (MQ=255)

                                                                                                                                      atgatTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGaa                                                                                                                                                                                                                                                                          <  7:1374433/149‑1 (MQ=255)

                                                                                                                                       tgatTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAg                                                                                                                                                                                                                                                                         <  8:961659/149‑1 (MQ=255)

                                                                                                                                          ttACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGa                                                                                                                                                                                                                                                                      >  5:647642/1‑149 (MQ=255)

                                                                                                                                            aCAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGc                                                                                                                                                                                                                                                                    <  8:776769/149‑1 (MQ=255)

                                                                                                                                               gCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCggg                                                                                                                                                                                                                                                                 <  6:1044339/149‑1 (MQ=255)

                                                                                                                                                    cACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTAAATTACCGCCTCTTGTCTCCCTCCGGTTNCCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGCTTAGCCACCATAAACATCAAAATCGAAGTACTTTNTCGCTAATTTCCCGTAAGTAACGTCAGCGCGCATTTCGGCAAAGGNTNTTTTAAGTGCTTCGAGCAGTTCGtttt                                                                                                                               >  3:34126/1‑272 (MQ=255)

                                                                                                                                                     aCTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttct                                                                                                                                                                                                                                                           <  8:718093/149‑1 (MQ=255)

                                                                                                                                                     aCTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttct                                                                                                                                                                                                                                                           <  6:413471/149‑1 (MQ=255)

                                                                                                                                                     aCTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttct                                                                                                                                                                                                                                                           <  7:511808/149‑1 (MQ=255)

                                                                                                                                                        aCAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCtttctttc                                                                                                                                                                                                                                                        >  6:214145/1‑149 (MQ=255)

                                                                                                                                                             gCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGa                                                                                                                                                                                                                                                   >  8:344773/1‑149 (MQ=255)

                                                                                                                                                              cTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGAt                                                                                                                                                                                                                                                  <  5:459847/149‑1 (MQ=255)

                                                                                                                                                               tGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATg                                                                                                                                                                                                                                                 <  5:1382921/149‑1 (MQ=255)

                                                                                                                                                                gATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATgg                                                                                                                                                                                                                                                >  8:1271988/1‑149 (MQ=255)

                                                                                                                                                                gATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATgg                                                                                                                                                                                                                                                >  6:1434954/1‑149 (MQ=255)

                                                                                                                                                                 aTAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAACTCGAAGTAATTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCCCATTTCGTCACAGGCTTTGTTCAGTGCTTCGCACCGTTCGTTCTCTTCTTTACGCAg                                                                                                                  >  4:81960/1‑274 (MQ=255)

                                                                                                                                                                 aTAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGATTTCTTTCACCGATTGCCACTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTCCCGTCAGCGCGCATTTCGGCAACGGCTTTGTTCATTGCTTCGAGCAGTTCTTTATCTTCTTTACGCAg                                                                                                                  >  2:283109/1‑274 (MQ=255)

                                                                                                                                                                   aGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGa                                                                                                                                                                                                                                             <  7:763397/149‑1 (MQ=255)

                                                                                                                                                                     ccAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGTGGGCTTTCTTTCACCGATGGCGAtt                                                                                                                                                                                                                                           >  5:580450/1‑149 (MQ=255)

                                                                                                                                                                     ccAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGAtt                                                                                                                                                                                                                                           >  8:497947/1‑149 (MQ=255)

                                                                                                                                                                      cAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTa                                                                                                                                                                                                                                          >  5:1230335/1‑149 (MQ=255)

                                                                                                                                                                      cAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTa                                                                                                                                                                                                                                          >  7:554925/1‑149 (MQ=255)

                                                                                                                                                                          tCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGcca                                                                                                                                                                                                                                      <  6:270774/149‑1 (MQ=255)

                                                                                                                                                                          tCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGcca                                                                                                                                                                                                                                      >  5:587562/1‑149 (MQ=255)

                                                                                                                                                                          tCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGcca                                                                                                                                                                                                                                      >  5:758241/1‑149 (MQ=255)

                                                                                                                                                                                 tGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTa                                                                                                 >  2:62736/1‑275 (MQ=255)

                                                                                                                                                                                     gAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATc                                                                                                                                                                                                                           >  8:56909/1‑149 (MQ=255)

                                                                                                                                                                                             ccTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGa                                                                                                                                                                                                                   >  6:750083/1‑149 (MQ=255)

                                                                                                                                                                                              cTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                                                                                                                                                  <  7:459589/149‑1 (MQ=255)

                                                                                                                                                                                              cTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                                                                                                                                                  <  6:1121465/149‑1 (MQ=255)

                                                                                                                                                                                              cTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                                                                                                                                                  <  6:77981/149‑1 (MQ=255)

                                                                                                                                                                                               tGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctt                                                                                                                                                                                                                                                              >  7:959139/1‑104 (MQ=255)

                                                                                                                                                                                                gTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGt                                                                                                                                                                                                                >  7:1357297/1‑149 (MQ=255)

                                                                                                                                                                                                   aaaTAACACCTGAAAACCCAAACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTAc                                                                                                                                                                                                              <  8:1374288/148‑1 (MQ=255)

                                                                                                                                                                                                    aaTAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACtt                                                                                                                                                                                                            <  7:534173/149‑1 (MQ=255)

                                                                                                                                                                                                    aaTAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGCGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACtt                                                                                                                                                                                                            >  5:258394/1‑149 (MQ=255)

                                                                                                                                                                                                         cacCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCCAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACAGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGTCCCATTCCGGTCCCTACCCCAAACCGTTTTTCATCtttc                                                                          >  4:232286/1‑273 (MQ=255)

                                                                                                                                                                                                            cTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAggg                                                                                                                                                                                                                                                                       >  5:262356/1‑82 (MQ=255)

                                                                                                                                                                                                            cTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                                                                                                                                                  <  6:1293783/135‑1 (MQ=255)

                                                                                                                                                                                                            cTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                                                                                                                                                  >  5:1293783/1‑135 (MQ=255)

                                                                                                                                                                                                             tGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTaa                                                                                                                                                                                                   >  7:977406/1‑149 (MQ=255)

                                                                                                                                                                                                             tGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTaa                                                                                                                                                                                                   >  5:809409/1‑149 (MQ=255)

                                                                                                                                                                                                                  aCCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATg                                                                                                                                                                                                                                                 <  7:58211/98‑1 (MQ=255)

                                                                                                                                                                                                                    ccATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGAt                                                                                                                                                                                                                                                  <  6:933062/95‑1 (MQ=255)

                                                                                                                                                                                                                    ccATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGAt                                                                                                                                                                                                                                                  >  5:933062/1‑95 (MQ=255)

                                                                                                                                                                                                                    ccATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAAc                                                                                                                                                                                                                              >  5:850506/1‑115 (MQ=255)

                                                                                                                                                                                                                    ccATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCg                                                                                                                                                                                            >  8:1004632/1‑149 (MQ=255)

                                                                                                                                                                                                                     cATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                                                                                                                                                  <  6:472356/126‑1 (MQ=255)

                                                                                                                                                                                                                       tACAACCTGCTGCCTGTCCTGTCGTGTGTGGTAAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCg                                                           <  1:34126/275‑1 (MQ=255)

                                                                                                                                                                                                                          aacaTGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTa                                                                                                                                                                                      >  8:682931/1‑149 (MQ=255)

                                                                                                                                                                                                                                  gcCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAgc                                                                                                                                                                              <  8:554925/149‑1 (MQ=255)

                                                                                                                                                                                                                                  gcCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTGTCATCTTTACCGAACGGGCCCCCGGATTTggactc                                                >  2:305253/1‑270 (MQ=255)

                                                                                                                                                                                                                                  gcCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATC‑AAATCGaa                                                                                                                                                                                                                  <  8:1154481/112‑1 (MQ=255)

                                                                                                                                                                                                                                  gcCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATC‑AAATCGaa                                                                                                                                                                                                                  >  7:1154481/1‑112 (MQ=255)

                                                                                                                                                                                                                                       tcctgtcGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCAGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCa                                                                                                                                                                         >  5:1165006/1‑149 (MQ=255)

                                                                                                                                                                                                                                         ctgtcGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCAtt                                                                                                                                                                       >  7:1417664/1‑149 (MQ=255)

                                                                                                                                                                                                                                           gtcGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTc                                                                                                                                                                     <  8:281896/149‑1 (MQ=255)

                                                                                                                                                                                                                                           gtcGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTc                                                                                                                                                                     <  7:210911/149‑1 (MQ=255)

                                                                                                                                                                                                                                            tcGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCg                                                                                                                                                                    <  8:612986/149‑1 (MQ=255)

                                                                                                                                                                                                                                             cGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCgg                                                                                                                                                                   <  5:992438/149‑1 (MQ=255)

                                                                                                                                                                                                                                                 tgtgGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGGATTCTTTACCGACAGtttg                                 >  4:109772/1‑275 (MQ=255)

                                                                                                                                                                                                                                                 tgtgGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCAAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCCAAAGgtt                                  >  4:302405/1‑274 (MQ=255)

                                                                                                                                                                                                                                                 tgtgGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCAGGCCGGTCCCTACGCCAAACCGTTTTTCATCTTTAAACACCGGGCCACGGACTTTTGAATCtttcacaaaaggtt                                  >  4:243575/1‑263 (MQ=255)

                                                                                                                                                                                                                                                  gtggtgAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAg                                                                                                                                                              <  5:214145/149‑1 (MQ=255)

                                                                                                                                                                                                                                                   tggtgAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAgg                                                                                                                                                             >  7:333703/1‑149 (MQ=255)

                                                                                                                                                                                                                                                   tggtgAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAgg                                                                                                                                                             >  5:778538/1‑149 (MQ=255)

                                                                                                                                                                                                                                                     gtgAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAg                                                                                                                                                              <  8:1139083/146‑1 (MQ=255)

                                                                                                                                                                                                                                                     gtgAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAg                                                                                                                                                              >  7:1139083/1‑146 (MQ=255)

                                                                                                                                                                                                                                                     gtgAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCt                                                                                                                                                           <  5:791730/149‑1 (MQ=255)

                                                                                                                                                                                                                                                       gAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCttt                                                                                                                                                         <  5:662694/149‑1 (MQ=255)

                                                                                                                                                                                                                                                        aaTTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCgg                                                                                                                                                                   >  6:647611/1‑138 (MQ=255)

                                                                                                                                                                                                                                                         aTTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGt                                                                                                                                                       <  8:225100/149‑1 (MQ=255)

                                                                                                                                                                                                                                                          ttACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGtt                                                                                                                                                      >  5:1313981/1‑149 (MQ=255)

                                                                                                                                                                                                                                                          ttACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGtt                                                                                                                                                      >  6:669982/1‑149 (MQ=255)

                                                                                                                                                                                                                                                           tACCGCCTCTTGGCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACAATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGGACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTGAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGaaa                       <  2:219041/275‑1 (MQ=255)

                                                                                                                                                                                                                                                             ccGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAg                                                                                                                                                   <  6:1357732/149‑1 (MQ=255)

                                                                                                                                                                                                                                                              cGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTtctc                                                                                                                                                                                             >  7:1049388/1‑106 (MQ=255)

                                                                                                                                                                                                                                                              cGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGt                                                                                                                                                  >  6:1230534/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                 tcgtttttcccccccGTTTACCNNAAGGTNGTGGGCTTTTTTTCACCGTTGGTGATTATCCGCCATAAACAGCAGAATCTAAGGACTTTTTCGCTAATTTCTCGTAAGTCCCTTCATCGCGCTTTTCGGCAAGCGCTTTGTTCAGTGCTTCCCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCg                 <  3:482/266‑1 (MQ=255)

                                                                                                                                                                                                                                                                 ctctTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCt                                                                                                                                               >  7:1010720/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                   ctTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTc                                                                                                                                             >  6:812818/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                     ccgctcccttcGGGTAACCGGAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCg                                                                                                                                                                                            <  7:567796/97‑1 (MQ=38)

                                                                                                                                                                                                                                                                     ccgctcccttcGGGTAACCGGAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTcgc                                                                                                                                           >  7:616807/4‑149 (MQ=255)

                                                                                                                                                                                                                                                                      cgctcccttcGGGTAACCGGAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAAtt                                                                                                                                                                                                 <  6:1326271/92‑1 (MQ=38)

                                                                                                                                                                                                                                                                      cgctcccttcGGGTAACCGGAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTcgcg                                                                                                                                          >  7:1203352/3‑149 (MQ=255)

                                                                                                                                                                                                                                                                         tcccttcGGGTAACCGGAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGt                                                                                                                                                       >  6:61710/1‑133 (MQ=255)

                                                                                                                                                                                                                                                                         tcccttcGGGTAACCGGAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGt                                                                                                                                                       <  5:61710/133‑1 (MQ=255)

                                                                                                                                                                                                                                                                           cctccGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGtt                                                                                                                                     >  7:482250/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                             tccGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGAcc     <  2:109772/275‑1 (MQ=255)

                                                                                                                                                                                                                                                                             tccGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGAcc     <  4:62736/275‑1 (MQ=255)

                                                                                                                                                                                                                                                                              gcGGTTACCCGAAGGGAGCAGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCCAAAGCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGTCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGt                                                                                                                                  <  6:1058829/148‑1 (MQ=255)

                                                                                                                                                                                                                                                                              ccGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGt                                                                                                                                  <  6:84119/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                                 gTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATACACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCCTAAGTACCGTCAGCGCGCATTTCGGCAATGGCTTTGTTCAGTGCTTCGCGCTGTTCGTTATCTTCTTTACGCATGCCCATTCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACATCCAGGCCCCCGAAGTTGTACTCTTTACAGACCGGTTGTTTGAGGCAACATTCCCTGGcgaccacctga  >  4:185861/1‑264 (MQ=255)

CAATCTGCAAACCGTAGAAAATCAGCAGCATCAACACTAAGTCCGGCACGCCACGAATCAGCGTGGTGTACCCTTCGAAAATCAGCCCCGAAAGCCGATTTTGCGAGAGCTTACCGCCAGCGCCAATTAAACCGATGATTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCA  >  NC_000913/2425694‑2426240

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