BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,522,558	G→T	5.1%	A177A (GCG→GCT)	yrfG →	GMP/IMP nucleotidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,522,558	0	G	T	5.1%	1243.8 / 30.3	323	A177A (GCG→GCT)	yrfG	GMP/IMP nucleotidase
Reads supporting (aligned to +/- strand): ref base G (158/148); new base T (9/8); total (167/156)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 348 total aligned reads displayed.

GAGCATACCGGTCTGGACGCACACCTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATCCCCTCGCTAATGTGAAGGAGACGCCA  >  CP009273/3522409‑3522707

gAGCATACCGGTCTGGACGCACACCTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACagcgaagcg                                                                                                                                                       <  2:4660254/150‑1 (MQ=255)

gAGCATACCGGTCTGGACGCACACCTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACagcgaagcg                                                                                                                                                       <  1:2753786/150‑1 (MQ=255)

  gCATACCGGTCTGGACGCACACCTTGATTTATTACTGTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGAt                                                                                                                                                     <  1:2120845/150‑1 (MQ=255)

      aCCGGTCTGGACGCACACCTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATtctc                                                                                                                                                 <  2:1955340/150‑1 (MQ=255)

       ccGGTCTGGACGCACACCTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCg                                                                                                                                                <  2:1201437/150‑1 (MQ=255)

       ccGGTCTGGACGCACACCTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCg                                                                                                                                                <  2:779188/150‑1 (MQ=255)

       ccGGTCTGGACGCACACCTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCg                                                                                                                                                <  2:1254270/150‑1 (MQ=255)

         gctCTGGACGCACACCTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGAt                                                                                                                                              >  1:3948175/3‑150 (MQ=255)

          gTCTGGACGCACACCTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGAtg                                                                                                                                             <  2:2581570/150‑1 (MQ=255)

            cTGGACGCACACCTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCTATTCTCGAtgct                                                                                                                                           <  2:294495/150‑1 (MQ=255)

            cTGGACGCACACCTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGAtgct                                                                                                                                           <  2:515618/150‑1 (MQ=255)

             tGGACGCACACCTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGAtgctg                                                                                                                                          <  1:3120260/150‑1 (MQ=255)

                       ccTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCTATTCTCGATGCTGCCGCGCAAtt                                                                                                                                >  1:2636995/1‑150 (MQ=255)

                       ccTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAAtt                                                                                                                                <  1:4138311/150‑1 (MQ=255)

                        cTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAAttt                                                                                                                               <  2:1737995/150‑1 (MQ=255)

                        cTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAAttt                                                                                                                               >  1:2252117/1‑150 (MQ=255)

                         ttGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTg                                                                                                                              <  2:2216148/150‑1 (MQ=255)

                         ttGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTg                                                                                                                              >  2:1431799/1‑150 (MQ=255)

                          tGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTgg                                                                                                                             <  1:1487757/150‑1 (MQ=255)

                              ttattaCTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTAtt                                                                                                                         <  2:853521/150‑1 (MQ=255)

                               tattaCTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTc                                                                                                                        <  2:376530/150‑1 (MQ=255)

                                 ttaCTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGt                                                                                                                      >  2:1007329/1‑150 (MQ=255)

                                 ttaCTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGt                                                                                                                      >  2:4548974/1‑150 (MQ=255)

                                  taCTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGtt                                                                                                                     >  2:2244021/1‑150 (MQ=255)

                                  taCTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGtt                                                                                                                     <  1:4620657/150‑1 (MQ=255)

                                   aCTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTa                                                                                                                    >  1:805876/1‑150 (MQ=255)

                                   aCTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTa                                                                                                                    <  2:4450420/150‑1 (MQ=255)

                                   aCTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTa                                                                                                                    >  1:4312732/1‑150 (MQ=255)

                                    cTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTAc                                                                                                                   <  1:3581823/150‑1 (MQ=255)

                                       tCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGc                                                                                                                >  2:96939/1‑150 (MQ=255)

                                        ccacccacACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCACTTTGGTATTCGTTACTGcc                                                                                                               >  2:2658796/1‑150 (MQ=255)

                                        ccacccacACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGcc                                                                                                               >  2:3131904/1‑150 (MQ=255)

                                          acccacACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTAATGCCTc                                                                                                             >  1:3209390/1‑150 (MQ=255)

                                           cccacACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCg                                                                                                            >  1:3846185/1‑150 (MQ=255)

                                            ccacACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCgg                                                                                                           <  2:945227/150‑1 (MQ=255)

                                             cacacaTTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGc                                                                                                          >  2:385214/1‑150 (MQ=255)

                                                acaTTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTg                                                                                                       >  2:260461/1‑150 (MQ=255)

                                                   tttGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACt                                                                                                    <  2:3066233/150‑1 (MQ=255)

                                                     tGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACtca                                                                                                  >  2:4117401/1‑148 (MQ=255)

                                                         tATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCt                                                                                              <  1:3741332/150‑1 (MQ=255)

                                                         tATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCt                                                                                              <  1:1093927/150‑1 (MQ=255)

                                                           tCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCTATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGa                                                                                            >  1:4668015/1‑150 (MQ=255)

                                                             cGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCTATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGAtt                                                                                          >  2:2055432/1‑150 (MQ=255)

                                                             cGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGAtt                                                                                          >  1:2916062/1‑150 (MQ=255)

                                                             cGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGAtt                                                                                          >  1:220128/1‑150 (MQ=255)

                                                              gAAAGAGGATCAGCGGGTATGGCATGCGGGGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTc                                                                                         <  2:2847989/150‑1 (MQ=255)

                                                               aaaGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTcc                                                                                        >  1:883834/1‑150 (MQ=255)

                                                                aaGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCg                                                                                       >  2:241451/1‑150 (MQ=255)

                                                                  gagGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCggg                                                                                     <  2:2720557/150‑1 (MQ=255)

                                                                   agGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTCATCCTGATTCCgggc                                                                                    >  2:3596849/1‑149 (MQ=255)

                                                                     gATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGAtt                                                                                  <  1:550368/150‑1 (MQ=255)

                                                                     gATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGAtt                                                                                  <  2:4170531/150‑1 (MQ=255)

                                                                     gATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGAtt                                                                                  >  1:2132557/1‑150 (MQ=255)

                                                                      aTCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCTATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTg                                                                                 <  2:553479/150‑1 (MQ=255)

                                                                       tCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCTATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGc                                                                                >  1:1139830/1‑150 (MQ=255)

                                                                        cAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGcc                                                                               <  2:3539789/150‑1 (MQ=255)

                                                                         aGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCg                                                                              <  2:2679503/150‑1 (MQ=255)

                                                                            ggTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCTATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCgaga                                                                           >  1:1026096/1‑150 (MQ=255)

                                                                            ggTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCgaga                                                                           >  2:4523147/1‑150 (MQ=255)

                                                                             gTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGaa                                                                          >  1:3913548/1‑150 (MQ=255)

                                                                             gTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGaa                                                                          >  1:168277/1‑150 (MQ=255)

                                                                               tATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAAc                                                                        >  1:4211035/1‑150 (MQ=255)

                                                                                aTGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTCATCCTGATTCCGGGATTGCCGAGAAACa                                                                       >  2:3628870/1‑150 (MQ=255)

                                                                                aTGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACa                                                                       >  2:1838893/1‑150 (MQ=255)

                                                                                aTGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACa                                                                       >  2:3186735/1‑150 (MQ=255)

                                                                                         ggTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACg                                                              >  1:3466007/1‑150 (MQ=255)

                                                                                         ggTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACg                                                              >  2:3319208/1‑150 (MQ=255)

                                                                                            ggCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCa                                                           <  2:1010528/150‑1 (MQ=255)

                                                                                             gCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCAt                                                          >  1:1918955/1‑150 (MQ=255)

                                                                                              ccGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATc                                                         >  1:1761739/1‑150 (MQ=255)

                                                                                                gAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCg                                                       >  2:2183048/1‑150 (MQ=255)

                                                                                                   gCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCa                                                    >  2:715733/1‑150 (MQ=255)

                                                                                                     tACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACt                                                  <  1:1070372/150‑1 (MQ=255)

                                                                                                     tACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACt                                                  <  1:127481/150‑1 (MQ=255)

                                                                                                        gggTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGaa                                               <  1:2823120/150‑1 (MQ=255)

                                                                                                         ggTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAAt                                              >  2:4204015/1‑150 (MQ=255)

                                                                                                         ggTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAAt                                              >  1:4589848/1‑150 (MQ=255)

                                                                                                          gTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATg                                             <  2:3749980/150‑1 (MQ=255)

                                                                                                            ctgaaagctgaaagAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGAc                                           >  2:448767/1‑150 (MQ=255)

                                                                                                             tgaaagctgaaagAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACt                                          <  2:4314573/150‑1 (MQ=255)

                                                                                                               aaagctgaaagAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTAc                                        >  2:3413262/1‑150 (MQ=255)

                                                                                                               aaagctgaaagAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTAc                                        >  1:147398/1‑150 (MQ=255)

                                                                                                               aaagctgaaagAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTAc                                        >  2:1453454/1‑150 (MQ=255)

                                                                                                                aagctgaaagAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTAcc                                       <  2:3012563/150‑1 (MQ=255)

                                                                                                                 agctgaaagAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTAccg                                      <  1:124713/150‑1 (MQ=255)

                                                                                                                  gctgaaagAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTAccgc                                     <  1:1387602/150‑1 (MQ=255)

                                                                                                                   ctgaaagAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTAccgcc                                    >  2:1703920/1‑150 (MQ=255)

                                                                                                                      aaagAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCgccgcc                                 >  1:2428378/1‑150 (MQ=255)

                                                                                                                       aagAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCt                                <  1:3150911/150‑1 (MQ=255)

                                                                                                                        agAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTg                               <  2:3521115/150‑1 (MQ=255)

                                                                                                                                 tttATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATCCCCTCg                      >  1:4426357/1‑150 (MQ=255)

                                                                                                                                 tttATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATCCCCTCg                      >  1:2041509/1‑150 (MQ=255)

                                                                                                                                     ttGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATCCCCTCGCTaa                  >  1:2523326/1‑150 (MQ=255)

                                                                                                                                      tgatgaCAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATCCCCTCGCTAAt                 >  1:981396/1‑150 (MQ=255)

                                                                                                                                        atgaCAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATCCCCTCGCTAAtgt               <  2:697177/150‑1 (MQ=255)

                                                                                                                                        atgaCAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATCCCCTCGCTAAtgt               <  1:3236578/150‑1 (MQ=255)

                                                                                                                                         tgaCAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATCCCCTCGCTAAtgtg              >  1:3681036/1‑150 (MQ=255)

                                                                                                                                           aCAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATCCCCTCGCTAATGTGaa            <  1:2474608/150‑1 (MQ=255)

                                                                                                                                                gaagcgaTTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATCCCCTCGCTAATGTGAAGgaga       >  1:3379778/1‑150 (MQ=255)

                                                                                                                                                     gaTTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATCCCCTCGCTAATGTGAAGGAGACGCCa  >  2:211437/1‑150 (MQ=255)

GAGCATACCGGTCTGGACGCACACCTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATCCCCTCGCTAATGTGAAGGAGACGCCA  >  CP009273/3522409‑3522707

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