BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,889,967	T→A	12.2%	A421A (GCA→GCT)	ygcS ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,889,967	0	T	A	12.2%	782.4 / 87.8	255	A421A (GCA→GCT)	ygcS	putative transporter
Reads supporting (aligned to +/- strand): ref base T (97/127); new base A (14/17); total (111/144)
Fisher's exact test for biased strand distribution p-value = 8.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

Only 100 of 256 total aligned reads displayed.

TCCGGTCCGCTTTTACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTACTTACCGCCGCGCCAAGTCGACTCATGGCAGTGGCAAAACCGAC  >  CP009273/2889820‑2890116

tCCGGTCCGCTTTTACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGAGc                                                                                                                                                     <  1:3552118/150‑1 (MQ=255)

 ccGGTCCGCTTTTACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGcc                                                                                                                                                    <  1:1109041/150‑1 (MQ=255)

 ccGGTCCGCTTTTACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGcc                                                                                                                                                    <  2:2264561/150‑1 (MQ=255)

    gTCCGCTTTTACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGAGCCCAt                                                                                                                                                 >  2:4004728/1‑150 (MQ=255)

      ccGCTTTTACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAg                                                                                                                                               <  2:2563151/150‑1 (MQ=255)

       cGCTTTTACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGc                                                                                                                                              >  2:1122036/1‑150 (MQ=255)

           tttACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAgg                                                                                                                                          <  2:3833985/150‑1 (MQ=255)

             tACCCTGAGCCGCCGCCTGCCACTCTGTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACAGTTCCCGCCGCCACCAGCGGGAGGGCTTTCGTTTCTGGTGCCCAGAGCCAgggc                                                                                                                                        >  2:2274207/1‑148 (MQ=255)

                ccTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACCTTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTcaca                                                                                                                                     <  2:3107693/150‑1 (MQ=255)

                 cTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGAGCCCATAGCCAGGTCacaa                                                                                                                                    >  2:2698606/1‑150 (MQ=255)

                      ccgccgccTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACaaaa                                                                                                                               >  2:4408365/1‑150 (MQ=255)

                          cgccTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGAGCCCATAGCCAGGTCACAACAAaaccaa                                                                                                                           >  2:1076492/1‑150 (MQ=255)

                           gccTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAaaccaac                                                                                                                          <  2:2677476/150‑1 (MQ=255)

                           gccTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGAGCCCATAGCCAGGTCACAACAAaaccaac                                                                                                                          <  2:1036677/150‑1 (MQ=255)

                            ccTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAaaccaacc                                                                                                                         >  2:2465885/1‑150 (MQ=255)

                                 cACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACcaacaa                                                                                                                    >  2:65861/1‑150 (MQ=255)

                                 cACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGAGCCCATAGCCAGGTCACAACAAAACCAACcaacaa                                                                                                                    <  1:1614195/150‑1 (MQ=255)

                                  aCTCTTTGTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCaacaac                                                                                                                   <  1:2937037/150‑1 (MQ=255)

                                        tttCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTc                                                                                                             <  1:4480529/150‑1 (MQ=255)

                                        tttCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTc                                                                                                             <  1:1027139/150‑1 (MQ=255)

                                        tttCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGAGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTc                                                                                                             <  2:2923728/150‑1 (MQ=255)

                                         ttCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCg                                                                                                            <  1:3189055/150‑1 (MQ=255)

                                           ccGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGcc                                                                                                          <  1:3950826/150‑1 (MQ=255)

                                           ccGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGcc                                                                                                          <  2:3829829/150‑1 (MQ=255)

                                            cGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCa                                                                                                         <  2:3912171/150‑1 (MQ=255)

                                            cGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCa                                                                                                         <  1:302774/150‑1 (MQ=255)

                                             gCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAg                                                                                                        <  1:2978648/150‑1 (MQ=255)

                                             gCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAg                                                                                                        >  1:615847/1‑150 (MQ=255)

                                              ccAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAgg                                                                                                       <  2:1545451/150‑1 (MQ=255)

                                                aGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAg                                                                                                     <  1:1500970/150‑1 (MQ=255)

                                                aGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAg                                                                                                     <  1:500659/150‑1 (MQ=255)

                                                aGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGAGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAg                                                                                                     >  2:752722/1‑150 (MQ=255)

                                                 gCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGc                                                                                                    <  2:4069674/150‑1 (MQ=255)

                                                 gCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGc                                                                                                    <  1:1856646/150‑1 (MQ=255)

                                                   aTTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCaa                                                                                                  >  1:3980906/1‑150 (MQ=255)

                                                     tCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATa                                                                                                <  1:2476657/150‑1 (MQ=255)

                                                     tCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATa                                                                                                <  1:1528746/150‑1 (MQ=255)

                                                       cGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAg                                                                                              <  1:3784807/150‑1 (MQ=255)

                                                       cGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAg                                                                                              >  2:4181801/1‑150 (MQ=255)

                                                        ggCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAgg                                                                                             <  1:760041/150‑1 (MQ=255)

                                                         gCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGt                                                                                            <  2:2104386/150‑1 (MQ=255)

                                                            tCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGAGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGAc                                                                                         <  1:3120555/150‑1 (MQ=255)

                                                               ggTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTg                                                                                      <  1:2361093/150‑1 (MQ=255)

                                                                gTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGc                                                                                     >  2:3126834/1‑150 (MQ=255)

                                                                 ttCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCa                                                                                    >  2:1295055/1‑150 (MQ=255)

                                                                  tCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCAt                                                                                   <  1:1916178/150‑1 (MQ=255)

                                                                   cGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGAGCCCATAGCCAGGTCACAACAAAACCAACCAACAACCCTGTCGCCAGGAGCAATAAGGTGACTTGCAtt                                                                                  <  2:4139317/150‑1 (MQ=255)

                                                                    ggCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTc                                                                                 >  1:2679631/1‑150 (MQ=255)

                                                                    ggCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGAGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTc                                                                                 >  2:1358180/1‑150 (MQ=255)

                                                                     gCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTcc                                                                                <  1:1234663/150‑1 (MQ=255)

                                                                      cTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTccc                                                                               <  1:3653059/150‑1 (MQ=255)

                                                                       tttAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTcccc                                                                              <  2:3752126/150‑1 (MQ=255)

                                                                        ttAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCa                                                                             <  1:570129/150‑1 (MQ=255)

                                                                        ttAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCa                                                                             >  1:3425336/1‑150 (MQ=255)

                                                                        ttAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGAGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCa                                                                             >  2:3444626/1‑150 (MQ=255)

                                                                        ttAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGAGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCa                                                                             >  1:283582/1‑150 (MQ=255)

                                                                          aaaCGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACt                                                                           <  2:421657/150‑1 (MQ=255)

                                                                           aaCGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTg                                                                          >  1:4128104/1‑150 (MQ=255)

                                                                           aaCGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTg                                                                          <  2:3670811/150‑1 (MQ=255)

                                                                              gCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTgcgc                                                                       >  1:3472423/1‑150 (MQ=255)

                                                                                      aaTGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCAccc                                                               >  2:3569042/1‑150 (MQ=255)

                                                                                       aTGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCa                                                              <  1:4251845/150‑1 (MQ=255)

                                                                                        tGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCAc                                                             <  1:512900/150‑1 (MQ=255)

                                                                                        tGTTCATTCGCACCTCCGACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTGTCTGGAGCCCATAGCCAGGTCACAACAAAACCAACCAACACCACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCAc                                                             <  2:1215908/150‑1 (MQ=255)

                                                                                         gTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACg                                                            >  2:1800703/1‑150 (MQ=255)

                                                                                             aTTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGgcag                                                        >  2:1145174/1‑150 (MQ=255)

                                                                                               tCGCACCTCCTACAGTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGgcagca                                                      <  1:3369255/150‑1 (MQ=255)

                                                                                                 gCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAgg                                                    >  2:801804/1‑150 (MQ=255)

                                                                                                   aCCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGcc                                                  <  1:1751741/150‑1 (MQ=255)

                                                                                                   aCCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGcc                                                  >  1:2703425/1‑150 (MQ=255)

                                                                                                    cctcctACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGAGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCa                                                 <  1:2257378/150‑1 (MQ=255)

                                                                                                       cctACATGTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTa                                              <  1:4289308/150‑1 (MQ=255)

                                                                                                        ccaCATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTAc                                             <  2:1621408/148‑1 (MQ=255)

                                                                                                         tACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTACt                                            >  2:669159/1‑150 (MQ=255)

                                                                                                          aCATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTACtt                                           <  2:1787119/150‑1 (MQ=255)

                                                                                                          aCATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTACtt                                           <  1:262474/150‑1 (MQ=255)

                                                                                                             tttCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTACTTAcc                                        <  2:2222119/150‑1 (MQ=255)

                                                                                                             tttCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTACTTAcc                                        <  2:4493676/150‑1 (MQ=255)

                                                                                                               tCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTACTTAccgc                                      >  2:1006988/1‑150 (MQ=255)

                                                                                                                 ccgccgCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTACTTAccgccg                                    >  1:2057923/1‑150 (MQ=255)

                                                                                                                      gccACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCGGGCCAGTACTTACCGCCGCGCCa                               <  1:933001/150‑1 (MQ=255)

                                                                                                                      gccACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTACTTACCGCCGCGCCa                               <  1:4440382/150‑1 (MQ=255)

                                                                                                                          ccaGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTACTTACCGCCGCGCCAAGTc                           <  1:1758081/150‑1 (MQ=255)

                                                                                                                          ccaGCGGGAGGGCTTTAGTGTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTACTTACCGCCGCGCCAAGTc                           <  1:182642/150‑1 (MQ=255)

                                                                                                                               gggAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTACTTACCGCCGCGCCAAGTCGACTc                      >  1:1328646/1‑150 (MQ=255)

                                                                                                                                 gAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTACTTACCGCCGCGCCAAGTCGACTCAt                    <  1:2079696/150‑1 (MQ=255)

                                                                                                                                 gAGGGCTTTAGTTTCTGGAGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTACTTACCGCCGCGCCAAGTCGACTCAt                    >  1:3500369/1‑150 (MQ=255)

                                                                                                                                  aGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTACTTACCGCCGCGCCAAGTCGACTCATg                   >  2:2425488/1‑150 (MQ=255)

                                                                                                                                  aGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTACTTACCGCCGCGCCAAGTCGACTCATg                   >  1:3807730/1‑150 (MQ=255)

                                                                                                                                  aGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTACTTACCGCCGCGCCAAGTCGACTCATg                   <  1:1230243/150‑1 (MQ=255)

                                                                                                                                   gggCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTACTTACCGCCGCGCCAAGTCGACTCATgg                  >  2:2495373/1‑150 (MQ=255)

                                                                                                                                   gggCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTACTTACCGCCGCGCCAAGTCGACTCATgg                  >  2:4091533/1‑150 (MQ=255)

                                                                                                                                    ggCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCTGTACTTACCGCCGCGCCAAGTCGACTCATGGc                 <  1:4408365/150‑1 (MQ=255)

                                                                                                                                    ggCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTACTTACCGCCGCGCCAAGTCGACTCATGGc                 <  1:2132502/150‑1 (MQ=255)

                                                                                                                                     gCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTACTTACCGCCGCGCCAAGTCGACTCATGGCa                >  2:1191441/1‑150 (MQ=255)

                                                                                                                                        ttAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTACTTACCGCCGCGCCAAGTCGACTCATGGCAGTg             <  2:1682298/150‑1 (MQ=255)

                                                                                                                                         tAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTACTTACCGCCGCGCCAAGTCGACTCATGGCAGTgg            >  2:2069719/1‑150 (MQ=255)

                                                                                                                                            tttCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTACTTACCGCCGCGCCAAGTCGACTCATGGCAGTGGCaa         >  2:4481657/1‑150 (MQ=255)

                                                                                                                                                tGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTACTTACCGCCGCGCCAAGTCGACTCATGGCAGTGGCAAAAcc     >  2:2439676/1‑150 (MQ=255)

                                                                                                                                                 ggTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTACTTACCGCCGCGCCAAGTCGACTCATGGCAGTGGCAAAACCg    <  1:2747095/150‑1 (MQ=255)

TCCGGTCCGCTTTTACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTACTTACCGCCGCGCCAAGTCGACTCATGGCAGTGGCAAAACCGAC  >  CP009273/2889820‑2890116

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: `‑`