BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,286,313	T→A	5.2%	intergenic (+26/+77)	proL → / ← yejO	tRNA‑Pro/pseudogene, autotransporter outer membrane homology;putative transport; Not classified; putative ATP‑binding component of a transport system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,286,313	0	T	A	5.2%	586.6 / 15.7	175	intergenic (+26/+77)	proL/yejO	tRNA‑Pro/pseudogene, autotransporter outer membrane homology;putative transport; Not classified; putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand): ref base T (81/84); new base A (2/7); total (84/91)
Fisher's exact test for biased strand distribution p-value = 1.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.12e-02

Only 100 of 180 total aligned reads displayed.

ACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGCTGACTGACAAATGACCGCTTAACGTCGGGGTAAACGATGACCTTATACCAGCCTGATAAACGCCACTGGTTCCAGCCACATCATTATTAAATTTGCCATCGTCATTAACTTTCACCTGGTT  >  NC_000913/2286049‑2286578

                aaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGTTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAACATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTGTATATCTGCtat                                                                                                                                                                                                                                                 >  2:263604/1‑275 (MQ=255)

                aaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTTTCGGGGGTCGGAGGTTCAAATCCTCCTGGGCCGAACAAAAATCCCAAGAAAAAACCAACCCTTACGTTTGGCTTTTTTATAACTGCaat                                                                                                                                                                                                                                                 >  4:387481/1‑275 (MQ=255)

                                                 cAAGTGATGAGAGAAGAAAAGCGTTTTAACGCAAACTGGTTAATTAAAAATCAGTTAGCGAAAATAATTAATGGCAAGCGGGGGGGAAAGAGGTAGTAGTAGCGACCACGGGTGGGCACGTAGCGCAGCCTGGAAG‑CACACCGTGAATGGGGGTCGGGGGTCGGAAGTTCAAATCCGCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGc                                                                                                                                                                                                                <  4:211915/275‑1 (MQ=255)

                                                                      cGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAg                                                                                                                                                                                           <  2:508498/275‑1 (MQ=255)

                                                                                   aCTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAAc                                                                                                                                                                                                                                    >  1:58003/1‑221 (MQ=255)

                                                                                                       ttAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATaaa                                                                                                                                                                                                                                     >  4:428815/1‑200 (MQ=255)

                                                                                                           cGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGa                                                                                                                                                                                                                                 <  2:428815/200‑1 (MQ=255)

                                                                                                           cGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGa                                                                                                                                                                                                                           <  1:114303/206‑1 (MQ=255)

                                                                                                                    ttACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCtt                                                                                                                                                                                                                                                                           >  5:253563/1‑149 (MQ=255)

                                                                                                                    ttACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCtt                                                                                                                                                                                                                                                                           >  7:351742/1‑149 (MQ=255)

                                                                                                                    ttACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCtt                                                                                                                                                                                                                                                                           >  6:122426/1‑149 (MQ=255)

                                                                                                                    ttACTTGCAATCGGTGTGG‑AAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTa                                                                                                                                                                                                                                                                          <  7:237745/149‑1 (MQ=255)

                                                                                                                          gccaTCGGTGTTGAAAACGGGAATTATAGGCGGCAACAGTCGGCAACGAGCGCCACCTTGTTGCCGCACCG‑CAAGGGGTGGCCGGGGGCGGAAGTTCAAATCCTTTCGTGGCCACCAAAAAACCCAAGAAAAAACCAACCCTTACggtt                                                                                                                                                                                                                                                                     <  8:383194/146‑1 (MQ=255)

                                                                                                                          gCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggtt                                                                                                                                                                                                                                                                     >  7:13099/1‑149 (MQ=255)

                                                                                                                           cAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTgcgc                                                                                                                                       >  2:587528/1‑274 (MQ=255)

                                                                                                                             aTCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggttggt                                                                                                                                                                                                                                                                  <  8:1024502/149‑1 (MQ=255)

                                                                                                                             aTCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggtt                                                                                                                                                                                                                                                                     <  8:474509/146‑1 (MQ=255)

                                                                                                                              tCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGggttggtt                                                                                                                                                                                                                                                                 <  8:771779/149‑1 (MQ=255)

                                                                                                                              tCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggttggtt                                                                                                                                                                                                                                                                 >  8:187806/1‑149 (MQ=255)

                                                                                                                               cGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttt                                                                                                                                                                                                                                                                <  5:212630/149‑1 (MQ=255)

                                                                                                                                ggTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGtttt                                                                                                                                                                                                                                                               <  6:994846/149‑1 (MQ=255)

                                                                                                                                  tgtgGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGtttttt                                                                                                                                                                                                                                                             <  5:1013926/149‑1 (MQ=255)

                                                                                                                                   gtgGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttttc                                                                                                                                                                                                                                                            <  6:192210/149‑2 (MQ=255)

                                                                                                                                      gAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTtata                                                                                                                                                                                                                                                         >  5:143162/1‑149 (MQ=255)

                                                                                                                                             gtagtaTTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaa                                                                                                                                                                                                                                                  <  7:462627/149‑1 (MQ=255)

                                                                                                                                                 taTTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaattaa                                                                                                                                                                                                                                              >  5:931835/1‑149 (MQ=255)

                                                                                                                                                 taTTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaactaa                                                                                                                                                                                                                                              <  8:523246/149‑1 (MQ=255)

                                                                                                                                                 taTTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTTCACGGCGCTTCTACCCCTGCGCCATTGCCATAGctgctgtctctt                                                                                                                >  1:91104/1‑266 (MQ=255)

                                                                                                                                                 taTTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGATGGGTTTTTTATATCTGCaattaa                                                                                                                                                                                                                                              >  6:226141/1‑149 (MQ=255)

                                                                                                                                                  aTTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaattaat                                                                                                                                                                                                                                             <  8:299143/149‑1 (MQ=255)

                                                                                                                                                   ttAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaattaatt                                                                                                                                                                                                                                            >  8:97698/1‑149 (MQ=255)

                                                                                                                                                    tAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACTCTTACGGTTGGTTTTTTTATATCTGCAATTAATTc                                                                                                                                                                                                                                           <  6:568598/149‑1 (MQ=255)

                                                                                                                                                    tAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTc                                                                                                                                                                                                                                           <  6:748835/149‑1 (MQ=255)

                                                                                                                                                    tAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTc                                                                                                                                                                                                                                           <  8:485018/149‑1 (MQ=255)

                                                                                                                                                     agcagcCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCg                                                                                                                                                                                                                                          >  8:834061/1‑149 (MQ=255)

                                                                                                                                                      gcagcCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGa                                                                                                                                                                                                                                         >  8:883513/1‑149 (MQ=255)

                                                                                                                                                        agcCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCACCTGAGTATTCCACGGCGATTCTACACCTTCCCCATTCCCATAGCTGACAGACAAATGAcgcc                                                                                                         >  4:32214/1‑272 (MQ=255)

                                                                                                                                                          ccACGAGTCGGCACGTAGCGCAGCCTGGTAG‑TGCACCTTCATGGGGTGTCGGGGGTCGGAGGCTCAACTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCCTACGGTTGGTTTTTTTATATCTGCTATTAATTCGATAAACAGACCGTGACACATCACAGCCGGCTTATTTTCTGTTATCAGAACGTCTAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCCACATTGCCATAGCTGACCAACAAATGACCACt                                                                                                        >  2:169042/1‑274 (MQ=255)

                                                                                                                                                          ccACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATaaa                                                                                                                                                                                                                                     <  6:253563/149‑1 (MQ=255)

                                                                                                                                                                   ggCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTg                                                                                                                                                                                                                            >  6:937587/1‑149 (MQ=255)

                                                                                                                                                                         tAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACAt                                                                                                                                                                                                                      <  6:734649/149‑1 (MQ=255)

                                                                                                                                                                         tAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACAt                                                                                                                                                                                                                      <  5:1084682/149‑1 (MQ=255)

                                                                                                                                                                          aGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGCTGACTGACAAATGACCGCTTAACGTCCGGGTAAAc                                                                                        >  2:349982/1‑274 (MQ=255)

                                                                                                                                                                          aGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACAAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACCCATCACCGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGTTTCTACCCCTGCCACATTCCCGCCGTTGACTGACAAACGACCGCCTACCGACGCGGAGAACg                                                                                       >  4:633107/1‑275 (MQ=255)

                                                                                                                                                                             ggacccTGGTTG‑AGCACCGTCATGGTGTGTCGTGGGTTGGAGGTTCAACTCCTCTCGTGGCGATCAAAAATCCCAAGAAAATACCAACCCTTACGTTTGTGTGTTTTATATCTGCAATTAATTCTATAAACAGACCGTGACACAACACAGCCTGTTTATTTTCTGTTATCGGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGAGCAATTGCCATAGCTGACTGACAAATGACCGCTTAACGTCGGGGTAAACGAt                                                                                     <  4:263604/270‑1 (MQ=255)

                                                                                                                                                                             gcAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATcaca                                                                                                                                                                                                                  >  5:89960/1‑149 (MQ=255)

                                                                                                                                                                               aGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTATCGTGACGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCAACCACGCCTGAGTATTCCACGGCGAATCTACCCCTGCGCCATTCCCATTGCTGCCTGACAAATTACCGCCTAACGGCGGGGTaaaggcagcc                                                                                  >  4:431312/1‑268 (MQ=255)

                                                                                                                                                                                gCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGCTGACTGACAAATGACCGCTTAACGTCGGGGTAAACGATGAcc                                                                                 <  2:387481/275‑1 (MQ=255)

                                                                                                                                                                                 ccTGGTAG‑CGCACCGTCGTGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGCTGACTGACAAATGACCGCTTAACGTCGGGGTAACCGATGACCt                                                                                >  1:164889/1‑275 (MQ=255)

                                                                                                                                                                                 ccTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGTAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCt                                                                                                                                                                                                              <  6:305310/149‑1 (MQ=255)

                                                                                                                                                                                 ccTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCg                                                                                                                                                                                                                              >  4:710384/1‑133 (MQ=255)

                                                                                                                                                                                 ccTGGTAG‑AGGACCGGCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGCTGACTGACAAATGACCGCTTAACGTCGGGGTAAACGATGAcc                                                                                 <  2:193994/274‑1 (MQ=255)

                                                                                                                                                                                   tGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGt                                                                                                                                                                                                            <  5:312661/149‑1 (MQ=255)

                                                                                                                                                                                    ggTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGtt                                                                                                                                                                                                           <  7:466152/149‑1 (MQ=255)

                                                                                                                                                                                     gTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGac                                                                                                                                                                                                                          <  2:710384/133‑1 (MQ=255)

                                                                                                                                                                                     gTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAAGTACAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACTCTTATGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGttg                                                                                                                                                                                                          >  6:788245/1‑148 (MQ=255)

                                                                                                                                                                                              ccGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTc                                                                                                                                                                                                    >  5:92246/1‑147 (MQ=255)

                                                                                                                                                                                               cGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGt                                                                                                                                                                                                 <  7:623456/149‑1 (MQ=255)

                                                                                                                                                                                                   aTGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATc                                                                                                                                                                                             <  7:834061/149‑1 (MQ=255)

                                                                                                                                                                                                   aTGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGCTGACTGACAAATGACCGCTTAACGTCGGGGTAAACGATGACCTTATACCAGCCTGATaaa                                                               <  4:587528/275‑1 (MQ=255)

                                                                                                                                                                                                        gtgtCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAAc                                                                                                                                                                                        <  7:950472/149‑1 (MQ=255)

                                                                                                                                                                                                         tgtCGGGGGTCGGAGGTTCAAATCCTCTGGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTATATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGCTGACTGACAAATGACCGCTTAACGTCGGGGTAAACGATGACCTTATACCAGCCTGATAAACGCCa                                                          <  2:319941/274‑1 (MQ=255)

                                                                                                                                                                                                         tgtCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGAAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGCTGACTGACAAATGACCGCTTAACGTCGGGGTAAACGATGACCTTATACCAGCCTGATAAACGCCAc                                                         <  2:633107/275‑1 (MQ=255)

                                                                                                                                                                                                               gggTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGa                                                                                                                                                                                          <  2:150692/140‑1 (MQ=255)

                                                                                                                                                                                                                  tCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACca                                                                                                                                                                              <  5:83879/149‑1 (MQ=255)

                                                                                                                                                                                                                    ggAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGa                                                                                                                                                                                                                           >  7:867652/1‑102 (MQ=255)

                                                                                                                                                                                                                     gAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAt                                                                                                                                                                                                        >  4:197119/1‑120 (MQ=255)

                                                                                                                                                                                                                      aGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACAcc                                                                                                                                                                          <  5:607396/149‑1 (MQ=255)

                                                                                                                                                                                                                       ggTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACAccc                                                                                                                                                                         >  8:28911/1‑149 (MQ=255)

                                                                                                                                                                                                                        gTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCg                                                                                                                                                                        >  5:65080/1‑149 (MQ=255)

                                                                                                                                                                                                                        gTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCg                                                                                                                                                                        >  8:552908/1‑149 (MQ=255)

                                                                                                                                                                                                                           cAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCt                                                                                                                                                                     >  5:732853/1‑149 (MQ=255)

                                                                                                                                                                                                                             aaTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGa                                                                                                                                                                   <  7:187806/149‑1 (MQ=255)

                                                                                                                                                                                                                                  tctcGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTAtt                                                                                                                                                              >  8:522431/1‑149 (MQ=255)

                                                                                                                                                                                                                                   ctcGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTc                                                                                                                                                             <  5:821823/149‑1 (MQ=255)

                                                                                                                                                                                                                                    tcGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGCATTcc                                                                                                                                                            >  5:388964/1‑149 (MQ=255)

                                                                                                                                                                                                                                        gCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACgg                                                                                                                                                        >  7:218994/1‑149 (MQ=255)

                                                                                                                                                                                                                                           gACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGACTATTCCACGGCGa                                                                                                                                                     >  8:684447/1‑149 (MQ=255)

                                                                                                                                                                                                                                             ccAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGAtt                                                                                                                                                   <  5:1024036/149‑1 (MQ=255)

                                                                                                                                                                                                                                              cAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTc                                                                                                                                                  <  6:89960/149‑1 (MQ=255)

                                                                                                                                                                                                                                               aaaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACgg                                                                                                                                                        >  8:759351/1‑142 (MQ=255)

                                                                                                                                                                                                                                               aaaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACgg                                                                                                                                                        <  7:759351/142‑1 (MQ=255)

                                                                                                                                                                                                                                                aaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTa                                                                                                                                                >  6:403281/1‑149 (MQ=255)

                                                                                                                                                                                                                                                 aaaTCGCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCTGCCTGAGTATTCCGCGGCGATTCTAc                                                                                                                                               >  8:670531/1‑149 (MQ=255)

                                                                                                                                                                                                                                                  aaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTAcc                                                                                                                                              >  5:55819/1‑149 (MQ=255)

                                                                                                                                                                                                                                                   aTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTAccc                                                                                                                                             <  7:552908/149‑1 (MQ=255)

                                                                                                                                                                                                                                                      ccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATc                                                                                                                                                                                             >  7:539171/1‑98 (MQ=255)

                                                                                                                                                                                                                                                       cAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGt                                                                                                                                                                                                                             <  5:980889/65‑1 (MQ=255)

                                                                                                                                                                                                                                                       cAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTc                                                                                                                                                             <  6:82067/129‑1 (MQ=255)

                                                                                                                                                                                                                                                         aaaaaaaaCCAACCCTTACGGTTGGTTTTTTTATTACTGCAATTAATTCGATAAACAGCCCCTGAAACATCACAGCCNTTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGCTGACTGACAAATGACCGCTTAACGTCGGGGTAAACGATGACCTTATACCAGCCTGATAAACGCCACTGGTTCCAGCCACATCATTATTAAATTTGCCATCGTCATTAACTTTCa         <  3:164889/273‑1 (MQ=255)

                                                                                                                                                                                                                                                          taaaaaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCTATAACCCTACCGTTACACATCACAGCCTGTTTATTTTCTGTTCTCAGAACGTCCATACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGCTGACTGACAAATGACCGCTTAACGTCGGGGTAAACGATGACCTTATACCAGCCTGATAAACGCCACTGGTTCCAGCCACATCATTATTAAATTTGCCATCGTCATTAACTTTCAc        <  1:224/274‑1 (MQ=255)

                                                                                                                                                                                                                                                            aaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGt                                                                                                                                                                                                            <  7:497040/77‑1 (MQ=37)

                                                                                                                                                                                                                                                             aaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTTTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCAtt                                                                                                                                   >  7:173208/1‑149 (MQ=255)

                                                                                                                                                                                                                                                             aaaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCAtt                                                                                                                                   >  6:641598/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                accaaccCATACGGTTGTTTTTTTTATATCTGCAATTAATTCGATAACCAGCCCGTGACACATCCCAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCAGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGCTGACTGACAAATGACCGCTTAACGTCGGGGTAAACGATGACCTTATACCAGCCTGATAAACGCCACTGGTTCCAGCCACATCATTATTAAATTTGCCATCGTCATTAACTTTCACCTGGtt  <  4:125236/275‑1 (MQ=255)

                                                                                                                                                                                                                                                                  caaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCAt                                                                                                                              >  6:692745/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                  caaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCAt                                                                                                                              <  6:554343/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                    accCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACAccc                                                                                                                                                                         >  2:608856/1‑104 (MQ=255)

                                                                                                                                                                                                                                                                       cTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGctg                                                                                                                         >  7:924662/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                        ttACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCt                                                                                                                                                                     <  4:608856/104‑1 (MQ=255)

ACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAG‑CGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGCTGACTGACAAATGACCGCTTAACGTCGGGGTAAACGATGACCTTATACCAGCCTGATAAACGCCACTGGTTCCAGCCACATCATTATTAAATTTGCCATCGTCATTAACTTTCACCTGGTT  >  NC_000913/2286049‑2286578

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