BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,432,676	C→A	9.3%	E65* (GAA→TAA)	yhdN ←	DUF1992 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,432,676	0	C	A	9.3%	1005.3 / 50.3	269	E65* (GAA→TAA)	yhdN	DUF1992 family protein
Reads supporting (aligned to +/- strand): ref base C (120/124); new base A (12/13); total (132/137)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

Only 100 of 280 total aligned reads displayed.

CTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCG  >  CP009273/3432527‑3432825

  tGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTtctc                                                                                                                                                     >  1:4454591/1‑150 (MQ=255)

   gTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTtctct                                                                                                                                                    >  2:2477336/1‑150 (MQ=255)

   gTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTtctct                                                                                                                                                    >  2:4383002/1‑150 (MQ=255)

     cAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTcc                                                                                                                                                  >  2:3962451/1‑150 (MQ=255)

       gCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGt                                                                                                                                                >  1:1945498/1‑150 (MQ=255)

         aTAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTg                                                                                                                                              >  1:1022440/1‑150 (MQ=255)

         aTAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTg                                                                                                                                              >  2:3569095/1‑150 (MQ=255)

          tAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGc                                                                                                                                             >  1:2573644/1‑150 (MQ=255)

             tCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCa                                                                                                                                          >  1:3126809/1‑150 (MQ=255)

                  gcgTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACTGATACCTTTGAGAATATCCCGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTc                                                                                                                                     >  2:34093/1‑150 (MQ=255)

                   cgTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCt                                                                                                                                    <  1:62033/150‑1 (MQ=255)

                     tAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTgg                                                                                                                                  >  2:795716/1‑150 (MQ=255)

                     tAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTgg                                                                                                                                  >  2:795728/1‑150 (MQ=255)

                         aaaTCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGt                                                                                                                              >  2:4549437/1‑150 (MQ=255)

                           aTCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTaa                                                                                                                            <  1:3652260/150‑1 (MQ=255)

                           aTCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTATCTCCGTTGCTCAAGTTCTGGCGGTaa                                                                                                                            <  2:4091005/150‑1 (MQ=255)

                            tCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAg                                                                                                                           >  2:812465/1‑150 (MQ=255)

                             cGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGc                                                                                                                          <  2:187816/150‑1 (MQ=255)

                              ggTACTCAATCCATCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCGTGATTTTGCGGATCATCGTGTCGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCa                                                                                                                         <  2:2109603/150‑1 (MQ=255)

                              ggTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCa                                                                                                                         <  2:2338758/150‑1 (MQ=255)

                               gTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCaa                                                                                                                        <  2:462107/150‑1 (MQ=255)

                               gTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCaa                                                                                                                        <  2:3504485/150‑1 (MQ=255)

                               gTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCaa                                                                                                                        >  1:3321252/1‑150 (MQ=255)

                                 aCTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAAcc                                                                                                                      <  2:4633314/150‑1 (MQ=255)

                                 aCTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAAcc                                                                                                                      <  2:2677169/150‑1 (MQ=255)

                                  cTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCg                                                                                                                     >  1:2672551/1‑150 (MQ=255)

                                   tCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCgg                                                                                                                    >  2:682475/1‑150 (MQ=255)

                                    cAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTATTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGc                                                                                                                   >  2:1325043/1‑150 (MQ=255)

                                       tCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCAtt                                                                                                                <  2:572671/150‑1 (MQ=255)

                                        ccAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCAttc                                                                                                               >  2:1732001/1‑150 (MQ=255)

                                        ccAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCAttc                                                                                                               >  2:1019730/1‑150 (MQ=255)

                                           gCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCAttcttc                                                                                                            >  1:1000986/1‑150 (MQ=255)

                                           gCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCAttcttc                                                                                                            >  1:994939/1‑150 (MQ=255)

                                            cTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCa                                                                                                           <  2:4341915/150‑1 (MQ=255)

                                             ttGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTATCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAg                                                                                                          <  2:542549/150‑1 (MQ=255)

                                                  gCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGc                                                                                                     >  2:1758436/1‑150 (MQ=255)

                                                        tCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATAcc                                                                                               <  2:3699366/150‑1 (MQ=255)

                                                          aaTTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATAcccc                                                                                             <  2:3814574/150‑1 (MQ=255)

                                                           aTTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCg                                                                                            >  2:3736915/1‑150 (MQ=255)

                                                              ccAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCAc                                                                                         >  2:1140120/1‑150 (MQ=255)

                                                                 gTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTa                                                                                      <  1:1818431/150‑1 (MQ=255)

                                                                   aaTGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAAt                                                                                    >  1:1230247/1‑150 (MQ=255)

                                                                    aTGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAAtt                                                                                   >  2:1247744/1‑150 (MQ=255)

                                                                    aTGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTATCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAAtt                                                                                   <  1:3889699/150‑1 (MQ=255)

                                                                     tGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATACAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATACTTCAGCAAGCGATACCCCGCACGTAATTc                                                                                  >  1:555470/1‑150 (MQ=255)

                                                                       aCAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCg                                                                                >  2:2500669/1‑150 (MQ=255)

                                                                       aCAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCg                                                                                >  2:2252781/1‑150 (MQ=255)

                                                                       aCAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTATCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCg                                                                                <  2:2010578/150‑1 (MQ=255)

                                                                        cAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCgg                                                                               >  1:1683242/1‑150 (MQ=255)

                                                                         aaTCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGt                                                                              <  2:3109964/150‑1 (MQ=255)

                                                                         aaTCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGt                                                                              <  2:2783510/150‑1 (MQ=255)

                                                                         aaTCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGt                                                                              >  1:829727/1‑150 (MQ=255)

                                                                              aCGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCAc                                                                         >  2:32262/1‑150 (MQ=255)

                                                                              aCGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCAc                                                                         <  1:4290432/150‑1 (MQ=255)

                                                                                 gCTTACCTCTTGATATTGCGGATCATCGTGAAGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCTGTGGCACGTg                                                                      <  1:4140125/150‑1 (MQ=255)

                                                                                 gCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTg                                                                      <  2:71285/150‑1 (MQ=255)

                                                                                    aaCCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTgaga                                                                   <  1:2315827/150‑1 (MQ=255)

                                                                                    aaCCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTgaga                                                                   >  2:1202828/1‑150 (MQ=255)

                                                                                      ccTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAat                                                                 <  2:2898489/150‑1 (MQ=255)

                                                                                      ccTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTATCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAat                                                                 <  2:1578317/150‑1 (MQ=255)

                                                                                       ctctTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAatc                                                                <  1:3255677/150‑1 (MQ=255)

                                                                                         ctTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAatcat                                                              >  2:1618697/1‑150 (MQ=255)

                                                                                          ttGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAatcatc                                                             >  2:4209717/1‑150 (MQ=255)

                                                                                          ttGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAatcatc                                                             <  2:1245735/150‑1 (MQ=255)

                                                                                           tGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAatcatca                                                            >  1:3624496/1‑150 (MQ=255)

                                                                                            gATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAatcatcat                                                           >  1:1609384/1‑150 (MQ=255)

                                                                                              tatTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATcc                                                         <  1:221973/150‑1 (MQ=255)

                                                                                                   cGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTAt                                                    <  1:1186686/150‑1 (MQ=255)

                                                                                                      atcatcGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCaa                                                 <  1:757155/150‑1 (MQ=255)

                                                                                                        catcGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATg                                               <  1:2254196/150‑1 (MQ=255)

                                                                                                        catcGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATg                                               >  2:1621997/1‑150 (MQ=255)

                                                                                                           cGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGGGAGAATCATCATCCAGTATCAATGGtt                                            >  1:22640/1‑150 (MQ=255)

                                                                                                               aCGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTcgcc                                        <  1:1477646/150‑1 (MQ=255)

                                                                                                                  gATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCt                                     <  2:3361261/150‑1 (MQ=255)

                                                                                                                  gATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCt                                     <  1:4081411/150‑1 (MQ=255)

                                                                                                                       cTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTg                                <  1:2066014/150‑1 (MQ=255)

                                                                                                                        tttGAGAATATCCAGAAGCTGAATTGCTTATCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGc                               >  1:3343/1‑150 (MQ=255)

                                                                                                                         ttGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGcc                              >  1:699788/1‑150 (MQ=255)

                                                                                                                         ttGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGcc                              <  1:2643928/150‑1 (MQ=255)

                                                                                                                           gagaATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAg                            >  1:2352926/1‑150 (MQ=255)

                                                                                                                           gagaATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAg                            <  1:1136134/150‑1 (MQ=255)

                                                                                                                            agaATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAgg                           >  2:1957071/1‑150 (MQ=255)

                                                                                                                              aaTATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGtt                         >  2:1927661/1‑150 (MQ=255)

                                                                                                                               atatCCAGAAGCTGAATTGCTTATCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTa                        <  1:1436514/150‑1 (MQ=255)

                                                                                                                                 atCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATc                      <  2:1686004/150‑1 (MQ=255)

                                                                                                                                  tCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCa                     >  1:3285717/1‑150 (MQ=255)

                                                                                                                                  tCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCa                     >  1:2904648/1‑150 (MQ=255)

                                                                                                                                   ccAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCaa                    <  2:1384160/150‑1 (MQ=255)

                                                                                                                                    cAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCaaa                   >  1:3925588/1‑150 (MQ=255)

                                                                                                                                     aGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAAc                  <  2:3001002/150‑1 (MQ=255)

                                                                                                                                     aGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAAc                  >  2:4205926/1‑150 (MQ=255)

                                                                                                                                         gCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCAc              <  2:607370/150‑1 (MQ=255)

                                                                                                                                          cTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCAcc             <  1:4479044/150‑1 (MQ=255)

                                                                                                                                            gAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCtt           >  1:915498/1‑150 (MQ=255)

                                                                                                                                            gAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCtt           <  1:3034985/150‑1 (MQ=255)

                                                                                                                                             aaTTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCttt          >  2:3823346/1‑150 (MQ=255)

                                                                                                                                                  cTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCtt     >  1:3753123/1‑150 (MQ=255)

                                                                                                                                                   ttCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTg    >  2:3996175/1‑150 (MQ=255)

                                                                                                                                                    tctctcCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTgc   <  2:3528287/150‑1 (MQ=255)

                                                                                                                                                     ctctcCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTgcg  <  2:181214/150‑1 (MQ=255)

CTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCG  >  CP009273/3432527‑3432825

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