BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,270,144	G→C	5.4%	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.4%	384.0 / 11.3	112	intergenic (+17/‑94)	aphA/yjbQ	acid phosphatase/phosphotransferase, class B, non‑specific/thiamine phosphate synthase
Reads supporting (aligned to +/- strand): ref base G (52/54); new base C (2/4); total (54/58)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.14e-01

Only 100 of 115 total aligned reads displayed.

TTTTCATATTCCTGCCACCAACATGAATCCGGTGATCTTTGCGGGCGATAAACCAGGGCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAACTGCGATC  >  NC_000913/4269875‑4270406

ttttCATATTCCTGCCACCAACATGAATCCGGTGATCTTTGCGGGCGATAAACCAGGGCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGct                                                                                                                                                                                                                                                                   <  2:373493/275‑1 (MQ=255)

                         attCCCGTGATCTTTGCGGGCGATAAACCGGTGCAAAATACCAAATCGCAATGGCTGCAGGATAACAATATCCGCATTTTTTCTGGCGATTCTGATAATGATATTACCGCCGCCCGCGATGTCGGCGCTCGTGGTCTCCGCATTCTGCGCGCCTCCAACTCTACCTACGAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGATGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATca                                                                                                                                                                                                                                          <  4:193746/273‑1 (MQ=255)

                                                      aGGGCAAAATACAAAATCGCAATGGCTGAAGGATAAAAATATCCGAATTTTTGATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCa                                                                                                                                                                                                             <  1:140915/275‑1 (MQ=255)

                                                        ggCAAAATACCAAATCGCAATGGCTGCAGGATAAAAATATCCGCATGGTTTATGGCGATTCTGATACTGATATTACCGCCCCACGCAATGTCGCCCCTCGTGCTATCCGCATTCTGCGCGCCTCCCACTCTACCTACCACCCCCTCCCACAAGCGGGTGCGGTTGGTGAGGCGGTTATCTTCAATTCAACATACTTACA‑CGCGCGGCGAGCGGGATGACCCCCCGAATCGTTTTTTGTTAAGCTACCCTTTTCCTTTTCCCTGTCTTGCTGCaaa                                                                                                                                                                                                           >  4:299313/1‑275 (MQ=255)

                                                        ggCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCaaa                                                                                                                                                                                                           <  2:351777/275‑1 (MQ=255)

                                                          aacaagataaaaACGTAATGGGTTCAGNNAAAAAATAGCCGATTTGTTTTGTGCGATTCNAAAAAGTCTCTTACGGCCGCAGGCCATGTTGGCGCTCGTGGTATCCCCATTCTGCTCGTCTCCAACTCTACCTANAAACCCTTGCCACATGCGGGTGCGTTTGGTGAAGATGTGTTCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTGTCGCTGTCTTGCTGCaaa                                                                                                                                                                                                           <  3:69658/265‑1 (MQ=255)

                                                             aagaCAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGAt                                                                                                                                                                                                      <  2:457951/272‑1 (MQ=255)

                                                                                              ttacgagtgttttATGGTGATTCTGGTAAGGATCTTACCGCCGCCCGCGATGGCGGGGCTCGGGGTAGCCGGATTCGGCGCGCCACCAGCGCCACCTACACACCCTGGCAACAAGCGGGTGCGTATGGTGAAGAAGGGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTGTTTCCTTTCGCTTTCTTGGTGCAAACTGATTAAGAGACTTTTATCAAGGGGCAGCACAtgtg                                                                                                                                                                      <  4:184402/265‑1 (MQ=255)

                                                                                               tCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGa                                                                                                                                                                                                                                                                       >  1:100206/1‑176 (MQ=255)

                                                                                                   aaTTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGct                                                                                                                                                                                                                                                                   <  3:100206/176‑1 (MQ=255)

                                                                                                                  ttCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCt                                                                                                                                                 <  2:183210/275‑1 (MQ=255)

                                                                                                                         aaTGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTg                                                                                                                                                                                                                                                                        >  5:115576/1‑149 (MQ=255)

                                                                                                                          aTGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGa                                                                                                                                                                                                                                                                       >  8:33860/1‑149 (MQ=255)

                                                                                                                              tatTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGct                                                                                                                                                                                                                                                                   <  7:23265/149‑1 (MQ=255)

                                                                                                                                 tACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCg                                                                                                                                  >  2:205562/1‑275 (MQ=255)

                                                                                                                                  aCCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCATCACGCTCTcc                                                                                                                                                                                                                                                               <  5:302755/149‑1 (MQ=255)

                                                                                                                                  aCCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCATCACGCTCTcc                                                                                                                                                                                                                                                               <  6:331684/149‑1 (MQ=255)

                                                                                                                                     gccgcACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCaaa                                                                                                                                                                                                           <  2:490440/198‑1 (MQ=255)

                                                                                                                                      ccgcACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTa                                                                                                                                                                                                                                              >  1:132653/1‑162 (MQ=255)

                                                                                                                                       cgcACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGGGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGtt                                                                                                                            >  2:270096/1‑275 (MQ=255)

                                                                                                                                         cACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGc                                                                                                                                                                                                                                                        >  8:260911/1‑149 (MQ=255)

                                                                                                                                          aCGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATca                                                                                                                                                                                                                                          <  3:132653/162‑1 (MQ=255)

                                                                                                                                                       gcTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTTAAGAGGTGATCGTCAATTCAGAATACTGACAGAGAG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATca                                                                                                                                                                                                                                          >  8:389088/1‑149 (MQ=255)

                                                                                                                                                         tCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATcaca                                                                                                                                                                                                                                        >  5:301717/1‑149 (MQ=255)

                                                                                                                                                         tCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATcaca                                                                                                                                                                                                                                        <  8:221132/149‑1 (MQ=255)

                                                                                                                                                                 tCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCATGCCACAAACGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCacac                                                                                                                                                                                                                                       <  7:409626/142‑1 (MQ=255)

                                                                                                                                                                 tCCGCATTCGGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCacac                                                                                                                                                                                                                                       >  8:409626/1‑142 (MQ=255)

                                                                                                                                                                    gCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATc                                                                                                                                                                                                                                <  5:32467/146‑1 (MQ=255)

                                                                                                                                                                    gCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATc                                                                                                                                                                                                                                >  6:32467/1‑146 (MQ=255)

                                                                                                                                                                     cATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCttt                                                                                                                                                                                                                            <  7:223924/149‑1 (MQ=255)

                                                                                                                                                                      aTTCTGAGCGCCTCCAACTCGACCTACAAACCCTTGCCACAAGCGGGGGCGTTTGGGGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGGGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAg                                                                                             <  4:270092/275‑1 (MQ=255)

                                                                                                                                                                         cTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCt                                                                                                                                                                                                                        <  6:115576/149‑1 (MQ=255)

                                                                                                                                                                                  tCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCttt                                                                                                                                                                                                                            >  5:156048/1‑136 (MQ=255)

                                                                                                                                                                                  tCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCt                                                                                                                                                                                                                    >  6:216275/1‑144 (MQ=255)

                                                                                                                                                                                           acctacAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGc                                                                                                                                                                                                                         >  4:509357/1‑130 (MQ=255)

                                                                                                                                                                                           acctacAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCcgcgc                                                                           >  2:370361/1‑272 (MQ=255)

                                                                                                                                                                                             ctacAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTa                                                                                                                                                                                                    >  5:237902/1‑149 (MQ=255)

                                                                                                                                                                                               acAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTc                                                                                                                                                                                                                     <  2:509362/130‑1 (MQ=255)

                                                                                                                                                                                               acAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGagca                                                                                                                                                                               <  4:206804/168‑1 (MQ=255)

                                                                                                                                                                                               acAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACa                                                                    <  2:299317/275‑1 (MQ=255)

                                                                                                                                                                                                cAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCaaa                                                                                                                                                                                                           >  6:317262/1‑139 (MQ=255)

                                                                                                                                                                                                cAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCaaa                                                                                                                                                                                                           <  5:317262/139‑1 (MQ=255)

                                                                                                                                                                                                cAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAaga                                                                                                                                                                                                 <  7:17181/149‑1 (MQ=255)

                                                                                                                                                                                                  aaCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAagaga                                                                                                                                                                                               <  8:207727/149‑1 (MQ=255)

                                                                                                                                                                                                  aaCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAagaga                                                                                                                                                                                               <  7:260911/149‑1 (MQ=255)

                                                                                                                                                                                                    cccTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGt                                                                                                                                                                                             >  8:398489/1‑149 (MQ=255)

                                                                                                                                                                                                      cTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTAt                                                                                                                                                                                                                                 >  6:450540/1‑111 (MQ=255)

                                                                                                                                                                                                      cTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTAt                                                                                                                                                                                                                                 <  5:450540/111‑1 (MQ=255)

                                                                                                                                                                                                      cTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGttt                                                                                                                                                                                           >  5:82448/1‑149 (MQ=255)

                                                                                                                                                                                                      cTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGttt                                                                                                                                                                                           >  8:296341/1‑149 (MQ=255)

                                                                                                                                                                                                        tGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGTTTAAGAGAGTTTTa                                                                                                                                                                                         <  7:33860/149‑1 (MQ=255)

                                                                                                                                                                                                        tGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTa                                                                                                                                                                                         <  8:454280/149‑1 (MQ=255)

                                                                                                                                                                                                        tGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTa                                                                                                                                                                                         <  6:301717/149‑1 (MQ=255)

                                                                                                                                                                                                            acaAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCaa                                                                                                                                                                                     >  6:6686/1‑149 (MQ=255)

                                                                                                                                                                                                            acaAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCaa                                                                                                                                                                                     >  6:252315/1‑149 (MQ=255)

                                                                                                                                                                                                               aGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTCATCCTTTCGCTGTCTTGCTGCATACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGACGAAATTCTGAATCAGCTGGCTGATATCCCGCGCGTTAACTTCGGCTCACTGCAtt                                                     >  2:506210/1‑273 (MQ=255)

                                                                                                                                                                                                                gCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGag                                                                                                                                                                                 <  5:211822/149‑1 (MQ=255)

                                                                                                                                                                                                                   ggTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGt                                                                                                                                                                                                                                                     >  1:15668/1‑78 (MQ=255)

                                                                                                                                                                                                                   ggTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGAtt                                                                                                                                                                                                     >  5:421809/1‑126 (MQ=255)

                                                                                                                                                                                                                   ggTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGAtt                                                                                                                                                                                                     <  6:421809/126‑1 (MQ=255)

                                                                                                                                                                                                                     tGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTTACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCCCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCCCGTTAACATCGGCTTACTGCATATGTtgc                                               >  1:139943/1‑274 (MQ=255)

                                                                                                                                                                                                                       cGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGttttt                                                                                                                                                                                                                                                 <  3:15668/78‑1 (MQ=255)

                                                                                                                                                                                                                          ttGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAAc                                         <  2:114587/275‑1 (MQ=255)

                                                                                                                                                                                                                                aaGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATc                                                                                                                                                                 >  5:257392/1‑149 (MQ=255)

                                                                                                                                                                                                                                        aTCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCACGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGCGTTTTATCAAGGAGCAGCACATGTGGTATCCAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGTCTGATATGCCGCGCGTTAACCGCCGCTTACTGCACCTGTTGCTGCAACATAACTCCGCccct                           >  4:462905/1‑275 (MQ=255)

                                                                                                                                                                                                                                           gTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTa                                                                                                                                                                                         <  7:225211/114‑1 (MQ=255)

                                                                                                                                                                                                                                           gTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTa                                                                                                                                                                                         >  8:225211/1‑114 (MQ=255)

                                                                                                                                                                                                                                           gTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTc                                                                                                                                                      >  7:374961/1‑149 (MQ=255)

                                                                                                                                                                                                                                            tCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTc                                                                                                                                                                                                                           <  8:442854/79‑1 (MQ=255)

                                                                                                                                                                                                                                            tCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTc                                                                                                                                                                                                                           >  7:442854/1‑79 (MQ=255)

                                                                                                                                                                                                                                            tCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCaaaa                                                                                                                                                             <  3:94736/141‑1 (MQ=255)

                                                                                                                                                                                                                                             cAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGa                       <  2:462910/274‑1 (MQ=255)

                                                                                                                                                                                                                                               aTTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGTTGTCTTGCTGCAAACTGATTAAGAGCGTTTTATCAAGGAGCAGCACATGTGGTCTCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTGACGGAAGAAATTCTGAAAAAGCTTGCCGATATGCCGCCCGTTCACATCGGCTTACTTCATCTGTTGCTGCAACATACCTCCGACCCCCTGaca                     >  1:3111/1‑274 (MQ=255)

                                                                                                                                                                                                                                               aTTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGc                                                                                                                                                  >  6:357013/1‑149 (MQ=255)

                                                                                                                                                                                                                                               aTTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTTTTGCTGCAACATACCCACGCCTCTCTgcaa                     >  4:64575/1‑271 (MQ=255)

                                                                                                                                                                                                                                                         cTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAAAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGTTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACAt                                                                   <  2:187600/218‑1 (MQ=255)

                                                                                                                                                                                                                                                          tGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCt                                   >  2:34540/1‑249 (MQ=255)

                                                                                                                                                                                                                                                          tGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGATTACTGCATCTGTTGCTGCAACATACACCCGCCTCTCTGACACTTAATGCGAAc         >  2:114170/1‑275 (MQ=255)

                                                                                                                                                                                                                                                              agagCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAAtctc                                                                                                                                   >  8:111492/1‑149 (MQ=255)

                                                                                                                                                                                                                                                              agagCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGc                               <  4:34540/249‑1 (MQ=255)

                                                                                                                                                                                                                                                               gagCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAACTGCg     >  2:396729/1‑274 (MQ=255)

                                                                                                                                                                                                                                                               gagCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTATTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCg                                                                                                                                  <  7:389088/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                ttcG‑GGAGAGCATCACGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGt                                                                                                                                 >  5:123961/3‑149 (MQ=255)

                                                                                                                                                                                                                                                                agCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGt                                                                                                                                 <  7:111492/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                 gCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTg                                                                                                                                <  6:82448/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                  cG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTAAATAAACTGGCTGATATGCCGCGCGCTAACATCGGCTTACTACTTCTGTTGCTGCAAAATACCTCCGCCGCCCTGACACTGAATgaaacctcggatc  >  4:437408/1‑263 (MQ=255)

                                                                                                                                                                                                                                                                  cG‑GGAGAGCATCACGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTgg                                                                                                                               <  6:123961/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                   g‑ggAGAGCATCACGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAaga                                                                                                                                                                                                 <  7:301334/82‑1 (MQ=39)

                                                                                                                                                                                                                                                                     ggAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGt                                                                                                                             >  5:123457/1‑149 (MQ=255)

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

                                                                                                                                                                                                                                                                       agagCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGttt                                                                                                                           <  5:357013/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                       agagCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATga             >  2:197172/1‑259 (MQ=255)

                                                                                                                                                                                                                                                                        gagCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGc                                                                                                                                                                                                                         >  2:12082/1‑54 (MQ=255)

                                                                                                                                                                                                                                                                         agCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCa                                                                                                                         <  7:168888/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                           cGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTtgc                                                                                                                                                                                                                 >  8:230467/1‑59 (MQ=255)

                                                                                                                                                                                                                                                                           cGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTtgc                                                                                                                                                                                                                 <  7:230467/59‑1 (MQ=255)

                                                                                                                                                                                                                                                                           cGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATCAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTg                                                                                                                                >  1:18467/1‑140 (MQ=255)

                                                                                                                                                                                                                                                                            gTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTc                                                                                                                                                                                                                     <  4:12082/54‑1 (MQ=255)

                                                                                                                                                                                                                                                                              gATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTgg                                                                                                                    >  5:396195/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                              gATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTgg                                                                                                                    <  7:403771/149‑1 (MQ=255)

TTTTCATATTCCTGCCACCAACATGAATCCGGTGATCTTTGCGGGCGATAAACCAGGGCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCG‑GGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAACTGCGATC  >  NC_000913/4269875‑4270406

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