BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,377,792	A→T	6.3%	intergenic (‑70/+19)	blc ← / ← ampC	outer membrane lipoprotein cell division and growth lipocalin/penicillin‑binding protein; beta‑lactamase, intrinsically weak

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,377,792	0	A	T	6.3%	547.4 / 19.3	159	intergenic (‑70/+19)	blc/ampC	outer membrane lipoprotein cell division and growth lipocalin/penicillin‑binding protein; beta‑lactamase, intrinsically weak
Reads supporting (aligned to +/- strand): ref base A (65/83); new base T (4/6); total (70/89)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01

Only 100 of 159 total aligned reads displayed.

TGTTGCGGTGACTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATCCGCCGGTCGCCCCTGTTTTATGTACCCATGATGCGCGTACTGCAGGAGTTGGGGGCGTAATCGCTTTTACGGGGCGTGCTGCCAGTGCAATTTTATTGTCA  >  NC_000913/4377534‑4378036

tGTTGCGGTGACTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGaa                                                                                                                                                                                                                                       <  2:518171/275‑1 (MQ=255)

        tGACTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGCATTTTAATg                                                                                                                                                                                                                               >  2:705893/1‑275 (MQ=255)

           cTTTTTCCAGTCCACGTTAAAAGCGGGTGTCAAAACGGGCAATCCGATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTag                                                                                                                                                                                                                            <  4:595601/275‑1 (MQ=255)

           cTTGTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACGGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTag                                                                                                                                                                                                                            <  2:195163/275‑1 (MQ=255)

                       accattCAAAGCGCGGATAAAAACGGGCAACTTCAACCCCGGAACCAAGATAACGTTTGGGGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGCGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAAt                                                                                                                                                                                                                <  4:107205/271‑1 (MQ=255)

                                    gTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACGGTAGAGCGTTAAGAATCTGCCAGgcggcg                                                                                                                                                                                                   >  2:539540/1‑275 (MQ=255)

                                                                 tACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGAACCCGAt                                                                                                                                                                                                                                           >  1:168469/1‑206 (MQ=255)

                                                                   ccGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACAAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGa                                                                                                                                                                                                                                            <  2:93184/203‑1 (MQ=255)

                                                                       gaTAACGTTTTGGGTCGAAATTATTTACCACGGTAACGCCACGGGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCGGTTGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGc                                                                                                                                                                <  4:78550/275‑1 (MQ=255)

                                                                        aTAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTGTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGcc                                                                                                                                                               <  2:139025/275‑1 (MQ=255)

                                                                                 ggCGTCGAAATTATTTACTACGGTCACGCCACGCGGCTGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACg                                                                                                                                                      <  2:522770/275‑1 (MQ=255)

                                                                                  gCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTCACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACg                                                                                                                                                      <  4:452898/274‑1 (MQ=255)

                                                                                                                cgcgGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAAt                                                                                                                                                                                                                                                      <  6:819571/148‑1 (MQ=255)

                                                                                                                agcggaggcgTAGTGGAACTGCATGCAACGACCAGAAATGCAGCGGTCGGTACGGGAACGAGAGGGGGCAGGCGCATAAATGTTTCCTTACTTGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGGTGGAATTTTACTGTAGGGCGTAAAGAATCTGCCAGGCGGCGTCGACTCTCGCCGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATGC‑CCAGCTCTTTTTCTGGAATAAACGCGa                                                                                                                      <  4:538325/274‑1 (MQ=255)

                                                                                                                 gcggcggcgTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAAtt                                                                                                                                                                                                                                                     <  5:197800/148‑1 (MQ=255)

                                                                                                                 gcggcggcgTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAAtt                                                                                                                                                                                                                                                     <  6:604764/148‑1 (MQ=255)

                                                                                                                 gcggcggcgTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTc                                                                                                                                                                                                                                                    <  7:182188/149‑1 (MQ=255)

                                                                                                                  cggcggcgTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTTTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCg                                                                                                                                                                                                                                                   <  8:64138/149‑1 (MQ=255)

                                                                                                                  cggcggcgTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTc                                                                                                                                                                                                                                                    <  5:326594/148‑1 (MQ=255)

                                                                                                                     aggcgTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTATGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGa                                                                                                                                                                                                                                                 <  6:710989/147‑1 (MQ=255)

                                                                                                                      ggcgTAGGAGAACTGCAGGCAACGACCAGAAATCCAGCTGTCGCTGCGGCAACGCGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTTAGTGTAGATGACAGCAAGGCAAAGCGGAGAAAACGTCCGAAACTTCGGACCCGATGGAATTTTACTGTAGAGCGTTACGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCCTCCTCCCG‑TACACTCTCTCTTTTTCTGGAATCAAcgcgcaatcg                                                                                                                 >  4:456941/1‑268 (MQ=255)

                                                                                                                       gcgTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAcct                                                                                                                                                                                                                                              <  7:238513/149‑2 (MQ=255)

                                                                                                                       gcgTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAcct                                                                                                                                                                                                                                              <  8:236694/149‑2 (MQ=255)

                                                                                                                        cgTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAcct                                                                                                                                                                                                                                              <  6:134353/148‑2 (MQ=255)

                                                                                                                         gTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCg                                                                                                                                                                                                                                             >  5:246962/1‑148 (MQ=255)

                                                                                                                         gTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCg                                                                                                                                                                                                                                             <  6:556312/148‑1 (MQ=255)

                                                                                                                          tAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgttt                                                                                                                                                                  <  2:663517/222‑1 (MQ=255)

                                                                                                                               gaACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGaa                                                                                                                                                                                                                                       >  5:43083/1‑148 (MQ=255)

                                                                                                                                 aCTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAAGGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGAAAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAtt                                                                                                                                                                                                                                     <  6:705785/148‑1 (MQ=255)

                                                                                                                                  cTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAtttt                                                                                                                                                                                                                                   >  8:600583/1‑149 (MQ=255)

                                                                                                                                   tGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTa                                                                                                                                                                                                                                  >  8:176482/1‑149 (MQ=255)

                                                                                                                                    gCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTa                                                                                                                                                                                                                                  >  6:885979/1‑148 (MQ=255)

                                                                                                                                     cAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTAc                                                                                                                                                                                                                                 <  5:792847/148‑1 (MQ=255)

                                                                                                                                       ggCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTCTTTTTCTGGaa                                                                                                                               >  2:676841/1‑243 (MQ=255)

                                                                                                                                         cAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGAt                                                                                                                                                                                                                                           <  5:571956/134‑1 (MQ=255)

                                                                                                                                         cAACGAACAGAAATGCCGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTGTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATCCGCCGGTCGc                                                                                             <  2:60759/275‑1 (MQ=255)

                                                                                                                                          aaCGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTaga                                                                                                                                                                                                                           >  8:549109/1‑149 (MQ=255)

                                                                                                                                          aaCGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTaga                                                                                                                                                                                                                           >  7:322847/1‑149 (MQ=255)

                                                                                                                                          aaCGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGTTAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAGTCCGCCGGTCGc                                                                                             >  2:113413/1‑274 (MQ=255)

                                                                                                                                            cGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATCCGCCGGTCGccc                                                                                           <  2:75360/274‑1 (MQ=255)

                                                                                                                                            cGACCAGAAATGCAGCTGGCGCGGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATCCGCCGGTCGccc                                                                                           <  2:48096/274‑1 (MQ=255)

                                                                                                                                             gACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCg                                                                                                                                                                                                                        >  8:396183/1‑149 (MQ=255)

                                                                                                                                               ccAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGtt                                                                                                                                                                                                                      <  7:396183/149‑1 (MQ=255)

                                                                                                                                                 aGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTa                                                                                                                                                                                                                     >  6:780522/1‑148 (MQ=255)

                                                                                                                                                      tGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAAt                                                                                                                                                                                                                <  5:756026/148‑1 (MQ=255)

                                                                                                                                                       gCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCATGATAC‑CCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATCCGCCGGTCGCCCCTGTTTTATg                                                                                 <  2:150840/273‑1 (MQ=255)

                                                                                                                                                         aGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTg                                                                                                                                                                                                             <  5:885979/148‑1 (MQ=255)

                                                                                                                                                          gCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCTAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATCCGCCGGTCGCCCCTGTTTTATGTa                                                                               >  2:228208/1‑272 (MQ=255)

                                                                                                                                                          gCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATCCGCCGGTCGCCCCTGTTTTATGTAccc                                                                            <  2:370640/275‑1 (MQ=255)

                                                                                                                                                            tGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCa                                                                                                                                                                                                          <  6:246962/148‑1 (MQ=255)

                                                                                                                                                             gTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCTGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATCCGCCGGTCGCCCCTGATTTATGTACCCatg                                                                         >  2:483918/1‑275 (MQ=255)

                                                                                                                                                                 cTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAggcgg                                                                                                                                                                                                     >  5:275563/1‑148 (MQ=255)

                                                                                                                                                                  tGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAggcggc                                                                                                                                                                                                    >  5:72499/1‑148 (MQ=255)

                                                                                                                                                                        aaCGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAcctt                                                                                                                                                                                                                                             >  7:364616/1‑99 (MQ=255)

                                                                                                                                                                         aCGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGAct                                                                                                                                                                                             >  6:39098/1‑148 (MQ=255)

                                                                                                                                                                           gagaGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAcctt                                                                                                                                                                                                                                             >  8:340638/1‑96 (MQ=255)

                                                                                                                                                                           gagaGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAcctt                                                                                                                                                                                                                                             <  7:340638/98‑3 (MQ=255)

                                                                                                                                                                              agGGAGTAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCCCt                                                                                                                                                                                       <  7:600583/149‑1 (MQ=255)

                                                                                                                                                                                     aGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTg                                                                                                                                                                                 <  5:126344/148‑1 (MQ=255)

                                                                                                                                                                                       gcgcATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGa                                                                                                                                                                              >  8:736014/1‑149 (MQ=255)

                                                                                                                                                                                         gcATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGAt                                                                                                                                                                             <  5:780522/148‑1 (MQ=255)

                                                                                                                                                                                            tAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTa                                                                                                                                                                                                                     >  6:425152/1‑105 (MQ=255)

                                                                                                                                                                                            tAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGc                                                                                                                                                                                                            >  7:315342/1‑114 (MQ=255)

                                                                                                                                                                                            tAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGc                                                                                                                                                                                                            <  8:315342/114‑1 (MQ=255)

                                                                                                                                                                                             aaaattgtCCTTACGGGGTTTTTTCTAAGGTGAGATGAAAGAAAGGAAAACCGGAGACAAGGTCGGAAAATTCGGGCCCGATGGAAGTTTACGGTAGAGTGTTAAGAATCTGCCAGTCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATCCGCCGGTCGCCCCTGTTTTATGTACCCATGATGCGCGTACTGCAGGAGTTGGGGGCGTAAt                                          <  4:539538/267‑1 (MQ=255)

                                                                                                                                                                                                 gTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgt                                                                                                                                                                    <  8:160921/149‑1 (MQ=255)

                                                                                                                                                                                                   ttCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgttt                                                                                                                                                                  <  8:391590/149‑1 (MQ=255)

                                                                                                                                                                                                   ttCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgtt                                                                                                                                                                   <  6:93847/148‑1 (MQ=255)

                                                                                                                                                                                                       ttACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGcc                                                                                                                                                               <  5:543168/148‑1 (MQ=255)

                                                                                                                                                                                                         aCTGGTTTTTTTCTAAGTGTAGATGACAGCAGGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATccgc                                                                                                    >  1:191682/1‑204 (MQ=255)

                                                                                                                                                                                                         aCTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCCCGATAC‑CCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATCCGCCGGTCGCCCCTGTTTTATGTACCCATGATTCCGCTACTGCAGGAGTTGGGGGCGTATTCTGTTTTActgct                             >  4:213390/1‑271 (MQ=255)

                                                                                                                                                                                                                tttttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACg                                                                                                                                                      <  6:72777/148‑1 (MQ=255)

                                                                                                                                                                                                                      taagtttaaagtaCAGCAAGGAAAAGCGGAGAAAAGGGCCGAAAATTCGGCCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCATGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATCCGCCGGTCGCCCCTGTTTTATGTACCCATGATGCGCGTACTGCTGGAGTTGGGGGCGTAATCGCTTTTACGGGGCGTGCTGCCAGTg                <  4:69043/263‑1 (MQ=255)

                                                                                                                                                                                                                          tgtAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGActc                                                                                                                                                                                            <  7:255524/100‑1 (MQ=255)

                                                                                                                                                                                                                          tgtAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGctc                                                                                                                                          >  8:83463/1‑149 (MQ=255)

                                                                                                                                                                                                                            tAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTCtt                                                                                                                                        >  8:580640/1‑149 (MQ=255)

                                                                                                                                                                                                                                 gACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTCTTTTTc                                                                                                                                    <  6:752337/148‑1 (MQ=255)

                                                                                                                                                                                                                                 gACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTCTTTTTc                                                                                                                                    <  5:695294/148‑1 (MQ=255)

                                                                                                                                                                                                                                  aCAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTCTTTTTCt                                                                                                                                   <  6:171805/148‑1 (MQ=255)

                                                                                                                                                                                                                                  aCAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATCCGCCGGTCGCCCCTGTTTTATGTACCCATGATGCGCGTACTGCAGGAGTTGGGGGCGTAATCGCtttt                                   >  2:637897/1‑244 (MQ=255)

                                                                                                                                                                                                                                  aCAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATCTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATCCGCCGGTCGCCCCTGTTTTATGTACCCATGATGCGCGTACTGCAGGAGTTGGGGGCGTAATCGATTTTCCGGGGCGTGCTGCCACTGCAATTTTATTGt    >  4:411173/1‑275 (MQ=255)

                                                                                                                                                                                                                                    aGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACt                                                                                                                                                                                                                                >  2:125299/1‑54 (MQ=255)

                                                                                                                                                                                                                                    aGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGActc                                                                                                                                                                                            >  2:43738/1‑90 (MQ=255)

                                                                                                                                                                                                                                    aGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTNTTTTTCTGGAATAAACGCGACATAGCTACCCAATCCGCCgggc                                                                                               >  3:196660/1‑180 (MQ=255)

                                                                                                                                                                                                                                    aGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTCTTTTTCTgg                                                                                                                                 >  6:550993/1‑148 (MQ=255)

                                                                                                                                                                                                                                    aGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATCCGCCGGTCGCCCCTGTTTTATGTACCCATGATGCGCGTACTGCAGGAGTTGGGGGCGTAATCGCTTTTACGGGGCGTGCTGCCAGTGCAATTTTATTGTCa  <  2:342077/275‑1 (MQ=255)

                                                                                                                                                                                                                                    aGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTATAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTAGCTGGATTGGGATAGTTTTTGTTTGCCAGCTTCACGATAC‑CCAGCTCTTTTTATGGAATAAACGCGACTAAGCTACCAAATCCGCCGGTc                                                                                               >  4:429324/1‑182 (MQ=255)

                                                                                                                                                                                                                                    aGCAAAGAAAAGCGGAGAAAATGTCCGAAAAGTGCGACCCGATGGAATTTTACTGTAGAGGGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATCCGCCGGTCGCCCCTGTTTTATGTACCCATGATGCGCGTACTGCAGGAGTTGGGGGCGTAATCGCTTTTACGGGGCGTGCTGCCAGTGCAATTTTATTGTCa  <  4:167735/275‑1 (MQ=255)

                                                                                                                                                                                                                                      cAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTCTTTTTCTGGAATaaa                                                                                                                           >  2:789790/1‑152 (MQ=255)

                                                                                                                                                                                                                                        aGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCa                                                                                                                                                        <  4:353824/122‑1 (MQ=255)

                                                                                                                                                                                                                                        aGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTCTTTTTCTGGAATAAAcgcg                                                                                                                       <  2:426445/154‑1 (MQ=255)

                                                                                                                                                                                                                                          gAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTCTTTTTCTGGAATAAAcgcg                                                                                                                       <  4:789788/152‑1 (MQ=255)

                                                                                                                                                                                                                                                  gagaAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATCCGCCGGTCGccc                                                                                           <  2:779748/172‑1 (MQ=255)

                                                                                                                                                                                                                                                    gaAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTCtttt                                                                                                                                      >  6:874103/1‑127 (MQ=255)

                                                                                                                                                                                                                                                    gaAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTCtttt                                                                                                                                      <  5:874103/127‑1 (MQ=255)

                                                                                                                                                                                                                                                     atcgggTCCGAATTTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTCTTTTTCTGGAATAAACGCGACATAGc                                                                                                                >  5:133070/5‑148 (MQ=255)

                                                                                                                                                                                                                                                     aaaaGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGcc                                                                                                                                                               >  4:763270/1‑102 (MQ=255)

                                                                                                                                                                                                                                                     aaaaGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTCTTTTTCTGGAATAAACGCGACATAGCt                                                                                                               <  7:580640/149‑1 (MQ=255)

                                                                                                                                                                                                                                                     aaaaGGACCGAAAATTCGGACCCGCTGGAATTCTACTGTAGAGTGTTAAGAATTTGCCAGGCGGCGCCGACTCTCGCTGGATTGGGATAGTTTTTGCTTGCCAGCATCACGACAC‑CCAGCTCTTTTTTTGGAATAAACGCGACATAGCAACCAAACCCGCCggac                                                                                               >  4:132339/1‑163 (MQ=255)

                                                                                                                                                                                                                                                         ggTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTCTTTTTCTGGAATAAACGCGACATAGCTAcc                                                                                                            >  5:160415/1‑148 (MQ=255)

TGTTGCGGTGACTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAC‑CCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATCCGCCGGTCGCCCCTGTTTTATGTACCCATGATGCGCGTACTGCAGGAGTTGGGGGCGTAATCGCTTTTACGGGGCGTGCTGCCAGTGCAATTTTATTGTCA  >  NC_000913/4377534‑4378036

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