BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,334,705	C→T	7.8%	intergenic (+24/‑204)	ispB → / → sfsB	octaprenyl diphosphate synthase/malPQ operon transcriptional activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,334,705	0	C	T	7.8%	625.8 / 29.3	203	intergenic (+24/‑204)	ispB/sfsB	octaprenyl diphosphate synthase/malPQ operon transcriptional activator
Reads supporting (aligned to +/- strand): ref base C (93/94); new base T (8/8); total (101/102)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.19e-02

Only 100 of 208 total aligned reads displayed.

ATGCATCATGGCACACCAGAACAGGCACAGATGATCCGTACCGCCATCGAACAGGGTAACGGTCGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGTTCTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGC  >  NC_000913/3334439‑3334960

acagagtcgggCACAACAGAACCGGCACAGATGATCCGGAGCGGCATCGAACAGGGTAACGGGCTGCATCGTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGAGTGCTTGAATGGACGCGTCAGCGGGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTGACAGGGGCTCCC‑GGACACCCCGTGGCGAGAAGCACTCATCGGCGTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGa                                                                                                                                                                                                                                                         <  4:399295/266‑1 (MQ=255)

      catGGCACACCAGAACAGGCACAGATGATCCGTACCGCCATCGAACAGGGTAACGGTCGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTCCAGGTGCTCCC‑GGACGCCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACCTCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCCCAGC‑GCGCGGGCTGGTTc                                                                                                                                                                                                                                                   >  2:485803/1‑275 (MQ=255)

        tGGCACACCAGAACAGGCACAGATGATCCGTACCGCCATCGAACAGGGTAACGGTCGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATGCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCCGTaa                                                                                                                                                                                                                                             >  1:59531/1‑279 (MQ=255)

                          acggatgatCCGTACCGCCGTCGAACAGGGTACCGGTCGCCAGCTTCTGGAACCGTGTCTGTAAGCAATGGACGCTTGTGGAGCTTTTGAATGGACGCGTCAGCGTGCCGAGGCAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTGGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCGGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAac                                                                                                                                                                                                                                <  4:22063/271‑1 (MQ=255)

                             gatgatCCGTACCGCCATCGAACAGGGTAACGGTCGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCACC‑CCGCGGGATGGTTCCCGTAAGTTCCATAAAAacccct                                                                                                                                                                                                                            >  4:261247/1‑271 (MQ=255)

                              atgatCCGTACCGCCATCGAACAGGGTAACGGTCGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATaaa                                                                                                                                                                                                                                   <  4:485555/267‑1 (MQ=255)

                                                           cGGTTGTCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGGGGATCGCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAggg                                                                                                                                                                                              <  2:542533/275‑1 (MQ=255)

                                                               cGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAg                                                                                                                                                                                                                      <  2:312990/247‑1 (MQ=255)

                                                                gCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                                                                                                                                         <  2:594867/275‑1 (MQ=255)

                                                                gCCATCTTCTGCAACCCGTTCTGGAAGGCATTAACGATTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGACTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                                                                                                                                         <  2:723967/275‑1 (MQ=255)

                                                                     cttttGGAAACGGTTCTGGAATCAATGAACGCATATGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATc                                                                                                                                                                                     <  2:386520/274‑1 (MQ=255)

                                                                                   tCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGaa                                                                                                                                                                      <  2:555961/275‑1 (MQ=255)

                                                                                                          gATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTcac                                                                                                                                                               <  2:243078/259‑1 (MQ=255)

                                                                                                                ttGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTaa                                                                                                                                                                                                                                             >  2:700074/1‑175 (MQ=255)

                                                                                                                    aTGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTc                                                                                                                                                                                                                                         <  4:700074/175‑1 (MQ=255)

                                                                                                                          gcgTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgc‑gcg                                                                                                                                                                                                                                                              >  8:471896/1‑148 (MQ=255)

                                                                                                                            gTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGc‑gcgcg                                                                                                                                                                                                                                                            <  7:496230/148‑1 (MQ=255)

                                                                                                                             tCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCgg                                                                                                                                                                                                                                                           >  8:36718/1‑148 (MQ=255)

                                                                                                                                 cGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATgg                                                                                                                                                                                                                                                      <  6:506702/149‑1 (MQ=255)

                                                                                                                                   tGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCCCATCCCGCGCGCT‑GCGCGGGATgagg                                                                                                                                                                                                                                                    <  5:706950/149‑4 (MQ=255)

                                                                                                                                   tGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCCCATCCCGCGCGCT‑GCGCGGGATgagg                                                                                                                                                                                                                                                    <  6:270828/149‑4 (MQ=255)

                                                                                                                                      cGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGTACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCa                                                                                                                                                                                                                                                 <  5:311467/149‑1 (MQ=255)

                                                                                                                                       gAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAg                                                                                                                                                                                                                                                >  5:901322/1‑149 (MQ=255)

                                                                                                                                        aGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAg                                                                                                                                                                                                                                                >  8:834724/1‑148 (MQ=255)

                                                                                                                                          gaagaaGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGT‑GCGCGGGATGGTTCCAGTa                                                                                                                                                                                                                                              <  7:359335/148‑1 (MQ=255)

                                                                                                                                           aagaagCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTaa                                                                                                                                                                                                                                             <  8:878139/148‑1 (MQ=255)

                                                                                                                                            agaagCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCTCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGt                                                                                                                                                                                                                                           >  6:291550/1‑149 (MQ=255)

                                                                                                                                               agCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCT‑GCGCGGGa                                                                                                                                                                                                                                                         <  8:309939/132‑1 (MQ=255)

                                                                                                                                               agCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCT‑GCGCGGGa                                                                                                                                                                                                                                                         >  7:309939/1‑132 (MQ=255)

                                                                                                                                                gCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCa                                                                                                                                                                                                                                       >  6:806923/1‑149 (MQ=255)

                                                                                                                                                gCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCa                                                                                                                                                                                                                                       >  5:488448/1‑149 (MQ=255)

                                                                                                                                                  aGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATa                                                                                                                                                                                                                                     >  5:133570/1‑149 (MQ=255)

                                                                                                                                                      aatgcccgcGCAGCGTTACAGGTGCTCCCGGGACA‑CCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCGGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATGCTCGACag                                                                                                    <  4:658523/266‑1 (MQ=255)

                                                                                                                                                          ccATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAacac                                                                                                                                                                                                                              >  8:515797/1‑148 (MQ=255)

                                                                                                                                                          ccATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACt                                                                                                                                                                                                                             >  6:722959/1‑149 (MQ=255)

                                                                                                                                                              cGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCGTAAAAACACTTAtt                                                                                                                                                                                                                         <  5:806923/149‑1 (MQ=255)

                                                                                                                                                              cGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTAt                                                                                                                                                                                                                          <  7:118918/148‑1 (MQ=255)

                                                                                                                                                              cGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTg                                                                                                                                                                                                        <  3:40982/166‑1 (MQ=255)

                                                                                                                                                              cGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTg                                                                                                                                                                                                        >  1:40982/1‑166 (MQ=255)

                                                                                                                                                                         aGGTGCTCCC‑GGACACCCCTTGGCGTGAAGCACTCATCGGCCGCGAGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGGAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTAGTCAGctctgcaa                                                                                                                                                                                                              <  6:287090/149‑5 (MQ=255)

                                                                                                                                                                             gCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGTTCTGAATGAAAGCGa                                                                            <  4:560127/275‑1 (MQ=255)

                                                                                                                                                                              cTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCt                                                                                                                                                                                                         <  6:14337/149‑1 (MQ=255)

                                                                                                                                                                               tCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTg                                                                                                                                                                                                        >  5:560281/1‑149 (MQ=255)

                                                                                                                                                                                ccc‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTg                                                                                                                                                                                                        >  7:378836/1‑148 (MQ=255)

                                                                                                                                                                                ccc‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGa                                                                                                                                                                                                       >  6:98040/1‑149 (MQ=255)

                                                                                                                                                                                    ggACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                                                                                                                                         <  2:624211/160‑1 (MQ=255)

                                                                                                                                                                                         ccccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAgg                                                                                                                                                                                               <  6:901322/149‑1 (MQ=255)

                                                                                                                                                                                         ccccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCCCATCCCGCGCGCT‑GCGCGGGATg                                                                                                                                                                                                                                                       >  7:452617/1‑93 (MQ=39)

                                                                                                                                                                                          cccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATa                                                                                                                                                                                                 <  5:640646/146‑1 (MQ=255)

                                                                                                                                                                                          cccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATa                                                                                                                                                                                                 >  6:640646/1‑146 (MQ=255)

                                                                                                                                                                                           ccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATaa                                                                                                                                                                                                              >  7:435576/1‑132 (MQ=255)

                                                                                                                                                                                           ccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAggg                                                                                                                                                                                              <  8:435903/148‑1 (MQ=255)

                                                                                                                                                                                           ccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                                                                                                                                             <  5:793653/149‑1 (MQ=255)

                                                                                                                                                                                                gCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                                                                                                                                         >  7:258727/1‑148 (MQ=255)

                                                                                                                                                                                                   agaAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCAt                                                                                                                                                                                      >  7:607222/1‑148 (MQ=255)

                                                                                                                                                                                                    gaAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCtata                                                                                                                                                                                                               >  2:154797/1‑122 (MQ=255)

                                                                                                                                                                                                    gaAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGATCTGAATGAAAGCGACTAACTTAAGGAGTGAGGaaaa                                                      >  2:458405/1‑275 (MQ=255)

                                                                                                                                                                                                     aaGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCtt                                                                                                                                                                                   >  5:890233/1‑149 (MQ=255)

                                                                                                                                                                                                          cTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtg                                                                                                                                                                                           >  8:273444/1‑136 (MQ=255)

                                                                                                                                                                                                          cTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTtatat                                                                                                                                                                              >  6:564404/1‑149 (MQ=255)

                                                                                                                                                                                                          cTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGATCTTAATGAAAGCGACTAACTTAAGGGTTGCGGCAAAATGGa                                                 >  4:154533/1‑274 (MQ=255)

                                                                                                                                                                                                           tCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATg                                                                                                                                                                                                                                                       >  7:807704/1‑75 (MQ=255)

                                                                                                                                                                                                             aTCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATAtt                                                                                                                                                                           >  6:19775/1‑149 (MQ=255)

                                                                                                                                                                                                             aTCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATAtt                                                                                                                                                                           <  5:512057/149‑1 (MQ=255)

                                                                                                                                                                                                             aTCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATAtt                                                                                                                                                                           >  5:954547/1‑149 (MQ=255)

                                                                                                                                                                                                               cGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGtt                                                                                                                                                                                                                                                    <  7:573665/74‑1 (MQ=255)

                                                                                                                                                                                                                 gCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTg                                                                                                                                                                        >  7:1028327/1‑148 (MQ=255)

                                                                                                                                                                                                                      gcgcACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATAtt                                                                                                                                                                  <  5:246422/149‑1 (MQ=255)

                                                                                                                                                                                                                      gcgcACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATAtt                                                                                                                                                                  >  5:546069/1‑149 (MQ=255)

                                                                                                                                                                                                                        tcacaTCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCTGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTca                                                                                                                                                                <  5:1051342/148‑1 (MQ=255)

                                                                                                                                                                                                                        gcACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTc                                                                                                                                                                 <  7:870610/148‑1 (MQ=255)

                                                                                                                                                                                                                          acaTCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATa                                                                                                                                                                                                 <  6:430595/114‑1 (MQ=255)

                                                                                                                                                                                                                          acaTCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATa                                                                                                                                                                                                 >  5:430595/1‑114 (MQ=255)

                                                                                                                                                                                                                              cGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACActc                                                                                                                                                           <  8:10734/148‑1 (MQ=255)

                                                                                                                                                                                                                                cTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCttt                                                                                                                                                        >  5:362188/1‑149 (MQ=255)

                                                                                                                                                                                                                                cTGTTCAACGCGATCGTTAATCACCTCACCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACATGAAAAATAGGGATGTGCATCTTTATATAGTCTGAATATTCACACTCTTTCCAGGAACTTTTTAGAGAAATAGGCCATCAGGAGTATAGTGATGCGCGACAGAAGAAGGGTTATGAATGACAGCGACAAACTTAAGGAGTGAGGAAAAAGGGAAAGtaatttaattacatggactc                          >  4:605522/1‑265 (MQ=255)

                                                                                                                                                                                                                                  gTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTAc                                                                                                                                                      <  6:853978/149‑1 (MQ=255)

                                                                                                                                                                                                                                   ttCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACa                                                                                                                                                     >  6:1009612/1‑149 (MQ=255)

                                                                                                                                                                                                                                            aTCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGaa                                                                                                                                                 >  5:170250/1‑144 (MQ=255)

                                                                                                                                                                                                                                             tCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                                                                                                                                         >  4:379560/1‑103 (MQ=255)

                                                                                                                                                                                                                                             tCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGTTCTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGAcatcat             >  2:320874/1‑275 (MQ=255)

                                                                                                                                                                                                                                                      ccctcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATa                                                                                                                                  <  6:939939/149‑1 (MQ=255)

                                                                                                                                                                                                                                                       cctcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATa                                                                                                                                  >  8:624646/1‑148 (MQ=255)

                                                                                                                                                                                                                                                        ctcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAgg                                                                                                                                <  6:546069/149‑1 (MQ=255)

                                                                                                                                                                                                                                                         tcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                                                                                                                                         >  1:119154/1‑91 (MQ=255)

                                                                                                                                                                                                                                                         tcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                                                                                                                                         <  3:119154/91‑1 (MQ=255)

                                                                                                                                                                                                                                                         tcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAg                                                                                                                                 <  1:24522/147‑1 (MQ=255)

                                                                                                                                                                                                                                                         tcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAg                                                                                                                                 >  3:24522/1‑147 (MQ=255)

                                                                                                                                                                                                                                                         tcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGTTCTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTgc  >  2:687800/1‑274 (MQ=255)

                                                                                                                                                                                                                                                           ccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCa                                                                                                                             <  6:667993/149‑1 (MQ=255)

                                                                                                                                                                                                                                                              tcATCCCGCGCAGC‑GCGCGGGATTGTTCAAGTAAGTTCCATAAAAACACTTATTCATCTCTATAACCCCTGAAAAATAGGGATGTGTTTTTTTATATATTCTGAATATTCACAGTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCa                                                                                                                          <  6:321561/149‑1 (MQ=255)

                                                                                                                                                                                                                                                               cATCCCGCGC‑GCTGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGaa                                                                                                                                                                                                      >  5:1042222/1‑72 (MQ=38)

                                                                                                                                                                                                                                                               cATCCCGCGC‑GCTGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTaga                                                                                                                                        <  5:666342/134‑1 (MQ=255)

                                                                                                                                                                                                                                                                aTCCCGCGCGCT‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAg                                                                                                                         >  8:179060/1‑148 (MQ=255)

                                                                                                                                                                                                                                                                aTCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAAc                                                                                                                                                >  2:768705/1‑125 (MQ=255)

                                                                                                                                                                                                                                                                 tCCCGCGCGCT‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATAtt                                                                                                                                                                           >  8:71556/1‑97 (MQ=255)

                                                                                                                                                                                                                                                                 tCCCGCGCGCT‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAgg                                                                                                                        >  8:421212/1‑148 (MQ=255)

                                                                                                                                                                                                                                                                 tCCCGCGCGCT‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAgg                                                                                                                        >  7:656755/1‑148 (MQ=255)

                                                                                                                                                                                                                                                                 tCCCGCGCAGC‑GCTCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAgg                                                                                                                        <  8:1028327/148‑1 (MQ=255)

                                                                                                                                                                                                                                                                    cgcgCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACtttt                                                                                                                                            <  4:768705/125‑1 (MQ=255)

ATGCATCATGGCACACCAGAACAGGCACAGATGATCCGTACCGCCATCGAACAGGGTAACGGTCGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCC‑GGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGTTCTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGC  >  NC_000913/3334439‑3334960

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