BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,033,430	G→T	13.9%	D73Y (GAC→TAC)	fldB →	flavodoxin 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,033,430	0	G	T	13.9%	1102.0 / 112.2	337	D73Y (GAC→TAC)	fldB	flavodoxin 2
Reads supporting (aligned to +/- strand): ref base G (139/151); new base T (16/31); total (155/182)
Fisher's exact test for biased strand distribution p-value = 8.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 354 total aligned reads displayed.

ATTATCGGCCCAGAACTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCTACTGGCCAACGGAA  >  CP009273/3033283‑3033579

attatCGGCCCAGAACTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTctac                                                                                                                                                     <  2:4417303/150‑1 (MQ=255)

attatCGGCCCAGAACTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTcgac                                                                                                                                                     >  2:905760/1‑150 (MQ=255)

attatCGGCCCAGAACTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTcgac                                                                                                                                                     >  1:3859255/1‑150 (MQ=255)

 ttatCGGCCCAGAACTGGTGACCGTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTcgacg                                                                                                                                                    <  1:634113/150‑1 (MQ=255)

     cGGCCCAGAACTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGAccn                                                                                                                                                <  2:4220010/150‑2 (MQ=255)

      ggCCCAGAACTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTg                                                                                                                                               >  1:5085624/1‑150 (MQ=255)

       gCCCAGAACTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGa                                                                                                                                              <  2:3280220/150‑1 (MQ=255)

       gCCCAGAACTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGa                                                                                                                                              <  1:384224/150‑1 (MQ=255)

        cccAGAACTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGaa                                                                                                                                             >  2:3322310/1‑150 (MQ=255)

        cccAGAACTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGaa                                                                                                                                             <  1:823743/150‑1 (MQ=255)

              aCTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGa                                                                                                                                       <  2:3764941/150‑1 (MQ=255)

               cTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGaa                                                                                                                                      <  2:2785857/150‑1 (MQ=255)

                tGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAg                                                                                                                                     >  1:1290278/1‑150 (MQ=255)

                 ggTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAgg                                                                                                                                    >  2:3836628/1‑150 (MQ=255)

                   tGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTa                                                                                                                                  <  2:569751/150‑1 (MQ=255)

                   tGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTa                                                                                                                                  >  2:4497399/1‑150 (MQ=255)

                    gACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTaa                                                                                                                                 >  1:24535/1‑150 (MQ=255)

                    gACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTaa                                                                                                                                 <  2:4865424/150‑1 (MQ=255)

                    gACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCAGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTaa                                                                                                                                 >  1:2262750/1‑150 (MQ=255)

                     aCCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTaaa                                                                                                                                >  1:3266692/1‑150 (MQ=255)

                      ccTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTaaaa                                                                                                                               <  1:4916158/150‑1 (MQ=255)

                      ccTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTaaaa                                                                                                                               <  1:4007898/150‑1 (MQ=255)

                        ttACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAAtt                                                                                                                             <  1:5174598/150‑1 (MQ=255)

                           cATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAAttgtt                                                                                                                          >  1:1907579/1‑150 (MQ=255)

                             tAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTgc                                                                                                                        <  1:5535191/150‑1 (MQ=255)

                             tAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTgc                                                                                                                        <  2:1589255/150‑1 (MQ=255)

                             tAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTgc                                                                                                                        >  2:4482805/1‑150 (MQ=255)

                                 cTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTg                                                                                                                    >  2:5164944/1‑150 (MQ=255)

                                 cTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTg                                                                                                                    >  2:3370124/1‑150 (MQ=255)

                                   cAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTa                                                                                                                  <  2:2296745/150‑1 (MQ=255)

                                           acTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTg                                                                                                          <  1:2938434/150‑1 (MQ=255)

                                            cTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTgg                                                                                                         >  1:3578887/1‑150 (MQ=255)

                                              ccccGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCg                                                                                                       <  1:1648011/150‑1 (MQ=255)

                                               cccGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGACGCCGTCTGGGATCAGCTCTCCGACCTGAACCGTGAAGGTCAAATTGTTGCGCTGTATGGGCttttgct                                                                                                      >  2:2476786/1‑145 (MQ=255)

                                               cccGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGa                                                                                                      >  1:746099/1‑150 (MQ=255)

                                               cccGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTAGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGa                                                                                                      >  1:2143643/1‑150 (MQ=255)

                                                ccGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGAt                                                                                                     <  1:3369785/150‑1 (MQ=255)

                                                 cGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATc                                                                                                    >  1:3140228/1‑150 (MQ=255)

                                                  gAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCa                                                                                                   >  2:755475/1‑150 (MQ=255)

                                                       tAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTgg                                                                                              <  1:1965623/150‑1 (MQ=255)

                                                       tAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTgg                                                                                              >  2:3075626/1‑150 (MQ=255)

                                                       tAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTgg                                                                                              >  1:422463/1‑150 (MQ=255)

                                                        aaTGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCACCTggg                                                                                             >  2:5505320/1‑150 (MQ=255)

                                                          tGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGAt                                                                                           <  1:4326482/150‑1 (MQ=255)

                                                            gAGCAGTACGTTGTGCTCATTCTGGGTATCCAGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATAc                                                                                         <  2:1937953/150‑1 (MQ=255)

                                                               cAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGc                                                                                      <  1:4138102/150‑1 (MQ=255)

                                                               cAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGc                                                                                      <  1:2727005/150‑1 (MQ=255)

                                                                aGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCg                                                                                     <  1:4538283/150‑1 (MQ=255)

                                                                aGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCg                                                                                     >  2:4023784/1‑150 (MQ=255)

                                                                aGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCg                                                                                     <  2:2029257/150‑1 (MQ=255)

                                                                aGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCg                                                                                     >  1:4071204/1‑150 (MQ=255)

                                                                aGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCg                                                                                     >  1:5616977/1‑150 (MQ=255)

                                                                  tACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAg                                                                                   <  1:3787335/150‑1 (MQ=255)

                                                                     gATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTgg                                                                                >  1:3439596/1‑150 (MQ=255)

                                                                       tgtgCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGtt                                                                              <  1:3928311/150‑1 (MQ=255)

                                                                        gtgCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTc                                                                             <  2:1879913/150‑1 (MQ=255)

                                                                           cTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTc                                                                          <  2:4178199/150‑1 (MQ=255)

                                                                           cTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTc                                                                          <  1:4275904/150‑1 (MQ=255)

                                                                           cTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTc                                                                          <  2:1731552/150‑1 (MQ=255)

                                                                              aTTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTCAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGAt                                                                       >  2:5011489/1‑150 (MQ=255)

                                                                                 cTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATgcg                                                                    <  1:3543416/150‑1 (MQ=255)

                                                                                 cTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATgcg                                                                    <  2:2062627/150‑1 (MQ=255)

                                                                                   gggTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCt                                                                  <  2:1979445/150‑1 (MQ=255)

                                                                                   gggTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCt                                                                  <  2:5557047/150‑1 (MQ=255)

                                                                                    ggTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTc                                                                 <  1:57820/150‑1 (MQ=255)

                                                                                      tATCCCGACCTGGGATTTTGGTGAAATCCAGGACGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCgg                                                               >  2:1446082/1‑150 (MQ=255)

                                                                                       aTCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGt                                                              <  1:3628941/150‑1 (MQ=255)

                                                                                         cccGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTAt                                                            >  1:3022269/1‑150 (MQ=255)

                                                                                            gACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTAtgct                                                         >  1:5392893/1‑150 (MQ=255)

                                                                                                tGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCAt                                                     >  1:1516879/1‑150 (MQ=255)

                                                                                                tGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCAt                                                     >  1:951977/1‑150 (MQ=255)

                                                                                                 gggATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATg                                                    >  2:2540455/1‑150 (MQ=255)

                                                                                                  ggATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGa                                                   <  2:2409507/150‑1 (MQ=255)

                                                                                                  ggATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGa                                                   >  2:694231/1‑150 (MQ=255)

                                                                                                    aTTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACa                                                 >  2:424879/1‑150 (MQ=255)

                                                                                                    aTTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACa                                                 >  1:4360917/1‑150 (MQ=255)

                                                                                                     ttttGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACaa                                                >  1:4038108/1‑150 (MQ=255)

                                                                                                             aaaTCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGa                                        <  1:260566/150‑1 (MQ=255)

                                                                                                                    ggAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAgg                                 >  1:104987/1‑150 (MQ=255)

                                                                                                                     gAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGc                                >  2:3199366/1‑150 (MQ=255)

                                                                                                                     gAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGc                                <  2:2917789/150‑1 (MQ=255)

                                                                                                                      aaGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCg                               >  2:4666231/1‑150 (MQ=255)

                                                                                                                           tGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAg                          >  2:3739307/1‑150 (MQ=255)

                                                                                                                            gggAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGt                         <  1:1762888/150‑1 (MQ=255)

                                                                                                                            gggAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGt                         >  1:2419329/1‑150 (MQ=255)

                                                                                                                             ggAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTTAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGtt                        <  2:1192712/150‑1 (MQ=255)

                                                                                                                              gAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTtc                       <  2:2324738/150‑1 (MQ=255)

                                                                                                                               aaGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTtcg                      <  2:435630/150‑1 (MQ=255)

                                                                                                                                 gCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTtcgtc                    <  2:1845042/150‑1 (MQ=255)

                                                                                                                                 gCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTtcgtc                    <  2:1666933/150‑1 (MQ=255)

                                                                                                                                 gCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTtcgtc                    <  2:1324162/150‑1 (MQ=255)

                                                                                                                                      cTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCTa               >  1:2090115/1‑150 (MQ=255)

                                                                                                                                        gggATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCTACt             <  2:2505684/150‑1 (MQ=255)

                                                                                                                                          gATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCTACTgg           >  1:2801788/1‑150 (MQ=255)

                                                                                                                                          gATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCTACTgg           >  2:4548512/1‑150 (MQ=255)

                                                                                                                                           aTCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCTACTGGc          >  2:4614704/1‑150 (MQ=255)

                                                                                                                                              aGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCTACTGGCCaa       >  2:5071248/1‑150 (MQ=255)

                                                                                                                                                cTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCTACTGGCCAACg     >  2:220637/1‑150 (MQ=255)

                                                                                                                                                cTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCTACTGGCCAACg     >  1:5337777/1‑150 (MQ=255)

                                                                                                                                                 tCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCTACTGGCCAACgg    <  1:473116/150‑1 (MQ=255)

ATTATCGGCCCAGAACTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCTACTGGCCAACGGAA  >  CP009273/3033283‑3033579

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