BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,815,043	C→A	7.5%	R86S (CGT→AGT)	gmk →	guanylate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,815,043	0	C	A	7.5%	932.0 / 33.0	251	R86S (CGT→AGT)	gmk	guanylate kinase
Reads supporting (aligned to +/- strand): ref base C (102/130); new base A (14/5); total (116/135)
Fisher's exact test for biased strand distribution p-value = 1.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 262 total aligned reads displayed.

CTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGACCGCCGTCTA  >  CP009273/3814894‑3815192

cTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTcgc                                                                                                                                                       >  2:3846760/1‑150 (MQ=255)

 tGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTcgcg                                                                                                                                                      >  1:3118741/1‑150 (MQ=255)

     tCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAg                                                                                                                                                  <  1:3112973/150‑1 (MQ=255)

     tCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAg                                                                                                                                                  <  2:188788/150‑1 (MQ=255)

       acacacCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGc                                                                                                                                                <  1:3347299/150‑1 (MQ=255)

          caccacacGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCAt                                                                                                                                             >  2:232277/1‑150 (MQ=255)

            ccacacGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCATTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTg                                                                                                                                           <  1:4492859/150‑1 (MQ=255)

                   ccAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGt                                                                                                                                    >  1:1247172/1‑150 (MQ=255)

                    cAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTa                                                                                                                                   >  2:1472082/1‑150 (MQ=255)

                     aaCCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTAc                                                                                                                                  <  1:1465691/150‑1 (MQ=255)

                      aCCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACt                                                                                                                                 <  2:1861151/150‑1 (MQ=255)

                        cgcgTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTgg                                                                                                                               <  2:691831/150‑1 (MQ=255)

                         gcgTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTGGc                                                                                                                              <  2:1650486/150‑1 (MQ=255)

                          cgTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCg                                                                                                                             >  1:325052/1‑150 (MQ=255)

                          cgTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCg                                                                                                                             >  1:326475/1‑150 (MQ=255)

                          cgTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCg                                                                                                                             >  1:3343379/1‑150 (MQ=255)

                           gTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGa                                                                                                                            <  2:2940003/150‑1 (MQ=255)

                               tGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCgg                                                                                                                        >  1:240553/1‑150 (MQ=255)

                                 gTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGgtg                                                                                                                      <  2:3728881/150‑1 (MQ=255)

                                   gAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGAAGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTc                                                                                                                    <  2:4401661/150‑1 (MQ=255)

                                    aaGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCg                                                                                                                   <  2:310114/150‑1 (MQ=255)

                                    aaGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCg                                                                                                                   >  1:510846/1‑150 (MQ=255)

                                     aGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGa                                                                                                                  >  2:1702782/1‑150 (MQ=255)

                                      gTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGAt                                                                                                                 >  1:1411893/1‑150 (MQ=255)

                                         cACGGTGAACATTATTTCTTTGTTCATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGtt                                                                                                              >  1:2994274/1‑150 (MQ=255)

                                           cGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGtttt                                                                                                            >  2:2101764/1‑150 (MQ=255)

                                               gAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTtctc                                                                                                        >  1:1242534/1‑150 (MQ=255)

                                                  cATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGat                                                                                                     <  1:2675810/150‑1 (MQ=255)

                                                  cATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGat                                                                                                     <  1:2675248/150‑1 (MQ=255)

                                                  cATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGat                                                                                                     >  1:688270/1‑150 (MQ=255)

                                                   attattTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGata                                                                                                    >  1:1768373/1‑150 (MQ=255)

                                                      attTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCg                                                                                                 >  2:3288728/1‑150 (MQ=255)

                                                      attTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCg                                                                                                 <  1:3508444/150‑1 (MQ=255)

                                                      attTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCg                                                                                                 <  1:563386/150‑1 (MQ=255)

                                                         tCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACt                                                                                              >  2:3485132/1‑150 (MQ=255)

                                                         tCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACt                                                                                              <  2:713657/150‑1 (MQ=255)

                                                          cTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTg                                                                                             <  2:1063448/150‑1 (MQ=255)

                                                          cTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTg                                                                                             <  1:2905975/150‑1 (MQ=255)

                                                            ttGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGc                                                                                           <  1:1779581/150‑1 (MQ=255)

                                                             tGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTACCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCa                                                                                          >  2:3825461/1‑150 (MQ=255)

                                                               ttAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAgg                                                                                        >  1:4310239/1‑150 (MQ=255)

                                                                    cATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGgcgcg                                                                                   >  2:149254/1‑150 (MQ=255)

                                                                    cATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGgcgcg                                                                                   >  2:1957241/1‑150 (MQ=255)

                                                                    cATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGgcgcg                                                                                   >  2:2449395/1‑150 (MQ=255)

                                                                     atgatgAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGgcgcgc                                                                                  >  1:3382075/1‑150 (MQ=255)

                                                                     atgatgAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGACACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGgcgcgc                                                                                  >  2:3974493/1‑150 (MQ=255)

                                                                      tgatgaATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCgca                                                                                 >  2:4232507/1‑150 (MQ=255)

                                                                          gaATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCaa                                                                             >  1:2319540/1‑150 (MQ=255)

                                                                          gaATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCaa                                                                             >  1:4378658/1‑150 (MQ=255)

                                                                               tAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCg                                                                        >  2:2599412/1‑150 (MQ=255)

                                                                               tAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCg                                                                        >  2:3218048/1‑150 (MQ=255)

                                                                               tAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCg                                                                        <  1:2271211/150‑1 (MQ=255)

                                                                                 aaGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGcc                                                                      >  1:1661241/1‑150 (MQ=255)

                                                                                   gAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCag                                                                    <  1:3812678/150‑1 (MQ=255)

                                                                                    aaaTGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCaga                                                                   <  2:3118800/150‑1 (MQ=255)

                                                                                    aaaTGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCaga                                                                   >  2:4181667/1‑150 (MQ=255)

                                                                                       tGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGAGCGCGCAGCAAATTCGCCagaaga                                                                >  2:111146/1‑150 (MQ=255)

                                                                                        gATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGAt                                                               <  2:492535/150‑1 (MQ=255)

                                                                                         aTTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATg                                                              <  1:501781/150‑1 (MQ=255)

                                                                                         aTTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATg                                                              >  2:553278/1‑150 (MQ=255)

                                                                                             gCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGc                                                          >  2:4122454/1‑150 (MQ=255)

                                                                                             gCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGc                                                          >  1:4061091/1‑150 (MQ=255)

                                                                                               agagATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCAc                                                        >  2:1534977/1‑150 (MQ=255)

                                                                                               agagATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCAc                                                        >  2:3564911/1‑150 (MQ=255)

                                                                                                      cGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGa                                                 <  2:3293865/150‑1 (MQ=255)

                                                                                                      cGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGa                                                 <  2:2687188/150‑1 (MQ=255)

                                                                                                          ccTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTAt                                             <  1:3507713/150‑1 (MQ=255)

                                                                                                          ccTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTAt                                             <  1:3532983/150‑1 (MQ=255)

                                                                                                           cTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATc                                            >  1:2195738/1‑150 (MQ=255)

                                                                                                              gAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCttt                                         >  2:256086/1‑150 (MQ=255)

                                                                                                              gAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGAGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCttt                                         <  2:3688341/150‑1 (MQ=255)

                                                                                                               aaCACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTa                                        <  1:3564911/150‑1 (MQ=255)

                                                                                                                acacGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTAt                                       <  1:1047516/150‑1 (MQ=255)

                                                                                                                 cacGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTAtt                                      >  1:1753550/1‑150 (MQ=255)

                                                                                                                   cGCAGAAGTTTTTGGTAATTACTATGACACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTAtttt                                    <  1:3974493/150‑1 (MQ=255)

                                                                                                                    gCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTa                                   >  1:127575/1‑150 (MQ=255)

                                                                                                                     cAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTAc                                  >  2:4386020/1‑150 (MQ=255)

                                                                                                                      aGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTAcc                                 >  1:2108645/1‑150 (MQ=255)

                                                                                                                            ttttGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTc                           <  1:3233546/150‑1 (MQ=255)

                                                                                                                            ttttGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTc                           <  1:956126/150‑1 (MQ=255)

                                                                                                                                     ttACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGa                  <  2:1475049/150‑1 (MQ=255)

                                                                                                                                     ttACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGa                  >  2:3758926/1‑150 (MQ=255)

                                                                                                                                     ttACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGa                  >  2:3758969/1‑150 (MQ=255)

                                                                                                                                     ttACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGa                  <  1:4146885/150‑1 (MQ=255)

                                                                                                                                      tACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGaa                 >  2:51747/1‑150 (MQ=255)

                                                                                                                                       aCTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAAc                >  2:3005617/1‑150 (MQ=255)

                                                                                                                                           tGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGa            >  1:3424542/1‑150 (MQ=255)

                                                                                                                                            ggCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGAc           <  2:1899168/150‑1 (MQ=255)

                                                                                                                                            ggCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGAc           <  2:1899069/150‑1 (MQ=255)

                                                                                                                                             gCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGAcc          >  1:2550913/1‑150 (MQ=255)

                                                                                                                                             gCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGAcc          >  1:4225792/1‑150 (MQ=255)

                                                                                                                                              cACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGAccg         <  1:2364338/150‑1 (MQ=255)

                                                                                                                                              cACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGAccg         <  2:1157076/150‑1 (MQ=255)

                                                                                                                                              cACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGAccg         <  2:1156785/150‑1 (MQ=255)

                                                                                                                                              cACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGAccg         <  2:1083040/150‑1 (MQ=255)

                                                                                                                                                cTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGAccgcc       <  1:3003600/150‑1 (MQ=255)

                                                                                                                                                cTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGAccgcc       <  1:3059045/150‑1 (MQ=255)

                                                                                                                                                    gcgTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGACCGCCGTCt   <  1:3846744/150‑1 (MQ=255)

                                                                                                                                                    gcgTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGACCGCCGTCt   <  2:80465/150‑1 (MQ=255)

                                                                                                                                                     cgTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGACCGCCGTCTa  >  2:2145899/1‑150 (MQ=255)

CTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATGCGTTCCTCGAACACGCAGAAGTTTTTGGTAATTACTATGGCACTTCGCGTGAGGCCATTGAGCAAGTACTGGCGACCGGTGTCGATGTTTTTCTCGATATCGACTGGCAGGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGACCGCCGTCTA  >  CP009273/3814894‑3815192

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