BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,946,246	G→A	noncoding (2543/2904 nt)	rrlC →	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,946,246	0	G	A	96.0%	134.5 / ‑2.5	77	noncoding (2543/2904 nt)	rrlC	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base G (2/1); new base A (33/40); total (35/42)
Fisher's exact test for biased strand distribution p-value = 5.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCTGGTCGGACATCAGGAGGTTAGTGCAATGGCATAAGCCAGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGT‑TCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑GGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGC  >  NC_000913/3946001‑3946474

aTCCTGGTCGGACATCAGGAGGTTAGTGCAATGGCATAAGCCAGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGT‑TCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATAACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTAt                                                                                                                                                                                                                                  >  2:14970/1‑251 (MQ=14)

  ccTGGTCGGACATCAGGAGGTTAGTGCAATGGCATAAGCCAGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGT‑TCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATg                                                                                                                                                                                                                                 <  1:259589/250‑1 (MQ=14)

     ggTCGGACATCAGGAGGTTAGTGCAATGGCATAAGCCAGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGT‑TCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGc                                                                                                                                                                                                                               <  1:25409/249‑1 (MQ=14)

        cgtccATCAGGAGGTTAGGGCAAGGGTATAAGCGAGCTTGACGGCGAGCGGGACGGCGCGAGCAGGGGCGAAAGCAGGTCATAGTGATCCGG‑GGTGTCTGAATGGAAGGGCCATCGCTCAACGGATAAAATGTCCTTCGTGGATAATTCGCGGATACCGCCCCAGAGTTCATATCGTCGGCGGTGTGTGGTACCGCGATGGCGAGTCCGCTCATCCGGGAGCTGAAGTAGGTCCCA‑A‑TGGttggtctggt                                                                                                                                                                                                                           <  2:75858/246‑10 (MQ=2)

        cGGACATTAGGAGGTTAGTGCAATGGGATAAGCCAGCTTGACTGCGAGCGCGACGGCGCGAGCTGGTGCGAAAGCAGGTCAGAGTGATCCGGTGGT‑TCTGAAGGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGTATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGtt                                                                                                                                                                                                                           <  1:462696/250‑1 (MQ=11)

           aCATCAGGAGGTTAGTGCAATGGCATAAGCCAGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGT‑TCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGcc                                                                                                                                                                                                                       <  1:114523/251‑1 (MQ=14)

             aTCAGGAGGTTAGTGCAATGACATAAGCCAGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGT‑TCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCAt                                                                                                                                                                                                                     <  2:238397/251‑1 (MQ=14)

                gggaggTTAGTGCAATGGCATAAGCCAGCTTGACTGCGAGCGTGACGGCGCGAGCCGGTGCGAAAGCAGGTCATAGTGATCCGGTGGT‑TCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTAATATCGCCGGCGGTGTTTGGCACCTCGATGTCAGCTCATCACATCCTGGAGATGAAGTAGGTCCCACA‑GGG‑ATGGCTGttacccctttt                                                                                                                                                                                                                  >  2:110000/2‑242 (MQ=11)

                           gCAATGGCATAAGCCAGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGT‑TCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTAcgc                                                                                                                                                                                                       <  2:225921/251‑1 (MQ=14)

                                ggCATAAGCCAGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGT‑TCTGAATGGAAGGGCCAGCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGc                                                                                                                                                                                                   <  2:5660/250‑1 (MQ=14)

                                    tAAGCCCGCTTGAATGCGAGCGTGCCGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGT‑TCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCAACTCGCAGTCGGCTCATCACATCATGGAGCTGAAGGCGGGTCCT‑A‑GGGGATGGCTGTTCGCCATTTAAAGTGGGAAGAGAGCTgggg                                                                                                                                                                                              >  2:162045/1‑250 (MQ=11)

                                        ccAGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGT‑TCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGttt                                                                                                                                                                                            >  1:101256/1‑249 (MQ=14)

                                            cTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGT‑TCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcg                                                                                                                                                                                      >  1:231025/1‑251 (MQ=14)

                                               gACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGT‑TCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgtc                                                                                                                                                                                    <  2:101256/250‑1 (MQ=14)

                                                aCTGCGAGCGTGACGGCGCGCGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGT‑TCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgtcg‑t                                                                                                                                                                                 >  1:148790/1‑251 (MQ=14)

                                                aCTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGT‑TCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGAGAGCTGGGTTTAGAAcgtcgg                                                                                                                                                                                  >  1:423176/1‑250 (MQ=14)

                                                   gCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGT‑TCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑Tgag                                                                                                                                                                              <  1:341234/251‑1 (MQ=14)

                                                      aGCGTGACGACGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGT‑TCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCCCCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGGCG‑TGGGAc                                                                                                                                                                            >  2:428602/1‑250 (MQ=9)

                                                         gTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGT‑TCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGGAAGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGtt                                                                                                                                                                        >  2:349884/1‑251 (MQ=14)

                                                            aCGGCGCGAGCAGGAGCGAAAGCAGGTCATAGTGATCCGGTGGT‑TCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCgg                                                                                                                                                                     >  1:409316/1‑251 (MQ=14)

                                                                     ggaGGTGCGAAAGCAGGTCATAGTGAGCCGGTGGT‑GCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGGACGCCGGGGATAAAAGGCTGATACCGCCCAAGAGTTCATAGCGACGGCGGTGTTTGGCACCGCGATGTCGGCTCATCACATCCGGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCGGTTCGCCATGTAAAGTGGTACGCGAGCTGGGTGTAGAACGGCG‑TGAGACAGTTCGGGCCCTagct                                                                                                                                                            <  2:166260/249‑4 (MQ=11)

                                                                     gCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGT‑TCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCt                                                                                                                                                            <  1:14970/251‑1 (MQ=14)

                                                                        ggTGCGAAATCATGTCATAGTGATCCGGTGGT‑TCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGTTACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCTCATCCTGGAGCTGAAGTAGGCCCCA‑A‑AGGTATGTCTGTTCGCCAtt                                                                                                                                                                                                                    <  2:393809/193‑1 (MQ=11)

                                                                        ggTGCGAAAGCAGGTCATAGTGATCCGGTGGT‑TCTGAATGGAAGGACCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCaat                                                                                                                                                                                                                    >  1:393809/1‑191 (MQ=12)

                                                                               aaGCAGGTCATAGTGATCCGGTGGT‑TCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGgcg                                                                                                                                                  >  1:145225/1‑251 (MQ=14)

                                                                                aGCATGTCATAGGGATCCGGCGGT‑TCTGAATGGAAGGGACATCGCTCAACGGATAAAAGGTACTACGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGGTTGGCACCTCGATGCCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCG‑A‑AGGAATGGAAGTTCGCCATTTTAAGTGGGACGCAAGCTGGGTTTTGACAGTCGGTG‑GACCGTTCGGACCCTATCTGCCGTTGgcgc                                                                                                                                                 >  2:121183/1‑251 (MQ=11)

                                                                                 gCAGGTCATAGTGATCCGGTGGT‑TCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCt                                                                                                                                                <  2:156039/251‑1 (MQ=14)

                                                                                 gCAGGTCATAGTGATCCGGTGGT‑TCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGGACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCt                                                                                                                                                >  2:133915/1‑251 (MQ=14)

                                                                                  cAGGTCAGAGTGATCCGGTGGT‑TCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTg                                                                                                                                               >  1:254273/1‑251 (MQ=14)

                                                                                    ggTCATAGTGATCCGGCGGT‑TCTGAATGAACGGGACATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCCAGAGTTCATATCGACGGAGGTGTTTGGCACCTCGATGTCGGCTCATCAAATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTAAGGCTGTTTGCCATTTCAAGGGGTACGCGAGCTGGGTTTAGCACCTTG‑TGAGCCAGTTTGGTCCCTATCTGCCGTGGGCGCTggg                                                                                                                                             >  2:296303/1‑250 (MQ=11)

                                                                                     gTCATAGTGATCCGGTGGT‑TCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGgag                                                                                                                                            <  1:5823/251‑1 (MQ=14)

                                                                                               tCCGGTGGT‑TCTGAATGGAAAGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCCTATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGAAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGGGGTACGCGAGCTGGGTTTAGAACGTAG‑TGAGGCAGTTCGGTCCCTATCTGCCGTGGGCCCTGGAGAACTGCgggg                                                                                                                                  >  2:276370/1‑251 (MQ=11)

                                                                                                   gtggt‑TCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGgct                                                                                                                              <  1:76794/251‑1 (MQ=14)

                                                                                                    tggt‑TCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCCTCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGAATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTACCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGgctg                                                                                                                             >  1:425291/1‑251 (MQ=9)

                                                                                                          cTGAATGGAAGGGCCATCGCTCAACGGATAAACGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAACGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCt                                                                                                                        >  2:206976/1‑251 (MQ=14)

                                                                                                            gAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAg                                                                                                                      <  2:353730/251‑1 (MQ=14)

                                                                                                                    gggCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACgag                                                                                                                >  1:248363/1‑249 (MQ=14)

                                                                                                                       ccATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGAc                                                                                                           >  1:463168/1‑251 (MQ=14)

                                                                                                                       ccATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGAc                                                                                                           <  1:350200/251‑1 (MQ=14)

                                                                                                                           cGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTTGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTTGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGGGGGCGCTGGAGAACTGAGGTGGTCTGCTCCTAGTACTAGATGACCGGa                                                                                                       <  2:231025/251‑1 (MQ=11)

                                                                                                                           cGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGa                                                                                                       <  1:228901/251‑1 (MQ=14)

                                                                                                                               cAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTgg                                                                                                   <  1:25923/251‑1 (MQ=14)

                                                                                                                               cAACGGATAAAAGGTACTCAGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACGTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTg                                                                                                    <  1:86342/250‑1 (MQ=14)

                                                                                                                                 aCGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGAc                                                                                                 <  1:124818/251‑1 (MQ=14)

                                                                                                                                    gATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTGTGCCATTTAAAGGGGTACGCGAACTGGGTTTATAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGGGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGGGGACGc                                                                                               >  2:248001/1‑250 (MQ=11)

                                                                                                                                      tAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCAACTCGCTGTCGGCTCATAACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATc                                                                                            >  2:451376/1‑251 (MQ=9)

                                                                                                                                       aaaaGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATc                                                                                            <  1:133915/250‑1 (MQ=14)

                                                                                                                                         aaGGTACTTCGGGACTAACAGGCGGATACCGCCCAAGAGTCCATATCGACGGCGGTGTTTGGCTTCGCGATGTCGGCTAATCACATCCTGGAGTTGAAGTAGGTCCCA‑AAGGGT‑TGGCTGTTCGCGATTTAAAGTGTTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCCATCTGCCGTTGGCGCGTGAGAACTGAGTGGGGCTGCTCCTAGTACGAGTGGACCTCAGTGGACGCATctct                                                                                         <  2:19213/251‑4 (MQ=11)

                                                                                                                                              aCTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTaaa                                                                                                                                                                                                                <  1:252209/128‑1 (MQ=18)

                                                                                                                                              aCTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATACTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCCTTTaaa                                                                                                                                                                                                                >  2:252209/1‑128 (MQ=255)

                                                                                                                                                tCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCGCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGGTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCAGGGGGCGCTGGAGAAATGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTc                                                                                  >  1:450287/1‑251 (MQ=255)

                                                                                                                                                    gggATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCggg                                                                               <  1:429750/250‑1 (MQ=14)

                                                                                                                                                    gggATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGACGGCTCATCACAACCTGGAGCTGAAGTAGGTCCCC‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGACGGGGGCGCTGGAGAACTGGGGGGGGCTGCCCCTAGTACGGGAGGACCGGAGTGGACGCATCACTGGGGTTCgggg                                                                              >  2:429750/1‑250 (MQ=11)

                                                                                                                                                         aacaggaagATACCGTCCAAGAGTACCTATCGACAGAGGTGTATGGACCCTAAATGTCGGCTTCTCACATCATGGAGCCGAAGTAGGTCCCC‑A‑AGGCATCGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGCGACCGTTCGGTCCCTATCTGCCGTGGGCGCTGGGGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGAAGCATTACTGGTGTTCGGGTTGTc                                                                          >  2:47066/9‑250 (MQ=2)

                                                                                                                                                                tGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAAt                                                                  >  1:332613/1‑251 (MQ=14)

                                                                                                                                                                 gATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATg                                                                 <  1:403654/251‑1 (MQ=14)

                                                                                                                                                                      cGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGGCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACGGGTGTTCGGGTTGTCATGCCAATGGCACt                                                            <  2:218401/251‑1 (MQ=14)

                                                                                                                                                                         ccAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGcc                                                         <  1:428602/251‑1 (MQ=14)

                                                                                                                                                                          cAAGAGTTCATATCGACGGCGGTGTTTGGCCCCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGccc                                                        <  2:253093/251‑1 (MQ=14)

                                                                                                                                                                              agTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTCGTACGAGGGGACCGGAGTGGACGCATCACGGGTGGTCGGGTTGTCATGcc                                                                     >  2:155272/1‑234 (MQ=255)

                                                                                                                                                                              agTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGcc                                                                     <  1:155272/234‑1 (MQ=14)

                                                                                                                                                                               gTTCATATCGACGGCGGTGTTTGGCACCTCGATGGCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTGCGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAg                                                   <  1:396259/251‑1 (MQ=14)

                                                                                                                                                                                  cATATCGACGGGGGTGGTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGGTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTa                                                <  1:251460/251‑1 (MQ=9)

                                                                                                                                                                                     atCGACGGCGGTGTTTGGCTCCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGGAGCTAAAt                                             <  1:349884/251‑1 (MQ=14)

                                                                                                                                                                                     atCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTATCTaaa                                              <  1:163155/250‑1 (MQ=14)

                                                                                                                                                                                          cggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTCGGTCCCA‑A‑AGGTATGGCTGTTCGAAATTTAAAGTGGTACGCCAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGGCCCCATCTGCCGTGGGGGGTGGAGACATGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGGGGACGCCTCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGGAGCTAAATGCggg                                        >  1:263946/1‑250 (MQ=11)

                                                                                                                                                                                                       cACCTCTATGTCGGCTCATCACATCATGGAGCTGCAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGATCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGAGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAAGGCGGACGAGATAAGTGCt                           >  2:211266/1‑251 (MQ=11)

                                                                                                                                                                                                       cACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGCGATAAGTGCt                           >  2:181651/1‑251 (MQ=14)

                                                                                                                                                                                                            cGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCTCTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAg                      <  1:426208/251‑1 (MQ=14)

                                                                                                                                                                                                            cGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGaa                        <  1:144243/249‑1 (MQ=14)

                                                                                                                                                                                                             gATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGgtgt                                                                                    <  2:155746/188‑1 (MQ=17)

                                                                                                                                                                                                             gATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGgtgt                                                                                    >  1:155746/1‑188 (MQ=17)

                                                                                                                                                                                                                gTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGCAAGCATc                  >  1:264158/1‑251 (MQ=14)

                                                                                                                                                                                                                   ggCTCAGCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGGGCTGAAAGCATCTaa               <  2:391994/251‑1 (MQ=14)

                                                                                                                                                                                                                        atcaCATCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACg          <  2:9519/251‑1 (MQ=14)

                                                                                                                                                                                                                           acaTCCTGGAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGaa        >  1:297664/1‑250 (MQ=14)

                                                                                                                                                                                                                                  ggAGCTGAAGTAGGTCCCA‑A‑AGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGTGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGc  <  1:296303/249‑1 (MQ=14)

ATCCTGGTCGGACATCAGGAGGTTAGTGCAATGGCATAAGCCAGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGT‑TCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCA‑A‑GGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCG‑TGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGC  >  NC_000913/3946001‑3946474

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`