BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,698,589	C→A	5.8%	A198S (GCA→TCA) ‡	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,589	0	C	A	5.8%	788.6 / 23.9	221	A198S (GCA→TCA) ‡	malI	transcriptional repressor of Mal regulon
Reads supporting (aligned to +/- strand): ref base C (100/108); new base A (5/8); total (105/116)
Fisher's exact test for biased strand distribution p-value = 5.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 222 total aligned reads displayed.

CGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCACTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCTGCAGCTCCGGCAATTACCACACCATCGACACCCTGATTCAGTAACAGTGAAAACCGCTGTGCCAGC  >  NC_000913/1698316‑1698857

cGGTCAATGCCGCCCGCCGAACCTGTTGCGCAAAGTAACGATCTGTTCGGCTTTCCCCGCTTTGCCGCCCTGCTTTGGGCAACCCAAACCATCCCCCCATCGCATGAGTTTCTTTATGTCACACCACGTCACTGATGGTCGGTTTATGACTTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCa                                                                                                                                                                                                                                                                             <  4:106087/275‑1 (MQ=255)

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

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

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

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

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

                                                                                            gCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCACTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTGTAGTAGAGTTTCACAATAGACACCAACCCGTCCTGCACGGGTTAATGAGGAACTTTGCCCTCCAAGCCAGGCGAGCCGCTGACGCCCATTGCGAATGAAAGGCCCCGTCAACAACTGTGCAGCCTTCa                                                                                                                                                                                  >  4:29725/1‑274 (MQ=255)

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

                                                                                                                                      aTGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTGCAGAACGGGttgg                                                                                                                                                                                                                                                                     <  7:1419397/149‑3 (MQ=255)

                                                                                                                                       tGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAAt                                                                                                                                                                                                                                                                    <  8:1237352/149‑1 (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                                                                                                                                                                                                                                                              >  7:1182692/1‑149 (MQ=255)

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

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

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

                                                                                                                                                   gACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGc                                                                                                                                                                                                                                                        <  8:1222608/149‑1 (MQ=255)

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

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

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

                                                                                                                                                              gcCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCccagcca                                                                                                                                                                                                                                             >  5:1355832/1‑149 (MQ=255)

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

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

                                                                                                                                                                         ttCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCt                                                                                                                                                                                                                                  >  5:1090893/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)

                                                                                                                                                                                      ttCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACttt                                                                                                                                                                                                                                                          >  5:1305092/1‑112 (MQ=255)

                                                                                                                                                                                      ttCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTGCAGAACGGGttgggg                                                                                                                                                                                                                                                                   >  7:517163/1‑99 (MQ=255)

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

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

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

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

                                                                                                                                                                                                   tGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTcc                                                                                                                                                                                                        <  7:322271/149‑1 (MQ=255)

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

                                                                                                                                                                                                      aCTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTc                                                                                                                                                                                                     >  5:307355/1‑149 (MQ=255)

                                                                                                                                                                                                      aCTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTc                                                                                                                                                                                                     >  8:1313672/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)

                                                                                                                                                                                                          cAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACttt                                                                                                                                                                                                                                                          >  7:1147262/1‑92 (MQ=255)

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

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

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

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

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

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

                                                                                                                                                                                                                ccAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACtgtg                                                                                                                                                                                           <  6:1077158/149‑1 (MQ=255)

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

                                                                                                                                                                                                                  aaTCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCa                                                                                                                                                                                         <  5:508418/149‑1 (MQ=255)

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

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

                                                                                                                                                                                                                          tGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCAt                                                                                                                                                                                 <  7:533703/149‑1 (MQ=255)

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

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

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

                                                                                                                                                                                                                             aaCGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCAtgtt                                                                                                                                                                              <  6:900021/149‑1 (MQ=255)

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

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

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

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

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

                                                                                                                                                                                                                                      ccAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCg                                                                                                                                                                     >  7:698875/1‑149 (MQ=255)

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

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

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

                                                                                                                                                                                                                                        aaaTTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATACGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGaa                                                                                                                                                                   <  8:101789/149‑1 (MQ=255)

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

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

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

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

                                                                                                                                                                                                                                                  agagTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCACATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAAc                                                                                                                                                         <  7:30347/149‑1 (MQ=255)

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

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

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

                                                                                                                                                                                                                                                         cacaATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGt                                                                                                                                                               <  5:1133167/136‑1 (MQ=255)

                                                                                                                                                                                                                                                          acaATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCaa                                                                                                                                                 <  7:1313672/149‑1 (MQ=255)

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

                                                                                                                                                                                                                                                            aaTAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGa                                                                                                                                               >  8:1181063/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)

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

                                                                                                                                                                                                                                                                 tcattaaCCCGTGCAGAACGGGTTAATGAGGAACTTTGCCCTCCTAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACt                                                                                                                                          >  5:423270/6‑149 (MQ=255)

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

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

                                                                                                                                                                                                                                                                   cccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTcaac                                                                                                                                                                                                  <  5:143116/91‑1 (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)

                                                                                                                                                                                                                                                                            ttCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAg                                                                                                                               <  5:849147/149‑1 (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` `< 12 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`