BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,563,539	T→C	6.6%	intergenic (+21/+248)	glpD → / ← yzgL	sn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding/pseudogene, periplasmic solute binding protein homology

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,563,539	0	T	C	6.6%	517.3 / 17.6	152	intergenic (+21/+248)	glpD/yzgL	sn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding/pseudogene, periplasmic solute binding protein homology
Reads supporting (aligned to +/- strand): ref base T (63/79); new base C (5/5); total (68/84)
Fisher's exact test for biased strand distribution p-value = 7.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.04e-01

Only 100 of 159 total aligned reads displayed.

TCGCACTTACGGCAGCAACAGCGAGCTGCTGCTCGGCAATGCGGGAACGGTAAGCGATCTCGGGGAAGATTTCGGTCATGAGTTCTACGAAGCGGAGCTGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAG  >  NC_000913/3563266‑3563807

tCGCACTTACGGCAGCAACAGCGAGCTGCTGCTCGGCAATGCGGGAACGGTAAGCGATCTCGGGGAAGATTTCGGTCATGAGTTCTACGAAGCGGAGCTGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGt                                                                                                                                                                                                                                                                              >  2:342740/1‑274 (MQ=255)

gcGCACTTACGGCAGCAACGGCGAGCTGCTGGTCGGCAATGCGGGAACGGTAAGCGATCTCGGGGAAGATTTCGGTCATGAGTTCTACGAAGCGGAGCTGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGt                                                                                                                                                                                                                                                                              <  4:285620/273‑1 (MQ=255)

  gCACTTACGGCAGCAACAGCGAGCTGCTGCTCGGCAATGCGGGAACGGTAAGCGATCTCGGGGAAGATTTCGGTCATGAGTTCTACGAAGCGGAGCTGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCCCACAACAAGGCATGTGGCTAAATGCGGATCAACAAGCCCGTGTGAGTCAGTGGCTGGTGGAGTATACGCCGCATAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCa                                                                                                                                                                                                                                                                            >  2:526225/1‑274 (MQ=255)

  gCACTTACGGCAGCAACAGCGAGCTGCTGCTCGGCAATGCGGGAACGGTAAGCGATCTCGGGGAAGATTTCGGTCATGAGTTCTACGAAGCGGAGCTGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGGAAATTAAGGTAAGGGGATAAGGTAAg                                                                                                                                                                                                                                                                           >  4:467940/1‑275 (MQ=255)

   cACTTACGGCAGCAACAGCGAGCTGCTGCTCGGCAATGCGGGAACGGTAAGCGATCTCGGGGAAGATTTCGGTCATGAGTTCTACGAAGCGGAGCTGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTc                                                                                                                                                                                                                                                                             >  2:517021/1‑272 (MQ=255)

        aCGGCAGCAACAGCGAGCTGCTGCTCGGCAATGCGGGAACGGTAAGCGATCTCGGGGAAGATTTCGGTCATGAGTTCTACGAAGCGGAGCTGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACACGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGCTATTGCTGGCGTCGTAAATTAACGTCACGCCATCAGGTCAGAtctc                                                                                                                                                                                                                                                                      >  2:10842/1‑274 (MQ=255)

                             tgctCGGCAATGCGGGAACGGTAAGCGATCTCGGGGAAGATTTCGGTCATGAGTTCTACGAAGCGGAGCTGAAATATCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCTGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATTAGGTCAGATTTCAATCTGGCCTGAGACtgatga                                                                                                                                                                                                                                                 >  2:100805/1‑274 (MQ=255)

                             gtctCGGCAAGGCGGGAACGGTAAGCGATTTCGGGGAAGATTTCGGGCAAGAGTTCTACGAAGCGGAGCTGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGAc                                                                                                                                                                                                                                                <  4:322350/273‑1 (MQ=255)

                                       tGCGGGAACGGTAAGGGATTTCGGGGAAGATTTCGGTCATGAGTTTTACGAAGCGGAGCTGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACaaa                                                                                                                                                                                                                                       <  2:80888/274‑1 (MQ=255)

                                        gCGGGAACGGTAAGCGATCTCGGGGAAGATTTCGGTCATGAGTTCTACGAAGCGGAGCTGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACaaaa                                                                                                                                                                                                                                      <  2:223577/274‑1 (MQ=255)

                                                   aaGCGATCTCGGAGAAGATTTCGGTCATGAGTTCTACGAAGCGGAGCTGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATgcg                                                                                                                                                                                                                          <  2:265906/275‑1 (MQ=255)

                                                    gaggagacaggggAAGATTTCGGTTAGGATGTCTACGAAGCTGAGCTGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCTACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATgcg                                                                                                                                                                                                                          <  1:55955/265‑1 (MQ=255)

                                                    aGCGATCTCGGGGAAGATTTCGGTCATGAGTTCTACGAAGCGGAGCTGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTAAGATTTCAATCTGGCCTGAGACGGATGACAAACACAAAACTGCCGGATgcc                                                                                                                                                                                                                          >  4:191112/1‑273 (MQ=255)

                                                         tctcGGGGAAGAATTCGGGCAGGAGTTCTCCGAAGCGGAGCTGACATAACTGGTGGATCACGAATGGGTCCGCCGGGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACg                                                                                                                                                                                                                     <  4:367038/274‑1 (MQ=255)

                                                                      ttCGGTCATGAGTGCTACGCAGCGGAGCGGAAATACCTGGTTGATCCCGAATGGGTCCGCCGCTCCGACGACGGCCTGTGGCTTCGCACGAAACAAGGCATGTGGCTAAATTCGGATCAACAATCTCGGGTGAGTCAGTGGCTTGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGAGTTCAATCTGGCCTGAGACTGGAGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTAc                                                                                                                                                                                                       <  4:93559/275‑1 (MQ=255)

                                                                       tCGGTCATTAGTTCTACGAAGCGGAGCTGAAATACCTGGTGGATAACGCATGGTTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTAc                                                                                                                                                                                                       <  1:132446/274‑1 (MQ=255)

                                                                       tCGGTCATGAGTTCTACGAAGCGGAGCTGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTAc                                                                                                                                                                                                       <  2:382308/274‑1 (MQ=255)

                                                                       tCGGTACTGAGTTCGCCAAAGCGGAGCTGAAATACCTGGGGGATAACGAATGGGTCCGCCGAGCCGACGACGCCCTGGGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCGCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTAc                                                                                                                                                                                                       <  4:517020/274‑1 (MQ=255)

                                                                           tcgtggggtacaCGGAGCGGAGGTGAAATAGCGGGTGGGTAACGAAGGGGACCGCCGCGCCGGCGACGCCCTGTGGCGTCGCAACAAACAAGGCATGTGGCTAAATGCGGACAAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTAGACGCAGCAGAGGTTATCGCTGGCGTCGTAAAGTAACGTAAGGTGATCAGGTCAGATTGCAAGCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGt                                                                                                                                                                                                  <  4:129371/264‑1 (MQ=255)

                                                                               gAGTTCTACGAAGCGGAGCTGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTTAGTGGCTGGTGGAGTATACGTATCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTTATTAGGTCAGATTTCACTCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTAt                                                                                                                                                                                              >  1:11131/1‑275 (MQ=255)

                                                                                  tgctccGAAGAGGAGCGGAAATACCTGGTGGATCACGAATGGGTCGGCCGCGCCGACGACGCCCGGTGGCGGCGCAGAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGGAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGc                                                                                                                                                                                            <  4:403709/269‑1 (MQ=255)

                                                                                   tCTACGAAGCGGGGCTGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCTGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCa                                                                                                                                                                                           <  2:2415/274‑1 (MQ=255)

                                                                                          aGCGGAGCTGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATAtt                                                                                                                                                                                    <  2:175490/274‑1 (MQ=255)

                                                                                                  tGAAATACCTGGTGGATAACGAATTGGTCGGCCTCGGCGACGACGCCCTGTGGTGTCGCACAAAGCCAGGCATGTGGCTAAATTCGGATCAACGATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTACGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCa                                                                                                                                                                           <  1:160278/275‑1 (MQ=39)

                                                                                                    aaaTACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATg                                                                                                                                                                         <  2:191113/275‑1 (MQ=255)

                                                                                                                      agaggggTCCGCCGGGCCCTAGACGCCCTGGGGGGTCGCGAAAAACAGGGCATGTGGGTAAATGCGGATTAACAATTNGGGGTGAGTCAGTGGCTGGTGGAGTATACGCAGCGGGGGTTATCGCTGGAGTCGGAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATaa                                                                                                                                                       <  3:124568/271‑1 (MQ=37)

                                                                                                                          gggggtccggcgcagACGGCGCCAGGGNNCGTGGCAAAAAGCAAGGCATGTGGCTAAAGNGGGATCAACAATTTGGTGTGGGGCGGGGGCTGGTGGAGTGTACGCAGCAGAGGTTATGGCTGGCGTCGTAAATTNACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGa                                                                                                                                                     <  3:67932/264‑1 (MQ=37)

                                                                                                                                                ccTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCt                                                                                                                                                                                                                                                            >  7:1110824/1‑148 (MQ=255)

                                                                                                                                                 cTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTgg                                                                                                                                                                                                                                                               >  7:240535/1‑144 (MQ=255)

                                                                                                                                                 cTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTgg                                                                                                                                                                                                                                                               <  8:240535/144‑1 (MQ=255)

                                                                                                                                                          cGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACtgatg                                                                                                                                                                                                                                                  <  8:1110824/148‑1 (MQ=255)

                                                                                                                                                              caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAg                                                                                                                                                                                                                                                                           <  7:401224/119‑1 (MQ=255)

                                                                                                                                                              caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAg                                                                                                                                                                                                                                                                           >  8:401224/1‑119 (MQ=255)

                                                                                                                                                              caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGa                                                                                                                                                                                                                                                                          <  6:807694/120‑1 (MQ=255)

                                                                                                                                                              caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGAt                                                                                                                                                                                                                                                                         >  6:137934/1‑121 (MQ=255)

                                                                                                                                                              caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATc                                                                                                                                                                                                                                                                  >  7:128038/1‑128 (MQ=255)

                                                                                                                                                              caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATc                                                                                                                                                                                                                                                                  <  8:128038/128‑1 (MQ=255)

                                                                                                                                                              caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACtg                                                                                                                                                                                                                                                     >  8:906686/1‑141 (MQ=255)

                                                                                                                                                              caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACtg                                                                                                                                                                                                                                                     <  7:906686/141‑1 (MQ=255)

                                                                                                                                                              caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGCCAGATTGAAATCTg                                                                                                                                                                                                                                                                >  7:176137/1‑130 (MQ=255)

                                                                                                                                                                aaaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGa                                                                                                                                                                                                                                                                          >  5:117463/1‑118 (MQ=255)

                                                                                                                                                                aaaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGa                                                                                                                                                                                                                                                                          <  6:117463/118‑1 (MQ=255)

                                                                                                                                                                  aCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGt                                                                                                                                                                                                                                                                              <  3:52115/112‑1 (MQ=255)

                                                                                                                                                                  aCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACtga                                                                                                                                                                                                                                                    <  6:15833/138‑1 (MQ=255)

                                                                                                                                                                  aCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACtga                                                                                                                                                                                                                                                    >  5:15833/1‑138 (MQ=255)

                                                                                                                                                                                   aTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATgcgc                                                                                                                                                                                                                         >  6:536833/1‑148 (MQ=255)

                                                                                                                                                                                             acaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAAt                                                                                                                                                                                                                                                                   <  8:528077/96‑1 (MQ=255)

                                                                                                                                                                                             acaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAAt                                                                                                                                                                                                                                                                   >  7:528077/1‑96 (MQ=255)

                                                                                                                                                                                             acaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTgg                                                                                                                                                                                                                                                               >  6:351168/1‑100 (MQ=255)

                                                                                                                                                                                             acaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTgg                                                                                                                                                                                                                                                               <  5:351168/100‑1 (MQ=255)

                                                                                                                                                                                             acaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGa                                                                                                                                                                                                                              >  7:249933/1‑133 (MQ=255)

                                                                                                                                                                                             acaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGa                                                                                                                                                                                                                              <  8:249933/133‑1 (MQ=255)

                                                                                                                                                                                                aTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTc                                                                                                                                                                                                                                                                      >  5:34165/1‑90 (MQ=255)

                                                                                                                                                                                                aTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTc                                                                                                                                                                                                                                                                      <  6:34165/90‑1 (MQ=255)

                                                                                                                                                                                                      tgtgAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGNTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACNTTCACGTCGAATCCGGCATGAACAAAGAGCACTTTGTTAAAAATCNANACTACTTTTTAATTTAAGGAAATaacac                                                                       >  3:33493/1‑271 (MQ=255)

                                                                                                                                                                                                      tgtgAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCg                                                                        >  2:43472/1‑274 (MQ=255)

                                                                                                                                                                                                         gAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTgg                                                                                                                                                                                                   <  5:33956/148‑1 (MQ=255)

                                                                                                                                                                                                          aGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGt                                                                                                                                                                                                  <  5:536833/148‑1 (MQ=255)

                                                                                                                                                                                                           gTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCt                                                                                                                                                                                                         <  5:143794/140‑1 (MQ=255)

                                                                                                                                                                                                            tCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAac                                                                                                                                                                                                                                            <  6:560566/104‑1 (MQ=255)

                                                                                                                                                                                                            tCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCCCATCCGGCCTTCACAAAGCTCACTTTGTCAAAAATCTAACATACTTTTTAATTCCGGGAAATACCGCAAAGc                                                                  >  4:366545/1‑274 (MQ=255)

                                                                                                                                                                                                             cAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAgg                                                                                                                                                                                                            >  7:190559/1‑135 (MQ=255)

                                                                                                                                                                                                             cAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAgg                                                                                                                                                                                                            <  8:190559/135‑1 (MQ=255)

                                                                                                                                                                                                               gTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGc                                                                                                                                                                                                           >  4:467860/1‑134 (MQ=255)

                                                                                                                                                                                                                  gCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTa                                                                                                                                                                                               >  5:611114/1‑143 (MQ=255)

                                                                                                                                                                                                                  gCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTa                                                                                                                                                                                               <  6:611114/143‑1 (MQ=255)

                                                                                                                                                                                                                   cTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTAc                                                                                                                                                                                                       <  2:467861/134‑1 (MQ=255)

                                                                                                                                                                                                                       tggAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGcc                                                                                                                                                                                                                                 <  7:455430/104‑1 (MQ=255)

                                                                                                                                                                                                                           gTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATAcc                                                   <  1:58355/274‑1 (MQ=255)

                                                                                                                                                                                                                            tataCGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCaa                                                                                                                                                                                                                                                                    >  5:929044/1‑64 (MQ=255)

                                                                                                                                                                                                                                    gcagAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGCCAGATTGAAATCTga                                                                                                                                                                                                                                                               <  7:897701/61‑2 (MQ=37)

                                                                                                                                                                                                                                    gcagAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGCCAGATTGAAATCTga                                                                                                                                                                                                                                                               >  8:897701/1‑60 (MQ=37)

                                                                                                                                                                                                                                       gagGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAAtt                                                                                                                                                                               <  5:609220/138‑1 (MQ=255)

                                                                                                                                                                                                                                        agGTTATCGCTGGCGTCGGAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCtt                                                                                                                                                                    <  8:142344/148‑1 (MQ=255)

                                                                                                                                                                                                                                            tATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTgag                                                                                                                                                                                                                                                         >  8:404051/1‑59 (MQ=255)

                                                                                                                                                                                                                                            tATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTa                                                                                                                                                                                               >  6:292666/1‑117 (MQ=255)

                                                                                                                                                                                                                                            tATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAg                                 <  4:43471/275‑1 (MQ=255)

                                                                                                                                                                                                                                             aTCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAACCACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAg                                 >  2:105122/1‑274 (MQ=255)

                                                                                                                                                                                                                                             aTCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAAttt                                                                                                                                                                              >  2:495273/1‑133 (MQ=255)

                                                                                                                                                                                                                                             aTCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTc                                                                                                                                           >  2:470055/1‑168 (MQ=255)

                                                                                                                                                                                                                                              tCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGc                                                                                                                                                              >  5:517258/1‑148 (MQ=255)

                                                                                                                                                                                                                                                 cTGGCGTCGGAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCAt                                                                                                                                                                          <  4:495272/133‑1 (MQ=255)

                                                                                                                                                                                                                                                    gCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAAtt                                                                                                                                                                               >  1:18295/1‑125 (MQ=255)

                                                                                                                                                                                                                                                    gCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCATTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGcc                                                                 >  2:332909/1‑235 (MQ=255)

                                                                                                                                                                                                                                                              ttAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCa                                                                                                                                                                                                              >  6:1075819/1‑84 (MQ=255)

                                                                                                                                                                                                                                                              ttAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCa                                                                                                                                                                                                              <  5:1075819/84‑1 (MQ=255)

                                                                                                                                                                                                                                                                      aGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACaaaa                                                                                                                                                                                                                                      >  4:392424/1‑52 (MQ=255)

                                                                                                                                                                                                                                                                        gTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATAtt                                                                                                                                                                                    >  6:51148/1‑100 (MQ=255)

                                                                                                                                                                                                                                                                         tGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTa                                                                                                                                                                                               >  7:684809/1‑88 (MQ=255)

                                                                                                                                                                                                                                                                         tGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTa                                                                                                                                                                                               <  8:684809/88‑1 (MQ=255)

                                                                                                                                                                                                                                                                         ggtacAGGTCCGATTTCAACCTTGCCTGGGACGGATGACAAACACAAAATGTCCTGATGGCTTACTTTTAGCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAATGCGCACTTTGTCAAAAATCTACCCTACTTTTTAATTTAGGGAATTACCGCAAAGCCAACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGaaaaa     <  4:234362/270‑1 (MQ=255)

                                                                                                                                                                                                                                                                            tCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAgg                                                                                                                                                                                                            >  5:76229/1‑72 (MQ=255)

                                                                                                                                                                                                                                                                            tCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAgg                                                                                                                                                                                                            <  6:76229/72‑1 (MQ=255)

                                                                                                                                                                                                                                                                               ggTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCa                                                                                                                                                                                           >  8:977423/1‑86 (MQ=255)

                                                                                                                                                                                                                                                                               ggTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCa                                                                                                                                                                                           <  7:977423/86‑1 (MQ=255)

                                                                                                                                                                                                                                                                               ggTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACg                                                             <  3:135080/212‑1 (MQ=255)

                                                                                                                                                                                                                                                                                gTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCa                                                                                                                                                                           <  8:634165/101‑1 (MQ=255)

                                                                                                                                                                                                                                                                                gTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCa                                                                                                                                                                           >  7:634165/1‑101 (MQ=255)

                                                                                                                                                                                                                                                                                 tCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAg                                                                                                                                                                                                             <  7:108252/66‑1 (MQ=255)

                                                                                                                                                                                                                                                                                 tCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAg                                                                                                                                                                                                             >  8:108252/1‑66 (MQ=255)

TCGCACTTACGGCAGCAACAGCGAGCTGCTGCTCGGCAATGCGGGAACGGTAAGCGATCTCGGGGAAGATTTCGGTCATGAGTTCTACGAAGCGGAGCTGAAATACCTGGTGGATCACGAATGGGTCCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAG  >  NC_000913/3563266‑3563807

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