BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	773,553	T→C	5.1%	L171P (CTG→CCG) ‡	cydB →	cytochrome d terminal oxidase, subunit II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	773,553	0	T	C	5.1%	732.3 / 17.3	217	L171P (CTG→CCG) ‡	cydB	cytochrome d terminal oxidase, subunit II
Reads supporting (aligned to +/- strand): ref base T (102/103); new base C (6/5); total (109/108)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

Only 100 of 229 total aligned reads displayed.

CTGCGTTCTCCGGCTTCTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCACTGGCTGGCGTATGGGTGATGTACGGTATCGATGGTTATGTCGTGAAATCGACAATGGACCATTACGCAGCCTCTAACCCACTGAATAAAGAAGTGGTTCGTGAAGCTGGCGCATGGCTGGTTAACTTCAACAACA  >  NC_000913/773283‑773825

                    gtggtgATGATCCTCGTGCGTGCTTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAGGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGgat                                                                                                                                                                                                                                                          <  2:81013/275‑1 (MQ=255)

                    gtgGCGTTGATCATAGTGCTGGCGTCTTGGTGCTTCCGTCCGGTCGGTTTTGACTGCCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGgat                                                                                                                                                                                                                                                          <  2:192877/275‑1 (MQ=255)

                    ggggcgtggttCCTCGTGCTGTCGTTTTTGTTCTTCCGTCCGGTCGGTTTTGTATCCCGCTCCAAGTTTGACGAAACCCGCGGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCAGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGgat                                                                                                                                                                                                                                                          <  4:13652/265‑1 (MQ=255)

                                                ttttcttcCGTCCGGTCGGTTGTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTg                                                                                                                                                                                                                              <  2:13988/273‑1 (MQ=255)

                                                                                tgcAAGATTGAAGAAGGCCTCTGGCGTACAATGTGGTACTGGGGTCTTGTCATTGGTATCTTCGTTCCGCCGCTGGTAGTTGGTCTAGCGTTCGGTAACCTGTTGCAGGGGTTACCGTTTCACGTTGATGAATATCTGCGTCTTTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCg                                                                                                                                                                                              <  1:48764/273‑1 (MQ=255)

                                                                                                                gtgGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgtggt                                                                                                                                                                                                                                                                      <  4:295395/171‑1 (MQ=255)

                                                                                                                           gCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCctgctg                                                                                                                                                                                                                                                                                 >  7:730372/1‑149 (MQ=255)

                                                                                                                           gCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCctgctg                                                                                                                                                                                                                                                                                 >  5:1034745/1‑149 (MQ=255)

                                                                                                                             aTCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGc                                                                                                                                                                                                                                                                               >  5:1075393/1‑149 (MQ=255)

                                                                                                                              tCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTGTCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTg                                                                                                                                                 >  4:68524/1‑274 (MQ=255)

                                                                                                                                 tCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGc                                                                                                                                                                                                                                                                           <  6:679731/149‑1 (MQ=255)

                                                                                                                                       gTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgtggtg                                                                                                                                                                                                                                                                     >  5:434305/1‑149 (MQ=255)

                                                                                                                                       gTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgtgg                                                                                                                                                                                                                                                                       >  6:1140500/1‑147 (MQ=255)

                                                                                                                                       gTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgtgg                                                                                                                                                                                                                                                                       <  5:1140500/147‑1 (MQ=255)

                                                                                                                                            ttCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGt                                                                                                                                                                                                                                                                <  6:242025/149‑1 (MQ=255)

                                                                                                                                              cGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAg                                                                                                                                                                                                                                                              <  6:150540/149‑1 (MQ=255)

                                                                                                                                              cGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAg                                                                                                                                                                                                                                                              <  5:759627/149‑1 (MQ=255)

                                                                                                                                                ttCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAggg                                                                                                                                                                                                                                                            >  5:1188627/1‑149 (MQ=255)

                                                                                                                                                 tCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCACCGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCACTGGCTGGCGttt                                                                                                                             >  2:64178/1‑273 (MQ=255)

                                                                                                                                                  ccgccgCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGgat                                                                                                                                                                                                                                                          <  5:335049/149‑1 (MQ=255)

                                                                                                                                                     ccgcTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGgatgat                                                                                                                                                                                                                                                       <  6:1220201/149‑1 (MQ=255)

                                                                                                                                                     ccgcTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGgatgat                                                                                                                                                                                                                                                       <  7:590193/149‑1 (MQ=255)

                                                                                                                                                     ccgcTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGgatgat                                                                                                                                                                                                                                                       <  6:1365367/149‑1 (MQ=255)

                                                                                                                                                         gggTAATTGGGGGAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCGGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCACTGGCTGGCGTATGGGTGATg                                                                                                                     <  4:122244/274‑1 (MQ=255)

                                                                                                                                                          ggTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTa                                                                                                                                                                                                                                                  >  7:795874/1‑149 (MQ=255)

                                                                                                                                                           gTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTAc                                                                                                                                                                                                                                                 <  8:753106/149‑1 (MQ=255)

                                                                                                                                                           gTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCCAGCAGGccgaac                                                                                                                                                                                                                                                                      >  8:150057/1‑123 (MQ=255)

                                                                                                                                                             aaTTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCACTGGCTGGCGTATGGGTGATGTACg                                                                                                                 <  2:323570/275‑1 (MQ=255)

                                                                                                                                                              aTTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCa                                                                                                                                                                                                                                              <  8:849882/149‑1 (MQ=255)

                                                                                                                                                              aTTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCa                                                                                                                                                                                                                                              <  5:30324/149‑1 (MQ=255)

                                                                                                                                                              aTTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCa                                                                                                                                                                                                                                              <  5:422426/149‑1 (MQ=255)

                                                                                                                                                                      aGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCCAGCAGGccgaac                                                                                                                                                                                                                                                                      >  7:202791/1‑112 (MQ=255)

                                                                                                                                                                         gTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgtggt                                                                                                                                                                                                                                                                      >  5:240351/1‑114 (MQ=255)

                                                                                                                                                                         gTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCt                                                                                                                                                                                                                                   >  6:985380/1‑149 (MQ=255)

                                                                                                                                                                         gTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCt                                                                                                                                                                                                                                   >  7:883249/1‑149 (MQ=255)

                                                                                                                                                                         gTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGAGAAGGGATGATCATTACTCAGGGCGCa                                                                                                                                                                                                                                       >  6:332112/1‑145 (MQ=255)

                                                                                                                                                                          ttCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgtg                                                                                                                                                                                                                                                                        >  8:1269354/1‑111 (MQ=255)

                                                                                                                                                                             ggTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGTCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACACAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCGGCGTACCCGTGCAACGGCTCAGGTGGATGCGCGGGTGACAACGGGCTGTTTCGCCCTGGCTGGCGTATGGGTGGTGTACGGTAACGAGGAGTATGt                                                                                                 >  4:246878/1‑275 (MQ=255)

                                                                                                                                                                              gTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGTTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTg                                                                                                                                                                                                                              >  6:353207/1‑149 (MQ=255)

                                                                                                                                                                               tAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGGAGGGATGATCATTACTCAGGGCGCAACCTATCTGc                                                                                                                                                                                                                             >  7:1123265/1‑149 (MQ=255)

                                                                                                                                                                                  ccTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCaaa                                                                                                                                                                                                                          >  6:1260578/1‑149 (MQ=255)

                                                                                                                                                                                   cTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAAt                                                                                                                                                                                                                         >  8:208657/1‑149 (MQ=255)

                                                                                                                                                                                   cTGTTGCAGGGCGTACCGTTCAACGTGGAGGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCTGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAAt                                                                                                                                                                                                                         <  8:379828/149‑1 (MQ=255)

                                                                                                                                                                                    tGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATc                                                                                                                                                                                                                                >  6:239999/1‑141 (MQ=255)

                                                                                                                                                                                     gTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGa                                                                                                                                                                                                                                                                    <  7:376370/104‑1 (MQ=255)

                                                                                                                                                                                     gTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGc                                                                                                                                                                                                                       <  5:285093/149‑1 (MQ=255)

                                                                                                                                                                                      ttGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCAc                                                                                                                                                                                                   <  4:188401/168‑1 (MQ=255)

                                                                                                                                                                                         cAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTAc                                                                                                                                                                                                                   >  5:881597/1‑149 (MQ=255)

                                                                                                                                                                                          aGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTCGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTAcc                                                                                                                                                                                                                  >  7:1259141/1‑149 (MQ=255)

                                                                                                                                                                                                tACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTTGTGAGCGTAGGGATGATCTTTACTCAGGGCGCACCCTATCTGCAACTGCGTACCGTGGGCGATCTGCACCTGCGTACCCGTGCAACGGTTCAGGTTGCTGCGCAGGTGCCTCTGGTCTGTTTCGACATGGCTGGACTATGGGTGATCTACGGTTGCGATGGTTATGTCGTATACTCGACATTGGa                                                                               >  4:55218/1‑274 (MQ=255)

                                                                                                                                                                                                tACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGc                                                                                                                                                                                                            >  8:312906/1‑149 (MQ=255)

                                                                                                                                                                                                    gTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAAc                                                                                                                                                                                                        >  5:918064/1‑149 (MQ=255)

                                                                                                                                                                                                     ttCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACt                                                                                                                                                                                                       <  6:3935/149‑1 (MQ=255)

                                                                                                                                                                                                       cAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAg                                                                                                                                                                                                                                                              >  7:501272/1‑92 (MQ=255)

                                                                                                                                                                                                          cGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCa                                                                                                                                                                                                    <  2:163088/147‑1 (MQ=255)

                                                                                                                                                                                                           gTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCt                                                                                                                                                                                                 >  7:666417/1‑149 (MQ=255)

                                                                                                                                                                                                             tgatgaATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGc                                                                                                                                                                                               <  7:312906/149‑1 (MQ=255)

                                                                                                                                                                                                             tgatgaATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGc                                                                                                                                                                                               <  6:795265/149‑1 (MQ=255)

                                                                                                                                                                                                               atgaATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGt                                                                                                                                                                                             >  5:180720/1‑149 (MQ=255)

                                                                                                                                                                                                                 gaATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTAc                                                                                                                                                                                           <  8:348031/149‑1 (MQ=255)

                                                                                                                                                                                                                 gaATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTAc                                                                                                                                                                                           >  8:498149/1‑149 (MQ=255)

                                                                                                                                                                                                                  aaTATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTAcc                                                                                                                                                                                          >  8:1140921/1‑149 (MQ=255)

                                                                                                                                                                                                                         gCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCACTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCag                                                                                                                                                                                   <  6:729470/149‑2 (MQ=255)

                                                                                                                                                                                                                           ggttGTACTGCACCGGTAACTTCATGCAGTTGCTTAACCCGTTCGGGCTGCTGGCAGGCGTGGTAAGCGGAGGGATGATCATTACGCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGAACCTGCGTACCCGTGCAACGGCTCAAGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCACTGGCTGGCGTATGGGTGATGTACGGTATCGATGGTTATGTCGTGAAATCGACAATGGACCATTACGCAGCCTCTAACCCACTGaa                                                    <  2:55218/271‑1 (MQ=255)

                                                                                                                                                                                                                             cTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGc                                                                                                                                                                               <  8:730372/149‑1 (MQ=255)

                                                                                                                                                                                                                              tGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCt                                                                                                                                                                              >  6:765590/1‑149 (MQ=255)

                                                                                                                                                                                                                              tGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCt                                                                                                                                                                              >  8:831095/1‑149 (MQ=255)

                                                                                                                                                                                                                               gTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGaa                                                                                                                                                                                                         >  7:157672/1‑121 (MQ=255)

                                                                                                                                                                                                                                tactacACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCa                                                                                                                                                                            >  8:1065114/1‑149 (MQ=255)

                                                                                                                                                                                                                                   tacACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGt                                                                                                                                                                         >  8:928860/1‑149 (MQ=255)

                                                                                                                                                                                                                                        cGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTa                                                                                                                                                                                                                                                  <  6:487740/71‑1 (MQ=255)

                                                                                                                                                                                                                                        cGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTa                                                                                                                                                                                                                                                  >  5:487740/1‑71 (MQ=255)

                                                                                                                                                                                                                                        cGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTg                                                                                                                                                                    >  8:15796/1‑149 (MQ=255)

                                                                                                                                                                                                                                            aaCTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCt                                                                                                                                                                >  6:1024435/1‑149 (MQ=255)

                                                                                                                                                                                                                                            aaCTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCACTGGCTGGCGTATGGGTGATGTACGGTATCGATGGTTATGTCGTGAAATCGACAATGGACCATTACGCa                                                                     <  4:14282/240‑1 (MQ=255)

                                                                                                                                                                                                                                              cttcttCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTgg                                                                                                                                                              <  5:985380/149‑1 (MQ=255)

                                                                                                                                                                                                                                               ttcttcCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGt                                                                                                                                                             <  6:884063/149‑1 (MQ=255)

                                                                                                                                                                                                                                               ttcttcCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGt                                                                                                                                                             <  7:241938/149‑1 (MQ=255)

                                                                                                                                                                                                                                               ttcttcCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGt                                                                                                                                                             <  8:239133/149‑1 (MQ=255)

                                                                                                                                                                                                                                               ttcttcCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACCCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCCACGGCTCAGGTGGCTgcgc                                                                                                                                                                 >  5:1438806/1‑145 (MQ=255)

                                                                                                                                                                                                                                                tcttcCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTg                                                                                                                                                            >  5:1371749/1‑149 (MQ=255)

                                                                                                                                                                                                                                                tcttcCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTg                                                                                                                                                            <  6:434305/149‑1 (MQ=255)

                                                                                                                                                                                                                                                  ttcCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGac                                                                                                                                                          >  7:940650/1‑149 (MQ=255)

                                                                                                                                                                                                                                                  ttcCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGac                                                                                                                                                          >  6:704905/1‑149 (MQ=255)

                                                                                                                                                                                                                                                      aGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGaca                                                                                                                                                         >  6:1327603/1‑146 (MQ=255)

                                                                                                                                                                                                                                                       gTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTgg                                                                                                                                                     >  5:1262/1‑149 (MQ=255)

                                                                                                                                                                                                                                                        ttGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGt                                                                                                                                                    <  5:1472680/149‑1 (MQ=255)

                                                                                                                                                                                                                                                           cTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTg                                                                                                                                                 <  6:296567/149‑1 (MQ=255)

                                                                                                                                                                                                                                                              aaCCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGttt                                                                                                                                              <  5:582512/149‑1 (MQ=255)

                                                                                                                                                                                                                                                               aCCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTc                                                                                                                                             >  7:149723/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                 ccGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGc                                                                                                                                           >  5:201690/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                   gTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAAt                                                                                                                                                                                                                         <  8:601549/69‑1 (MQ=255)

                                                                                                                                                                                                                                                                   gTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGCGGCCGAATTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCAc                                                                                                                                         >  6:1055327/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                     cacgCCTGCCAGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGc                                                                                                                                                                                                                       <  5:1129173/66‑1 (MQ=255)

                                                                                                                                                                                                                                                                     cacgCCTGCCAGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGc                                                                                                                                                                                                                       >  6:1129173/4‑69 (MQ=255)

                                                                                                                                                                                                                                                                     cacgCCTGCCAGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGt                                                                                                                                                                         <  5:662538/112‑1 (MQ=255)

                                                                                                                                                                                                                                                                     cacgCCTGCCAGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCActg                                                                                                                                       >  7:1292501/4‑149 (MQ=255)

                                                                                                                                                                                                                                                                     cacgCCTGCCAGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCActg                                                                                                                                       >  6:473349/4‑149 (MQ=255)

                                                                                                                                                                                                                                                                          ctgctgGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCActggctgg                                                                                                                                  >  5:1282700/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                           tgctgGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGATGCGCTGGTGACACTGGTCTGTTTCGCACtggctggc                                                                                                                                 >  7:1452823/1‑149 (MQ=255)

CTGCGTTCTCCGGCTTCTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCACTGGCTGGCGTATGGGTGATGTACGGTATCGATGGTTATGTCGTGAAATCGACAATGGACCATTACGCAGCCTCTAACCCACTGAATAAAGAAGTGGTTCGTGAAGCTGGCGCATGGCTGGTTAACTTCAACAACA  >  NC_000913/773283‑773825

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