BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,145,313	C→A	100%	S683Y (TCC→TAC)	yegE →	putative diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,145,313	0	C	A	98.8%	933.5 / 1.2	251	S683Y (TCC→TAC)	yegE	putative diguanylate cyclase
Reads supporting (aligned to +/- strand): ref base C (3/1); new base A (116/131); total (119/132)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 253 total aligned reads displayed.

TGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACTCGGTGGGGATGAATTTGGTCTGCTGCTACCAGACTGT  >  NC_000913/2145076‑2145581

tGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTTATTCAGGACGTTACCGAATCACGCAAAAGTCTGCGCCAGCTGAGCTAAAGCCCCTACCATGATGACCTGACGCATCTCGCCAACCGACCCAGt                                                                                                                                                                                                                                          >  4:557415/1‑275 (MQ=255)

            ttGGCGACAACTGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAaaca                                                                                                                                                                                                                               <  2:357299/274‑1 (MQ=255)

              ggCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCatga                                                                                                                                                                                                                                                                         >  2:105248/1‑230 (MQ=255)

              ggCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAaacaac                                                                                                                                                                                                                             >  2:513991/1‑274 (MQ=255)

                            aTTAAGGGGGAAAAGTTTCAGTGCCGGAACCGCCCGTTCCGCGACCGGGCAAGATGTGTTGTGGCCCTGTCGGAGGGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGAATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCtgctgc                                                                                                                                                                                                              <  4:450504/275‑1 (MQ=255)

                             gtAATGGAGAACATTTACAGGGCCGCAACCTCACGTTCAGCGATCGAGCAAGATGTGGGGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCtgctgc                                                                                                                                                                                                              <  4:86121/273‑1 (MQ=255)

                                 tGGAGAACATTTACAGTTCTGTAACCTCCCGGCCCGCGTTCGAGCAAGATGTTGTGTTGCTCGGTCGGGGCGGCGGCAGTTACGAAGGGCATTACAGTATGACGCCGTGAAGTACTCTTGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTAACGAATCACGCAAAATGCTGCGCCAGCTGATCTACAGCGCCTACCATGATGCACTGACGAATCTCGCCAACCGCGCCAGTTTTGAGACACAACTGCGTATCCTGCTGCAAACg                                                                                                                                                                                                         <  4:311863/275‑1 (MQ=255)

                                     gaACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGAAGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCACATGATGCCCTGACGCATCTCGCCAACCGCGCCAGTTTTTAGANACAATNNTGTATCCTGCCGCacacaaggaaaa                                                                                                                                                                                                  >  3:15054/1‑270 (MQ=255)

                                                 tGCCCACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGGGCGGCGTCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTGCTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTac                                                                                                                                                                                             <  2:677354/270‑1 (MQ=255)

                                                    cGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACa                                                                                                                                                                                 <  1:120396/279‑1 (MQ=255)

                                                     gACACCTCACGTCCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTg                                                                                                                                                                                                                                                                  >  2:317567/1‑198 (MQ=255)

                                                       cacCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGAc                                                                                                                                                                                  <  2:292226/275‑1 (MQ=255)

                                                             aCGTTCCGCGATCGAGCAAGATGTGGTGCTCCCCTGTCGGAGCGGCAGCAGTTACGACGTGCATTACAGTATTACGCCGTTACGAACTCTGGCCGGCCGCAATATTGGTTCGCCTCTGGTGTTTCAGGACGTTCCCGAACCACGCAAAATGCCGCGACAGCCGAGTTACACCGCCAACCCTCATTCACCGACGCTTCCCCCTAACGCCGCCATCTTTTCGACACAACTCCCTTACCTTCAGCAAACGCTTAAACAGATCCCATCACCgcaaccc                                                                                                                                                                              >  4:679693/1‑268 (MQ=255)

                                                                                      gtgtTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGGCGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCtt                                                                                                                                                   >  2:24299/1‑275 (MQ=255)

                                                                                      gtgtTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATg                                                                                                                                                                               >  2:109887/1‑247 (MQ=255)

                                                                                        tttGCACTGTCGGAGCGGCGGTAGTTACGACTTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAg                                                                                                                                              <  1:103165/277‑1 (MQ=255)

                                                                                             aCTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCat                                                                                                                                                                                                                                                                           >  7:11923/1‑149 (MQ=255)

                                                                                             aCTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCat                                                                                                                                                                                                                                                                           <  6:1243484/149‑1 (MQ=255)

                                                                                                 tCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGANCTACAGCGCCTACCATGATGCACTGACGCATCTCCCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATNCTGCCGCCAACG‑GTAAACAGTACACATCAGCGACACGCCCAGGTGNTNCTCCACTTTGCTCCCtttcacaccgttaaat                                                                                                                                     >  3:78202/1‑265 (MQ=255)

                                                                                                 acggacaggaggaaGTTACGACGTGTATTACAGTATTAAGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGTCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTg                                                                                                                                         <  2:555030/261‑1 (MQ=255)

                                                                                                    gAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCAc                                                                                                                                                                                                                                                                    <  8:1424871/149‑1 (MQ=255)

                                                                                                     aGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTgag                                                                                                                                                                                                                                    <  2:80824/180‑1 (MQ=255)

                                                                                                          cggcAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGc                                                                                                                                                                                                                                                              <  6:1249991/149‑1 (MQ=255)

                                                                                                                 tACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGc                                                                                                                                                                                                                                                       <  7:37688/149‑1 (MQ=255)

                                                                                                                   ccacgTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGc                                                                                                                                                                                                                                                              <  8:1076440/138‑1 (MQ=255)

                                                                                                                    gacgTGCATTACAGTATTACGCCGTTACGTACTCTGGCCGGCAGTAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCTCGCAAACTGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCAGCTCGCCaa                                                                                                                                                                                                                                                    >  8:1037866/1‑149 (MQ=255)

                                                                                                                       gTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAAAGCTGCACCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAAcca                                                                                                                                                                                                                                                 >  6:320072/1‑148 (MQ=255)

                                                                                                                          cATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCgcgc                                                                                                                                                                                                                                              <  6:595443/149‑1 (MQ=255)

                                                                                                                          cATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGAGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCTTGCTGCaaa                                                                                                                                                                                                           >  4:38975/1‑184 (MQ=255)

                                                                                                                            ttACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCa                                                                                                                                                                                                                                            >  8:1400578/1‑149 (MQ=255)

                                                                                                                                 gTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACAGCGCCAGtttt                                                                                                                                                                                                                                       <  7:1091206/149‑1 (MQ=255)

                                                                                                                                  tATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTg                                                                                                                                                                                                                                      >  7:201409/1‑149 (MQ=255)

                                                                                                                                  tATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTg                                                                                                                                                                                                                                      >  8:1026764/1‑149 (MQ=255)

                                                                                                                                   aTTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTga                                                                                                                                                                                                                                     <  8:201767/149‑1 (MQ=255)

                                                                                                                                   aTTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTga                                                                                                                                                                                                                                     <  7:1071131/149‑1 (MQ=255)

                                                                                                                                   aTTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTga                                                                                                                                                                                                                                     <  7:795674/149‑1 (MQ=255)

                                                                                                                                   aTTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTga                                                                                                                                                                                                                                     <  6:1527631/149‑1 (MQ=255)

                                                                                                                                    ttACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTgag                                                                                                                                                                                                                                    <  7:1026764/149‑1 (MQ=255)

                                                                                                                                       cGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGaaa                                                                                                                                                                                                                                 >  8:670792/1‑149 (MQ=255)

                                                                                                                                           gTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAaacaac                                                                                                                                                                                                                             <  8:11923/149‑1 (MQ=255)

                                                                                                                                           gTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAaacaac                                                                                                                                                                                                                             <  6:414556/149‑1 (MQ=255)

                                                                                                                                            ttAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACt                                                                                                                                                                                                                            >  8:1458387/1‑149 (MQ=255)

                                                                                                                                              aaGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGc                                                                                                                                                                                                                          <  8:472619/149‑1 (MQ=255)

                                                                                                                                                 tACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTa                                                                                                                                                                                                                       >  8:799037/1‑149 (MQ=255)

                                                                                                                                                    tctGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCatgat                                                                                                                                                                                                                                                                        >  7:397998/1‑97 (MQ=255)

                                                                                                                                                    tctGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATcc                                                                                                                                                                                                                    >  6:853382/1‑149 (MQ=255)

                                                                                                                                                    tctGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATcc                                                                                                                                                                                                                    <  5:179199/149‑1 (MQ=255)

                                                                                                                                                    tctGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATcc                                                                                                                                                                                                                    >  6:1366990/1‑149 (MQ=255)

                                                                                                                                                      tGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGCCCTGgtg                                                                                                                                                                       >  2:346223/1‑191 (MQ=255)

                                                                                                                                                        gACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCctgct                                                                                                                                                                                                                <  7:1429101/149‑1 (MQ=255)

                                                                                                                                                          cGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGCCCTGGTGTTTa                                                                                                                                                                   <  4:346219/191‑1 (MQ=255)

                                                                                                                                                             cagcaATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGc                                                                            <  2:575951/275‑1 (MQ=255)

                                                                                                                                                                    atTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTaaa                                                                                                                                                                                                   >  7:1276988/1‑149 (MQ=255)

                                                                                                                                                                    atTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGAt                                                                                                                                                        <  2:448840/192‑1 (MQ=255)

                                                                                                                                                                       ggttcggttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAg                                                                                                                                                                                                >  7:90842/1‑149 (MQ=255)

                                                                                                                                                                         ttcggttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTa                                                                                                                                                                                              >  6:974793/1‑149 (MQ=255)

                                                                                                                                                                         ttcggttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAAAG‑GTAAACAGTa                                                                                                                                                                                              <  7:1458387/149‑1 (MQ=255)

                                                                                                                                                                          tcggttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTAtgctg                                                                 >  2:232556/1‑273 (MQ=255)

                                                                                                                                                                           cggttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCa                                                                                                                                                                                                                                                                            <  8:1575053/70‑1 (MQ=255)

                                                                                                                                                                           cggttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACTCTTTGCTGCGCTCACTGGCGTCGTTAATGCTGCGTAtgctg                                                                 >  2:480484/1‑272 (MQ=255)

                                                                                                                                                                              ttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTa                                                                                                                                                                                              <  8:772836/144‑1 (MQ=255)

                                                                                                                                                                              ttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTa                                                                                                                                                                                              >  7:772836/1‑144 (MQ=255)

                                                                                                                                                                              ttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACAt                                                                                                                                                                                         >  5:1224722/1‑149 (MQ=255)

                                                                                                                                                                              ttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCAGGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACAt                                                                                                                                                                                         >  8:117787/1‑149 (MQ=255)

                                                                                                                                                                                  ggTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGc                                                                                                                                                                                     >  5:466604/1‑149 (MQ=255)

                                                                                                                                                                                   gTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTaaa                                                                                                                                                                                                   >  3:114022/1‑134 (MQ=255)

                                                                                                                                                                                     gATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGAc                                                                                                                                                                                  >  6:514090/1‑149 (MQ=255)

                                                                                                                                                                                      aTTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACa                                                                                                                                                                                 >  8:321797/1‑149 (MQ=255)

                                                                                                                                                                                         cAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGc                                                                                                                                                                              <  6:719158/149‑1 (MQ=255)

                                                                                                                                                                                           ggACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGccc                                                                                                                                                                            >  6:76087/1‑149 (MQ=255)

                                                                                                                                                                                           ggACGTTACCGAATCACGCAAAATGCTGCCCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAATCG‑GTCACCAGTACACATCAACGACATGCCCTGGTGTTTTACGCTCTTGATCGCTTTAAACCGGTGAATGACATCGCCCGGCATCCGGCGGGTTACGCTTTGCTGCGCGCACTCGCGCACTTCATGCTGAGTATCCTGCGTTCCAGctccgtcttgt                                              >  4:115518/1‑265 (MQ=255)

                                                                                                                                                                                            gACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGCCCt                                                                                                                                                                           <  5:100809/149‑1 (MQ=255)

                                                                                                                                                                                                 tACCGAATCACGCAAAATGCTGCGCCATCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGCCCTGgtgt                                                                                                                                                                      <  5:1352827/149‑1 (MQ=255)

                                                                                                                                                                                                  aCCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACAt                                                                                                                                                                                >  2:186015/1‑138 (MQ=255)

                                                                                                                                                                                                   ccGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGc                                                                                                                                                                              <  6:435166/139‑1 (MQ=255)

                                                                                                                                                                                                   ccGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGCCCTGGTGttt                                                                                                                                                                    >  7:161204/1‑149 (MQ=255)

                                                                                                                                                                                                    cGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGCCCTGGTGTTTa                                                                                                                                                                   <  8:314389/149‑1 (MQ=255)

                                                                                                                                                                                                    cGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGCCCTGGTGTTTa                                                                                                                                                                   >  5:1478262/1‑149 (MQ=255)

                                                                                                                                                                                                    cGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGCCCTGGTGTTTa                                                                                                                                                                   <  5:824161/149‑1 (MQ=255)

                                                                                                                                                                                                      aaTCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACg                                                                                                                                                                                                         <  2:767921/110‑1 (MQ=255)

                                                                                                                                                                                                      aaTCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGccc                                                                                                                                                                            <  4:186012/138‑1 (MQ=255)

                                                                                                                                                                                                          aCGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATc                                                                                                                                                             <  7:321797/149‑1 (MQ=255)

                                                                                                                                                                                                           cGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCatg                                                                                                                                                                                                                                                                          >  1:257183/1‑40 (MQ=255)

                                                                                                                                                                                                           cGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCatg                                                                                                                                                                                                                                                                          <  3:257183/40‑1 (MQ=255)

                                                                                                                                                                                                              aaaaTGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCa                                                                                                                                                                                       >  5:564520/1‑119 (MQ=255)

                                                                                                                                                                                                               aaaTGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTaaa                                                                                                                                                                                                   >  6:2868/1‑106 (MQ=255)

                                                                                                                                                                                                               aaaTGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGAc                                                                                                                                                                                  <  2:709844/123‑1 (MQ=255)

                                                                                                                                                                                                               aaaTGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAg                                                                                                                                              <  2:192545/159‑1 (MQ=255)

                                                                                                                                                                                                                aaTGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATc                                                                                                                                                       >  5:1606448/1‑149 (MQ=255)

                                                                                                                                                                                                                    ctgcGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCtt                                                                                                                                                   <  8:1276988/149‑1 (MQ=255)

                                                                                                                                                                                                                    ctgcGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCtt                                                                                                                                                   <  5:1694810/149‑1 (MQ=255)

                                                                                                                                                                                                                        gcCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTaaa                                                                                                                                               >  7:664863/1‑149 (MQ=255)

                                                                                                                                                                                                                          cAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATc                                                                                                                                                                                                                     <  4:714712/78‑1 (MQ=255)

                                                                                                                                                                                                                           aGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCg                                                                                                                                            <  8:418540/149‑1 (MQ=255)

                                                                                                                                                                                                                             cTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGt                                                                                                                                          <  6:1478262/149‑1 (MQ=255)

                                                                                                                                                                                                                              tGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTg                                                                                                                                         >  7:271489/1‑149 (MQ=255)

                                                                                                                                                                                                                                    aCAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATCCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGAc                                                                                                                                   <  5:321568/149‑1 (MQ=255)

                                                                                                                                                                                                                                    aCAGCGCCTACCATGATGCACTGACGCATCTCGCAAACCGCGCCAGTTTTGAGAAACAACTGCGTGTCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACTCGGTGGGGATGAATTTGGTCTGCTGCTACCAGACTGt  <  1:203917/278‑1 (MQ=255)

                                                                                                                                                                                                                                     cAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACa                                                                                                                                  <  7:742834/149‑1 (MQ=255)

                                                                                                                                                                                                                                     cAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACa                                                                                                                                  <  7:5093/149‑1 (MQ=255)

TGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACG‑GTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACTCGGTGGGGATGAATTTGGTCTGCTGCTACCAGACTGT  >  NC_000913/2145076‑2145581

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