BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,127,251	C→T	6.9%	intergenic (+26/+35)	rpmE → / ← yiiX	50S ribosomal subunit protein L31/putative lipid binding hydrolase, DUF830 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,127,251	0	C	T	6.9%	485.7 / 16.6	157	intergenic (+26/+35)	rpmE/yiiX	50S ribosomal subunit protein L31/putative lipid binding hydrolase, DUF830 family protein
Reads supporting (aligned to +/- strand): ref base C (72/73); new base T (5/6); total (77/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 164 total aligned reads displayed.

TTAAACCGAGGTTTTCCCATGAAAAAAGATATTCACCCGAAATACGAAGAAATTACTGCTAGCTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCG  >  NC_000913/4126995‑4127509

ttAAAACGAGGTTTTCCCATGAAAAAAGATATTCACCCGAAATACGAAGAAATTACTGCTAGCTGCTCTTGCGGTAACGTAAAGAAAATCCGTTCCACCGTTGGTCATGACCTGAACCTCGACGTGGGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGc                                                                                                                                                                                                                                                  <  1:15700/275‑1 (MQ=255)

            ttttCCATGACAAAAGATATTCACCCGAAATACGAAGAAATTAGTGCTCGCTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTgg                                                                                                                                                                                                                                       <  2:329468/274‑1 (MQ=255)

            tttCCCATTAAAAAAGATATTCACCAGAAATACGAAGAAATTACTGCTAGCTGCGCTTGCGGTAACGTAATGAAAATCCGCTCAACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTgg                                                                                                                                                                                                                                       <  2:141149/274‑1 (MQ=255)

                                            cGAAGAAATTACTGCTAGCTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGTGGCCATTCTTTg                                                                                                                                                                                                      >  2:147906/1‑275 (MQ=255)

                                                       cTGCTAGCTGCTCTTGGGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCgtggt                                                                                                                                                                                           <  4:190880/275‑1 (MQ=255)

                                                                           aaCGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATGACCACGa                                                                                                                                                                                                                         >  2:236060/1‑225 (MQ=255)

                                                                               tAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGGGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAAc                                                                                                                                                                                                                     <  4:236060/225‑1 (MQ=255)

                                                                                                  ctttGGTCATGACCAGAACCTCTCCTTGTCCAGCAAGGGCCGCCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGtt                                                                                                                                                <  4:147906/273‑1 (MQ=255)

                                                                                                            gACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGcgccgc                                                                                                                                                                                                                                                                    >  7:40467/1‑149 (MQ=255)

                                                                                                               cTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACg                                                                                                                                                                                                                                                                 >  5:118400/1‑149 (MQ=255)

                                                                                                                tGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACgg                                                                                                                                                                                                                                                                <  5:524802/149‑1 (MQ=255)

                                                                                                                     cTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCtt                                                                                                                                                                                                                                                           >  8:1322086/1‑149 (MQ=255)

                                                                                                                     cTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCtt                                                                                                                                                                                                                                                           >  6:1431133/1‑149 (MQ=255)

                                                                                                                     cTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCtt                                                                                                                                                                                                                                                           >  7:818452/1‑149 (MQ=255)

                                                                                                                       cGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGTGCGGCGCtttt                                                                                                                                                                                                                                                         <  7:803456/149‑1 (MQ=255)

                                                                                                                       cGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGAc                                                                                                                                                                                                 >  2:268853/1‑205 (MQ=255)

                                                                                                                        gACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTTAACAAGCATTTGAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCttttt                                                                                                                                                                                                                                                        >  6:1054329/1‑149 (MQ=255)

                                                                                                                         aCGTGTGCAGCAAGTGCCACCCGGTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCtttttt                                                                                                                                                                                                                                                       <  8:1213419/149‑1 (MQ=255)

                                                                                                                             gtgCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGCGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTCtgtg                                                                                                                                                                                                                                                   <  6:395956/149‑1 (MQ=255)

                                                                                                                              tgCAGCATGTTCGACCCGGTTTTCACTCACAAACAGCGTGATGCTACTAACCTTTCCCGTGTTCACATCTTCAACAACCGTTTAAACATCCCGGGCAGCAAATAAGGCCGGACGTAAGAAAAAAGCGCCGCACTGCGATTGTTTTTTGCCTGAATTTTGGTTGTATTACCCCGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGACTTGTGGCCCACCGAACACGGCTTGTGGCGGGACCACCGTGCACTCAACCGGAATAtttctaccggaaa                                                                                                                    >  4:155237/1‑264 (MQ=255)

                                                                                                                              tgCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGc                                                                                                                                                                                                                                                  <  6:790086/149‑1 (MQ=255)

                                                                                                                                 agcaAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTg                                                                                                                                                                                                                                               <  5:1250943/149‑1 (MQ=255)

                                                                                                                                    aaGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAAt                                                                                                                                                                                                                                            >  7:391210/1‑149 (MQ=255)

                                                                                                                                    aaGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAAt                                                                                                                                                                                                                                            <  8:818452/149‑1 (MQ=255)

                                                                                                                                        gCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTg                                                                                                                                                                                                                                        <  5:1431133/149‑1 (MQ=255)

                                                                                                                                          cACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGt                                                                                                                                                                                                                                      <  7:1151058/149‑1 (MQ=255)

                                                                                                                                            cccGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCtttt                                                                                                                                                                                                                                                         <  5:1247518/128‑1 (MQ=255)

                                                                                                                                            cccGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCtttt                                                                                                                                                                                                                                                         >  6:1247518/1‑128 (MQ=255)

                                                                                                                                             gcGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACATGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGt                                                                                                                                                                                                                                   >  6:850946/2‑149 (MQ=255)

                                                                                                                                               gTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTAt                                                                                                                                                                                                                                 <  8:1122814/149‑1 (MQ=255)

                                                                                                                                                 tcttcACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGc                                                                                                                                                                                                    >  2:213713/1‑176 (MQ=255)

                                                                                                                                                     cACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCAc                                                                                                                                                                                                                           <  6:1036624/149‑1 (MQ=255)

                                                                                                                                                        tGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGaa                                                                                                                                                                                                                        >  8:1438203/1‑149 (MQ=255)

                                                                                                                                                          gCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGaaaa                                                                                                                                                                                                                      >  8:1018046/1‑149 (MQ=255)

                                                                                                                                                          gCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGc                                                                                                                                                                                                    >  1:23670/1‑167 (MQ=255)

                                                                                                                                                            aaaCAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACa                                                                                                                                                                                                                    >  7:990599/1‑149 (MQ=255)

                                                                                                                                                            aaaCAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACa                                                                                                                                                                                                                    <  8:1312668/149‑1 (MQ=255)

                                                                                                                                                              aCAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGAcc                                                                                                                                                                                                <  3:23670/167‑1 (MQ=255)

                                                                                                                                                                     gATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCAtt                                                                                                                                                                                                           <  8:447543/149‑1 (MQ=255)

                                                                                                                                                                        gTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCtt                                                                                                                                                                                                        >  7:770186/1‑149 (MQ=255)

                                                                                                                                                                        gTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCtt                                                                                                                                                                                                        <  8:391210/149‑1 (MQ=255)

                                                                                                                                                                             tACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCg                                                                                                                                                                                                   >  7:131962/1‑149 (MQ=255)

                                                                                                                                                                             tACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCg                                                                                                                                                                                                   >  5:1444103/1‑149 (MQ=255)

                                                                                                                                                                               ccGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGAc                                                                                                                                                                                                 <  5:1007689/149‑1 (MQ=255)

                                                                                                                                                                                  gTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGt                                                                                                                                                                                                                                      >  6:1381251/1‑109 (MQ=255)

                                                                                                                                                                                  gTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGt                                                                                                                                                                                                                                      <  5:1381251/109‑1 (MQ=255)

                                                                                                                                                                                   tGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCgtg                                                                                                                                                                                             <  5:23590/149‑1 (MQ=255)

                                                                                                                                                                                       cGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGa                                                                                                                                                                                         <  7:611225/149‑1 (MQ=255)

                                                                                                                                                                                         tgtTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTc                                                                                                                                                                                                                                     <  5:323875/103‑1 (MQ=255)

                                                                                                                                                                                         tgtTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTc                                                                                                                                                                                                                                     >  6:323875/1‑103 (MQ=255)

                                                                                                                                                                                         tgtTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGACTTGTNGCGCCTCGACAACGGCTTGTTGCGAGAACACCGTTTCCTCAATCGGACNNTTCt                                                                                                                             >  4:351036/1‑207 (MQ=255)

                                                                                                                                                                                         tgtTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCGGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGt                                                                                                                                                                                       >  5:1102886/1‑149 (MQ=255)

                                                                                                                                                                                            tGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAATCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGACGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCtgctgc                                                      <  2:28137/275‑1 (MQ=255)

                                                                                                                                                                                            tGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCtgctgc                                                      <  2:49903/275‑1 (MQ=255)

                                                                                                                                                                                               ccGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGgc                                                                                                                                                                                 <  6:50559/149‑1 (MQ=255)

                                                                                                                                                                                                 gCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGgcgc                                                                                                                                                                               >  7:1901/1‑149 (MQ=255)

                                                                                                                                                                                                   ttCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTtgt                                                                                                                                                                                    >  8:1342222/1‑142 (MQ=255)

                                                                                                                                                                                                   ttCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCAt                                                                                                                                                                             <  8:476690/149‑1 (MQ=255)

                                                                                                                                                                                                   ttCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCAt                                                                                                                                                                             >  5:502261/1‑149 (MQ=255)

                                                                                                                                                                                                   ttCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTcc                                                                                                                                             >  4:283528/1‑181 (MQ=255)

                                                                                                                                                                                                    tCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATc                                                                                                                                                                            <  8:47900/149‑1 (MQ=255)

                                                                                                                                                                                                      aacaaGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGa                                                                                                                                                                          >  8:392825/1‑149 (MQ=255)

                                                                                                                                                                                                       acaaGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGaa                                                                                                                                                                         >  7:258324/1‑149 (MQ=255)

                                                                                                                                                                                                       acaaGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCaa                                                                                                                                         <  2:283528/181‑1 (MQ=255)

                                                                                                                                                                                                          aGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGa                                                                                                                                                                                         <  8:921939/130‑1 (MQ=255)

                                                                                                                                                                                                          aGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGAATTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAAc                                                                                                                                                                      <  7:974261/149‑1 (MQ=255)

                                                                                                                                                                                                              tttCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCt                                                                                                                                                                  >  7:881289/1‑149 (MQ=255)

                                                                                                                                                                                                                   aCATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTTGCGCATCGAAAACGGCTTGTgg                                                                                                                                                             >  5:1271246/1‑149 (MQ=255)

                                                                                                                                                                                                                    cagcccGGGCAGCAAATAAGGCCGGATTTCCGAAAAAAGCGCCGGACGGCGCTTTTTTTTTGCCTGAAGTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGcgcg                              <  4:280592/272‑1 (MQ=255)

                                                                                                                                                                                                                    cATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGc                                                                                                                                                            <  6:502261/149‑1 (MQ=255)

                                                                                                                                                                                                                    cATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGc                                                                                                                                                            >  8:633032/1‑149 (MQ=255)

                                                                                                                                                                                                                      tCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCga                                                                                                                                                          <  5:850946/149‑1 (MQ=255)

                                                                                                                                                                                                                      tCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCga                                                                                                                                                          <  8:68351/149‑1 (MQ=255)

                                                                                                                                                                                                                       cccGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCgag                                                                                                                                                         <  5:594158/149‑1 (MQ=255)

                                                                                                                                                                                                                       cccGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCgag                                                                                                                                                         >  5:1401490/1‑149 (MQ=255)

                                                                                                                                                                                                                         cGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGac                                                                                                                                                       <  5:505337/149‑1 (MQ=255)

                                                                                                                                                                                                                                 aaaTAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGttt                                                                                                                                               >  6:1020602/1‑149 (MQ=255)

                                                                                                                                                                                                                                     aaGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTc                                                                                                                                           >  8:224899/1‑149 (MQ=255)

                                                                                                                                                                                                                                     aaGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTc                                                                                                                                           >  6:685124/1‑149 (MQ=255)

                                                                                                                                                                                                                                       gACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCgag                                                                                                                                                         <  6:20085/133‑1 (MQ=255)

                                                                                                                                                                                                                                        aCCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAg                                                                                                                                        >  7:1353260/1‑149 (MQ=255)

                                                                                                                                                                                                                                          cGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCg                                                                                                                                      >  7:47293/1‑149 (MQ=255)

                                                                                                                                                                                                                                               ttCCGCACAAAGCGCCGCACTGCGCTTTTTTTGTGCCTGAATTTTTGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACtt     <  4:224980/273‑1 (MQ=255)

                                                                                                                                                                                                                                                gcccaaaaaaGCGTCGTACTGCGCTTTTTTTGGTCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCCTTCTTTGGCGACCGTGGTGAGTTGTGGCGGAGCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACCCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCg  <  2:81581/271‑1 (MQ=255)

                                                                                                                                                                                                                                                 acaaaaaaaGCGCCGTGCGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATAtt                                                                                                                               >  7:201551/4‑149 (MQ=255)

                                                                                                                                                                                                                                                   gAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCt                                                                                                                             >  6:1215868/1‑149 (MQ=255)

                                                                                                                                                                                                                                                   gAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCt                                                                                                                             <  6:52404/149‑1 (MQ=255)

                                                                                                                                                                                                                                                    aaaaaaGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCtt                                                                                                                            <  5:670779/149‑1 (MQ=255)

                                                                                                                                                                                                                                                    aaaaaaGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCtt                                                                                                                            <  5:1215868/149‑1 (MQ=255)

                                                                                                                                                                                                                                                    aaaaaaGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCtt                                                                                                                            >  5:77952/1‑149 (MQ=255)

                                                                                                                                                                                                                                                    aaaaaaGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCtt                                                                                                                            <  8:1901/149‑1 (MQ=255)

                                                                                                                                                                                                                                                    aaaaaaGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCtt                                                                                                                            <  6:447845/149‑1 (MQ=255)

                                                                                                                                                                                                                                                      aaaaGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGc                                                                                                                          <  7:1388335/149‑1 (MQ=255)

                                                                                                                                                                                                                                                       aaaGCGCCGTGCGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGcc                                                                                                                         >  7:1033929/1‑149 (MQ=255)

                                                                                                                                                                                                                                                       aaaGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGcc                                                                                                                         <  8:990599/149‑1 (MQ=255)

                                                                                                                                                                                                                                                        aaGCGCCGTGCGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCATTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCg                                                                                                                        >  8:571243/1‑149 (MQ=255)

                                                                                                                                                                                                                                                        aaGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTgg                                                                                                                                                                                                     >  8:596523/1‑72 (MQ=255)

                                                                                                                                                                                                                                                        aaGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTgg                                                                                                                                                                                                     <  7:596523/72‑1 (MQ=255)

                                                                                                                                                                                                                                                           cgccgcACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCGTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTaa                                                                                                                     <  8:1251443/149‑1 (MQ=255)

                                                                                                                                                                                                                                                            gccgcACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCGATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAAc                                                                                                                    >  8:1019623/1‑149 (MQ=255)

TTAAACCGAGGTTTTCCCATGAAAAAAGATATTCACCCGAAATACGAAGAAATTACTGCTAGCTGCTCTTGCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCG  >  NC_000913/4126995‑4127509

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