BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,177,172	T→A	6.6%	intergenic (+44/‑186)	tufB → / → secE	translation elongation factor EF‑Tu 2/preprotein translocase membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,177,172	0	T	A	6.6%	895.0 / 39.0	257	intergenic (+44/‑186)	tufB/secE	translation elongation factor EF‑Tu 2/preprotein translocase membrane subunit
Reads supporting (aligned to +/- strand): ref base T (119/121); new base A (7/10); total (126/131)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.36e-03

Only 100 of 264 total aligned reads displayed.

CGTTCTTCAAAGGCTACCGTCCGCAGTTCTACTTCCGTACTACTGACGTGACTGGTACCATCGAACTGCCGGAAGGCGTAGAGATGGTAATGCCGGGCGACAACATCAAAATGGTTGTTACCCTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCGGGGTGAAAATGTTTGTAGAAAACTTCTGACAGGTTGGTTTATGAGTGCGAATACCGAAGCTCAAGGAAGCGGGCGCGGCCTGGAAGCGATGAAGTGGGTCGTTGTGGTGGCA  >  NC_000913/4176908‑4177429

                                    gTACTACTGACGTGACTGGTACCATCGAACTGCCGGAAGGCGTAGAGATGGTAATGCCGGGCGTCAACATCAAAATGGTTGTTACCCTGATCCACCAGATCGCGATGGACGACGGTCTGCGTTTCGCAGTCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATACCCGATAACATTTGCCGCAATGCGCACTAAACGGGCATCATTTGATGCCCTTTTTGCACGCTTTCTTACCAGCACCTGGCTCATCCGt                                                                                                                                                                                                                     >  4:464911/1‑275 (MQ=18)

                                            gactggACTGGTACCATCGACCTGCCGGAAGGCGTAGAGATTGTAATGCCGGGCGACAACATCAAAATGGTTGTTACCCTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCt                                                                                                                                                                                                             <  2:537817/270‑1 (MQ=255)

                                            gaagTGACGTGTACCATCGAATGGCAGGAAGGCGTAGGAATGGTAATGCCGGGCAACAACATCAAAAGGCGTGTTACCCTGATGCACCCGATCGCGATGGACGACAGACTGCGTTTCGCAACCCGTGAAGGTGGCCGAACCGTTGGCGGAGGCGTTCTAGGAAAAGCACTGAGCTAATTGCCGATAACACTTGACGCAATGCGAACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGGTTTCGTACCAGATCCTGGCTCATCAGTGATTTTCt                                                                                                                                                                                                             <  4:546773/272‑1 (MQ=255)

                                            gACGTGACTGGTACCATCGAACTGCCGGAAGGCGTAGAGATGGTAATGCCGGGCGACAACATCAAAATGGTTGTTACCCTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCt                                                                                                                                                                                                             <  2:65085/275‑1 (MQ=255)

                                                  aCTGGTACCATCGAACTGCCGGAAGGCGTAGAGATGGTAATGCCGGGCGCCAACATCAAAATGGTTGTTACCCTGATCCACCCGATAGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGCAGGCGGCCGTACCGTTTGCGCGGGCGTTGTAGCAAAAGTTCGGAGCTAATTGCCGATAACATTTGAAGCAATGCGCACTAAAAAGGCATCATTTGATGCCCTTTTTGCCTCCTTTCGTACCAGAACCTGGCCCATCAGGGATTTTCttttttc                                                                                                                                                                                                       >  4:492459/1‑271 (MQ=21)

                                                  aCTGGTACCATCGAACTGCCGGAAGGCGTAGAGATGGTAATGCCGGGCGACAACATCAAAATGGTTGTTACCCTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGAGGCACTGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCGTTCGTACCAGAACCTGCCTCATAAGTGATTTTCTTTGTCATata                                                                                                                                                                                                   >  1:20133/1‑277 (MQ=255)

                                                                  tGCCGGAAGGCGTAGAGATGGTAATGCCGGGCGACAACATCAAAATGGTTGTTACCCTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGAc                                                                                                                                                                                       <  2:594685/275‑1 (MQ=255)

                                                                     cGGAAGGCGTAGAGATGGTAATGCCGGGCGACAACATCAAAATGGTTGTTACCCTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCa                                                                                                                                                                                                                                                   <  4:134387/212‑1 (MQ=255)

                                                                              tAGAGATGGTAATGCCGGGCGACAACATCAAAATGGTTGTTACCCTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTg                                                                                                                                                                            <  2:465259/274‑1 (MQ=255)

                                                                               tgagaTGGTAATGTCGGGCGACAACATCAAAATGGTGGTTACCCTGATCCACCCGACCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAg                                                                                                                                                                          <  4:233632/274‑1 (MQ=255)

                                                                                gagaTGGTAATGCCGGGCGACAACATCAAAATGGTTGTTACCCTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAgg                                                                                                                                                                         <  2:624538/275‑1 (MQ=255)

                                                                                      gTAATGCCGGGCGACAACATCAAAATGGTTGTTACCCTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGNCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCCCGTCTTCGTACCAGAACCNGGCCCATACGTGTATTTCTTTTGTATAATTNNTCCTTAGACAGGCCCCGtttatgtggttaa                                                                                                                                                                  >  3:115061/1‑265 (MQ=18)

                                                                                          tGCCGGGCGACAACATCAAAATGGTTGTTACCCTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGAc                                                                                                                                                                                       <  2:234852/251‑1 (MQ=255)

                                                                                            ccGGGCGACAACATCAAAATGGTTGTTACCCTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGaa                                                                                                                                                           <  1:115061/277‑1 (MQ=255)

                                                                                                            aaaTGGTTGTTACCCGGCGCCACGCCTGCGGGCTGGACGACGGTTTGCGTTTCGCAATCCGTGAAGGCTGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGAACTAAAAGGGCATCAGTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGt                                                                                                                                             <  4:597981/275‑1 (MQ=255)

                                                                                                                     ttACCCTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAttt                                                                                                                                                                                                                                                                  <  5:6268/149‑1 (MQ=255)

                                                                                                                     ttACCCTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAttt                                                                                                                                                                                                                                                                  <  7:227531/149‑1 (MQ=255)

                                                                                                                      tACCCTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTg                                                                                                                                                                                                                                                                 >  5:330723/1‑149 (MQ=255)

                                                                                                                        cccTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGATGCAATGCGCACTAAAAGGGCATCATTTGAt                                                                                                                                                                                                                                                               <  6:1264833/149‑1 (MQ=255)

                                                                                                                         ccTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGa                                                                                                                                                                                                                                                                <  6:1177334/147‑1 (MQ=255)

                                                                                                                         ccTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATg                                                                                                                                                                                                                                                              >  6:677283/1‑149 (MQ=255)

                                                                                                                         ccTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATg                                                                                                                                                                                                                                                              >  8:1296991/1‑149 (MQ=255)

                                                                                                                         ccTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATg                                                                                                                                                                                                                                                              >  8:1326014/1‑149 (MQ=255)

                                                                                                                         ccTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATg                                                                                                                                                                                                                                                              >  5:855625/1‑149 (MQ=255)

                                                                                                                            gATCCACCCTATCGGGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGc                                                                                                                              <  2:448662/274‑1 (MQ=255)

                                                                                                                            gATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGccc                                                                                                                                                                                                                                                           <  7:1156429/149‑1 (MQ=255)

                                                                                                                            gATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGccc                                                                                                                                                                                                                                                           <  7:17987/149‑1 (MQ=255)

                                                                                                                               ccACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCttt                                                                                                                                                                                                                                                        <  8:765363/149‑1 (MQ=255)

                                                                                                                                cACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCtttt                                                                                                                                                                                                                                                       >  6:355392/1‑149 (MQ=255)

                                                                                                                                 aCCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCttttt                                                                                                                                                                                                                                                      <  7:1282526/149‑1 (MQ=255)

                                                                                                                                     gATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGTGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCAc                                                                                                                                                                                                                                                  >  5:1401264/1‑149 (MQ=255)

                                                                                                                                       tCGCGATGAACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGc                                                                                                                                                                                                                                                <  6:1401264/149‑1 (MQ=255)

                                                                                                                                           gATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTc                                                                                                                                                                                                                                            <  7:861659/149‑1 (MQ=255)

                                                                                                                                            aTGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCg                                                                                                                                                                                                                                           <  7:420173/149‑1 (MQ=255)

                                                                                                                                              ggACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTa                                                                                                                                                                                                                                         >  6:314683/1‑149 (MQ=255)

                                                                                                                                                     ggTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAAc                                                                                                                                                                                                                                  <  5:355392/149‑1 (MQ=255)

                                                                                                                                                      gTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAAcc                                                                                                                                                                                                                                 <  5:695142/149‑1 (MQ=255)

                                                                                                                                                      gTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAAcc                                                                                                                                                                                                                                 >  7:705658/1‑149 (MQ=255)

                                                                                                                                                       tCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTc                                                                                                                                                                                                                                            <  7:1577162/137‑1 (MQ=255)

                                                                                                                                                       tCTGCGTTTCGCAATCCGTAAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCt                                                                                                                                                                                                                                <  7:1149628/149‑1 (MQ=255)

                                                                                                                                                        cTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTg                                                                                                                                                                                                                               <  7:1208576/149‑1 (MQ=255)

                                                                                                                                                             tttCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGc                                                                                                                                                                                                                                                    >  8:454868/1‑123 (MQ=255)

                                                                                                                                                              ttCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCtcat                                                                                                                                                                                                                         <  5:572119/149‑1 (MQ=255)

                                                                                                                                                                 gCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAg                                                                                                                                                                                                                      >  7:944920/1‑149 (MQ=255)

                                                                                                                                                                  cAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGt                                                                                                                                                                                                                     <  6:878344/149‑1 (MQ=255)

                                                                                                                                                                       cGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCGGGGTGAAAATGTTTGTAGAAAACTTCTGACAg                                                                                  <  2:278771/275‑1 (MQ=255)

                                                                                                                                                                         tGAAGGTGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTc                                                                                                                                                                                                              <  5:900939/149‑1 (MQ=255)

                                                                                                                                                                         tGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCt                                                                                                                                                                                                                                                          >  8:646155/1‑105 (MQ=255)

                                                                                                                                                                           aaGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCtt                                                                                                                                                                                                            >  5:441694/1‑149 (MQ=255)

                                                                                                                                                                           aaGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCtt                                                                                                                                                                                                            <  7:1448124/149‑1 (MQ=255)

                                                                                                                                                                           aaGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCtt                                                                                                                                                                                                            <  5:1197960/149‑1 (MQ=255)

                                                                                                                                                                              gcggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGt                                                                                                                                                                                                         >  8:309242/1‑149 (MQ=255)

                                                                                                                                                                               cggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTc                                                                                                                                                                                                        <  8:501078/149‑1 (MQ=255)

                                                                                                                                                                                 gcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGt                                                                                                                                                                              >  4:397471/1‑173 (MQ=255)

                                                                                                                                                                                 gcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATc                                                                                                                                                               >  4:112232/1‑188 (MQ=255)

                                                                                                                                                                                     tACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGct                                                                                                                                                                                  <  4:773298/165‑1 (MQ=255)

                                                                                                                                                                                     tACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAg                                                                                                                                                                          <  2:397475/173‑1 (MQ=255)

                                                                                                                                                                                       ccGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAAcc                                                                                                                                                                                                                                 <  5:1147619/116‑1 (MQ=255)

                                                                                                                                                                                        cGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCAtt                                                                                                                                                                                               <  6:1493928/149‑1 (MQ=255)

                                                                                                                                                                                         gTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCtc                                                                                                                                                                                                                           >  7:629752/1‑120 (MQ=255)

                                                                                                                                                                                          ttGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGc                                                                                                                                                                                             >  6:297098/1‑149 (MQ=255)

                                                                                                                                                                                             gcgcGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTga                                                                                                                                                                                          >  6:129431/1‑149 (MQ=255)

                                                                                                                                                                                                cgGGGGTGGTAGCAAAAGTACGGGGCTAATGGCCGTAAACATTTGAAGAAAGTGGGAATAAAAGGGGATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGGAATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGGAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCGGGGTGAAAATGTTTGTAGAAAACTTCTGACAGGTTGGTTTATGAGTGCGAATACCg                                                          <  4:415885/274‑1 (MQ=255)

                                                                                                                                                                                                   gCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCttttt                                                                                                                                                                                                                                                      >  7:1341295/1‑83 (MQ=255)

                                                                                                                                                                                                         tAGCCAAAGTTCTGTGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCCTCAGTGATTTTCTTTGTCACAATCATTGCTGAGACAGGCTCTGt                                                                                                                                                                              <  6:163293/149‑1 (MQ=255)

                                                                                                                                                                                                                  ttCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGt                                                                                                                                                                     >  8:386305/1‑149 (MQ=255)

                                                                                                                                                                                                                  ttCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGt                                                                                                                                                                     <  6:906974/149‑1 (MQ=255)

                                                                                                                                                                                                                   tCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGTGATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCGGGGTGAAAATGTTTGTAGAAAACTTCTGACAGGTTGGTTTATGAGTGCGAATACCGAAGCTCAAGGAAGCGGgcg                                       <  2:510141/274‑1 (MQ=255)

                                                                                                                                                                                                                   tCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCGGGGTGAAAATGTTTGTAGAAAACTTCTGACAGGTTGGTTTATGAGTGCGAATACCGAAGCTCAAGGAAGCGGgcg                                       <  2:492463/274‑1 (MQ=255)

                                                                                                                                                                                                                           aaTTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCGGGGTGAAAATGTTTGTAGAAAACTTCTGACAGGTTGGTTTATGAGTGCGAATACCGAAGCTCAAGGAAGCGGGCGCGGCCTGGa                              <  2:355378/275‑1 (MQ=255)

                                                                                                                                                                                                                            aTTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCt                                                                                                                                                                                                             >  4:331725/1‑99 (MQ=255)

                                                                                                                                                                                                                            aTTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCGGGGTGAAAATGTTTGTAGAAAACTTCTAACAGGTTGGTTTATGAGTGCGAAAACCGAAGCCCAAGAAAGCGGGCGCGGCCGGGaa                             >  4:276779/1‑275 (MQ=255)

                                                                                                                                                                                                                              tGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGaaaa                                                                                                                                                         <  8:849745/149‑1 (MQ=255)

                                                                                                                                                                                                                                ccGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGt                                                                                                                                                                                                         <  2:331729/99‑1 (MQ=255)

                                                                                                                                                                                                                                ccGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGa                                                                                                                                                            <  4:472566/144‑1 (MQ=255)

                                                                                                                                                                                                                                ccGATAACATTGGACGCAATGCGAACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCGGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCGGGGTGAAAATGTTTGTAGAAAACTTCTGACAGGTTGGTTTATGAGTGCGAATACCGAAGCTCAAGGAAGCGGGCGCGGCCTGGAAGCGATGaa                     <  1:267349/279‑1 (MQ=255)

                                                                                                                                                                                                                                  gATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTa                                                                                                                                                     >  6:400219/1‑149 (MQ=255)

                                                                                                                                                                                                                                    tAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAAt                                                                                                                                                   >  7:1601484/1‑149 (MQ=255)

                                                                                                                                                                                                                                        aTTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAg                                                                                                                                                        <  8:939424/140‑1 (MQ=255)

                                                                                                                                                                                                                                        aTTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATAcgc                                                                                                                                               >  7:325364/1‑149 (MQ=255)

                                                                                                                                                                                                                                         tttGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTg                                                                                                                                                                                                                               >  5:1463869/1‑68 (MQ=255)

                                                                                                                                                                                                                                         tttGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATAcgcg                                                                                                                                              >  7:468891/1‑149 (MQ=255)

                                                                                                                                                                                                                                          ttGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGt                                                                                                                                             >  8:1460724/1‑149 (MQ=255)

                                                                                                                                                                                                                                          ttGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGt                                                                                                                                             >  8:1537329/1‑149 (MQ=255)

                                                                                                                                                                                                                                            gACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGttt                                                                                                                                           >  5:51328/1‑149 (MQ=255)

                                                                                                                                                                                                                                             aCGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCACAATCAtt                                                                                                                                                                                               <  6:786250/96‑1 (MQ=255)

                                                                                                                                                                                                                                              cGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGct                                                                                                                                                                                  <  2:181118/108‑1 (MQ=255)

                                                                                                                                                                                                                                              cGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCGGGGTGAAAATGTTTGTAGAAAACTTCTGACAGGTTGGTTTATGAGTGCGAATACCGAAGCTCAAGGAAGCg                                           >  4:525546/1‑243 (MQ=255)

                                                                                                                                                                                                                                                 aaTGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGAtt                                                                                                                                      >  8:260346/1‑149 (MQ=255)

                                                                                                                                                                                                                                                  aTGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGAttt                                                                                                                                     <  5:1283422/149‑1 (MQ=255)

                                                                                                                                                                                                                                                  aTGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGAttt                                                                                                                                     <  5:1674571/149‑1 (MQ=255)

                                                                                                                                                                                                                                                       cACTAAAAGGGCATCATTTTATTGCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCGGGGTGAAAATGTTTGTAGAAAACTTCTGACAGGTTGGTTTATGAGTGCGAATACCGAAGCTCAAGGAAGCGGGCGCGGCCTGGAAGCGATGAAGTGGGTCGTTGTGGTGGCa  <  2:687358/275‑1 (MQ=255)

                                                                                                                                                                                                                                                       cACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCGGGGTGAAAATGTTTGTAGAAAACTTCTGACAGGTTGGCTTATGAGTGCGAATACCGAAGCTCAAGGAAGCGGGCGCGGCCTGGAAGCGATGAAGTGGGTCGTTgtgg       <  2:130538/270‑1 (MQ=255)

                                                                                                                                                                                                                                                             aaGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTcg                                                                                                                          >  8:294494/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               gggCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCt                                                                                                                                                                                                             <  1:252609/64‑1 (MQ=255)

                                                                                                                                                                                                                                                                   atcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATcg                                                                                                                    >  6:225040/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                    tcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGAtt                                                                                                                                                                                                                 >  2:482067/1‑55 (MQ=255)

                                                                                                                                                                                                                                                                       tttGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCggg                                                                                                                <  6:1123924/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                       tttGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCggg                                                                                                                >  5:1099494/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                       tttGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCggg                                                                                                                <  5:297098/149‑1 (MQ=255)

CGTTCTTCAAAGGCTACCGTCCGCAGTTCTACTTCCGTACTACTGACGTGACTGGTACCATCGAACTGCCGGAAGGCGTAGAGATGGTAATGCCGGGCGACAACATCAAAATGGTTGTTACCCTGATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCGGGGTGAAAATGTTTGTAGAAAACTTCTGACAGGTTGGTTTATGAGTGCGAATACCGAAGCTCAAGGAAGCGGGCGCGGCCTGGAAGCGATGAAGTGGGTCGTTGTGGTGGCA  >  NC_000913/4176908‑4177429

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