BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,971,151	T→C	5.4%	intergenic (+18/‑120)	ygdQ → / → ygdR	UPF0053 family inner membrane protein/DUF903 family verified lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,971,151	0	T	C	5.4%	804.7 / 22.9	240	intergenic (+18/‑120)	ygdQ/ygdR	UPF0053 family inner membrane protein/DUF903 family verified lipoprotein
Reads supporting (aligned to +/- strand): ref base T (115/112); new base C (4/9); total (119/121)
Fisher's exact test for biased strand distribution p-value = 2.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-02

Only 100 of 241 total aligned reads displayed.

TTTATTATGATGGCAGCCGTGGTGATTGCCGTAGGCGTGATGATGTTCGCCGCGCGCTCGATTGGTGATTTTGTCGAACGCCATCCTTCGGTAAAAATGCTGGCGCTCTCTTTCCTGATTCTGGTGGGCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATAATTTCTGCAGTCGGACTGGCGTTTGCTGTTTCCGGGTGTTCCAGTGATTACGTCATGGCGACCAAAGATGGCCGTATGATTTTGACCGATGGAAAACCTGAAATTGATGATGATACCGGGCTG  >  NC_000913/2970882‑2971414

tttATTATGATGGCAGCCGTGGTGATTGCCGTAGGCGTGATGATGTTCGCCGCGCGCTCGATTGGTGATTTTGTCGAACGCCATCCTTCGGTAAAAATGCTGGCGCTCTCTTTCCTGATTCTGGTGGGCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCggg                                                                                                                                                                                                                                                                    <  2:520190/275‑1 (MQ=255)

          tagCAGCCGTGGTGATTGCCGTAGGCGTGATGATGTTCGCCGCGCGCTCGATTGGTGATTTTGTCGAACGCCATCCTTCGGTAAAAATGCTGGCGCTCTCTTTCCTGATTCTGGTGGGCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCACCAAAAAGAATCCGCCCTGATTCTCCTTACGCTTCCCTGCGGGAGGGcttcc                                                                                                                                                                                                                                                           >  4:268636/3‑270 (MQ=255)

            gCAGCCGTGGTGATTGTCGTAGGCGGGTTGATTTTGGCCGCGCGCTCGAGTTGGGATTTTGGCGAACGGCATACTTCGGTAAAAATGCTGGCGCTCTCTTTCCTGATTCTGGTGGGCTTTACCCTGATTCTGGAAAGTTTCGACATCAAAGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTactca                                                                                                                                                                                                                                                        <  4:330709/275‑1 (MQ=255)

            gCAGCCGTGGTGATTGCCGTAGGCGTGATGATGTTCGCCGCGCGCTCGATTGGTGATTTTGTCGAACGCCATCCTTCGGTAAAAATGCTGGCGCTCTCTTTCCTGATTCTGGTGGGCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTactca                                                                                                                                                                                                                                                        <  2:351840/275‑1 (MQ=255)

              aGCCGTGGTGATTGCCGTAGGCGTGATGATGTTCGCCGCGCGCTCGATTGGTGATTTTGTCGAACGCCATCCTTCGGTAAAAATGCTGGCGCTCTCTTTCCTGATTCTGGTGGGCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCggg                                                                                                                                                                                                                                                                    <  2:268636/261‑1 (MQ=255)

              aGCCGTGGTGATTGCCGTAGGCGTGATGATGTTCGCCGCGCGCTCGATTGGTGATTTTGTCGAACGCCATCCTTCGGTAAAAATGCTGGCGCTCTCTTTCCTGATTCTGGTGGGCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCggg                                                                                                                                                                                                                                                                    <  2:422091/261‑1 (MQ=255)

                                aGGCGTGATGATGTTCGCCGCGCGCTCGATTGGTGAGTTTGTCGAACGCCATCCTTCGGTAAAAATGCTGGCGCTCTCTTTCCTGATTCTGGTGGGCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGtt                                                                                                                                                                                                                                    <  2:74282/275‑1 (MQ=255)

                                       atgatgTTCGCCGCGCGCTCGATTGGTGATTTTGTCGAACGCCATCCTTCGGTAAAAATGCTGGCGCTCTCTTTCCTGATTCTGGTGGGCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACNTCTACTTCGCGATGTTCTTCTTTATTTCGGTTGACAGCCACAATTTGCTTCGCAACAAAAAGAATCCGCCCTGANTCTCCGTACGCTCTCCTGCGGGCGGGCGTCNNCCTACCATTTACCCTCCCACTTCAATAtt                                                                                                                                                                                                                             >  3:89208/1‑275 (MQ=255)

                                               cgccgcGCGCTCGATTGGTGATTTTGTCGAACGCCATCCTTCGGTAAAAATGCTGGCGCTCTCTTTCCTGATTCTGGTGGGCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCAGTACGCTCTCCTGCGGGAGGGCGCACTCACCTCTTTTACCCTCACAGTTTAAGATTTTCCTGc                                                                                                                                                                                                                      >  2:292638/1‑274 (MQ=255)

                                                   gcgcgcTCGATTGGTGATTTTGTCGAACGCCATCCTTCGGTAAAAATGCTGGCGCTCTCTTTCCTGATTCTGGTGGGCTTTACCCTGATTCTGGAAAGTTTCGACAGCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTc                                                                                                                                                                                                                  <  1:89208/274‑1 (MQ=255)

                                                    cgcgcTCGATTGGTGATTTTGTCGAACGCCATCCTTCGGTAAAAATGCTGGCGCTCTCTTTCCTGATTCTGGTGGGCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCa                                                                                                                                                                                                                 <  4:237768/274‑1 (MQ=255)

                                                                                                     ggCGCTCTCTTTCCTGATTCTGGTGGGCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCt                                                                                                                                                                                                                     >  2:205621/1‑221 (MQ=255)

                                                                                                         ctctctTTCCTGATTCTGGTGGGCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCa                                                                                                                                                                                                                 <  4:205621/221‑1 (MQ=255)

                                                                                                                      ttCTGGTGGGCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCt                                                                                                                                                                                                                        >  4:217002/1‑201 (MQ=255)

                                                                                                                          ggtggGCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCtt                                                                                                                                                                                                                    <  2:217002/201‑1 (MQ=255)

                                                                                                                                 tttACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAgg                                                                                                                                                                                                                                                                 <  6:393740/149‑1 (MQ=255)

                                                                                                                                 tttACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAgg                                                                                                                                                                                                                                                                 >  7:802309/1‑149 (MQ=255)

                                                                                                                                 tttACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAgg                                                                                                                                                                                                                                                                 <  7:328545/149‑1 (MQ=255)

                                                                                                                                       cTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTac                                                                                                                                                                                                                                                           >  7:602209/1‑149 (MQ=255)

                                                                                                                                       cTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTac                                                                                                                                                                                                                                                           <  5:1020009/149‑1 (MQ=255)

                                                                                                                                        tGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTact                                                                                                                                                                                                                                                          >  5:338310/1‑149 (MQ=255)

                                                                                                                                        tGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTact                                                                                                                                                                                                                                                          >  6:304093/1‑149 (MQ=255)

                                                                                                                                        tGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTact                                                                                                                                                                                                                                                          >  7:623142/1‑149 (MQ=255)

                                                                                                                                           ttCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTactcac                                                                                                                                                                                                                                                       <  7:890755/149‑1 (MQ=255)

                                                                                                                                              tGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAgg                                                                                                                                                                                                                                                                 <  6:894247/136‑1 (MQ=255)

                                                                                                                                              tGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTacgga                                                                                                                                                                                                                                                        <  5:588581/145‑4 (MQ=255)

                                                                                                                                                gAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCggg                                                                                                                                                                                                                                                                    >  2:239615/1‑131 (MQ=255)

                                                                                                                                                      ttCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTAccc                                                                                                                                                                                                                                            >  5:769148/1‑149 (MQ=255)

                                                                                                                                                       tCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCt                                                                                                                                                                                                                                           <  7:972891/149‑1 (MQ=255)

                                                                                                                                                               cACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGttt                                                                                                                                                                                                                                   >  5:62887/1‑149 (MQ=255)

                                                                                                                                                                     ccGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGAtt                                                                                                                                                                                                                             >  6:472422/1‑149 (MQ=255)

                                                                                                                                                                      cGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGAttt                                                                                                                                                                                                                            >  7:869363/1‑149 (MQ=255)

                                                                                                                                                                      cGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGAttt                                                                                                                                                                                                                            >  5:886453/1‑149 (MQ=255)

                                                                                                                                                                      cGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGAttt                                                                                                                                                                                                                            >  7:775118/1‑149 (MQ=255)

                                                                                                                                                                       gAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGAtttt                                                                                                                                                                                                                           >  6:556217/1‑149 (MQ=255)

                                                                                                                                                                       gAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGAtttt                                                                                                                                                                                                                           >  7:518142/1‑149 (MQ=255)

                                                                                                                                                                        aaaGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTc                                                                                                                                                                                                                          >  5:873673/1‑149 (MQ=255)

                                                                                                                                                                        aaaGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTc                                                                                                                                                                                                                          <  6:122976/149‑1 (MQ=255)

                                                                                                                                                                         aaGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTTAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATAATTTCTGCAGTCGGACTGGCGTTTGCTGTTTcc                                                                                            <  2:596098/274‑1 (MQ=255)

                                                                                                                                                                               aCATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCttt                                                                                                                                                                                                                   >  2:181618/1‑149 (MQ=255)

                                                                                                                                                                                  tCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCaa                                                                                                                                                                                                                <  7:833174/149‑1 (MQ=255)

                                                                                                                                                                                  tCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCaa                                                                                                                                                                                                                >  8:267258/1‑149 (MQ=255)

                                                                                                                                                                                  tCTACTTCGCGATGTTCTTCTCTATTGCGGGTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCaa                                                                                                                                                                                                                >  6:966737/1‑149 (MQ=255)

                                                                                                                                                                                   cTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAg                                                                                                                                                                                                               <  4:181618/149‑1 (MQ=255)

                                                                                                                                                                                     aCTTCGCGATGTTCTTATCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGat                                                                                                                                                                                                             >  6:368718/1‑149 (MQ=255)

                                                                                                                                                                                      cTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTg                                                                                                                                                                                                                       >  4:576887/1‑138 (MQ=255)

                                                                                                                                                                                       ttCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGatat                                                                                                                                                                                                           >  7:557730/1‑149 (MQ=255)

                                                                                                                                                                                        tCGCGATGTGCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGatata                                                                                                                                                                                                          <  5:472422/149‑1 (MQ=255)

                                                                                                                                                                                         cgcgATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCCCAGTTTAAGATTTTCCTGCTTTCAAGatatat                                                                                                                                                                                                         <  6:1012640/149‑1 (MQ=255)

                                                                                                                                                                                             aTGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACg                                                                                                                                                                                                     >  8:386665/1‑149 (MQ=255)

                                                                                                                                                                                                  cttctCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGt                                                                                                                                                                                                <  8:620869/149‑1 (MQ=255)

                                                                                                                                                                                                   ttctCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGtt                                                                                                                                                                                               <  5:936070/149‑1 (MQ=255)

                                                                                                                                                                                                    tctcTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGttt                                                                                                                                                                                              >  7:259366/1‑149 (MQ=255)

                                                                                                                                                                                                        tATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCt                                                                                                                                                                                                                        >  4:289521/1‑119 (MQ=255)

                                                                                                                                                                                                        tATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATaa                                                                                                                                                                                          >  6:185195/1‑149 (MQ=255)

                                                                                                                                                                                                        tATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATaa                                                                                                                                                                                          >  8:412127/1‑149 (MQ=255)

                                                                                                                                                                                                             cGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTAcc                                                                                                                                                                                                                                             >  2:662835/1‑93 (MQ=255)

                                                                                                                                                                                                              ggTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGAc                                                                                                                                                                                    <  7:404826/149‑1 (MQ=255)

                                                                                                                                                                                                                 tGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTca                                                                                                                                                                                                                                         <  4:662835/93‑1 (MQ=255)

                                                                                                                                                                                                                  gAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGt                                                                                                                                                                                                                                     >  1:39046/1‑96 (MQ=255)

                                                                                                                                                                                                                   aaaGCCGCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATAATTTCTGCAGTCGGACTGGCGTTTGCTGTTTCCGGGTGTTCCAGTGATTACGTCATGGCGACCAAAGATGGCCGTa                                                 <  1:151736/275‑1 (MQ=255)

                                                                                                                                                                                                                     aGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATc                                                                                                                                                                             <  7:708760/149‑1 (MQ=255)

                                                                                                                                                                                                                       ccTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAt                                                                                                                                                                           >  7:1038142/1‑149 (MQ=255)

                                                                                                                                                                                                                         tCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtag                                                                                                                                                                         <  8:869363/149‑1 (MQ=255)

                                                                                                                                                                                                                         tCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtag                                                                                                                                                                         <  8:802309/149‑1 (MQ=255)

                                                                                                                                                                                                                          cAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtagg                                                                                                                                                                        <  6:612974/149‑1 (MQ=255)

                                                                                                                                                                                                                           aaCTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAAc                                                                                                                                                                                                      <  8:349559/118‑1 (MQ=255)

                                                                                                                                                                                                                           aaCTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtaggt                                                                                                                                                                       >  6:2767/1‑149 (MQ=255)

                                                                                                                                                                                                                              ttGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATc                                                                                                                                                                             <  8:496029/140‑1 (MQ=255)

                                                                                                                                                                                                                                    cGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAAt                                                                                                                                                              <  8:518142/149‑1 (MQ=255)

                                                                                                                                                                                                                                      caacaaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGa                                                                                                                                                            >  7:290931/1‑149 (MQ=255)

                                                                                                                                                                                                                                       aacaaAAAGAATCCGCTCTGATACTCCGTACTCTCTCCTGCGGGAGGGCGTACTCGCTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAg                                                                                                                                                           <  6:852063/149‑1 (MQ=255)

                                                                                                                                                                                                                                       aacaaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAg                                                                                                                                                           <  7:386665/149‑1 (MQ=255)

                                                                                                                                                                                                                                       aacaaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAg                                                                                                                                                           <  8:444661/149‑1 (MQ=255)

                                                                                                                                                                                                                                       aacaaAAAGAATCCGCTCTCATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGACAAATGAg                                                                                                                                                           >  7:635305/1‑149 (MQ=255)

                                                                                                                                                                                                                                        acaaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAgg                                                                                                                                                          <  5:385348/149‑1 (MQ=255)

                                                                                                                                                                                                                                            aaaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTaa                                                                                                                                                                                                                                 >  7:730261/1‑74 (MQ=255)

                                                                                                                                                                                                                                                     gCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATgca                                                                                                                                             >  8:250512/1‑147 (MQ=255)

                                                                                                                                                                                                                                                       tctGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGaaaa                                                                                                                                           >  6:786434/1‑149 (MQ=255)

                                                                                                                                                                                                                                                          gATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAAt                                                                                                                                        <  8:603999/149‑1 (MQ=255)

                                                                                                                                                                                                                                                             aCTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATggg                                                                                                                                     <  8:314028/149‑1 (MQ=255)

                                                                                                                                                                                                                                                               tCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATggg                                                                                                                                     >  5:822062/1‑147 (MQ=255)

                                                                                                                                                                                                                                                               tCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATggg                                                                                                                                     <  6:822062/147‑1 (MQ=255)

                                                                                                                                                                                                                                                                 cGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATaa                                                                                                                                                                                                       >  3:165531/1‑79 (MQ=255)

                                                                                                                                                                                                                                                                   tACGCTTTCCTGCGGGAGGGCGTACTCACTCCTTTTGCCCTCACAGTTCAAGATTTTCCTGCTTTCACGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCTAAATGAGGTTACCCGGATGAAAAAATGGGCAGTAATAATTTCTGCAGTCGGACTGGCGTTTGCTGTTTCCGGGTGTTCAAGTGTTTCCGTCATGGCTACCAAAGACGGCCGTTACATTTTGCCCGATGGAAAACTGGAAATTGATGATGTTACCGGGTTg  >  4:86935/1‑274 (MQ=255)

                                                                                                                                                                                                                                                                   tACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGa                                                                                                                                                                                     >  7:1052902/1‑95 (MQ=255)

                                                                                                                                                                                                                                                                   tACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGa                                                                                                                                                                                     <  8:1052902/95‑1 (MQ=255)

                                                                                                                                                                                                                                                                   tACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGaaaaaa                                                                                                                                         >  4:517641/1‑139 (MQ=255)

                                                                                                                                                                                                                                                                    aCGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGtaat                                                                                                                              <  8:775118/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                     cGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGc                                                                                                                                                                                                                      >  3:103923/1‑60 (MQ=255)

                                                                                                                                                                                                                                                                     cGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTc                                                                                                                                                                                                   <  1:165531/79‑1 (MQ=255)

                                                                                                                                                                                                                                                                      gCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGtaataa                                                                                                                            >  8:658428/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                        tctcCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATg                                                                                                                                                             <  7:482585/114‑1 (MQ=255)

                                                                                                                                                                                                                                                                        tctcCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATAAtt                                                                                                                          >  8:20019/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                         ctcCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTc                                                                                                                                                                                                                  <  1:103923/60‑1 (MQ=255)

                                                                                                                                                                                                                                                                          tcCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGaaaaaa                                                                                                                                         >  4:395911/1‑132 (MQ=255)

                                                                                                                                                                                                                                                                          tcCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATAATTTc                                                                                                                        <  8:74174/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                          tcCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATAATTTc                                                                                                                        >  6:461463/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                            cTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACta                                                                                                                                                                                  <  2:565514/89‑1 (MQ=255)

                                                                                                                                                                                                                                                                            cTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATAATTTCTg                                                                                                                      >  5:149548/1‑149 (MQ=255)

TTTATTATGATGGCAGCCGTGGTGATTGCCGTAGGCGTGATGATGTTCGCCGCGCGCTCGATTGGTGATTTTGTCGAACGCCATCCTTCGGTAAAAATGCTGGCGCTCTCTTTCCTGATTCTGGTGGGCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATAATTTCTGCAGTCGGACTGGCGTTTGCTGTTTCCGGGTGTTCCAGTGATTACGTCATGGCGACCAAAGATGGCCGTATGATTTTGACCGATGGAAAACCTGAAATTGATGATGATACCGGGCTG  >  NC_000913/2970882‑2971414

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