BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,270,144	G→C	5.5%	intergenic (+17/‑94)	aphA → / → yjbQ	acid phosphatase/phosphotransferase, class B, non‑specific/thiamine phosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,270,144	0	G	C	5.5%	573.6 / 19.5	165	intergenic (+17/‑94)	aphA/yjbQ	acid phosphatase/phosphotransferase, class B, non‑specific/thiamine phosphate synthase
Reads supporting (aligned to +/- strand): ref base G (83/72); new base C (5/4); total (89/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.54e-02

Only 100 of 165 total aligned reads displayed.

TTTTCATATTCCTGCCACCAACATGAATCCGGTGATCTTTGCGGGCGATAAACCAGGGCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAACTGCGATCCCACCGTA  >  NC_000913/4269875‑4270414

ttttCATATTCCTGCCACCAACATGAATCCGGTGATCTTTGCGGGCGATAAACCAGGGCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGct                                                                                                                                                                                                                                                                           <  4:56847/275‑1 (MQ=255)

ttttCATATTCCTGCCACCAACATGAATCCGGTGATCTTTGCGGGCGATAAACCAGGGCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGct                                                                                                                                                                                                                                                                           <  2:185563/275‑1 (MQ=255)

                         aaTCCGGGGATCTTGTCGGGCAATAAACCAGGGAAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAAGTTTTTATGGCGATGCTGATAAGGACATTACCGCCGCACGCGATGTCGGCTCTTGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTGTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATca                                                                                                                                                                                                                                                  <  1:86645/275‑1 (MQ=255)

                                                         ggaaaaaccAAAATCGTAATTGCTGCAGGATAAAAAGATCCGTATGTTTTATGGCGATTCGGTTTATTCTATTCCCGACGCACGCGAAGTCGGCGCTCGTGGTATCCGCATTCTGCGGGCCTCCAACTCGACCTTCAAGCCCTTGCGACAAGCGGGTGCGTTTGGTGAAGAGGTTATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGGGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCaaa                                                                                                                                                                                                                   <  4:656682/266‑1 (MQ=255)

                                                             aaTACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTgc                                                                                                                                                                                                               >  2:702185/1‑273 (MQ=255)

                                                                                               tCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGttttttt                                                                                                                                                                                                                                                       >  1:134082/1‑200 (MQ=255)

                                                                                               tCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCACGAATGCTGTTTTTTTa                                                                                                                                                                                                                                                      >  2:312935/1‑201 (MQ=255)

                                                                                                                             atatTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATg                                                                                                                                                                                                                                                                             >  6:658654/1‑148 (MQ=255)

                                                                                                                              tatTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATATTGACAGAGCGGGAGAGCGTGATGct                                                                                                                                                                                                                                                                           >  8:309828/1‑149 (MQ=255)

                                                                                                                              tatTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGc                                                                                                                                                                                                                                                                            >  6:121906/1‑148 (MQ=255)

                                                                                                                                ttACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGctct                                                                                                                                                                                                                                                                         <  7:333178/149‑1 (MQ=255)

                                                                                                                                 tACCGCCGCATGCGATGTCGGCGCTCGTGGTATCTGCATTCTGCGCGCCTCCAACTCTACTTTCAAACCCTTGCCACCAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCCGAATACTGNAATCGCGGTAGAGCGTGATGCTCNCCCGNATGCTTTTTTTTTAATCATACTTTTATCCTTTCGCTGACTTGCTGCAAACTGa                                                                                                                                                                                                               >  4:718522/1‑206 (MQ=255)

                                                                                                                                 tACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTATAAACCCTTGCCACAAGCGGTTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCtctc                                                                                                                                                                                                                                                                        <  7:556275/149‑1 (MQ=255)

                                                                                                                                 tACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGcac                                                                                                                                                                                   >  2:353889/1‑234 (MQ=255)

                                                                                                                                   ccgccgCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTcc                                                                                                                                                                                                                                                                       >  5:506423/1‑148 (MQ=255)

                                                                                                                                   ccgccgCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGAGTGATGCTCTcc                                                                                                                                                                                                                                                                       <  5:305616/148‑1 (MQ=255)

                                                                                                                                     gccgcACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACAtgt                                                                                                                                                                               <  4:353887/234‑1 (MQ=255)

                                                                                                                                      ccgcACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGt                                                                                                                                     >  1:132635/1‑275 (MQ=255)

                                                                                                                                          aCGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATg                                                                                                                                                                                                                                                                             <  4:426214/135‑1 (MQ=255)

                                                                                                                                              gATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGttt                                                                                                                                                                                                                                                           >  7:510637/1‑149 (MQ=255)

                                                                                                                                               aTGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGttt                                                                                                                                                                                                                                                           >  5:830752/1‑148 (MQ=255)

                                                                                                                                                 gTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGttttt                                                                                                                                                                                                                                                         <  5:658654/148‑1 (MQ=255)

                                                                                                                                                    ggCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTa                                                                                                                                                                                                                                                      <  6:119644/148‑1 (MQ=255)

                                                                                                                                                          cGTGGTATCGGCATTCTGCGCGCCTCCAACTCTCCCTACAAACCCTGGCCACAAACGGGTGCGTTTGGTGTATAGGACATCGTCAATTCAGAATACCGCCTGAGCGGGAGAGCGTGATGCTCTCCCGAATGCCGGTTTTTTAATCCCACCTTTACCCTTTCTCTGTCTTGCTGCAAACTGATTAAGCGAGTTTTATCAGGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTGACCGCCACATTCCGAGGGTTTCATTAGGAAAGGgaagaa                                                                                                                 >  4:264594/1‑275 (MQ=255)

                                                                                                                                                              gTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCttt                                                                                                                                                                                                                                           >  8:763300/1‑149 (MQ=255)

                                                                                                                                                                aTCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTAt                                                                                                                                                                                                                                         >  7:281682/1‑149 (MQ=255)

                                                                                                                                                                aTCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTAt                                                                                                                                                                                                                                         >  7:92299/1‑149 (MQ=255)

                                                                                                                                                                 tCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATc                                                                                                                                                                                                                                        >  8:80653/1‑149 (MQ=255)

                                                                                                                                                                    gCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCt                                                                                                                                                                                                                                      >  5:54144/1‑148 (MQ=255)

                                                                                                                                                                     cATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCCTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGATTc                                                                                                       >  1:231019/1‑274 (MQ=255)

                                                                                                                                                                      aTTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCttt                                                                                                                                                                                                                                    <  5:736893/148‑1 (MQ=255)

                                                                                                                                                                               gcCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTc                                                                                                                                                                                                                                   <  8:503810/140‑1 (MQ=255)

                                                                                                                                                                                ccTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTtg                                                                                                                                                                                                                          <  5:121906/148‑1 (MQ=255)

                                                                                                                                                                                   ccAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTtgctg                                                                                                                                                                                                                       >  6:265727/1‑148 (MQ=255)

                                                                                                                                                                                          tacctacAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATg                                                                                                                                                                                                                                                                             <  4:170379/87‑1 (MQ=255)

                                                                                                                                                                                           acctacAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAAt                                                                                                                                                                                                                                                    >  8:687878/1‑111 (MQ=255)

                                                                                                                                                                                           acctacAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCaa                                                                                                                                                                       >  4:59043/1‑188 (MQ=255)

                                                                                                                                                                                           acctacAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGtt                                                                                >  2:338462/1‑275 (MQ=255)

                                                                                                                                                                                           acctacAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGCTTTTTTAATCCCACCTTTATCCTTTCGCTGTCTTGCTGAAAACTGTTTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGGGGGTTTCATCTGGTAACGGAAGAAATTTTGAATCAGCTGGCTGATATGAcgctg                                                                                   >  2:290915/1‑270 (MQ=255)

                                                                                                                                                                                            cctacAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGtt                                                                                >  2:78820/1‑274 (MQ=255)

                                                                                                                                                                                               acAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTaa                                                                                                                                                                                                           >  6:849718/1‑148 (MQ=255)

                                                                                                                                                                                               acAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCa                                                                                                                                                                                              <  2:357522/161‑1 (MQ=255)

                                                                                                                                                                                                cAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAag                                                                                                                                                                                                          >  6:813958/1‑148 (MQ=255)

                                                                                                                                                                                                  aaCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCaa                                                                                                                                                                                             >  1:215032/1‑159 (MQ=255)

                                                                                                                                                                                                  aaCCCTTGCCACAAGCGGGTGAGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAagag                                                                                                                                                                                                        >  6:344728/1‑148 (MQ=255)

                                                                                                                                                                                                    cccTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGCATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGc                                                                       <  2:264595/275‑1 (MQ=255)

                                                                                                                                                                                                     ccTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGtt                                                                                                                                                                                                    <  7:763300/149‑1 (MQ=255)

                                                                                                                                                                                                       ttGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCacac                                                                                                                                                                                                                                               <  8:258119/104‑1 (MQ=255)

                                                                                                                                                                                                       ttGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCacac                                                                                                                                                                                                                                               >  7:258119/1‑104 (MQ=255)

                                                                                                                                                                                                       ttGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAagaga                                                                                                                                                                                                       <  8:51018/144‑1 (MQ=255)

                                                                                                                                                                                                       ttGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAagaga                                                                                                                                                                                                       >  7:51018/1‑144 (MQ=255)

                                                                                                                                                                                                       ttGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGtttt                                                                                                                                                                                                  <  7:80653/149‑1 (MQ=255)

                                                                                                                                                                                                         gCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACCCCTTTATTCTTTCGCTGTCTTGCTGCAAACTGATTTAGCGCGTTTTATCaa                                                                                                                                                                                             >  4:794017/1‑152 (MQ=255)

                                                                                                                                                                                                            acaAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCaaa                                                                                                                                                                                                                   <  2:671708/127‑1 (MQ=255)

                                                                                                                                                                                                            acaAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTaa                                                                                                                                                                                                           <  8:538100/135‑1 (MQ=255)

                                                                                                                                                                                                             caAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAg                                                                                                                                                                                            <  7:256157/149‑1 (MQ=255)

                                                                                                                                                                                                             caAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGag                                                                                                                                                                                         <  2:794019/152‑1 (MQ=255)

                                                                                                                                                                                                              aaGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAg                                                                                                                                                                                            >  5:727426/1‑148 (MQ=255)

                                                                                                                                                                                                               aGCGGGTGCTTTTGGTGAAGAGGTGATCGCCAATTCAGAATACTTACAGAGCGTGAGAGCGTGATGTCCTCCCGCATGCTGTTTTTTTAATCACACCTTTATCCTTTCGTTTCCTTGCTGCAAACTGATTAAGAGAGCTTTATCAAGGAGCATCACATGTGGTATTAAAAGACGCTCAGGCTTAGCCCAAAATTTCGTGGGTTTCATCTGCTAACTGTTGAATTTCTTAATCAGCTGGCTGATATTCCGCCGGTTAACTACGGTTTACTgaacct                                                            >  4:334900/1‑270 (MQ=255)

                                                                                                                                                                                                               aGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCa                                                                                                                                                             >  2:171135/1‑178 (MQ=255)

                                                                                                                                                                                                                   ggTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGagca                                                                                                                                                                                       >  6:554239/1‑148 (MQ=255)

                                                                                                                                                                                                                    gTGCGGTTGGTGAAGAGGTGATCGTCAATTCATAATACTGACAGATCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTtgc                                                       <  4:78819/275‑1 (MQ=255)

                                                                                                                                                                                                                        gTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGcacac                                                                                                                                                                                 >  7:78289/1‑148 (MQ=255)

                                                                                                                                                                                                                         tttGGTGAAGAGGTGATCGTCAATTTAGAATACTGACAGAGCGGGAGAGCATCACGCTCTcc                                                                                                                                                                                                                                                                       <  6:691334/62‑1 (MQ=25)

                                                                                                                                                                                                                         tttGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTaa                                                                                                                                                                                                                                                     >  5:714127/1‑80 (MQ=255)

                                                                                                                                                                                                                         tttGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTaa                                                                                                                                                                                                                                                     <  6:714127/80‑1 (MQ=255)

                                                                                                                                                                                                                         tttGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCaaa                                                                                                                                                                                                                   <  4:328082/114‑1 (MQ=255)

                                                                                                                                                                                                                         tttGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCATCACGCTCTcc                                                                                                                                                                                                                                                                       >  5:691334/1‑62 (MQ=37)

                                                                                                                                                                                                                           tGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACAtgt                                                                                                                                                                               <  5:813958/148‑1 (MQ=255)

                                                                                                                                                                                                                            ggTGAAGAGGTGATCGTCAATTCAGAATATTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTgg                                                                                                                                                                             <  7:309828/149‑1 (MQ=255)

                                                                                                                                                                                                                                aaGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATc                                                                                                                                                                         >  7:692863/1‑149 (MQ=255)

                                                                                                                                                                                                                                    ggTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGc                                       <  1:86230/275‑1 (MQ=255)

                                                                                                                                                                                                                                     gTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAg                                                                                                                                                                    >  8:383597/1‑149 (MQ=255)

                                                                                                                                                                                                                                     gTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAg                                                                                                                                                                    <  8:281682/149‑1 (MQ=255)

                                                                                                                                                                                                                                     gTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGc                                       <  2:673367/274‑1 (MQ=255)

                                                                                                                                                                                                                                      tGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAg                                                                                                                                                                    <  6:516357/148‑1 (MQ=255)

                                                                                                                                                                                                                                        aTCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCctct                                   >  4:117932/1‑275 (MQ=255)

                                                                                                                                                                                                                                        aTCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCCCGTTAACATCGGCTTACCGCATCTTTTTCTGCCAAATAcccccccccca                                   >  4:671707/1‑272 (MQ=255)

                                                                                                                                                                                                                                         tCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACg                                                                                                                                                                 >  5:790424/1‑148 (MQ=255)

                                                                                                                                                                                                                                            tCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTgg                                                                                                                                                                             <  7:578830/133‑1 (MQ=255)

                                                                                                                                                                                                                                                ttCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGc                                                                >  2:474887/1‑238 (MQ=255)

                                                                                                                                                                                                                                                 tCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCt                                                                                                                                                         >  6:461128/1‑148 (MQ=255)

                                                                                                                                                                                                                                                 tCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAAt                                                                                                        <  2:592373/197‑1 (MQ=255)

                                                                                                                                                                                                                                                          tGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAagag                                                                                                                                                                                                        >  6:520886/1‑92 (MQ=255)

                                                                                                                                                                                                                                                          tGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCt                                                                                                                              >  2:725328/1‑166 (MQ=255)

                                                                                                                                                                                                                                                            aCAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGAt                                                                                                                                                                                                              >  2:661125/1‑84 (MQ=255)

                                                                                                                                                                                                                                                            aCAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAagaga                                                                                                                                                                                                       >  3:186098/1‑91 (MQ=255)

                                                                                                                                                                                                                                                            aCAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAAtc                                                                                                                                             >  7:873526/1‑149 (MQ=255)

                                                                                                                                                                                                                                                              agagCGGGAGAGCNTGATNCTCTCCCGAATGCTGTTTTTTTANTCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTa                                                                                                                          <  4:725326/166‑1 (MQ=255)

                                                                                                                                                                                                                                                               gagCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCctct                                   >  2:585619/1‑252 (MQ=255)

                                                                                                                                                                                                                                                                agCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAag                                                                                                                                                                                                          <  4:661123/84‑1 (MQ=255)

                                                                                                                                                                                                                                                                agCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGttt                                                                                                                                                                                                   <  1:186099/91‑1 (MQ=255)

                                                                                                                                                                                                                                                                 gCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGt                                                                                                                                         <  6:675065/148‑1 (MQ=255)

                                                                                                                                                                                                                                                                  cGGGAGAGCATCACGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTa                                                                                                                                                                           <  5:678938/113‑1 (MQ=255)

                                                                                                                                                                                                                                                                     gagaGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGtt                                                                                                                                    <  7:383597/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                     gagaGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGt                                                                                                                                     >  5:881387/1‑148 (MQ=255)

                                                                                                                                                                                                                                                                     gagaGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGt                                                                                                                                     >  5:524496/1‑148 (MQ=255)

                                                                                                                                                                                                                                                                      agagCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAACTGCGATCCCAcc     >  2:460093/1‑275 (MQ=255)

                                                                                                                                                                                                                                                                          cGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCAt                                                                                                                                >  5:334284/1‑148 (MQ=255)

                                                                                                                                                                                                                                                                          cGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCGGCAACATACCTCCGCCTCTCTGACCCTTACTGAGAACTGCGATCCCACCGTa  >  2:147034/1‑274 (MQ=255)

TTTTCATATTCCTGCCACCAACATGAATCCGGTGATCTTTGCGGGCGATAAACCAGGGCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAACTGCGATCCCACCGTA  >  NC_000913/4269875‑4270414

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: `‑`