BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,344,033	A→G	5.5%	L105P (CTG→CCG)	yciT ←	global regulator of transcription; DeoR family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,344,033	0	A	G	5.5%	571.5 / 14.9	179	L105P (CTG→CCG)	yciT	global regulator of transcription; DeoR family
Reads supporting (aligned to +/- strand): ref base A (80/84); new base G (8/2); total (90/89)
Fisher's exact test for biased strand distribution p-value = 9.95e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.64e-02

Only 100 of 184 total aligned reads displayed.

TTTTCCAGCACGGCATTGACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATCGGTTACGCTAACCTGACCCTGGT  >  NC_000913/1343769‑1344306

ttttCCAGCACGGCATTGACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCgg                                                                                                                                                                                                                                                                           >  2:316554/1‑273 (MQ=255)

                                          cgACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTTCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAA‑ATGGCCAACCATACTTTCGCTTTTTTTCTGATCCACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCCTACCGTGATGATAGTGCACTTTTTCTTCCACTCGCCCAAAGACCGGGCGATTAGGGCATTTCGTCTGCCATTTTGGATAAACCTGGTTTc                                                                                                                                                                                                                               >  4:490927/1‑275 (MQ=255)

                                            aCCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGc                                                                                                                                                                                                                             <  4:455690/275‑1 (MQ=255)

                                                                   gccagccaTCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCa                                                                                                                                                                                                      >  1:81289/1‑275 (MQ=255)

                                                                                       aaTGGTTTGCTGAAAGGTACCTGTTGGATGCACGTGC‑TGTTAATAGGGCCAACCTTACTTTCGGTTTTTTTCTGATACACGCCACCGAGAAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCtt                                                                                                                                                                                  <  4:559498/275‑1 (MQ=255)

                                                                                         tGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCCCCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCCGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTACGCCGGGTTGAACCAGGGAACGGGCAACATCGGACAGTTCACACTTCa                                                                                                                                                                                >  4:75180/1‑275 (MQ=255)

                                                                                             ttGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGt                                                                                                                                                                            >  2:366716/1‑275 (MQ=255)

                                                                                              tGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACGCGGCGAGTTCACGCTTCAGCGta                                                                                                                                                                           >  2:396656/1‑275 (MQ=255)

                                                                                                        cACCTGTTGGATGCACTGGCGTGTCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGtt                                                                                                                                                                                                                                 >  2:591165/1‑211 (MQ=255)

                                                                                                              ttGGATGCACTGGCGTGTCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATa                                                                                                                                                                                                                                          >  2:451169/1‑196 (MQ=255)

                                                                                                                ggATGCACTGGCGTGTCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCATCGTATAGTTGCTCATCATatg                                                                                                                                                          >  1:12745/1‑272 (MQ=255)

                                                                                                                      aCTGGCGTGTCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCaga                                                                                                                                                                                                                                                                                 <  5:627200/149‑1 (MQ=255)

                                                                                                                        tGGCGTGTCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGt                                                                                                                                                                                                                                                                               >  5:120539/1‑149 (MQ=255)

                                                                                                                        tGGCGTGTCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGt                                                                                                                                                                                                                                                                               >  5:43672/1‑149 (MQ=255)

                                                                                                                        tGGCGTGTCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGAc                                                                                                                                                 >  4:61717/1‑275 (MQ=255)

                                                                                                                           cGTGTCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATa                                                                                                                                                                                                                                          >  2:523651/1‑183 (MQ=255)

                                                                                                                             tgtCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTACTCGCCCGGACTCTggg                                                                                                                                                                                                                                                                          <  7:827969/149‑1 (MQ=255)

                                                                                                                               tCAA‑AGTGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGa                                                                                                                                                                                                                                      <  4:523651/183‑1 (MQ=255)

                                                                                                                               ccAA‑AGGGCCCACCATATTTTCGCTTTTTTTCTGATACACGCCACTGAGCAAAATAACTTCTCAAGGTTCGTCTTTCAGCAAATGCGCGACGTAGCTGCTGACCGTGCTGATAGTGCCATTTTTCTTCTGCACGCCCCGAGTCTGGGCGAGTAGGGAATTGCTGCTGCCATTTTCGACAAAGATGGTTTCGCCCGGTTGAAACAGCGACGGGGCAAACTCGGCGAGTTCACGATTCCGAGTATTGTTGCACAACATCGGGGTTTCGagggaattc                                                                                                                                          >  4:593511/2‑267 (MQ=255)

                                                                                                                                 aa‑aGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAg                                                                                                                                                                                                                                                                      >  6:1042069/1‑149 (MQ=255)

                                                                                                                                    aGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTa                                                                                                                                                                                                                                                                    <  8:309565/149‑1 (MQ=255)

                                                                                                                                      ggCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAgg                                                                                                                                                                                                                                                                  <  5:261386/149‑1 (MQ=255)

                                                                                                                                      ggCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAgg                                                                                                                                                                                                                                                                  <  6:455523/149‑1 (MQ=255)

                                                                                                                                      ggCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTACTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAgg                                                                                                                                                                                                                                                                  <  8:6886/149‑1 (MQ=255)

                                                                                                                                       gCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAggg                                                                                                                                                                                                                                                                 >  5:319390/1‑149 (MQ=255)

                                                                                                                                        ccAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGc                                                                                                                                                                                                                                                                >  7:705749/1‑149 (MQ=255)

                                                                                                                                          aaCCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCAt                                                                                                                                                                                                                                                              >  5:866194/1‑149 (MQ=255)

                                                                                                                                           ccctttttttctttttttttCANNAACACNACCCCGGCAAAAATACCTGCACCAGGGGCGCGTTGCAGCAAATGGCCGATTTGGCTGCTGACCGTGATTATGGTGTCCTTTTGCTTATGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGa                                                                                                                               <  3:9082/262‑1 (MQ=255)

                                                                                                                                                  tttCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCa                                                                                                                                                                                                                                                                                   <  7:578471/120‑1 (MQ=255)

                                                                                                                                                  tttCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCa                                                                                                                        <  2:61717/275‑1 (MQ=255)

                                                                                                                                                       cTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGTTTTCGCCCGGTTGAACCAGCGACGCGCCAAACCCGGCGAGTTCACGCCTCAGCGTTTACTTGCTCAACATAGGGGTTTCGACCTCCTCAATATCAAGCGAAACACCAAAccc                                                                                                                   >  4:621565/1‑275 (MQ=255)

                                                                                                                                                        tttttttCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCAttt                                                                                                                                                                                                                                                >  5:68371/1‑149 (MQ=255)

                                                                                                                                                        tttttttCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATGGCTACTGCCAttt                                                                                                                                                                                                                                                >  6:850423/1‑149 (MQ=255)

                                                                                                                                                        tttttttCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGAAGATCGTGACATTTTACTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCAttt                                                                                                                                                                                                                                                >  6:789976/1‑149 (MQ=255)

                                                                                                                                                           ttttCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCg                                                                                                                                                                                                                                             >  8:458545/1‑149 (MQ=255)

                                                                                                                                                                tGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCgg                                                                                                                                                                                                                                                                           <  5:365187/114‑1 (MQ=255)

                                                                                                                                                                tGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCgg                                                                                                                                                                                                                                                                           >  6:365187/1‑114 (MQ=255)

                                                                                                                                                                  aTACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGt                                                                                                                                                                                                                                  >  4:328183/1‑153 (MQ=255)

                                                                                                                                                                   tACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGAt                                                                                                                                                                                                                                     <  5:445421/149‑1 (MQ=255)

                                                                                                                                                                    acacGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATg                                                                                                                                                                                                                                    <  6:909521/149‑1 (MQ=255)

                                                                                                                                                                     cacGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATCTCCGATAAAGATgg                                                                                                                                                                                                                                   <  7:82103/149‑1 (MQ=255)

                                                                                                                                                                      acGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGt                                                                                                                                                                                                                                  <  6:284090/149‑1 (MQ=255)

                                                                                                                                                                      acGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGt                                                                                                                                                                                                                                  >  8:19006/1‑149 (MQ=255)

                                                                                                                                                                         ccACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTc                                                                                                                                                                                                                               <  6:562139/149‑1 (MQ=255)

                                                                                                                                                                          cACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCg                                                                                                                                                                                                                              <  5:1042069/149‑1 (MQ=255)

                                                                                                                                                                           aCCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGc                                                                                                                                                                                                                             >  7:1017817/1‑149 (MQ=255)

                                                                                                                                                                            ccGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTc                                                                                                                                                                                                                               <  7:908600/146‑1 (MQ=255)

                                                                                                                                                                            ccGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTc                                                                                                                                                                                                                               >  8:908600/1‑146 (MQ=255)

                                                                                                                                                                            ccGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCCGCCATTTTCGATAAAGATGGTTTCGcc                                                                                                                                                                                                                            <  7:722156/149‑1 (MQ=255)

                                                                                                                                                                                   aaaTAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGa                                                                                                                                                                                                                     >  7:117744/1‑149 (MQ=255)

                                                                                                                                                                                   aaaTAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGa                                                                                                                                                                                                                     <  6:326785/149‑1 (MQ=255)

                                                                                                                                                                                   aaaTAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGa                                                                                                                                                                                                                     <  6:736022/149‑1 (MQ=255)

                                                                                                                                                                                    aaTAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGaa                                                                                                                                                                                                                    <  7:660658/149‑1 (MQ=255)

                                                                                                                                                                                    aaTAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGt                                                                                                                                                                       >  4:246584/1‑194 (MQ=255)

                                                                                                                                                                                       aaCTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCa                                                                                                                                                                                                                 <  6:796107/149‑1 (MQ=255)

                                                                                                                                                                                        aCTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCt                                                                                                                                                                   <  2:246584/194‑1 (MQ=255)

                                                                                                                                                                                         ctccacaAGGCGCGTTTTTCATCAAACGCGCGATGTAGCTGCTGACCGTTATGTTAGTGACATTTTTCTCCTGCCCGCCCCGGGTCCGGGCGCTTAGGGCATTGCTGCTGCCATTGTCGATAAAGATGTTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAgt                                                                                 <  4:546230/272‑1 (MQ=255)

                                                                                                                                                                                         cTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGc                                                                                                                                                                                                               <  5:765543/149‑1 (MQ=255)

                                                                                                                                                                                         cTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATAATACTGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTCACACCCTTTTTTGaag                                                                                 >  1:189710/1‑273 (MQ=255)

                                                                                                                                                                                             acaAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATgg                                                                                                                                                                                                                                   <  6:64037/125‑1 (MQ=255)

                                                                                                                                                                                              caAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCa                                                                                                                                                                                                                 >  4:213054/1‑142 (MQ=255)

                                                                                                                                                                                                    gcgTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTa                                                                                                                                                                                                                                                                    <  8:453015/85‑1 (MQ=255)

                                                                                                                                                                                                    gcgTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTa                                                                                                                                                                                                                                                                    >  7:453015/1‑85 (MQ=255)

                                                                                                                                                                                                    gcgTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGc                                                                                                                                                                                                                             >  8:225478/1‑124 (MQ=255)

                                                                                                                                                                                                    gcgTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGc                                                                                                                                                                                                                             <  7:225478/124‑1 (MQ=255)

                                                                                                                                                                                                    gcgTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCaaa                                                                                                                                                                                                    >  8:884334/1‑149 (MQ=255)

                                                                                                                                                                                                     cgTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATa                                                                                                                                                                                                                                          >  4:447907/1‑110 (MQ=255)

                                                                                                                                                                                                            cAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGa                                                                                                                                                                                            <  6:386968/149‑1 (MQ=255)

                                                                                                                                                                                                              gCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGAt                                                                                                                                                                                                                                           <  6:71839/100‑1 (MQ=255)

                                                                                                                                                                                                              gCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGAt                                                                                                                                                                                                                                           >  5:71839/1‑100 (MQ=255)

                                                                                                                                                                                                              gCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGt                                                                                                                                                                                          >  7:1219/1‑149 (MQ=255)

                                                                                                                                                                                                              gCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGt                                                                                                                                                                                          <  7:369016/149‑1 (MQ=255)

                                                                                                                                                                                                                aaaTGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTc                                                                                                                                                                                        >  5:157722/1‑149 (MQ=255)

                                                                                                                                                                                                                 aaTGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCa                                                                                                                                                                                       <  6:157267/149‑1 (MQ=255)

                                                                                                                                                                                                                 aaTGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCa                                                                                                                                                                                       <  8:1012600/149‑1 (MQ=255)

                                                                                                                                                                                                                             aGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGta                                                                                                                                                                           <  8:27724/149‑1 (MQ=255)

                                                                                                                                                                                                                               ctgctgACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGtata                                                                                                                                                                         <  7:458545/149‑1 (MQ=255)

                                                                                                                                                                                                                                     aCCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCt                                                                                                                                                                   <  8:335689/149‑1 (MQ=255)

                                                                                                                                                                                                                                        gTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTg                                                                           <  2:633351/234‑1 (MQ=255)

                                                                                                                                                                                                                                            tgatAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATa                                                                                                                                                            <  8:417523/149‑1 (MQ=255)

                                                                                                                                                                                                                                              atAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGTCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACg                                                                                                                                                          >  8:599966/1‑149 (MQ=255)

                                                                                                                                                                                                                                               tAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACgg                                                                                                                                                         <  8:1219/149‑1 (MQ=255)

                                                                                                                                                                                                                                                   gACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGttt                                                                                                                                                     <  6:43672/149‑1 (MQ=255)

                                                                                                                                                                                                                                                    aCATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTTCTGCCATTTTCGATAAAGATGGGTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTAAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATCgg                       <  2:621565/274‑1 (MQ=255)

                                                                                                                                                                                                                                                       tttttCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGc                                                                                                                                                                    <  6:514145/130‑1 (MQ=255)

                                                                                                                                                                                                                                                       tttttCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGAc                                                                                                                                                 <  6:179784/149‑1 (MQ=255)

                                                                                                                                                                                                                                                       tttttCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGAc                                                                                                                                                 <  8:117744/149‑1 (MQ=255)

                                                                                                                                                                                                                                                        ttttCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCGCGCTTCAGCGTATAGTTGCTCATCATAc                                                                                                                                                           <  2:132674/138‑1 (MQ=255)

                                                                                                                                                                                                                                                        ttttCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGa                                                                                   <  2:87550/210‑1 (MQ=255)

                                                                                                                                                                                                                                                            gccctacTCGCCCGGACTCTGGGCGAGCAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGtcat                                                                                                                                            >  6:880422/7‑149 (MQ=255)

                                                                                                                                                                                                                                                             ttctGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGtcatc                                                                                                                                           >  7:255267/1‑149 (MQ=255)

                                                                                                                                                                                                                                                             ttctGCTCGCCCAGAGTCC‑GGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGtcatca                                                                                                                                          >  8:858166/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               ctacTCGCCCGGACTCTGGGCGAGCAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCAc                                                                                                                                         >  6:281887/4‑149 (MQ=255)

                                                                                                                                                                                                                                                                tacTCGCCCGGACTCTGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACt                                                                                                                                        >  7:252985/3‑149 (MQ=255)

                                                                                                                                                                                                                                                                 gCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGt                                                                                                                                                                                                                                  >  2:616362/1‑58 (MQ=255)

                                                                                                                                                                                                                                                                     gCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCaa                                                                                                                                   <  5:1054005/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                       ccAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGc                                                                                                                                 <  7:394078/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                        cAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCg                                                                                                                                >  5:796258/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                        cAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATCGGTTACGCTAACCTGACCCTGGt  <  2:490927/275‑1 (MQ=255)

                                                                                                                                                                                                                                                                        cAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGAGAAAGATGGTTTCGCCCGGTTGAACCAGCGAGGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTACGACGTCATCACTATCAAGCg                                                                                                                                <  5:789976/149‑1 (MQ=255)

TTTTCCAGCACGGCATTGACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAA‑AGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATCGGTTACGCTAACCTGACCCTGGT  >  NC_000913/1343769‑1344306

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