BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,479,020	A→T	6.0%	intergenic (+131/+14)	rraB → / ← yjgM	protein inhibitor of RNase E/GNAT family putative N‑acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,479,020	0	A	T	6.0%	836.4 / 31.0	235	intergenic (+131/+14)	rraB/yjgM	protein inhibitor of RNase E/GNAT family putative N‑acetyltransferase
Reads supporting (aligned to +/- strand): ref base A (111/108); new base T (7/7); total (119/116)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.45e-01

Only 100 of 240 total aligned reads displayed.

GGTTGAACAACTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTCGCAAATATCCGAC  >  NC_000913/4478751‑4479286

                                  gACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGAGGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACtt                                                                                                                                                                                                                                     <  2:689342/275‑1 (MQ=255)

                                    cGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGACGATTTTGTCGATGAAGACGATGACAGAGTTCGCCACTAATTAAGCGTAGGGCATCAGGCCCAGGTCGTTCCTGCAACCTGGTGAGTTTGTAGCCCTGATAAGGACTTAATTTCACATCCGGCATAAACAACACGTACTGTGTTAATGTTGTAATGCCCGCGAGTTCAACCCAACGGAATTTTTTTAGACTTGTCGCaaactccaccattcc                                                                                                                                                                                                                                   >  4:400169/1‑262 (MQ=255)

                                                  gATGGGGCACTTACTTTGAAGATCCTAACGTCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGGATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGc                                                                                                                                                                                                                      <  2:389788/274‑1 (MQ=255)

                                                    tGGGGAACTTACTTTGAAGATCCAAAGGGGGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCgcgc                                                                                                                                                                                                                                                  <  4:637563/244‑1 (MQ=255)

                                                      gggCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGAGGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCa                                                                                                                                                                                                               <  1:287222/277‑1 (MQ=255)

                                                       gtgacttagatTGAAGGTCCGAAAGGCGGATATGGCGAGGATTAAGATTTTGGCGATGAAGACGACGGCGGGTTGTTCAACTAATTAAGCGTAGCGCGTCAGGCCCAGGTAGTTCCTTCAACCTGGGGAATGTGTAGCCCGGATAAGTCATTCATGTCACATCCGTCATAAAGAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCGCGTg                                                                                                                                                                                                                 <  4:29055/271‑1 (MQ=255)

                                                       ggCATTTACGTTGAAGATCCTAACGGCTACTATGTGGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCAACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTTGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTg                                                                                                                                                                                                                 <  4:509387/274‑1 (MQ=255)

                                                       ggCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACaa                                                                                                                                                                                                                                <  2:284341/259‑1 (MQ=255)

                                                                         ccTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCa                                                                                                                                                                                                                                                                    <  3:79406/205‑1 (MQ=255)

                                                                            gcggnngngggtGGGGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTNGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTg                                                                                                                                                                                                                 <  3:191206/242‑1 (MQ=255)

                                                                                 cGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGAAATGCTCa                                                                                                                                                                                      >  2:352151/1‑275 (MQ=255)

                                                                                   aaGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCCCAAAGcca                                                                                                                                                                                >  1:206165/1‑278 (MQ=255)

                                                                                         gCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCa                                                                                                                                                                               >  2:595244/1‑274 (MQ=255)

                                                                                             cgaTGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCa                                                                                                                                                                                      <  2:398345/263‑1 (MQ=255)

                                                                                                                         gACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCACCCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGCTAATGTTCGAATGCCCGCGAGTTGAaca                                                                                                                                                                                                                                                                            >  5:1337682/1‑149 (MQ=255)

                                                                                                                          aCGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAacac                                                                                                                                                                                                                                                                           >  5:689534/1‑149 (MQ=255)

                                                                                                                            ggAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACAcgc                                                                                                                                                                                                                                                                         >  7:1000313/1‑149 (MQ=255)

                                                                                                                              aGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCgg                                                                                                                                                                                                                                                                       <  5:466748/149‑1 (MQ=255)

                                                                                                                              aGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCgg                                                                                                                                                                                                                                                                       <  6:1160228/149‑1 (MQ=255)

                                                                                                                              aGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCgg                                                                                                                                                                                                                                                                       <  6:921466/149‑1 (MQ=255)

                                                                                                                               gTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTCCGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCggg                                                                                                                                                                                                                                                                      <  5:500491/149‑1 (MQ=255)

                                                                                                                                  cGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACa                                                                                                                                                                                                                          >  2:635358/1‑190 (MQ=255)

                                                                                                                                      aCTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAttttt                                                                                                                                                                                                                                                               <  5:977148/149‑1 (MQ=255)

                                                                                                                                         aaTTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTa                                                                                                                                                                                                                                                            <  5:153659/149‑1 (MQ=255)

                                                                                                                                         aaTTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTa                                                                                                                                                                                                                                                            <  7:239873/149‑1 (MQ=255)

                                                                                                                                               gCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCTTGCAGCCAAGCGCATAGTCGATATGCTCAAAACCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCTGTCCTTCCCAGAtcac                                                                                                                         >  4:132662/1‑271 (MQ=255)

                                                                                                                                                 gTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTcgc                                                                                                                                                                                                                                                    <  5:105057/149‑1 (MQ=255)

                                                                                                                                                  tAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTc                                                                                                                                                                                                                                    >  4:5350/1‑164 (MQ=255)

                                                                                                                                                  tAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCCAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATagcagc                                                                                                                      >  2:533510/1‑274 (MQ=255)

                                                                                                                                                    gcgcATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCgg                                                                                                                                                                                                                                                                       >  5:773067/1‑127 (MQ=255)

                                                                                                                                                    gcgcATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGtt                                                                                                                   <  2:109557/275‑1 (MQ=255)

                                                                                                                                                                  ggTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTca                                                                                                                                                                                                                                   <  7:1210977/149‑1 (MQ=255)

                                                                                                                                                                      gTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAAt                                                                                                                                                                                                                               >  5:1131418/1‑149 (MQ=255)

                                                                                                                                                                      gTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTCCACAAt                                                                                                                                                                                                                               >  7:724479/1‑149 (MQ=255)

                                                                                                                                                                      gTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTCCACAAt                                                                                                                                                                                                                               <  8:724479/149‑1 (MQ=255)

                                                                                                                                                                       ttCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATc                                                                                                                                                                                                                              <  5:1072737/149‑1 (MQ=255)

                                                                                                                                                                         ccTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGa                                                                                                                                                                                                                            <  8:1041073/149‑1 (MQ=255)

                                                                                                                                                                           tGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCa                                                                                                                                                                                                                                        >  8:1417326/1‑135 (MQ=255)

                                                                                                                                                                           tGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACa                                                                                                                                                                                                                          <  5:1559242/149‑1 (MQ=255)

                                                                                                                                                                             cAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATg                                                                                                                                                                                                                        >  8:524797/1‑149 (MQ=255)

                                                                                                                                                                                 gtGGGGAAATTGTAACCCGGTTAAAGCATGAATGTGCGATCCGGCAAAAACAAAACGTACTTTGTTAATGTACGCATGCCCGCAAGTTGAACACGCGGGCATTATTGTAGAGTTCGCGCAACATCCGAACTTCACACTCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCAt                                                                                       <  4:280963/273‑1 (MQ=255)

                                                                                                                                                                                  tGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAttttt                                                                                                                                                                                                                                                               <  6:1697159/105‑1 (MQ=255)

                                                                                                                                                                                  tGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAttttt                                                                                                                                                                                                                                                               >  5:1697159/1‑105 (MQ=255)

                                                                                                                                                                                  tGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTc                                                                                                                                                                                                                                    >  8:1320854/1‑132 (MQ=255)

                                                                                                                                                                                  tGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCg                                                                                                                                                                                                                   >  6:1014647/1‑149 (MQ=255)

                                                                                                                                                                                   ggTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGt                                                                                                                                                                                                                  >  5:1440097/1‑149 (MQ=255)

                                                                                                                                                                                    gTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACCCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTTCAGc                                                                                                                                                                                                             >  1:50844/1‑153 (MQ=255)

                                                                                                                                                                                      gAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCa                                                                                                                                                                                                               <  5:563452/149‑1 (MQ=255)

                                                                                                                                                                                         tttGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTa                                                                                                                                                                                                                                                            >  6:27427/1‑101 (MQ=255)

                                                                                                                                                                                                    gATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAgcg                                                                                                                                                                                                       <  7:841077/143‑1 (MQ=255)

                                                                                                                                                                                                    gATAAGGCATTCAGTTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTGCGAATGCACGCGAGTTTAACACGCGGGCATTTTGTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGtttt                                                                    <  2:221681/274‑1 (MQ=255)

                                                                                                                                                                                                     aTAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGc                                                                                                                                                                                                                      >  1:236816/1‑127 (MQ=255)

                                                                                                                                                                                                          gCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCgg                                                                                                                                                                                                                                                                       >  6:440623/1‑73 (MQ=38)

                                                                                                                                                                                                           cATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATg                                                                                                                                                                                          <  8:232355/149‑1 (MQ=255)

                                                                                                                                                                                                            aTTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGc                                                                                                                                                                                         >  8:577739/1‑149 (MQ=255)

                                                                                                                                                                                                                 tGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCaaa                                                                                                                                                                                    >  7:13088/1‑149 (MQ=255)

                                                                                                                                                                                                                  gTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAg                                                                                                                                                                                   <  5:828682/149‑1 (MQ=255)

                                                                                                                                                                                                                  gTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAg                                                                                                                                                                                   >  8:299114/1‑149 (MQ=255)

                                                                                                                                                                                                                   tCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGc                                                                                                                                                                                  <  8:1473834/149‑1 (MQ=255)

                                                                                                                                                                                                                   tCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGc                                                                                                                                                                                  >  5:404558/1‑149 (MQ=255)

                                                                                                                                                                                                                   tCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGc                                                                                                                                                                                  >  5:1380754/1‑149 (MQ=255)

                                                                                                                                                                                                                    cacaTCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGcc                                                                                                                                                                                 <  6:1440097/149‑1 (MQ=255)

                                                                                                                                                                                                                            gCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCACGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATcc                                                                                                                                                                                                                                           >  8:1026365/1‑83 (MQ=255)

                                                                                                                                                                                                                              aTAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCa                                                                                                                                                                       >  7:904671/1‑149 (MQ=255)

                                                                                                                                                                                                                                 aaCAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATcc                                                                                                                                                                                                                                           >  7:605619/1‑78 (MQ=255)

                                                                                                                                                                                                                                  aCAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATaaa                                                                                                                                                                   >  8:1604789/1‑149 (MQ=255)

                                                                                                                                                                                                                                  aCAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGgaa                                     <  2:132664/275‑1 (MQ=255)

                                                                                                                                                                                                                                         gtacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAAt                                                                                                                                                            <  7:1494219/149‑1 (MQ=255)

                                                                                                                                                                                                                                           acgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGACGAAAATAc                             >  2:553605/1‑274 (MQ=255)

                                                                                                                                                                                                                                                 ttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCt                                                                                                                                                    >  7:1507383/1‑149 (MQ=255)

                                                                                                                                                                                                                                                   aaTGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCttt                                                                                                                                                  <  6:1112621/149‑1 (MQ=255)

                                                                                                                                                                                                                                                      gTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGc                                                                                                                                                                                                                      >  7:516624/1‑78 (MQ=255)

                                                                                                                                                                                                                                                      gTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTaaa                                                                                                                                               >  6:353187/1‑149 (MQ=255)

                                                                                                                                                                                                                                                       ttCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTaaaa                                                                                                                                              >  7:1055502/1‑149 (MQ=255)

                                                                                                                                                                                                                                                       ttCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTaaaa                                                                                                                                              >  6:1245981/1‑149 (MQ=255)

                                                                                                                                                                                                                                                          gAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAg                                                                                                                                           <  7:1604789/149‑1 (MQ=255)

                                                                                                                                                                                                                                                           aaTGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGcc                                                                                                                                                                                                            <  2:44103/83‑1 (MQ=255)

                                                                                                                                                                                                                                                            aTGCCCGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCg                                                                                                                                         >  7:1564028/1‑149 (MQ=255)

                                                                                                                                                                                                                                                             tGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTACGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCGCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCAtctctc                                                                                                    <  2:136051/185‑1 (MQ=255)

                                                                                                                                                                                                                                                              gCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGt                                                                                                                                       <  8:470880/149‑1 (MQ=255)

                                                                                                                                                                                                                                                               cccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCAt                                                                                                                                                                                                    >  3:42459/1‑87 (MQ=255)

                                                                                                                                                                                                                                                               cccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTc                                                                                                                                      <  6:513781/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                ccGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCt                                                                                                                                                    >  7:358401/1‑134 (MQ=255)

                                                                                                                                                                                                                                                                ccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCAc                                                                                                                                                                                                                                       >  7:1245943/1‑51 (MQ=255)

                                                                                                                                                                                                                                                                ccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCAt                                                                                                                                                                                                    <  6:1502521/86‑1 (MQ=255)

                                                                                                                                                                                                                                                                ccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCAt                                                                                                                                                                                                    >  5:1502521/1‑86 (MQ=255)

                                                                                                                                                                                                                                                                ccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTc                                                                                                                                      >  7:1358664/1‑148 (MQ=255)

                                                                                                                                                                                                                                                                ccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCg                                                                                                                                     >  5:911814/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                 cgcgAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGATCAAAGCCCAAATGCTCATAAAGCGCAATGGATTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTCTTTTGCCAGCCCTTTGCCGCGGATCGCCGGAAAAAAATACATCATTTGCATTTCGCAAATATCCGac  >  1:174247/1‑279 (MQ=255)

                                                                                                                                                                                                                                                                    gaatttataacGCCGCCATTTTTTTAGGCTCCGTGCAACATCCGAACTTCACATCTGAACTGGCCCGTGCACCCAAGCGCATGGTCGATTGGCTCAAAGCCCATATGCGCATAAAGCGCAAGGGCTTCCTTTAAAAAAGCGGTCGTTTCGAGAGAGCAGCGTTTGAAACCCATCTCTCTCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAGAATACATCTTTTGCAGTTCGCAAATATCCGa   <  4:612300/266‑1 (MQ=255)

                                                                                                                                                                                                                                                                    gAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTc                                                                                                                                 <  8:1539413/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                     aGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTcc                                                                                                                                <  5:1223237/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                      gTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCa                                                                                                                               <  6:1559293/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                      gTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCa                                                                                                                               <  5:238721/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                       ttGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAg                                                                                                                              >  6:678725/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                          aaCACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGata                                                                                                                                                                                            >  7:45846/1‑84 (MQ=255)

                                                                                                                                                                                                                                                                          aaCACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTaaaaaa                                                                                                                                            >  6:810705/1‑132 (MQ=255)

                                                                                                                                                                                                                                                                          aaCACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATa                                                                                                                           >  8:1105383/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                          aaCACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATa                                                                                                                           >  8:804388/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                             acGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATagca                                                                                                                        <  6:1359650/149‑1 (MQ=255)

GGTTGAACAACTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTCGCAAATATCCGAC  >  NC_000913/4478751‑4479286

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