BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,933,029	A→T	12.6%	H375L (CAC→CTC)	pykA →	pyruvate kinase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,933,029	0	A	T	12.6%	924.6 / 104.1	325	H375L (CAC→CTC)	pykA	pyruvate kinase II
Reads supporting (aligned to +/- strand): ref base A (149/134); new base T (25/16); total (174/151)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

Only 100 of 344 total aligned reads displayed.

GTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGC  >  CP009273/1932880‑1933178

gTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAacca                                                                                                                                                       <  2:92284/150‑1 (MQ=255)

gTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAacca                                                                                                                                                       <  1:2152743/150‑1 (MQ=255)

    ccGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTg                                                                                                                                                   >  1:5073100/1‑150 (MQ=255)

    ccGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTg                                                                                                                                                   <  1:2731438/150‑1 (MQ=255)

     cGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTCTTGCCATGTCAGCAATGTACGCAGCTAACCTCCTGa                                                                                                                                                  >  1:2624978/1‑150 (MQ=255)

     cGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGa                                                                                                                                                  >  1:1880143/1‑150 (MQ=255)

      gTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGaa                                                                                                                                                 <  1:572995/150‑1 (MQ=255)

         aGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAgg                                                                                                                                              >  2:951844/1‑150 (MQ=255)

         aGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAgg                                                                                                                                              >  1:2118505/1‑150 (MQ=255)

          gAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGc                                                                                                                                             >  2:701974/1‑150 (MQ=255)

          gAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGc                                                                                                                                             >  1:1081538/1‑150 (MQ=255)

            aaCCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGt                                                                                                                                           >  1:3121160/1‑150 (MQ=255)

                gTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACg                                                                                                                                       <  1:2559896/150‑1 (MQ=255)

                   gCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCg                                                                                                                                    <  1:3052613/150‑1 (MQ=255)

                    cAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGa                                                                                                                                   >  2:742694/1‑150 (MQ=255)

                         aTGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGatcatc                                                                                                                              <  1:4956149/150‑1 (MQ=255)

                            gcgcgcGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCTCCTGAAAGGCGTTACGGCGATCATCAcc                                                                                                                           <  1:539288/150‑1 (MQ=255)

                                gcgTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGa                                                                                                                       <  1:4440165/150‑1 (MQ=255)

                                 cgTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGAc                                                                                                                      >  2:3704443/1‑150 (MQ=255)

                                  gTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGAcc                                                                                                                     >  1:3587149/1‑150 (MQ=255)

                                    ttGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGa                                                                                                                   <  2:2227019/150‑1 (MQ=255)

                                     tGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCTCCTGAAAGGCGTTACGGCGATCATCACCATGACCGaa                                                                                                                  >  2:5077406/1‑150 (MQ=255)

                                         tGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATCTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCgg                                                                                                              >  2:3998022/1‑150 (MQ=255)

                                           ggTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGt                                                                                                            >  2:3947824/1‑150 (MQ=255)

                                            gTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTc                                                                                                           <  1:5109731/150‑1 (MQ=255)

                                              gCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGt                                                                                                         <  2:4030924/150‑1 (MQ=255)

                                               cGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTa                                                                                                        >  2:180412/1‑150 (MQ=255)

                                                ggAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTAc                                                                                                       <  2:3186137/150‑1 (MQ=255)

                                                 gAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTAcc                                                                                                      <  1:1652909/150‑1 (MQ=255)

                                                  aaaaaaTCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACcg                                                                                                     <  2:445829/150‑1 (MQ=255)

                                                   aaaaaTCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAATCTATTGCCATGTCAGCTATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACcgc                                                                                                    <  1:1189939/150‑1 (MQ=255)

                                                   aaaaaTCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACcgc                                                                                                    <  1:2024366/150‑1 (MQ=255)

                                                   aaaaaTCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACcgc                                                                                                    <  1:4173774/150‑1 (MQ=255)

                                                    aaaaTCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACcgcg                                                                                                   <  2:1732804/150‑1 (MQ=255)

                                                      aaTCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCt                                                                                                 >  2:1409715/1‑150 (MQ=255)

                                                       aTCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCtg                                                                                                <  2:4172629/150‑1 (MQ=255)

                                                        tCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCtga                                                                                               <  1:1452439/150‑1 (MQ=255)

                                                           cGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCtgatga                                                                                            >  2:4274556/1‑150 (MQ=255)

                                                           cGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCtgatga                                                                                            >  1:3913265/1‑150 (MQ=255)

                                                           cGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCtgatga                                                                                            <  2:2777669/150‑1 (MQ=255)

                                                            gAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCTCCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGGCGTACCGCGCTGATGAc                                                                                           <  1:3164148/150‑1 (MQ=255)

                                                               cATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTc                                                                                        >  1:4740136/1‑150 (MQ=255)

                                                                 tCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTccc                                                                                      <  2:1487293/150‑1 (MQ=255)

                                                                 tCAACGTTTCTAAACACAGTCTTGAAGTTCAGTTCGAAAATGTGTAAGAAGCTATTGCCATTTCAGCATTGTACTCAGCTAACAACCTGTAAGGCGTTACGGCGATCATCAACATGACCGAATCGGTTCGTACCGCTCTGATGACCTcta                                                                                      <  2:1109902/150‑3 (MQ=255)

                                                                    aCGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTa                                                                                   <  2:432755/150‑1 (MQ=255)

                                                                    aCGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTa                                                                                   >  2:130459/1‑150 (MQ=255)

                                                                       tttCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCa                                                                                <  1:3455909/150‑1 (MQ=255)

                                                                            aaaCACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCTCCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGGCGTACCGCGCTGATGACCTCCCGTATCAGctct                                                                           <  1:1533685/150‑1 (MQ=255)

                                                                             aaCACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTg                                                                          <  2:4231593/150‑1 (MQ=255)

                                                                               cacCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGt                                                                        >  2:3895605/1‑150 (MQ=255)

                                                                                acCGTCTGGACGTTCAGTTCGACTATGTGGAAGATGCTATTGCCATGTCAGCTATGTTCGCAGCTTACCACCTGAAAGGCGTTACGGGTATCATGACCATGAGCGATTCTGGTCGTACGGCGCTGTTTACCTCCCGTATCAGCTCTGGTc                                                                       <  2:530787/150‑1 (MQ=255)

                                                                                acCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTc                                                                       >  1:3822185/1‑150 (MQ=255)

                                                                                 ccGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGGACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCt                                                                      <  2:4559856/150‑1 (MQ=255)

                                                                                   gTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCTCCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGGCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGc                                                                    >  2:3553079/1‑150 (MQ=255)

                                                                                   gTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGc                                                                    >  1:2941864/1‑150 (MQ=255)

                                                                                    tCTGGACGTTCAGTTCGACAATCTGGAAGAAGCTATTTCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCAATACCTAATCGGGTCGTAACGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGcc                                                                   >  1:2835361/1‑150 (MQ=255)

                                                                                       ggACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCTCCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAAt                                                                >  2:2054446/1‑150 (MQ=255)

                                                                                        gACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAAtt                                                               >  2:1830671/1‑150 (MQ=255)

                                                                                            ttCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCg                                                           <  1:1791994/150‑1 (MQ=255)

                                                                                            ttCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCg                                                           >  1:3678586/1‑150 (MQ=255)

                                                                                              cAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCTCCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGGCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGcc                                                         >  1:80009/1‑150 (MQ=255)

                                                                                              cAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGcc                                                         >  2:3735093/1‑150 (MQ=255)

                                                                                              cAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGcc                                                         >  1:1018080/1‑150 (MQ=255)

                                                                                              cAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGcc                                                         >  2:1154288/1‑150 (MQ=255)

                                                                                               aGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCa                                                        >  2:606966/1‑150 (MQ=255)

                                                                                                     actaTGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCTCCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGGCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTcgc                                                  <  2:3112892/147‑1 (MQ=255)

                                                                                                     aCAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTcgc                                                  >  1:4584801/1‑150 (MQ=255)

                                                                                                      cAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTcgcg                                                 >  2:1529647/1‑150 (MQ=255)

                                                                                                      cAATGTGGAAGAAGCTAGTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTcgcg                                                 <  2:3830842/150‑1 (MQ=255)

                                                                                                        aTGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGcc                                               >  2:3446124/1‑150 (MQ=255)

                                                                                                           tgGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATg                                            <  2:2941864/150‑1 (MQ=255)

                                                                                                           tgGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATg                                            >  2:4824345/1‑150 (MQ=255)

                                                                                                            ggAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGa                                           >  1:2301795/1‑150 (MQ=255)

                                                                                                               agaagCTATTGCCATGTCAGCAATGTACGCAGCTAACCTCCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGGCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACg                                        >  1:3534041/1‑150 (MQ=255)

                                                                                                                gaagCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGt                                       >  2:2529702/1‑150 (MQ=255)

                                                                                                                  agCTATTGCCATGTCAGCAATGTACGCAGCTAACCTCCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGGCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTAc                                     >  1:4712157/1‑150 (MQ=255)

                                                                                                                   gCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACg                                    <  1:2465191/150‑1 (MQ=255)

                                                                                                                     tATTGCCATGTCAGCTATGTACTCAGCTAACCACCTGAATGGCGTTACTGCGATCATCACCATGACCGAATCGGGTCTTACCGCGCTGATGTCCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCt                                  <  1:4862485/150‑1 (MQ=255)

                                                                                                                     tATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCt                                  <  2:1245255/150‑1 (MQ=255)

                                                                                                                     tATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCt                                  >  1:4993169/1‑150 (MQ=255)

                                                                                                                      aTTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTg                                 >  2:4955075/1‑150 (MQ=255)

                                                                                                                       ttGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGa                                >  1:4915573/1‑150 (MQ=255)

                                                                                                                         gCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAAc                              >  1:1331849/1‑150 (MQ=255)

                                                                                                                          ccATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAAcc                             >  2:2649468/1‑150 (MQ=255)

                                                                                                                                cAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTg                       <  1:4220093/150‑1 (MQ=255)

                                                                                                                                  gCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGct                     <  1:1337016/150‑1 (MQ=255)

                                                                                                                                   cAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGctc                    <  2:1880143/150‑1 (MQ=255)

                                                                                                                                    aaTGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGctct                   >  1:4463381/1‑150 (MQ=255)

                                                                                                                                    aaTGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGctct                   >  1:1995148/1‑150 (MQ=255)

                                                                                                                                      tGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGctctct                 <  1:4097112/150‑1 (MQ=255)

                                                                                                                                      tGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGctctct                 >  2:512699/1‑150 (MQ=255)

                                                                                                                                       gTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTa                <  2:4761612/150‑1 (MQ=255)

                                                                                                                                        tACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTAt               >  2:1545531/1‑150 (MQ=255)

                                                                                                                                            cAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTg           <  1:4825758/150‑1 (MQ=255)

                                                                                                                                               cTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCg        <  1:4361700/150‑1 (MQ=255)

                                                                                                                                                tAACCTCCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGt       >  1:3263566/1‑150 (MQ=255)

                                                                                                                                                tAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGt       >  1:2277599/1‑150 (MQ=255)

                                                                                                                                                 aaCCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGtt      >  1:972815/1‑150 (MQ=255)

                                                                                                                                                   ccaccTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTAc    >  2:4352765/1‑150 (MQ=255)

                                                                                                                                                   ccaccTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTAc    <  1:2258071/150‑1 (MQ=255)

GTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGC  >  CP009273/1932880‑1933178

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