BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,984,478	G→T	47.8%	P208Q (CCG→CAG)	yecA ←	UPF0149 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,984,478	0	G	T	47.8%	151.5 / inf	253	P208Q (CCG→CAG)	yecA	UPF0149 family protein
Reads supporting (aligned to +/- strand): ref base G (64/68); new base T (79/42); total (143/110)
Fisher's exact test for biased strand distribution p-value = 7.84e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 257 total aligned reads displayed.

TGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTCTGGCGACAT  >  CP009273/1984330‑1984626

 gCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACggg                                                                                                                                                    <  1:1054572/150‑1 (MQ=255)

     tttttGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAg                                                                                                                                                >  2:3172884/1‑150 (MQ=255)

        ttGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGAt                                                                                                                                             <  2:4047808/150‑1 (MQ=255)

                ggCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGc                                                                                                                                     >  2:442585/1‑150 (MQ=255)

                 gCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGcc                                                                                                                                    <  1:315453/150‑1 (MQ=255)

                   aGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGcccc                                                                                                                                  <  1:980171/150‑1 (MQ=255)

                    gggTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCg                                                                                                                                 <  2:3572192/150‑1 (MQ=255)

                    gggTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCg                                                                                                                                 <  2:2961331/150‑1 (MQ=255)

                     ggTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCgg                                                                                                                                <  1:1754234/150‑1 (MQ=255)

                      gttgttACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCACGGATCGTTACGCCCCGGt                                                                                                                               >  2:178517/1‑150 (MQ=255)

                         gttACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTtc                                                                                                                            <  2:1343711/150‑1 (MQ=255)

                           tACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTtctc                                                                                                                          >  1:4173592/1‑150 (MQ=255)

                                   gTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCtt                                                                                                                  >  1:1121695/1‑150 (MQ=255)

                                    tGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCttt                                                                                                                 >  2:539177/1‑150 (MQ=255)

                                     ggggTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTg                                                                                                                >  1:1256767/1‑150 (MQ=255)

                                     ggggTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTg                                                                                                                <  1:1968567/150‑1 (MQ=255)

                                         tGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCg                                                                                                            <  1:2294096/150‑1 (MQ=255)

                                          gATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCgg                                                                                                           >  1:1648999/1‑150 (MQ=255)

                                          gATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCgg                                                                                                           >  2:2075544/1‑150 (MQ=255)

                                            tCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGtt                                                                                                         <  2:2781212/150‑1 (MQ=255)

                                               gATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCt                                                                                                      >  1:3635397/1‑150 (MQ=255)

                                                aTAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTg                                                                                                     >  2:4071509/1‑150 (MQ=255)

                                                 tAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACAACCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGa                                                                                                    <  1:716385/150‑1 (MQ=255)

                                                    ccTCAAAGTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTATTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACa                                                                                                 <  1:4087773/150‑1 (MQ=255)

                                                         aaTTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCttt                                                                                            >  2:2602534/1‑150 (MQ=255)

                                                         aaTTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCttt                                                                                            >  1:250770/1‑150 (MQ=255)

                                                         aaTTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCttt                                                                                            <  2:1005130/150‑1 (MQ=255)

                                                           ttCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCttttt                                                                                          <  2:2996722/150‑1 (MQ=255)

                                                            tCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTc                                                                                         >  2:3886267/1‑150 (MQ=255)

                                                            tCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTc                                                                                         >  2:1607743/1‑150 (MQ=255)

                                                              cTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCt                                                                                       >  2:2227320/1‑150 (MQ=255)

                                                              cTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTGTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCt                                                                                       >  1:2968075/1‑150 (MQ=255)

                                                               ttATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTg                                                                                      >  1:483067/1‑150 (MQ=255)

                                                               ttATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTg                                                                                      <  1:265867/150‑1 (MQ=255)

                                                                tATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGc                                                                                     <  1:4139990/150‑1 (MQ=255)

                                                                tATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGc                                                                                     >  2:64372/1‑150 (MQ=255)

                                                                 aTTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCg                                                                                    >  1:3559897/1‑150 (MQ=255)

                                                                  ttGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCgg                                                                                   <  1:4041897/150‑1 (MQ=255)

                                                                   tGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGa                                                                                  >  2:817711/1‑150 (MQ=255)

                                                                    gggTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGAt                                                                                 <  1:532012/150‑1 (MQ=255)

                                                                     ggTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATg                                                                                <  2:2076399/150‑1 (MQ=255)

                                                                      gTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATgc                                                                               >  1:2862470/1‑150 (MQ=255)

                                                                         ccAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGcc                                                                            <  1:3509679/150‑1 (MQ=255)

                                                                           aGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCAt                                                                          <  1:3212427/150‑1 (MQ=255)

                                                                            gAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATc                                                                         <  1:2643786/150‑1 (MQ=255)

                                                                             aaTTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATcc                                                                        <  2:902554/150‑1 (MQ=255)

                                                                              aTTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCa                                                                       >  2:1786894/1‑150 (MQ=255)

                                                                                tAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGt                                                                     >  2:559279/1‑150 (MQ=255)

                                                                                 aaCGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTa                                                                    >  2:2893594/1‑150 (MQ=255)

                                                                                 aaCGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTa                                                                    <  2:1089108/150‑1 (MQ=255)

                                                                                  aCGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTAc                                                                   >  1:578785/1‑150 (MQ=255)

                                                                                    gCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGc                                                                 >  2:2341412/1‑150 (MQ=255)

                                                                                       gACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATg                                                              >  1:3209267/1‑150 (MQ=255)

                                                                                          acCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAg                                                           >  2:2648565/1‑150 (MQ=255)

                                                                                             cAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTCGATc                                                        >  1:3586795/1‑150 (MQ=255)

                                                                                             cAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATc                                                        >  1:1303893/1‑150 (MQ=255)

                                                                                             cAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATc                                                        <  1:3886707/150‑1 (MQ=255)

                                                                                             cAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATc                                                        >  1:1359604/1‑150 (MQ=255)

                                                                                              aaTTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCa                                                       >  2:751880/1‑150 (MQ=255)

                                                                                               aTTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCaa                                                      >  1:3238941/1‑150 (MQ=255)

                                                                                                tttGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAg                                                     >  2:598762/1‑150 (MQ=255)

                                                                                                tttGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAg                                                     <  1:950917/150‑1 (MQ=255)

                                                                                                 ttGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAgc                                                    >  2:1278134/1‑150 (MQ=255)

                                                                                                 ttGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAgc                                                    <  2:2549076/150‑1 (MQ=255)

                                                                                                  tGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAgcg                                                   >  2:890816/1‑150 (MQ=255)

                                                                                                    gCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGcgcc                                                 >  2:2471510/1‑150 (MQ=255)

                                                                                                        ctTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCa                                             >  1:3234971/1‑150 (MQ=255)

                                                                                                        ctTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCa                                             <  2:3389743/150‑1 (MQ=255)

                                                                                                         ttAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAg                                            <  2:1400133/150‑1 (MQ=255)

                                                                                                             tGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGa                                        >  1:225923/1‑150 (MQ=255)

                                                                                                             tGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGa                                        >  1:2899359/1‑150 (MQ=255)

                                                                                                              gcaggcagCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGaa                                       >  1:1378278/1‑150 (MQ=255)

                                                                                                               caggcagCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAAt                                      >  1:2577236/1‑150 (MQ=255)

                                                                                                                aggcagCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATg                                     <  1:2735044/150‑1 (MQ=255)

                                                                                                                 ggcagCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATgg                                    <  2:130791/150‑1 (MQ=255)

                                                                                                                    agcaCTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCAt                                 <  2:782048/150‑1 (MQ=255)

                                                                                                                    agcaCTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCAt                                 <  2:3196322/150‑1 (MQ=255)

                                                                                                                      caCTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCt                               <  1:2427035/150‑1 (MQ=255)

                                                                                                                       aCTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTa                              <  1:756613/150‑1 (MQ=255)

                                                                                                                           cTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGct                          <  1:206139/150‑1 (MQ=255)

                                                                                                                              aaaTTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGctctc                       >  2:2715954/1‑150 (MQ=255)

                                                                                                                              aaaTTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGctctc                       >  2:2105575/1‑150 (MQ=255)

                                                                                                                               aaTTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGctctct                      <  1:3730220/150‑1 (MQ=255)

                                                                                                                                 tttCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTc                    >  2:3278791/1‑150 (MQ=255)

                                                                                                                                 tttCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTc                    <  1:3101014/150‑1 (MQ=255)

                                                                                                                                  ttCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCa                   >  1:2122073/1‑150 (MQ=255)

                                                                                                                                   tCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCaa                  <  2:3664890/150‑1 (MQ=255)

                                                                                                                                      tACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACg               <  2:1936125/150‑1 (MQ=255)

                                                                                                                                       aCCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGc              >  2:1153011/1‑150 (MQ=255)

                                                                                                                                       aCCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGc              >  2:1462760/1‑150 (MQ=255)

                                                                                                                                         cACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCtt            >  1:3343819/1‑150 (MQ=255)

                                                                                                                                           cTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTCt          <  1:3898142/150‑1 (MQ=255)

                                                                                                                                           cTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTCt          <  1:670217/150‑1 (MQ=255)

                                                                                                                                           cTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTCt          <  1:634127/150‑1 (MQ=255)

                                                                                                                                              ccGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTCTGGc       <  2:2457075/150‑1 (MQ=255)

                                                                                                                                               cGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTCTGGCg      >  1:2077601/1‑150 (MQ=255)

                                                                                                                                               cGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTCTGGCg      <  2:2102794/150‑1 (MQ=255)

                                                                                                                                                gCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTCTGGCGa     <  1:2491811/150‑1 (MQ=255)

                                                                                                                                                   cGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTCTGGCGACAt  >  1:1474524/1‑150 (MQ=255)

                                                                                                                                                   cGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTCTGGCGACAt  <  2:136811/150‑1 (MQ=255)

TGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTCTGGCGACAT  >  CP009273/1984330‑1984626

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