BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,335,654	T→A	5.9%	intergenic (‑64/+40)	eptA ← / ← adiC	lipid A phosphoethanolamine transferase/arginine:agmatine antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,335,654	0	T	A	5.9%	823.8 / 30.8	237	intergenic (‑64/+40)	eptA/adiC	lipid A phosphoethanolamine transferase/arginine:agmatine antiporter
Reads supporting (aligned to +/- strand): ref base T (90/133); new base A (5/9); total (95/142)
Fisher's exact test for biased strand distribution p-value = 7.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.91e-02

Only 100 of 240 total aligned reads displayed.

AGAGCTTAGTGTCAGGACAATATTAATCACGCTGAAAGCGACGACCGGCATCGACAAGAAAACCAGTACGTTATGCAGCGAATCCAGCGGCAGCGCCTGCAACACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATGCAGTAGAGGAAGGCAATGGTAGTAACTGCCAGATATGCCGGGCGTGCTTTACCAAAGTGACCGTGTCCGAGCAGCAGTAACGCCGCA  >  NC_000913/4335381‑4335925

agagCTTAGTGTCAGGACAATATTAATCACGCTGAAAGCGACGACCGGCATCGACAAGAAAACCAGTACGTTATGCAGCGAATCCAGCGGCAGCGCCTGCAACACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCCACATCACGGTGTTTCCATCGAACAAAGTGCGCATCTGCTCTCAAAGCAAAATTAAGCCAACCTTAAttt                                                                                                                                                                                                                                                                                >  2:587407/1‑275 (MQ=255)

                                        acgacCGGCATCGACAAGAAAACCAGTACGTTATGCAGCGAATCCAGCGGCAGCGCCTGCAACACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAAAAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAAtt                                                                                                                                                                                                                                        >  2:154887/1‑275 (MQ=255)

                                                                        aTGCAGCGAATCCAGCGGCAGCGCCTGCAACACCTGTTTAAAAAAGGCAAGATTCCGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATAc                                                                                                                                                                                                        <  2:291703/275‑1 (MQ=255)

                                                                              cGAATCCAGCGGCAGCGCCTGCAACACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGttt                                                                                                                                                                                                  <  2:439251/275‑1 (MQ=255)

                                                                                aaTCCAGCGGCAGCGCCTGCAACACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTAt                                                                                                                                                                                                                           >  4:441540/1‑248 (MQ=255)

                                                                                                  gCAAAACCTGTTTAAAAAAGGCAAAATTCAGGCAAATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAAtt                                                                                                                                                                              <  4:391593/275‑1 (MQ=255)

                                                                                                                  aaaGGCAAGATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCg                                                                                                                                                              <  2:300076/275‑1 (MQ=255)

                                                                                                                   aaGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCgg                                                                                                                                                             <  1:94584/275‑1 (MQ=255)

                                                                                                                                    gATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAg                                                                                                                                                                                                                                                                          >  6:1053019/1‑149 (MQ=255)

                                                                                                                                       agagATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGtt                                                                                                                                                                                                                                                                       <  5:154213/149‑1 (MQ=255)

                                                                                                                                       agagATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGtt                                                                                                                                                                                                                                                                       >  7:175890/1‑149 (MQ=255)

                                                                                                                                        gagaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGAACTTAAGGTTg                                                                                                                                                                                                                                                                      <  5:338419/149‑1 (MQ=255)

                                                                                                                                        aagaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTg                                                                                                                                                                                                                                                                      <  6:791526/148‑1 (MQ=255)

                                                                                                                                          gaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGc                                                                                                                                                                                                                                                                    <  6:917715/149‑1 (MQ=255)

                                                                                                                                          gaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGAACTTAAGGTTGGc                                                                                                                                                                                                                                                                    <  6:1008589/149‑1 (MQ=255)

                                                                                                                                               aaaaaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTa                                                                                                                                                                                                                                                                            >  4:731882/1‑136 (MQ=255)

                                                                                                                                               aaaaaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGc                                                                                                                                                                                                                                                        >  2:691507/1‑156 (MQ=255)

                                                                                                                                                aaagagggCCGACAATAACCAGGCGAGTGAATTGTAAGAGGATCTTTTTAGTAGGAACTTCAAAATCACTGGGTTTCCATCGAACAAAGTGCGCATATGCTCGCGAAGCGAAATTAAGCCAACCTTAAGTGCTGAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAAGCTTTGCTTATGGGTGCATCTAAGGGATACGAGTTTTTATGTAGCCGGTTGTAATTCGGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGAc                                                                                                                                   <  4:489592/266‑1 (MQ=255)

                                                                                                                                                  aaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGt                                                                                                                                                                                                                                                            >  5:890556/1‑149 (MQ=255)

                                                                                                                                                    gCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGttt                                                                                                                                                                                                                                                          <  5:1004932/149‑1 (MQ=255)

                                                                                                                                                       gCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGATGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCTGGGTGACAAATGACCACATAAcc                                                                                                                          >  2:367816/1‑274 (MQ=255)

                                                                                                                                                         taacaaTAGCCAGGCGAGCAAATTCAAAGAGGGTCTTTTTAGTCGTCGCTTAAACACCTCGGTGTCTCCATAGAACAATGTTTGCATATGCTCACAAGACAAAATTAACCCACCCTTACGTTCTTAAGATTGGCATTTATGTTTGCTgg                                                                                                                                                                                                                                                     >  6:1104958/2‑149 (MQ=255)

                                                                                                                                                         caacaaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTgg                                                                                                                                                                                                                                                     >  5:181869/1‑149 (MQ=255)

                                                                                                                                                          aacaaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGa                                                                                                                                                                                                                                                    <  6:590466/149‑1 (MQ=255)

                                                                                                                                                               tAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAg                                                                                                                                                                                                                                               <  7:1064145/149‑1 (MQ=255)

                                                                                                                                                               tAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAg                                                                                                                                                                                                                                               <  5:423683/149‑1 (MQ=255)

                                                                                                                                                               tAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATATTTGCTGGATTAAg                                                                                                                                                                                                                                               <  7:82771/149‑1 (MQ=255)

                                                                                                                                                                  ccAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGt                                                                                                                                                                                                                       >  2:78273/1‑170 (MQ=255)

                                                                                                                                                                   cAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGa                                                                                                                                                                                                                                           >  8:184240/1‑149 (MQ=255)

                                                                                                                                                                    aGGCGAGTAAATTCAAAGAGGGTCTTTTTTGTAGGCGTTTCCACATCACGGTGTTTCCATCGAACAAATTTCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGTTCATATGATTGCTGGCTTAAGGGGaa                                                                                                                                                                                                                                          >  8:508840/1‑149 (MQ=255)

                                                                                                                                                                      gCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGAt                                                                                                                                                                                                                                                   >  1:24363/1‑138 (MQ=255)

                                                                                                                                                                      gCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCAt                                                                                                                                                                                                                   <  4:78273/170‑1 (MQ=255)

                                                                                                                                                                      gCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCGGCGCATACATGGCGGTGATGACCATCAGGGTGACAATGGCCCCCATAACCCCTTTCGCTCCGGAcc                                                                                                          >  4:605487/1‑275 (MQ=255)

                                                                                                                                                                          gTAAATTTACAGAGGGTCTTTTTAGTAGTCGCTTCAACATAACGGTTCTTCCATCGAACAAAGGGCGCATAGGCTCTCAAAGCAAAATTAAGCCAAACTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAg                                                                                                                                                                                                                                               <  3:24363/138‑1 (MQ=255)

                                                                                                                                                                            aaaTTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCtt                                                                                                                                                                                                                                  <  6:496323/149‑1 (MQ=255)

                                                                                                                                                                             aaTTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCttt                                                                                                                                                                                                                                 <  7:358414/149‑1 (MQ=255)

                                                                                                                                                                              aTTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTg                                                                                                                                                                                                                                <  6:328516/149‑1 (MQ=255)

                                                                                                                                                                               ttCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGc                                                                                                                                                                                                                               >  6:693410/1‑149 (MQ=255)

                                                                                                                                                                               ttCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTAt                                                                                                                                                                                                                           >  2:180865/1‑153 (MQ=255)

                                                                                                                                                                                     gagGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTg                                                                                                                                                                                                                         <  8:435378/149‑1 (MQ=255)

                                                                                                                                                                                            tttttAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATc                                                                                                                                                                                                                  >  8:371724/1‑149 (MQ=255)

                                                                                                                                                                                            tttttAGTAGGCGCTTCAACACCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTgg                                                                                                                                                                                                                                                     <  5:551001/114‑1 (MQ=255)

                                                                                                                                                                                              tttAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGa                                                                                                                                    <  4:138554/225‑1 (MQ=255)

                                                                                                                                                                                                tagtagTCGCTGCAACATCTCGGTGTTTCCATCGAACAAAGTGCGCATATGCTCTCAAAGCAAAATTAAGCCAACCTGAAGTTCTTAAGGTTGGCTGTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACTGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATGCAGTAGAgg                                                                                <  1:188097/275‑1 (MQ=255)

                                                                                                                                                                                                  gtagGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAggg                                                                                                                                                                                                            >  6:860927/1‑149 (MQ=255)

                                                                                                                                                                                                  gtagGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCa                                                                                             >  2:423927/1‑260 (MQ=255)

                                                                                                                                                                                                   tagGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTa                                                                                                                                                                                                                                                                            >  1:61686/1‑84 (MQ=255)

                                                                                                                                                                                                    agGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTACGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCNTATTGGTGCATCTAAGGGATATGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATTGCGGTGATNACCATCAGGGTGACAAATGACCACATAACCTCTTTCGATCCGGACNCNACCCCCGCCCATATGCAATaaaagcaa                                                                             >  3:69709/1‑267 (MQ=255)

                                                                                                                                                                                                      gcgcTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATAc                                                                                                                                                                                                        <  5:1028583/149‑1 (MQ=255)

                                                                                                                                                                                                      gcgcTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATg                                                                                         <  4:423927/260‑1 (MQ=255)

                                                                                                                                                                                                       cgcTTCAACATCACGGTGTTTGCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATGCAGTAGAGGAAGGCa                                                                          <  2:117320/274‑1 (MQ=255)

                                                                                                                                                                                                        gcTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAgcgc                                                                                                                                                                        <  4:110009/179‑1 (MQ=255)

                                                                                                                                                                                                          ttCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATa                                                                                                                                                                                                         <  4:292976/144‑1 (MQ=255)

                                                                                                                                                                                                          ttCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGt                                                                                                                                                                                                    >  6:203374/1‑149 (MQ=255)

                                                                                                                                                                                                           tCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGtt                                                                                                                                                                                                   <  7:380768/149‑1 (MQ=255)

                                                                                                                                                                                                              aCATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGttttt                                                                                                                                                                                                <  7:958800/149‑1 (MQ=255)

                                                                                                                                                                                                                aTCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTAt                                                                                                                                                                                                                                                              >  6:975531/1‑85 (MQ=255)

                                                                                                                                                                                                                aTCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTAt                                                                                                                                                                                                                                                              <  5:975531/85‑1 (MQ=255)

                                                                                                                                                                                                                  cACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGt                                                                                                                                                                                            >  5:1089445/1‑149 (MQ=255)

                                                                                                                                                                                                                   aCGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGc                                                                                                                                                                                                                                                        >  6:738703/1‑88 (MQ=255)

                                                                                                                                                                                                                   aCGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAAt                                                                                                                                                                                                                                     <  8:832984/107‑1 (MQ=255)

                                                                                                                                                                                                                     ggTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGc                                                                                                                                                                                         >  8:821925/1‑149 (MQ=255)

                                                                                                                                                                                                                      gtgtTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGt                                                                                                                                                                                                                       <  8:254397/118‑1 (MQ=255)

                                                                                                                                                                                                                      gtgtTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGt                                                                                                                                                                                                                       >  7:254397/1‑118 (MQ=255)

                                                                                                                                                                                                                         tttCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGttt                                                                                                                                                                                                  <  6:662849/136‑1 (MQ=255)

                                                                                                                                                                                                                          ttCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGtt                                                                                                                                                                                    <  8:77417/149‑1 (MQ=255)

                                                                                                                                                                                                                              aTCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTaa                                                                                                                                                                                <  5:735683/149‑1 (MQ=255)

                                                                                                                                                                                                                              aTCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTaa                                                                                                                                                                                >  7:714918/1‑149 (MQ=255)

                                                                                                                                                                                                                              aTCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTaa                                                                                                                                                                                <  6:563392/149‑1 (MQ=255)

                                                                                                                                                                                                                               tCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAgg                                                                                                                                                                                                                                                                         >  5:644573/1‑59 (MQ=255)

                                                                                                                                                                                                                               tCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGAGATACGGGTTTTTATGTAGCCGGTTGTAAt                                                                                                                                                                               <  5:812473/149‑1 (MQ=255)

                                                                                                                                                                                                                                 gAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATGCAGTAGAGGAAGGCAATGGTAGTAACTGCCAGATATGCCggg                                               >  2:452146/1‑275 (MQ=255)

                                                                                                                                                                                                                                    cAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGTGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAgc                                                                                                                                                                          <  8:714918/149‑1 (MQ=255)

                                                                                                                                                                                                                                         tGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATa                                                                                                                                                                     <  7:781283/149‑1 (MQ=255)

                                                                                                                                                                                                                                               tatGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGc                                                                                                                                                               <  7:576328/149‑1 (MQ=255)

                                                                                                                                                                                                                                                    tCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACAt                                                                                                                                                                  <  7:1021501/141‑1 (MQ=255)

                                                                                                                                                                                                                                                    tCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtga                                                                                                                                                          <  5:1053019/149‑1 (MQ=255)

                                                                                                                                                                                                                                                    tCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtga                                                                                                                                                          <  8:461374/149‑1 (MQ=255)

                                                                                                                                                                                                                                                     cGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTCGTGCATCCAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATTACCACATAACCTCTTTCGCTCCGGACCCCACCACGCCCCAGATGCCGCAGAGTAAGGCAAGGGTAGTAACTCCCAAATATCCCGGCCTTCCTTTACCaaacgcac                            >  4:22387/1‑269 (MQ=255)

                                                                                                                                                                                                                                                      gcaaagcaaaATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtgatg                                                                                                                                                        >  6:480500/1‑149 (MQ=255)

                                                                                                                                                                                                                                                       caaagcaaaATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtgatga                                                                                                                                                       <  5:682783/149‑1 (MQ=255)

                                                                                                                                                                                                                                                          aggaaaaTTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTTGTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCACCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCa                                                                                                                                                    <  8:739469/146‑1 (MQ=255)

                                                                                                                                                                                                                                                            caaaATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtgatga                                                                                                                                                       >  8:127662/1‑144 (MQ=255)

                                                                                                                                                                                                                                                            caaaATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATc                                                                                                                                                  <  8:303993/149‑1 (MQ=255)

                                                                                                                                                                                                                                                            caaaATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATGCAGTAGAGGAAGGCAATGGTAGTAACTGCCAGATATGCCGGGCGTGCTTTACCAAAGTGACCGTGTCCg                    <  2:409576/275‑1 (MQ=255)

                                                                                                                                                                                                                                                             aaaaTTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCa                                                                                                                                                 <  5:323950/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                aTTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAggg                                                                                                                                              >  7:196860/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                 ttAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtgatga                                                                                                                                                       <  6:475074/139‑1 (MQ=255)

                                                                                                                                                                                                                                                                 ttAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtgatga                                                                                                                                                       >  5:475074/1‑139 (MQ=255)

                                                                                                                                                                                                                                                                 ttAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGt                                                                                                                                             <  8:587283/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                 ttAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGt                                                                                                                                             <  6:1038802/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                  aaaGCCAACCTTAAGAACTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTg                                                                                                                                            >  6:255068/2‑149 (MQ=255)

                                                                                                                                                                                                                                                                    aGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGAc                                                                                                                                          >  5:611526/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                    aGCCAACCTTAAGTGCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATGCAGTAGAGGAAGGCAATGGTAGTAACTGCCAGATATGCCGGGCGTGCTTTACCAAAGTGACCGTGTCCGAGcagcag            <  2:605487/275‑1 (MQ=255)

                                                                                                                                                                                                                                                                       cAACCTTAAGAACTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACaaa                                                                                                                                       >  5:705607/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                        aaCCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACTAAt                                                                                                                                      <  5:949474/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                        aaCCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAAt                                                                                                                                      <  8:111192/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                        aaCCTTAAGAACTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGtt                                                                                                                                                                                                   >  7:115701/1‑88 (MQ=255)

                                                                                                                                                                                                                                                                             tAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACca                                                                                                                                 <  6:173972/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                              aaGTTCTTAATGTTGGCTGTTATGTTTGCTGGATTAAGGGGGATTAATCTTTTCTTATTGGTGCATCTAAGGGATACGGGTTTTGATGTAGCCGGTTGGAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATGCAGTAGAGGAAGGCAATGGTAGTAACTGCCAGATATGCCGGGCGTGCTTTACCAAAGTGACCGTGTCCGAGCAGCAGTAACGCCGca  <  2:22387/275‑1 (MQ=255)

AGAGCTTAGTGTCAGGACAATATTAATCACGCTGAAAGCGACGACCGGCATCGACAAGAAAACCAGTACGTTATGCAGCGAATCCAGCGGCAGCGCCTGCAACACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAACCTCTTTCGCTCCGGACCCCACCACGGCCCAGATGCAGTAGAGGAAGGCAATGGTAGTAACTGCCAGATATGCCGGGCGTGCTTTACCAAAGTGACCGTGTCCGAGCAGCAGTAACGCCGCA  >  NC_000913/4335381‑4335925

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