BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	92,096	A→T	A228A (GCA→GCT)	ftsI →	peptidoglycan DD‑transpeptidase FtsI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	92,096	0	A	T	97.2%	339.5 / 0.6	106	A228A (GCA→GCT)	ftsI	peptidoglycan DD‑transpeptidase FtsI
Reads supporting (aligned to +/- strand): ref base A (1/2); new base T (54/49); total (55/51)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.56e-01

Only 100 of 106 total aligned reads displayed.

CATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCAGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGTAACCGTACCATCACCGACGTGTTTGAACCGGGCTCAACGGT  >  NC_000913/91850‑92338

cATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTgcgc                                                                                                                                                                                                                                                >  1:748946/1‑251 (MQ=255)

 aTGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTACGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTAAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTAACGGGCAGCCGGGTGAGAGACGTGTTCGTACAGACCGCTATGGTCGCGTAATTGAAGATATATCTTCTACTGACAGCCAGGCTgcga                                                                                                                                                                                                                                                >  2:541387/1‑249 (MQ=255)

   ggCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGcac                                                                                                                                                                                                                                              <  1:301311/250‑1 (MQ=255)

     cGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAAc                                                                                                                                                                                                                                           >  1:670449/1‑251 (MQ=255)

     cGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAAc                                                                                                                                                                                                                                           >  1:807841/1‑251 (MQ=255)

      ggaatacacaaaaaaaCTGAATCTGCCAGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAAGTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAAcc                                                                                                                                                                                                                                          <  2:663867/241‑1 (MQ=255)

       gACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCt                                                                                                                                                                                                                                         <  1:101321/251‑1 (MQ=255)

        aCTACATCAAAAAACTGAAACTGCCGTGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGACGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAACCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTg                                                                                                                                                                                                                                        >  1:68345/1‑251 (MQ=255)

        aCTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGCGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTg                                                                                                                                                                                                                                        >  2:596526/1‑251 (MQ=255)

            cATCCAAAAACTGAAACAGCCGGGGCTTCATCTGCGCGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAACTGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGGAGATATTTCTTCTACTGACAGCCAGGCTGCGCAAAACCTggagc                                                                                                                                                                                                                                    >  1:283178/1‑248 (MQ=255)

            cATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTCTCCATCCGGCGAAGTGACTGCTAACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGCTTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGgcgc                                                                                                                                                                                                                                    >  1:990877/1‑251 (MQ=255)

              acAAAAAACTAAATCTGCCGGGGATTTTTCTGCGTGAAGATCCTCGCCGTTATTATCCCTCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCTCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAATATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTg                                                                                                                                                                                                                                  <  2:204593/250‑1 (MQ=255)

                  aaaaCTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTa                                                                                                                                                                                                                              <  1:614957/251‑1 (MQ=255)

                     actgaaactgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATtg                                                                                                                                                                                                                           <  1:635847/251‑1 (MQ=255)

                              gCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACg                                                                                                                                                                                                                    <  2:946612/249‑1 (MQ=255)

                              gCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGc                                                                                                                                                                                                                   <  1:336348/250‑1 (MQ=255)

                               ccgtggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGTGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCt                                                                                                                                                                                                                 <  1:751529/247‑1 (MQ=255)

                                cGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTg                                                                                                                                                                                                                >  1:858612/1‑251 (MQ=255)

                                cGGGGATCCAACTTCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCGGTACTAACGTCGATAGGCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGAGAAATGGCTTACCGGTCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTg                                                                                                                                                                                                                <  2:620080/251‑1 (MQ=255)

                                   ggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAg                                                                                                                                                                                                             <  1:205909/251‑1 (MQ=255)

                                       tCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGgcgc                                                                                                                                                                                                         >  1:882048/1‑251 (MQ=255)

                                       tCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGgcgc                                                                                                                                                                                                         >  1:113091/1‑251 (MQ=255)

                                        cATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTCTTGATGAACGCCTGCAGGCGCt                                                                                                                                                                                                        >  2:931001/1‑251 (MQ=255)

                                         ctCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTCATTGAAGATATTTCTTCTACTGCCAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTg                                                                                                                                                                                                       >  2:284438/2‑251 (MQ=255)

                                                      atctcGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGCCAGCCAGGCTGCGCACAACTTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAAc                                                                                                                                                                                          >  2:35885/2‑251 (MQ=255)

                                                           gCCGTTATTATCCATACGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATCAATGGCTTCCCGGGCAGCCGGGTGAGCGCATTGTGCGTAACGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGTAATCCAGGCTGCGCACAACCTGGCGCTGAGTATTGTTGAAAGCTTGCAGGCGCTGGTTTATGGCGCACTgtac                                                                                                                                                                                     >  2:980817/1‑248 (MQ=255)

                                                             cGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAACGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGaaca                                                                                                                                                                                    >  1:456270/1‑250 (MQ=255)

                                                                 attagccatcctgCGAAGGGACTGCTCACCTCCTTGGCTTTACGAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGACTGCAGGCGCTGGTTTATCGCGAACTGAACAAtgcg                                                                                                                                                                               <  2:266019/239‑1 (MQ=255)

                                                                 attagccatccgGCGAAGTGACTGCGCACCTCATTTGCGTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGaacaatgc                                                                                                                                                                                <  2:937058/242‑4 (MQ=255)

                                                                 attATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAAtgcg                                                                                                                                                                               <  1:377510/249‑1 (MQ=255)

                                                                       catccgGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGcc                                                                                                                                                                         >  1:79127/3‑251 (MQ=255)

                                                                        atccgGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCt                                                                                                                                                                        <  2:538368/250‑1 (MQ=255)

                                                                        atccgGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCt                                                                                                                                                                        <  1:931001/250‑1 (MQ=255)

                                                                        atccgGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCt                                                                                                                                                                        <  1:902127/250‑1 (MQ=255)

                                                                          ccgGCGAAGGTACGGCTCTCCTCAGTGGCGTTACTAACGGCGATAGTCAAGGGATTGAGGGCGTGGAGAAGAGTTTCGATAAATGGCGTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGTTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGAGGAACGCCTTCGGGCTCCGGTTTATCGCGAACTGAACAATGCGGTGGCCtt                                                                                                                                                                       <  2:44405/250‑1 (MQ=255)

                                                                           cgGCGAAGTGTCTGCTCGCCTCATTGGCTTTACTAACGTCGAGAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCGATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCTCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGACCAATGCGGTGGCCTTTa                                                                                                                                                                     <  2:716713/251‑1 (MQ=255)

                                                                                    gACTGCTCACCGCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCt                                                                                                                                                             <  1:374992/250‑1 (MQ=255)

                                                                                     aCTGCTCACCTCATTGGCTTTACTAACGTCGAAAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGa                                                                                                                                                           <  1:304631/251‑1 (MQ=255)

                                                                                           cACCTCATTGGCTTTACTAACGTCGATAATCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCt                                                                                                                                                       <  2:453786/249‑1 (MQ=255)

                                                                                           cACCTCAGTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAATAGTTTCGATAAATGGCATACCGTGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCCGTGGCCTTTAACAAGGCTGAAGCTg                                                                                                                                                      <  1:315509/250‑1 (MQ=255)

                                                                                            aCCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGCTGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTACCAAGGCTGAATCTGGt                                                                                                                                                    >  2:179605/1‑251 (MQ=255)

                                                                                            aCCTCATCGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCAGCGCACATCCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGAGGTGGCCTTTAACAAGGCTGAATCTGGt                                                                                                                                                    >  2:458298/1‑251 (MQ=255)

                                                                                                cattggCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAgcg                                                                                                                                                <  1:817059/247‑1 (MQ=255)

                                                                                                cattggCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAgcg                                                                                                                                                <  1:756702/247‑1 (MQ=255)

                                                                                                   tggCTTTACTAACGTCGATAGTCAAGGGATTTAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGcc                                                                                                                                              <  2:812363/249‑1 (MQ=255)

                                                                                                           cTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCtggtg                                                                                                                                     >  2:685207/1‑251 (MQ=255)

                                                                                                             aaCGTCGATAGTCAAGTGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCtggtgg                                                                                                                                    <  2:192247/250‑1 (MQ=255)

                                                                                                               cgtcgtTAGTAAAGGGACTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCtggtgggtg                                                                                                                                 >  2:927077/1‑248 (MQ=255)

                                                                                                                gTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGt                                                                                                                                <  2:271849/251‑1 (MQ=255)

                                                                                                                  ggagaGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCa                                                                                                                              <  2:506924/247‑1 (MQ=255)

                                                                                                                  cGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCa                                                                                                                              <  1:625342/251‑1 (MQ=255)

                                                                                                                       gTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCAGCGCACATCCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACAcc                                                                                                                         <  1:458298/251‑1 (MQ=255)

                                                                                                                               aTTGAGGACGTTGCGAAGAGTTTCGATAAATGGCTTACCTGGCAGCCGGGGGAGCGCCTTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGCTATTTATTCTACTGCCAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAGCGCCTGCAAGCGCTGGTTTCTCGCGTACTGTACCATGCGCGGGCCGTTAACAAAGATGAATCTGGTAGCGCCGTGCTGTTGAATGTCacaaccgcagaagt                                                                                                                 >  2:632540/1‑238 (MQ=255)

                                                                                                                                   ggggCGTTGAGAAGAGTTTCGAGAAATGGCTTACCGGGCAGCCGGGGGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTg                                                                                                             <  2:958820/250‑1 (MQ=255)

                                                                                                                                         ttGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATg                                                                                                       >  1:413851/1‑251 (MQ=255)

                                                                                                                                           gagaAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGc                                                                                                     <  1:851716/251‑1 (MQ=255)

                                                                                                                                              aagaGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTaa                                                                                                  <  1:422305/251‑1 (MQ=255)

                                                                                                                                                  gttgcGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGc                                                                                              <  2:113091/247‑1 (MQ=255)

                                                                                                                                                      cGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCCGGCGATGGCTAACAGCCCGt                                                                                          >  1:117507/1‑251 (MQ=255)

                                                                                                                                                          aaaTGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGACGTGCTGGCGATGGCTAACAGCCCGTCAt                                                                                       >  2:926521/1‑250 (MQ=255)

                                                                                                                                                              ggCTTTCCGGGCAGCCGGGTGAGAGCATTTTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTTGCCTGTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCTCGTCATACaa                                                                                   <  2:616999/250‑1 (MQ=255)

                                                                                                                                                               gCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGCTGTCCACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATAAAAcc                                                                                 >  2:938647/1‑251 (MQ=255)

                                                                                                                                                                  tACCGGGCAGCTGGGTAAGCGCATTGTGCGTAAAGTCCGCTATGGTCTCGTAATTGAAGATATTTCTTCTCCTGACAGCCAGGCTGCGCCCATCCTGGCGCTGAGTATTGATGAACGCCTGCCGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGGAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCt                                                                               <  2:569048/250‑1 (MQ=255)

                                                                                                                                                                     cGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATCAAACCCTAAc                                                                            >  2:896938/1‑250 (MQ=255)

                                                                                                                                                                          aGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTg                                                                      >  2:935593/1‑251 (MQ=255)

                                                                                                                                                                           gCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATTTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGa                                                                     <  2:855051/251‑1 (MQ=255)

                                                                                                                                                                             cGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGATGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTACCAATCTGAGc                                                                   >  2:214742/1‑251 (MQ=255)

                                                                                                                                                                              gggTGCGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCg                                                                  >  2:222259/1‑251 (MQ=255)

                                                                                                                                                                               ggTGAGCGCATTGGGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGACTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGAGGCCTTTAACAAGGCTGAATCTGGGAGAGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACCACCCTAACAATCAGAGCgg                                                                 >  2:39565/1‑251 (MQ=255)

                                                                                                                                                                                 tGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGc                                                                >  1:792553/1‑250 (MQ=255)

                                                                                                                                                                                   aGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACg                                                             <  1:364187/251‑1 (MQ=255)

                                                                                                                                                                                       cATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGa                                                         <  2:801657/251‑1 (MQ=255)

                                                                                                                                                                                         ttGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGCACGCCTGCAGGCGCTGGTTTCTCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGaa                                                        >  2:292579/1‑250 (MQ=255)

                                                                                                                                                                                              cGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGc                                                  <  1:318502/251‑1 (MQ=255)

                                                                                                                                                                                                 aaaGACCGCTATGGTCGCGAAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACACCTGCAGGCGCTGGTTTATCGCGAACTGACCAATGCGGTGGCCTTTAACAAGGCTGATTCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATAAAACCCTAACCATCTGAGCGGACCGCCGACAGAGGCggt                                               >  2:432160/1‑249 (MQ=255)

                                                                                                                                                                                                  aaGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGGTTATCGCGAACTGAACAATGCGGGGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCTATg                                              <  2:677010/251‑1 (MQ=255)

                                                                                                                                                                                                     aCCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGCATCTGGTAGCGCCGTGCTGGTGGAAGTCAACAACGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCAGAAAGAGGCGacacgt                                           >  2:873173/1‑246 (MQ=255)

                                                                                                                                                                                                     aCCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGt                                           <  2:152553/251‑1 (MQ=255)

                                                                                                                                                                                                      ccGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGt                                           >  1:862341/1‑250 (MQ=255)

                                                                                                                                                                                                       cGCTATGGTCGCGTAATTGAAGCTATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATCCCACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGTa                                          >  2:39966/1‑250 (MQ=255)

                                                                                                                                                                                                       cGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGTa                                          >  1:4068/1‑250 (MQ=255)

                                                                                                                                                                                                       cGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGGGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAAAGGCGATGCGTa                                          >  2:148370/1‑250 (MQ=255)

                                                                                                                                                                                                           aTGGTCGCGTAATTGAAGATATTTCTTCTACTTACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTAAATCTGGTAGCGCCGTGCTGGTGGCTGTCAACACCGGTGAAGTGCTGGCGACGGCTAACGGCCCGTACTAACACCCTAACAATATGAGCGGCACGCCGAAAGAGGCGATGCGTAAccat                                     >  2:749508/1‑249 (MQ=255)

                                                                                                                                                                                                           aTGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGTAACCg                                      <  1:677090/250‑1 (MQ=255)

                                                                                                                                                                                                           aTGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGTAACCGt                                     <  2:264007/251‑1 (MQ=255)

                                                                                                                                                                                                             ggTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCTCAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACTACCCTAACAATCTGATCGGCACGCCGAAAGAGGCGATGCGTAACCGTa                                    <  1:927077/250‑1 (MQ=255)

                                                                                                                                                                                                             ggTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGTAACCGTAc                                   <  1:187199/251‑1 (MQ=255)

                                                                                                                                                                                                                 gcgTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTTATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGTAACCGTACCATc                               >  2:389777/1‑251 (MQ=255)

                                                                                                                                                                                                                      aTTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACGACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTAACGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCACCACCGGTGAAGTGCTGGCGATGGATAAAAGCCCGTCATACAACCCTAACAATCTACGCGGCACGCCGAAAGAGACGATGCGTAACCGTACCATccacgc                          >  2:106616/1‑246 (MQ=255)

                                                                                                                                                                                                                          aaGATCTTTCTTCTCATGACTGCCAGGCTGCGCACCACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTACCAAGGCTGAATCTGGTAGCGCCTTGCTGGTTGATGTCAACACCAGTGCAGTGCTGGAGATGGCTACCAGCCCGTCATGCAACCCTAACAATCTGAGCGGAACGCCGGAAGAGGCGATGCGTAACCGCAcccaccccgccctt                      >  2:201204/1‑239 (MQ=255)

                                                                                                                                                                                                                          aaGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGTAACCGTACCATCACCGCCgtg                      >  1:467906/1‑251 (MQ=255)

                                                                                                                                                                                                                            gATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGTAACCGTACCATCACCGACGGGtt                    >  1:263789/1‑251 (MQ=255)

                                                                                                                                                                                                                              tatTTCTTCTACTGCCAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGCTGCGTAACCGTACCATCACCGACGTGTTTg                  >  2:948554/1‑251 (MQ=255)

                                                                                                                                                                                                                               atTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGCAAGAGGCGATGCGTAACCGTCCCATCACCTCCGTGTTTg                  >  1:108801/1‑250 (MQ=255)

                                                                                                                                                                                                                                tttCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGTAACCGTACCATCACCGACGTGTTTGa                 >  2:797868/1‑250 (MQ=255)

                                                                                                                                                                                                                                      ctACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGTAACCGTACCATCACCGACGTGTTTGAACCTGGc          >  1:990655/1‑251 (MQ=255)

                                                                                                                                                                                                                                      ctACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGTAACCGTACCATCACCGACGTGTTTGAACCTGGc          >  1:429801/1‑251 (MQ=255)

                                                                                                                                                                                                                                        aCTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGTAACCGTACCATCACCGACGTGTTTGAACCTGGc          <  1:613009/249‑1 (MQ=255)

                                                                                                                                                                                                                                             cAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGTAACCGTACCATCACCGACGTGTTTGAACCTGGCTCAACgg   >  1:984725/1‑251 (MQ=255)

                                                                                                                                                                                                                                              aGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGTAACCGTACCATCACCGACGTGTTTGAACCTGGCTCAACGGt  <  1:130337/251‑1 (MQ=255)

CATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCAGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGTAACCGTACCATCACCGACGTGTTTGAACCGGGCTCAACGGT  >  NC_000913/91850‑92338

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