BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,042,403	A→T	6.0%	intergenic (+110/+86)	ygfZ → / ← yqfA	iron‑sulfur cluster repair protein, plumbagin resistance/hemolysin III family HylIII inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,042,403	0	A	T	6.0%	640.6 / 23.6	182	intergenic (+110/+86)	ygfZ/yqfA	iron‑sulfur cluster repair protein, plumbagin resistance/hemolysin III family HylIII inner membrane protein
Reads supporting (aligned to +/- strand): ref base A (91/80); new base T (4/7); total (95/87)
Fisher's exact test for biased strand distribution p-value = 3.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-02

Only 100 of 185 total aligned reads displayed.

TACCGGTACGGTACTGGCTGCGGTAAAACTGGAAGATGGTCAGGTCGTGGTACAGGTCGTCATG‑AATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAA  >  NC_000913/3042137‑3042667

             ctgcgtgggggaaaaCTGGAAGATGGTATGGTGGTGGTA‑ATGTGGCGCTGCGATAATGTAATGGGACCGCGTAGGATCTGCCCCGTACTGGACTATTCGAATACATTGCAGATCGAGCCGCTGTCGTATTCGCTCGAAGAGTAAATCTTGTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGAGCGTCTTATCATGTCTGGGATTTGATGCCGTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCgg                                                                                                                                                                                                                                                     <  4:618260/264‑1 (MQ=255)

             cTGGCTGCGGTAAAACTGGAAGATGGTCAGGTCGTGGTACAGGTCGTCATG‑AATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCgg                                                                                                                                                                                                                                                     <  2:43687/275‑1 (MQ=255)

                gCTGCGGTAAAACTGGAAGATGGTCAGGTCGTGGTACAGGTCGTCATG‑AATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTtg                                                                                                                                                                                                                                               >  1:70474/1‑278 (MQ=255)

                     ggTAAAACTGGAAGATGGTCAGGTCGTGGTACAGGTCGTCATG‑AATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCTCCTGATGCGACGCTGGCGGGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCACGTAGAACGGAGACGGCGTGCACGCCGAATCCGGCCGTTGTGCTcc                                                                                                                                                                                                                                         >  1:113641/1‑279 (MQ=255)

                                aaGATGGTCAGGTCGTGGTACAGGTCGGCATG‑AATAACTATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCAGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGcc                                                                                                                                                                                                                                   <  2:302830/274‑1 (MQ=255)

                                                                          atGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATa                                                                                                                                                                                         <  2:97021/275‑1 (MQ=255)

                                                                          atGGAACCGGATAGCATCTGCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTAACATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCg                                                                                                                                                                                                                                                            <  2:307048/208‑1 (MQ=255)

                                                                           tGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATa                                                                                                                                                                                         <  2:488694/274‑1 (MQ=255)

                                                                                   gATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGCCGAGGCCCGTAAATCGGATAAGGCGTTCACGCCGAATTCGGCCGTTGTTCTCCGCTGCCTTAGGCGAAGCTGGCCGGTCTTATCAGGCCTACAAAGGCAAACCCCtttag                                                                                                                                                                                >  4:489616/1‑272 (MQ=255)

                                                                                    aTAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGc                                                                                                                                                                               <  2:732228/275‑1 (MQ=255)

                                                                                            ttCCGCGTACGCGACGATGCGAATACATTGCATATCGGGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGAAGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCaa                                                                                                                                                                       <  2:565523/275‑1 (MQ=255)

                                                                                             tCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATATGGCGTTCACGCCGAATCCCGCAGTTGTGCTCCGATGCCTGATGCGACGCTGACGCTTCTTACCAGCACTACAAAGGCATGCCCAttcaccccccccaa                                                                                                                                                                       >  4:401687/1‑262 (MQ=255)

                                                                                                                        tGCATATCGAGCCGCTGCCGTATTCGCTCGAGGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATa                                                                                                                                                                                                                                                                          <  7:924225/148‑1 (MQ=255)

                                                                                                                        tGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATaa                                                                                                                                                                                                                                                                         <  5:434218/149‑1 (MQ=255)

                                                                                                                        tGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGca                                                                                                                                           >  4:304547/1‑275 (MQ=255)

                                                                                                                         gCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAg                                                                                                                                                                                                                                                                        >  5:421234/1‑149 (MQ=255)

                                                                                                                           atatCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGc                                                                                                                                                                                                                                                                      <  6:261657/149‑1 (MQ=255)

                                                                                                                               cGAGCCGCTGCCGTATTCGCTCGGAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGtt                                                                                                                                                                                                                                                                   >  7:789123/1‑148 (MQ=255)

                                                                                                                                 aGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCa                                                                                                                                                                                                                                                                 >  7:895886/1‑148 (MQ=255)

                                                                                                                                   ccGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACg                                                                                                                                                                                                                                                               <  7:107531/148‑1 (MQ=255)

                                                                                                                                   ccGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACg                                                                                                                                                                                                                                                               <  7:63466/148‑1 (MQ=255)

                                                                                                                                   ccGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACg                                                                                                                                                                                                                                                               >  7:475128/1‑148 (MQ=255)

                                                                                                                                      cTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCg                                                                                                                                                                                                                                                            >  8:807173/1‑148 (MQ=255)

                                                                                                                                        gCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATTCGGCGTGAACGcctta                                                                                                                                                                                                                                                          <  7:557205/148‑4 (MQ=39)

                                                                                                                                         ccGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATc                                                                                                                                                                                                                                                        >  6:778214/1‑149 (MQ=255)

                                                                                                                                                 gCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATcc                                                                                                                                                                                                                                                       >  6:529423/1‑142 (MQ=255)

                                                                                                                                                 gCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATcc                                                                                                                                                                                                                                                       <  5:529423/142‑1 (MQ=255)

                                                                                                                                                 gCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGt                                                                                                                                                                                                                                                 >  8:176303/1‑148 (MQ=255)

                                                                                                                                                 gCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCACCAGACAGTGCCACACCCAACCGCCGAGCACGAAACCg                                                                                                                  >  4:350038/1‑275 (MQ=255)

                                                                                                                                                   tCGAAGAGTAAATCTCTTTATCGCATCAGGTATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAAAAAGTGGCACACACTAACGCCTAGCACGAAGCAGTGCCAg                                                                                                            >  1:180296/1‑279 (MQ=255)

                                                                                                                                                   tCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTtg                                                                                                                                                                                                                                               <  7:807173/148‑1 (MQ=255)

                                                                                                                                                    cGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTtgt                                                                                                                                                                                                                                              <  8:475954/148‑1 (MQ=255)

                                                                                                                                                     gAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGc                                                                                                                                                                                                                                            >  5:866400/1‑149 (MQ=255)

                                                                                                                                                     gAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGc                                                                                                                                                                                                                                            >  6:953092/1‑149 (MQ=255)

                                                                                                                                                      aaGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGc                                                                                                                                                                                                                                            >  8:1001253/1‑148 (MQ=255)

                                                                                                                                                       agagTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGTTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCa                                                                                                             <  2:138776/274‑1 (MQ=255)

                                                                                                                                                           tAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGAt                                                                                                                                                                                                                                      >  6:751360/1‑149 (MQ=255)

                                                                                                                                                             aaTCTCTTTTGCGCATCAGGTATTTATGGCCGGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGAGCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCa                                                                                                      <  2:591792/275‑1 (MQ=255)

                                                                                                                                                             aaTCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGc                                                                                                                                                                                                                                    >  6:460576/1‑149 (MQ=255)

                                                                                                                                                             aaTCTCTTTATCGCATCAGGAATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCa                                                                                                      <  2:489616/275‑1 (MQ=255)

                                                                                                                                                               tctcTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCt                                                                                                                                                                                                                                  >  6:220125/1‑149 (MQ=255)

                                                                                                                                                                 tcgttATCGCATCAGGAATTTATCGCCTGATGCGAGGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCACAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGt                                                                                                  <  2:256583/272‑1 (MQ=255)

                                                                                                                                                                    ttATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATTCGGCGTGAACGccttat                                                                                                                                                                                                                                                         <  7:705651/121‑5 (MQ=38)

                                                                                                                                                                    ttATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATTCGGCGTGAACGccttat                                                                                                                                                                                                                                                         >  8:705651/1‑117 (MQ=38)

                                                                                                                                                                       tCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGa                                                                                                                                                                                                                           <  7:465685/148‑1 (MQ=255)

                                                                                                                                                                              aGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATTCGGCGTGAACGccttatc                                                                                                                                                                                                                                                        <  5:528141/112‑6 (MQ=38)

                                                                                                                                                                               ggCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATTCGGCGTGAACGcctta                                                                                                                                                                                                                                                          >  8:748218/1‑106 (MQ=38)

                                                                                                                                                                                gCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGcgcg                                                                                                                                                                                                                 >  5:370988/1‑149 (MQ=255)

                                                                                                                                                                                  aTTTCTCGCCTGATGCGACGCTGGGGGTGCTTATCATGTTTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTg                                                                                 <  4:47785/275‑1 (MQ=255)

                                                                                                                                                                                     tATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCg                                                                                                                                                                                                                                        >  3:96207/1‑121 (MQ=255)

                                                                                                                                                                                     tATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCg                                                                                                                                                                                                                                        <  1:96207/121‑1 (MQ=255)

                                                                                                                                                                                              aTGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGCGCACGAAGCCGTGCCAGATGGCATGGCTGTATGGAATGCGTTTGCAGACGTAGAAAATCa                                                                 >  1:29749/1‑279 (MQ=255)

                                                                                                                                                                                                  gACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACaa                                                                                                                                                                                                >  8:1036669/1‑148 (MQ=255)

                                                                                                                                                                                                  gACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGc                                                                                                                                                                               <  4:752201/165‑1 (MQ=255)

                                                                                                                                                                                                   aCGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAg                                                                                                                                                                                              <  5:595999/149‑1 (MQ=255)

                                                                                                                                                                                                   aCGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAg                                                                                                                                                                                              <  5:220125/149‑1 (MQ=255)

                                                                                                                                                                                                    cGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAg                                                                                                                                                                                              <  8:815126/148‑1 (MQ=255)

                                                                                                                                                                                                      cTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGa                                                                                                                                                                                                                                       <  7:273961/105‑1 (MQ=255)

                                                                                                                                                                                                      cTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGa                                                                                                                                                                                                                                       >  8:273961/1‑105 (MQ=255)

                                                                                                                                                                                                          cgcgTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCATACAAAGGCATAc                                                                                                                                                                                        <  7:74984/148‑1 (MQ=255)

                                                                                                                                                                                                           gcgTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATAccc                                                                                                                                                                                      <  5:687709/149‑1 (MQ=255)

                                                                                                                                                                                                           gcgTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATAccc                                                                                                                                                                                      <  5:26637/149‑1 (MQ=255)

                                                                                                                                                                                                            cgTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATTCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATAccc                                                                                                                                                                                      >  7:741039/1‑148 (MQ=255)

                                                                                                                                                                                                            cgTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATAccc                                                                                                                                                                                      <  7:316525/148‑1 (MQ=255)

                                                                                                                                                                                                               cTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCAtt                                                                                                                                                                                   >  7:434923/1‑148 (MQ=255)

                                                                                                                                                                                                                ttATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTAc                                                                                                                                                                                 <  5:250067/149‑1 (MQ=255)

                                                                                                                                                                                                                   tCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCCAGAAGCCGTGACAGATGGCATGGTTGTATGGAATGCGTTTGCCGACAGAAAAAATCCCCCCGAACGAAAAAAc                                                 >  4:142273/1‑274 (MQ=39)

                                                                                                                                                                                                                    cATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGcc                                                                                                                                                                              >  7:393799/1‑148 (MQ=255)

                                                                                                                                                                                                                      tGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTg                                                                                                                                                                            <  7:208844/148‑1 (MQ=255)

                                                                                                                                                                                                                        tCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGAt                                                                                                          >  4:382450/1‑212 (MQ=255)

                                                                                                                                                                                                                        tCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAAGGCGGTTGGAGACGTAAAAAACCACCCGGAGCGAAAAAACCAcgcc                                           >  4:251875/1‑275 (MQ=255)

                                                                                                                                                                                                                        tCTGGGATTTGATGCCTTTATCAAAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGcc                                                                                                                                                                                                     >  2:642475/1‑121 (MQ=255)

                                                                                                                                                                                                                          tggtttttGTTTCCTTTTGCACACCGNNGACGGNCCCCGTAGGTGTGGTAAGGCGGTCACGCCAATGTCGGCTGTGTTGCTCCGATGCCGGCGGCGACGCTGGCGCTTCTTATCAGGCCTACAAATGCATCCCCATTAAGCCTCCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTAc                                     <  3:29749/274‑1 (MQ=255)

                                                                                                                                                                                                                            ggATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACa                                                                                                                                                                                                 <  4:642475/121‑1 (MQ=255)

                                                                                                                                                                                                                            ggATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAat                                                                                                                                                                      >  7:328833/1‑148 (MQ=255)

                                                                                                                                                                                                                            ggATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCa                                                                                                      <  2:382450/212‑1 (MQ=255)

                                                                                                                                                                                                                                 tGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCCTGGTTGTATGGCATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCg                                   >  2:223827/1‑274 (MQ=255)

                                                                                                                                                                                                                                   aTGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAAt                                                                                                                                                              >  5:636743/1‑149 (MQ=255)

                                                                                                                                                                                                                                      ccTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGaa                                                                                                                                                                                                                                                          <  6:984489/54‑1 (MQ=255)

                                                                                                                                                                                                                                      ccTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGaa                                                                                                                                                                                                                                                          >  5:984489/1‑54 (MQ=255)

                                                                                                                                                                                                                                       cTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGa                                                                                                                                                           <  8:789123/148‑1 (MQ=255)

                                                                                                                                                                                                                                       cTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGAt                                                                                                                                                          >  5:792816/1‑149 (MQ=255)

                                                                                                                                                                                                                                        tttATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATc                                                                                                                                                         <  6:792816/149‑1 (MQ=255)

                                                                                                                                                                                                                                         ttATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATc                                                                                                                                                         >  8:586031/1‑148 (MQ=255)

                                                                                                                                                                                                                                            tCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAg                                                                                                                                                            >  5:767968/1‑142 (MQ=255)

                                                                                                                                                                                                                                            tCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGcc                                                                                                                                                      >  8:143189/1‑148 (MQ=255)

                                                                                                                                                                                                                                                aacaGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCTGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGaaa                                                                                                                                                 <  6:421234/149‑1 (MQ=255)

                                                                                                                                                                                                                                                aacaGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACaa                                                                                                                                                                >  8:298278/1‑134 (MQ=255)

                                                                                                                                                                                                                                                aacaGGGACGAGGCCCGGAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGaaa                                                                                                                                                 >  5:737599/1‑149 (MQ=255)

                                                                                                                                                                                                                                                    gggACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGc                                                                                                                                                       >  8:456120/1‑139 (MQ=255)

                                                                                                                                                                                                                                                    gggACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTgg                                                                                                                                             <  5:968892/149‑1 (MQ=255)

                                                                                                                                                                                                                                                     ggACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTgg                                                                                                                                             <  8:644641/148‑1 (MQ=255)

                                                                                                                                                                                                                                                     ggACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTgg                                                                                                                                             >  7:298502/1‑148 (MQ=255)

                                                                                                                                                                                                                                                     ggACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGTTAACGCTGCccg              >  2:323216/1‑275 (MQ=255)

                                                                                                                                                                                                                                                          aGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGcacaca                                                                                                                                       <  6:262630/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                    aTCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGa                                                                                                                             <  6:332790/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                        gATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCa                                                                                                                          >  8:207319/1‑148 (MQ=255)

                                                                                                                                                                                                                                                                         aTAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCAc                                                                                                                         >  7:1011519/1‑148 (MQ=255)

                                                                                                                                                                                                                                                                          tAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGcac                                                                                                                                          <  6:185258/130‑1 (MQ=255)

                                                                                                                                                                                                                                                                           aaGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGaa                                                                                                                      >  5:34615/1‑149 (MQ=255)

TACCGGTACGGTACTGGCTGCGGTAAAACTGGAAGATGGTCAGGTCGTGGTACAGGTCGTCATG‑AATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAA  >  NC_000913/3042137‑3042667

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