BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,888,584	G→T	64.2%	I185I (ATC→ATA)	ygcQ ←	putative flavoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,888,584	0	G	T	64.2%	288.6 / inf	315	I185I (ATC→ATA)	ygcQ	putative flavoprotein
Reads supporting (aligned to +/- strand): ref base G (48/65); new base T (98/104); total (146/169)
Fisher's exact test for biased strand distribution p-value = 3.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 316 total aligned reads displayed.

CTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGACAATGTTCAGTTGTTGCATTCC  >  CP009273/2888436‑2888733

ctgcGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGa                                                                                                                                                      <  2:1544/150‑1 (MQ=255)

 tgcGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGGTCATCACTCGCGCCCGACTGTCGCCCACTTCGGCTCCCAGCTTTTCTGCCCGCATAGCTAt                                                                                                                                                     >  1:36292/1‑150 (MQ=255)

    gcCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTc                                                                                                                                                  >  1:3797667/1‑150 (MQ=255)

    gcCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTc                                                                                                                                                  <  2:369044/150‑1 (MQ=255)

     ccAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTcc                                                                                                                                                 >  1:1763760/1‑150 (MQ=255)

      cAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCt                                                                                                                                                <  1:1323551/150‑1 (MQ=255)

      cAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCt                                                                                                                                                <  1:3368320/150‑1 (MQ=255)

      cAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCt                                                                                                                                                <  1:4780867/150‑1 (MQ=255)

       aGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTg                                                                                                                                               >  2:3877846/1‑150 (MQ=255)

        gAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTgg                                                                                                                                              <  2:2773616/150‑1 (MQ=255)

        gAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTgg                                                                                                                                              >  2:1430046/1‑150 (MQ=255)

         aaGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGt                                                                                                                                             >  1:2560615/1‑150 (MQ=255)

         aaGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGt                                                                                                                                             <  2:2976541/150‑1 (MQ=255)

          aGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGtt                                                                                                                                            >  1:4866009/1‑150 (MQ=255)

               caacaacaaTACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTg                                                                                                                                       <  2:1071966/150‑1 (MQ=255)

                  caacaaTACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCct                                                                                                                                    <  1:2176338/150‑1 (MQ=255)

                     caaTACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCtctc                                                                                                                                 >  2:4648950/1‑150 (MQ=255)

                     caaTACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCtctc                                                                                                                                 <  2:2121696/150‑1 (MQ=255)

                      aaTACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTcc                                                                                                                                <  1:3944236/150‑1 (MQ=255)

                      aaTACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTcc                                                                                                                                <  1:1238501/150‑1 (MQ=255)

                       aTACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCg                                                                                                                               >  1:597252/1‑150 (MQ=255)

                         acacacCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGcc                                                                                                                             <  1:3865530/150‑1 (MQ=255)

                              ccTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGcc                                                                                                                        <  2:245565/150‑1 (MQ=255)

                                tCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGcccc                                                                                                                      >  2:1991078/1‑150 (MQ=255)

                                  aGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCccac                                                                                                                    >  2:1462390/1‑150 (MQ=255)

                                    gTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCccacca                                                                                                                  <  2:702535/150‑1 (MQ=255)

                                       ccAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCccaccacca                                                                                                               >  1:2489200/1‑150 (MQ=255)

                                        cagcagATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAg                                                                                                              <  2:3296923/150‑1 (MQ=255)

                                         agcagATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGc                                                                                                             >  2:2374185/1‑150 (MQ=255)

                                             gATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCAccc                                                                                                         <  1:2574597/150‑1 (MQ=255)

                                              aTGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCg                                                                                                        <  1:1076251/150‑1 (MQ=255)

                                              aTGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCg                                                                                                        <  1:4477100/150‑1 (MQ=255)

                                               tGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGt                                                                                                       <  2:3514880/150‑1 (MQ=255)

                                                 tCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCt                                                                                                     <  2:4371098/150‑1 (MQ=255)

                                                  ccGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTg                                                                                                    <  2:4470462/150‑1 (MQ=255)

                                                  ccGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTg                                                                                                    >  1:314552/1‑150 (MQ=255)

                                                    ggAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGc                                                                                                  >  1:4256357/1‑150 (MQ=255)

                                                      aaaTACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCt                                                                                                >  2:4259316/1‑150 (MQ=255)

                                                      aaaTACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCt                                                                                                <  1:1301661/150‑1 (MQ=255)

                                                       aaTACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTc                                                                                               <  1:141119/150‑1 (MQ=255)

                                                       aaTACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTGTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTc                                                                                               <  2:4061867/150‑1 (MQ=255)

                                                        aTACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCa                                                                                              >  1:4346604/1‑150 (MQ=255)

                                                        aTACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCa                                                                                              <  2:3167867/150‑1 (MQ=255)

                                                          aCCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGc                                                                                            >  1:2327041/1‑150 (MQ=255)

                                                             aaTCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTaa                                                                                         <  1:747536/150‑1 (MQ=255)

                                                             aaTCACTTGTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTaa                                                                                         <  2:4117417/150‑1 (MQ=255)

                                                               tCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATg                                                                                       >  2:544386/1‑150 (MQ=255)

                                                                cACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATgg                                                                                      <  2:3955505/150‑1 (MQ=255)

                                                                  cTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGt                                                                                    >  2:4842720/1‑150 (MQ=255)

                                                                   tttttCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTc                                                                                   <  2:1374772/150‑1 (MQ=255)

                                                                    ttttCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTcg                                                                                  <  1:2718625/150‑1 (MQ=255)

                                                                     tttCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTcgc                                                                                 >  1:2292549/1‑150 (MQ=255)

                                                                         aGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGt                                                                             >  1:935925/1‑150 (MQ=255)

                                                                         aGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGt                                                                             >  1:4409801/1‑150 (MQ=255)

                                                                         aGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGt                                                                             >  2:1654585/1‑150 (MQ=255)

                                                                         aGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGt                                                                             <  2:4001864/150‑1 (MQ=255)

                                                                            aTCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAAc                                                                          <  1:3553559/150‑1 (MQ=255)

                                                                             tCAAAGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACa                                                                         >  2:3645020/1‑150 (MQ=255)

                                                                              cAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTCACAt                                                                        >  1:4709151/1‑150 (MQ=255)

                                                                              cAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACAt                                                                        >  2:3213319/1‑150 (MQ=255)

                                                                                 cGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACAtttt                                                                     >  2:856672/1‑150 (MQ=255)

                                                                                 cGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACAtttt                                                                     <  2:2818026/150‑1 (MQ=255)

                                                                                  gCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACAttttt                                                                    <  1:4534297/150‑1 (MQ=255)

                                                                                  gCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACAttttt                                                                    >  2:1208666/1‑150 (MQ=255)

                                                                                    cACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCa                                                                  >  1:1470654/1‑150 (MQ=255)

                                                                                     aCCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAg                                                                 <  2:1018663/150‑1 (MQ=255)

                                                                                     aCCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAg                                                                 >  1:1013018/1‑150 (MQ=255)

                                                                                      ccGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGa                                                                <  1:1800004/150‑1 (MQ=255)

                                                                                        gTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATc                                                              <  1:3515586/150‑1 (MQ=255)

                                                                                          tcatcaCTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCct                                                            >  2:3305025/1‑150 (MQ=255)

                                                                                            atcaCTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCctct                                                          <  1:1991078/150‑1 (MQ=255)

                                                                                              caCTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCtgt                                                        <  2:3332273/150‑1 (MQ=255)

                                                                                               aCTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCtgtg                                                       <  2:903538/150‑1 (MQ=255)

                                                                                               aCTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCtgtg                                                       >  2:1760136/1‑150 (MQ=255)

                                                                                                cTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGc                                                      <  2:1611131/150‑1 (MQ=255)

                                                                                                cTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGc                                                      <  2:4023396/150‑1 (MQ=255)

                                                                                                    cgcCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAAc                                                  >  1:2782523/1‑150 (MQ=255)

                                                                                                            tGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGt                                          >  2:1916159/1‑150 (MQ=255)

                                                                                                              tAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTcc                                        >  2:3388900/1‑150 (MQ=255)

                                                                                                               aGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCg                                       <  2:2291546/150‑1 (MQ=255)

                                                                                                               aGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCg                                       >  2:2055327/1‑150 (MQ=255)

                                                                                                                gCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCgg                                      >  1:2672591/1‑150 (MQ=255)

                                                                                                                 cccACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGa                                     <  2:1765185/150‑1 (MQ=255)

                                                                                                                 cccACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGa                                     >  1:4459646/1‑150 (MQ=255)

                                                                                                                  ccACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGaa                                    >  1:505438/1‑150 (MQ=255)

                                                                                                                  ccACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGaa                                    >  2:4172320/1‑150 (MQ=255)

                                                                                                                       tCGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCa                               <  1:2152754/150‑1 (MQ=255)

                                                                                                                        cGGCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCAc                              <  1:3304913/150‑1 (MQ=255)

                                                                                                                         ggCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCAcc                             >  2:2667185/1‑150 (MQ=255)

                                                                                                                          gCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCAccc                            <  1:2250582/150‑1 (MQ=255)

                                                                                                                          gCTCCCAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCAccc                            <  1:1462390/150‑1 (MQ=255)

                                                                                                                              ccAGCTTTTCTGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGa                        <  1:1552964/150‑1 (MQ=255)

                                                                                                                                aGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACGTTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGACa                      >  1:137101/1‑150 (MQ=255)

                                                                                                                                  cTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGACAAt                    <  1:4648950/150‑1 (MQ=255)

                                                                                                                                        tGCCAGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGACAATGTTCAg              <  1:1797689/150‑1 (MQ=255)

                                                                                                                                            aGCATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGACAATGTTCAgttgt          >  1:4680008/1‑150 (MQ=255)

                                                                                                                                              cATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGACAATGTTCAGttgttg        >  2:4342940/1‑150 (MQ=255)

                                                                                                                                              cATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGACAATGTTCAGttgttg        >  2:1379195/1‑150 (MQ=255)

                                                                                                                                              cATAGCTATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGACAATGTTCAGttgttg        >  2:3940187/1‑150 (MQ=255)

                                                                                                                                                   cGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGACAATGTTCAGTTGTTGCATTc   <  1:4101571/150‑1 (MQ=255)

CTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCCGACTGTAGCCCACTTCGGCTCCCAGCTTTTCTGCCAGCATAGCGATTTCCTGGTTATCTGCCTCTCCGCCCTGCCCCACCACCAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGACAATGTTCAGTTGTTGCATTCC  >  CP009273/2888436‑2888733

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