BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,698,587	T→A	5.8%	A198A (GCA→GCT) ‡	malI ←	transcriptional repressor of Mal regulon

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,698,587	0	T	A	5.8%	792.8 / 28.3	223	A198A (GCA→GCT) ‡	malI	transcriptional repressor of Mal regulon
Reads supporting (aligned to +/- strand): ref base T (101/109); new base A (5/8); total (106/117)
Fisher's exact test for biased strand distribution p-value = 5.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.07e-01

Only 100 of 224 total aligned reads displayed.

CGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCACTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCTGCAGCTCCGGCAATTACCACACCATCGACACCCTGATTCAGTAACAGTGAAAACCGCTGTGCCAGC  >  NC_000913/1698316‑1698857

cGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCACTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGc                                                                                                                                                                                                                                                                              <  2:78245/274‑1 (MQ=255)

      agcccgcaAGCGAAACCAGTCTCTCAACGTCACGCTGGACTCCCCTTTCCCCGCGTTGGCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGAAATCGTTTCGTTATAGCAAACCACGGCACTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGt                                                                                                                                                                                                                                                                        <  4:275961/271‑1 (MQ=255)

                         ttGGTCAAAGTAAAGCTCGACTCGGCTGTCCCCTCGTTGCCGCCCGTCTTTGGGCAAAACAAACCATGCCCCCTTCGCAATAGTTTCGTTATAGCAAACCAGGGCACTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTGGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTc                                                                                                                                                                                                                                                    <  4:11655/275‑1 (MQ=255)

                                                 gCTTTCCCCGCTTGTCCGCCCTGCTTTGATCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCACTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGcc                                                                                                                                                                                                                            <  2:273989/275‑1 (MQ=255)

                                                           ttttGCCTCCCCGATTTGAACAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCACTGATGGTCGGGTTATGACGTAAAAGCGCCGTTATAGCTTCCTCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAAt                                                                                                                                                                                                                  <  1:38504/274‑1 (MQ=255)

                                                                tcgccCTGCTTTGGGCAAACCGAACCAGGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCACTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGTACAATAGCCCCCAACCCTTTCTGCAc                                                                                                                                                                                                                                                                            <  2:82653/211‑1 (MQ=255)

                                                                                                                                 cACTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACgg                                                                                                                                                                                                                                                                          <  5:280439/149‑1 (MQ=255)

                                                                                                                                 cACTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACgg                                                                                                                                                                                                                                                                          <  7:374501/149‑1 (MQ=255)

                                                                                                                                 cACTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCACAAGCCGTTCTGCACgg                                                                                                                                                                                                                                                                          <  6:932059/149‑1 (MQ=255)

                                                                                                                                    tGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGtt                                                                                                                                                                                                                                                                       >  8:416580/1‑149 (MQ=255)

                                                                                                                                       tGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAAt                                                                                                                                                                                                                                                                    <  6:1470598/149‑1 (MQ=255)

                                                                                                                                       tGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTGCAGAACGGGttggg                                                                                                                                                                                                                                                                    <  7:938073/149‑4 (MQ=255)

                                                                                                                                        ggTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTGCAGAACGGGttgggg                                                                                                                                                                                                                                                                   <  5:1029629/149‑5 (MQ=255)

                                                                                                                                         gTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGa                                                                                                                                                                                                                                                                  <  6:1258899/149‑1 (MQ=255)

                                                                                                                                             ggTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGaa                                                                                                                                                                                                                                                              >  5:641061/1‑149 (MQ=255)

                                                                                                                                             ggTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGaa                                                                                                                                                                                                                                                              >  7:1182692/1‑149 (MQ=255)

                                                                                                                                                tATGACGTAAAAGCGCCGTGATAGCTTGCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAAAACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAATCCGTTCTGCACGGGTTAATGAGGAACTTTGcc                                                                                                                                                                                                                                                       >  1:71755/1‑153 (MQ=255)

                                                                                                                                                tATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACtt                                                                                                                                                                                                                                                           >  8:1224498/1‑149 (MQ=255)

                                                                                                                                                 aTGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACttt                                                                                                                                                                                                                                                          <  6:407219/149‑1 (MQ=255)

                                                                                                                                                   gACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGc                                                                                                                                                                                                                                                        <  6:1052333/149‑1 (MQ=255)

                                                                                                                                                    aCGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGcc                                                                                                                                                                                                                                                       <  6:31899/149‑1 (MQ=255)

                                                                                                                                                      gTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCt                                                                                                                                                                                                                                                     >  7:1161032/1‑149 (MQ=255)

                                                                                                                                                      gTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTGCAGAACGGGttgggg                                                                                                                                                                                                                                                                   <  6:889234/135‑5 (MQ=255)

                                                                                                                                                           aGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCccag                                                                                                                                                                                                                                                <  6:1217917/149‑1 (MQ=255)

                                                                                                                                                               ccGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCccagccag                                                                                                                                                                                                                                            <  7:305573/149‑1 (MQ=255)

                                                                                                                                                               ccGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCccagccag                                                                                                                                                                                                                                            >  5:348663/1‑149 (MQ=255)

                                                                                                                                                                cGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCaacaac                                                                                                                                                                                               >  4:310769/1‑193 (MQ=255)

                                                                                                                                                                   gATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGa                                                                                                                                                                                                                                        <  6:1188135/149‑1 (MQ=255)

                                                                                                                                                                        cTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGc                                                                                                                                                                                                                                   <  5:1093586/149‑1 (MQ=255)

                                                                                                                                                                         ttCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCt                                                                                                                                                                                                                                  >  5:106047/1‑149 (MQ=255)

                                                                                                                                                                            cgcgGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATgaga                                                                                                                                                                                                              >  2:186277/1‑166 (MQ=255)

                                                                                                                                                                             gcgGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATGTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATg                                                                                                                                                                                                                              >  8:242476/1‑149 (MQ=255)

                                                                                                                                                                                gcttgcttCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGccc                                                                                                                                                                                                                           >  5:1391655/1‑149 (MQ=255)

                                                                                                                                                                                gcttgcttCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCt                                                                                                                                                                                                          <  4:186277/166‑1 (MQ=255)

                                                                                                                                                                                      ttCTGGCTGGAAGTGCACTCCATCACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTGCAGAACGGGttgggg                                                                                                                                                                                                                                                                   <  8:517163/103‑5 (MQ=255)

                                                                                                                                                                                      ttCTGGCTGGAAGTGCACTCCAGCACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCg                                                                                                                                                                                                                     >  5:900021/1‑149 (MQ=255)

                                                                                                                                                                                      ttCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCg                                                                                                                                                                                                                     >  7:677870/1‑149 (MQ=255)

                                                                                                                                                                                           gctggAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATga                                                                                                                                                                                                                <  5:45409/149‑1 (MQ=255)

                                                                                                                                                                                           gctggAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATga                                                                                                                                                                                                                <  8:1162342/149‑1 (MQ=255)

                                                                                                                                                                                             tggAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATgaga                                                                                                                                                                                                              <  5:1149247/149‑1 (MQ=255)

                                                                                                                                                                                                 aGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCt                                                                                                                                                                                                          >  5:197296/1‑149 (MQ=255)

                                                                                                                                                                                                 aGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCt                                                                                                                                                                                                          >  7:1266275/1‑149 (MQ=255)

                                                                                                                                                                                                      aCTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTc                                                                                                                                                                                                     >  8:427641/1‑149 (MQ=255)

                                                                                                                                                                                                       cTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTca                                                                                                                                                                                                    >  5:578473/1‑149 (MQ=255)

                                                                                                                                                                                                        tCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTcaa                                                                                                                                                                                                   >  8:541744/1‑149 (MQ=255)

                                                                                                                                                                                                         ccAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCa                                                                                                                                                                                                                                                                             <  4:163370/74‑1 (MQ=255)

                                                                                                                                                                                                         ccAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTcaac                                                                                                                                                                                                  <  7:167954/149‑1 (MQ=255)

                                                                                                                                                                                                          cAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCTCATTGCGAATGAGATGCTCCGTcaaca                                                                                                                                                                                                 >  8:1261153/1‑149 (MQ=255)

                                                                                                                                                                                                          cAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTcaaca                                                                                                                                                                                                 >  8:1137119/1‑149 (MQ=255)

                                                                                                                                                                                                            acacCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCaacaac                                                                                                                                                                                               <  6:1391655/149‑1 (MQ=255)

                                                                                                                                                                                                              acCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTc                                                                                                                                                                                                     <  2:364567/141‑1 (MQ=255)

                                                                                                                                                                                                              acCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACtg                                                                                                                                                                                             >  7:1351058/1‑149 (MQ=255)

                                                                                                                                                                                                               cccAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACtgt                                                                                                                                                                                            <  7:1456415/149‑1 (MQ=255)

                                                                                                                                                                                                                ccAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACtgtg                                                                                                                                                                                           >  8:229253/1‑149 (MQ=255)

                                                                                                                                                                                                                ccAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACtgtg                                                                                                                                                                                           >  6:415495/1‑149 (MQ=255)

                                                                                                                                                                                                                 cAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGaa                                                                                                                                                                                                                   <  7:42721/124‑1 (MQ=255)

                                                                                                                                                                                                                    tCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCATCCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGc                                                                                                                                                                                       <  8:1351058/149‑1 (MQ=255)

                                                                                                                                                                                                                      gCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCt                                                                                                                                                                                     <  5:762572/149‑1 (MQ=255)

                                                                                                                                                                                                                       cTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTg                                                                                                                                                                                    >  7:101789/1‑149 (MQ=255)

                                                                                                                                                                                                                        tATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGc                                                                                                                                                                                   >  8:965996/1‑149 (MQ=255)

                                                                                                                                                                                                                          tGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCAt                                                                                                                                                                                 >  6:197014/1‑149 (MQ=255)

                                                                                                                                                                                                                           gAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCTGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTGCGTCAACAACTGTGCAGCCTGCAtg                                                                                                                                                                                >  8:1077711/1‑149 (MQ=255)

                                                                                                                                                                                                                           gAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCAtg                                                                                                                                                                                >  8:763947/1‑149 (MQ=255)

                                                                                                                                                                                                                            aaaCGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCAtgt                                                                                                                                                                               >  6:1318137/1‑149 (MQ=255)

                                                                                                                                                                                                                                ggCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTc                                                                                                                                                                           >  5:1203544/1‑149 (MQ=255)

                                                                                                                                                                                                                                    ggCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCggg                                                                                                                                                                       >  8:937424/1‑149 (MQ=255)

                                                                                                                                                                                                                                    ggCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCggg                                                                                                                                                                       >  6:1200607/1‑149 (MQ=255)

                                                                                                                                                                                                                                    ggCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCggg                                                                                                                                                                       <  7:39641/149‑1 (MQ=255)

                                                                                                                                                                                                                                     gCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCa                                                                                                                                                                                  <  5:170663/137‑1 (MQ=255)

                                                                                                                                                                                                                                      ccAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATga                                                                                                                                                                                                                >  8:213150/1‑106 (MQ=255)

                                                                                                                                                                                                                                       cAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGa                                                                                                                                                                                                                                >  7:99924/1‑89 (MQ=255)

                                                                                                                                                                                                                                       cAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCAt                                                                                                                                                                                 >  6:36781/1‑136 (MQ=255)

                                                                                                                                                                                                                                       cAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCAt                                                                                                                                                                                 <  5:36781/136‑1 (MQ=255)

                                                                                                                                                                                                                                       cAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGa                                                                                                                                                                    >  5:862935/1‑149 (MQ=255)

                                                                                                                                                                                                                                         aaTTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATcc                                                                                                                                                                                                                                     >  7:1078600/1‑82 (MQ=255)

                                                                                                                                                                                                                                         aaTTTTAATAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGAAGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAAc                                                                                                                                                                  <  5:990928/149‑1 (MQ=255)

                                                                                                                                                                                                                                          aTTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAAcc                                                                                                                                                                 >  7:699418/1‑149 (MQ=255)

                                                                                                                                                                                                                                           ttttAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCg                                                                                                                                                                >  5:1040285/1‑149 (MQ=255)

                                                                                                                                                                                                                                             ttAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTa                                                                                                                                                              <  5:1318137/149‑1 (MQ=255)

                                                                                                                                                                                                                                             ttAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTa                                                                                                                                                              <  6:578473/149‑1 (MQ=255)

                                                                                                                                                                                                                                                 tagAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCaa                                                                                                                                                          <  8:1182692/149‑1 (MQ=255)

                                                                                                                                                                                                                                                     gTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTGAATGAGGAACTTTGCCCTCCCCGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCTGCAGCTCCGGCAATTACCACACCATCGACACCCTGATTCAGTa                         <  2:293810/274‑1 (MQ=255)

                                                                                                                                                                                                                                                      ttGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATg                                                                                                                                                                                                            <  5:41362/94‑1 (MQ=255)

                                                                                                                                                                                                                                                      ttGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTc                                                                                                                                                                                                         >  4:45454/1‑97 (MQ=255)

                                                                                                                                                                                                                                                       tGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAAcatcat                                                                                                                                                    <  7:229253/149‑1 (MQ=255)

                                                                                                                                                                                                                                                          acaATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCaa                                                                                                                                                 >  6:818084/1‑149 (MQ=255)

                                                                                                                                                                                                                                                            aaTAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGa                                                                                                                                               >  7:524531/1‑149 (MQ=255)

                                                                                                                                                                                                                                                            aaTAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGa                                                                                                                                               >  7:844682/1‑149 (MQ=255)

                                                                                                                                                                                                                                                              tAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATACCTGACACGGGAAGCGCAAATCACCGGCATAGCTTCTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCTGCAGCTCCGCCAATTCCCACACCATCGACCCCCTGATTCAGTCACAGTGa                  >  4:69696/1‑272 (MQ=255)

                                                                                                                                                                                                                                                               aGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATaa                                                                                                                                            <  7:1108934/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                  ccccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTcaac                                                                                                                                                                                                  >  5:693466/1‑92 (MQ=255)

                                                                                                                                                                                                                                                                  ccccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTcaac                                                                                                                                                                                                  <  6:693466/92‑1 (MQ=255)

                                                                                                                                                                                                                                                                   cccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTcaac                                                                                                                                                                                                  <  5:143116/91‑1 (MQ=255)

                                                                                                                                                                                                                                                                   cccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACatcatc                                                                                                                                                   >  5:127932/1‑138 (MQ=255)

                                                                                                                                                                                                                                                                    ttaaCCCGTGCAGAACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGc                                                                                                                                       >  8:810022/3‑149 (MQ=255)

                                                                                                                                                                                                                                                                      aaCCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTcaac                                                                                                                                                                                                  >  6:168019/1‑88 (MQ=255)

                                                                                                                                                                                                                                                                        cccGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGc                                                                                                                                                                                          <  8:1291942/94‑1 (MQ=255)

                                                                                                                                                                                                                                                                        cccGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGc                                                                                                                                                                                          >  7:1291942/1‑94 (MQ=255)

                                                                                                                                                                                                                                                                            ttCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGCTCCGCTGATGTCCATCGCGACTGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCCACATCATCAAGATAACTTGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCCAAGGTCATCGTTACGTCCTGCAGCTCCGGCAATTACCCCCCCATCCACACCCTGATTCAGTAACAGTGCAAACCGCCTTGCCCgc  >  4:59848/1‑274 (MQ=255)

                                                                                                                                                                                                                                                                            ttCTGCACGGGTGAATGAGGAACTTTGCCCTCCCAGCCAGGCTATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTTTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCTGCAGCTCCGGCAATTACCACACCATCGACACCCTGATTCAGTAACa                      <  2:48487/254‑1 (MQ=255)

CGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCACTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCTGCAGCTCCGGCAATTACCACACCATCGACACCCTGATTCAGTAACAGTGAAAACCGCTGTGCCAGC  >  NC_000913/1698316‑1698857

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