BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,224,419	C→A	5.8%	R725S (CGT→AGT)	ygjK →	alpha‑glucosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,224,419	0	C	A	5.8%	1282.2 / 34.5	326	R725S (CGT→AGT)	ygjK	alpha‑glucosidase
Reads supporting (aligned to +/- strand): ref base C (138/169); new base A (7/12); total (145/181)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 326 total aligned reads displayed.

TCGATCCTAAAGAGTTCAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATAT  >  CP009273/3224273‑3224567

tCGATCCTAAAGAGTTCAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGtg                                                                                                                                                   >  2:128974/1‑150 (MQ=255)

 cGATCCTAAAGAGTTCAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGtga                                                                                                                                                  >  1:2448765/1‑150 (MQ=255)

  gATCCTAAAGAGTTCAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGtgat                                                                                                                                                 >  1:3486236/1‑150 (MQ=255)

   aTCCTAAAGAGTTCAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGtgatg                                                                                                                                                >  2:3268935/1‑150 (MQ=255)

    tCCTAAAGAGTTCAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGtgatga                                                                                                                                               >  1:1042294/1‑150 (MQ=255)

         aaGAGTTCAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGccc                                                                                                                                          <  2:4862915/150‑1 (MQ=255)

          agagTTCAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCt                                                                                                                                         >  2:535634/1‑150 (MQ=255)

          agagTTCAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCt                                                                                                                                         <  2:3406295/150‑1 (MQ=255)

             gTTCAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATAGTGATGATGCCCTGaa                                                                                                                                      >  2:2678078/1‑150 (MQ=255)

              ttCAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAg                                                                                                                                     <  1:4802940/150‑1 (MQ=255)

               tCAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGc                                                                                                                                    >  1:625832/1‑150 (MQ=255)

               tCAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGc                                                                                                                                    <  2:3599478/150‑1 (MQ=255)

                cAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCt                                                                                                                                   <  1:4145354/150‑1 (MQ=255)

                 aaCACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTg                                                                                                                                  >  2:3683584/1‑150 (MQ=255)

                  acacCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTgg                                                                                                                                 >  2:940776/1‑150 (MQ=255)

                  acacCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTgg                                                                                                                                 <  1:300250/150‑1 (MQ=255)

                   cacCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGc                                                                                                                                >  2:668620/1‑150 (MQ=255)

                    acCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCg                                                                                                                               >  2:816022/1‑150 (MQ=255)

                       tttGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGAt                                                                                                                            <  1:315928/150‑1 (MQ=255)

                        ttGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATa                                                                                                                           <  2:276393/150‑1 (MQ=255)

                         tggcccGCTGGGAACGGCGGCGTTAACCCATCCGGCTTTTGGCGCTGATCTCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATAc                                                                                                                          <  2:2724476/147‑1 (MQ=255)

                         tGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATAc                                                                                                                          <  2:4234071/150‑1 (MQ=255)

                         tGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATAc                                                                                                                          <  2:469245/150‑1 (MQ=255)

                           tCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATAGTGATGATGCCCTGAAGCTGGCGGATACGt                                                                                                                        >  1:98457/1‑150 (MQ=255)

                            cccGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGtt                                                                                                                       >  2:2150889/1‑150 (MQ=255)

                              cGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGttct                                                                                                                     <  2:1511575/150‑1 (MQ=255)

                                 tGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTcc                                                                                                                  >  1:5117763/1‑150 (MQ=255)

                                  gggAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCg                                                                                                                 <  2:3375792/150‑1 (MQ=255)

                                      aCGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCAc                                                                                                             >  2:3037776/1‑150 (MQ=255)

                                       cggcggCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACg                                                                                                            >  2:4150418/1‑150 (MQ=255)

                                           ggcgTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCaa                                                                                                        >  2:2916556/1‑150 (MQ=255)

                                             cgTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGTGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAg                                                                                                      <  2:2173837/150‑1 (MQ=255)

                                               ttAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAggg                                                                                                    <  1:4998140/150‑1 (MQ=255)

                                                 aaCCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGtt                                                                                                  <  2:4954926/150‑1 (MQ=255)

                                                   ccAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTaa                                                                                                >  2:1680808/1‑150 (MQ=255)

                                                     aaTCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAAcc                                                                                              >  2:1153817/1‑150 (MQ=255)

                                                     aaTCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAAcc                                                                                              >  1:2695340/1‑150 (MQ=255)

                                                     aaTCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAAcc                                                                                              >  1:1498837/1‑150 (MQ=255)

                                                      aTCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAAccg                                                                                             <  2:2879413/150‑1 (MQ=255)

                                                       tCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAAccgc                                                                                            >  1:471813/1‑150 (MQ=255)

                                                         cGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAAccgccg                                                                                          <  2:3463204/150‑1 (MQ=255)

                                                         cGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAAccgccg                                                                                          >  2:1247340/1‑150 (MQ=255)

                                                         cGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAAccgccg                                                                                          <  2:4297750/150‑1 (MQ=255)

                                                          ggCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGa                                                                                         <  2:1107069/150‑1 (MQ=255)

                                                          ggCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGa                                                                                         >  1:4001065/1‑150 (MQ=255)

                                                          ggCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGa                                                                                         <  1:4045887/150‑1 (MQ=255)

                                                            cTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATg                                                                                       <  1:2596291/150‑1 (MQ=255)

                                                              tttGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTCACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGa                                                                                     >  2:2070516/1‑150 (MQ=255)

                                                               ttGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGAc                                                                                    <  1:3242426/150‑1 (MQ=255)

                                                               ttGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGAc                                                                                    <  1:3242277/150‑1 (MQ=255)

                                                                 ggCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGTAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCa                                                                                  <  1:4607527/150‑1 (MQ=255)

                                                                 ggCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCa                                                                                  >  1:298278/1‑150 (MQ=255)

                                                                       gATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAg                                                                            >  1:2227705/1‑150 (MQ=255)

                                                                          atCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGaa                                                                         <  1:1116628/150‑1 (MQ=255)

                                                                          atCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGaa                                                                         <  1:5028939/150‑1 (MQ=255)

                                                                             tACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAAt                                                                      <  1:769489/150‑1 (MQ=255)

                                                                              aCTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAAtt                                                                     <  1:871035/150‑1 (MQ=255)

                                                                               cTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTa                                                                    >  1:4243735/1‑150 (MQ=255)

                                                                                     tGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAAccc                                                              >  2:2259136/1‑150 (MQ=255)

                                                                                          gcgTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGa                                                         <  2:2429729/150‑1 (MQ=255)

                                                                                            gTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACa                                                       >  2:3122555/1‑150 (MQ=255)

                                                                                             tATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAg                                                      <  1:4963686/150‑1 (MQ=255)

                                                                                             tATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAg                                                      <  2:2003357/150‑1 (MQ=255)

                                                                                             tATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAg                                                      <  1:806920/150‑1 (MQ=255)

                                                                                                    ggATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGcaca                                               >  1:4663224/1‑150 (MQ=255)

                                                                                                      aTCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGc                                             >  2:5083364/1‑150 (MQ=255)

                                                                                                      aTCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGc                                             >  2:2272860/1‑150 (MQ=255)

                                                                                                       tCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCa                                            >  1:1524729/1‑150 (MQ=255)

                                                                                                       tCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCa                                            >  1:1699077/1‑150 (MQ=255)

                                                                                                       tCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCa                                            <  2:3809668/150‑1 (MQ=255)

                                                                                                        cAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCaa                                           <  1:2132930/150‑1 (MQ=255)

                                                                                                         aGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAg                                          <  1:3782423/150‑1 (MQ=255)

                                                                                                            tCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAACATTACAACCCGCTGACCGGCGCACAGCAAGgcg                                       >  2:2816475/1‑150 (MQ=255)

                                                                                                              tGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCa                                     <  1:254775/150‑1 (MQ=255)

                                                                                                               ggTTTGGTCTGAAAGGGATGGGGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCAc                                    <  1:1795505/150‑1 (MQ=255)

                                                                                                                  ttGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCaa                                 <  2:3756934/150‑1 (MQ=255)

                                                                                                                       cTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTtc                            >  1:4303557/1‑150 (MQ=255)

                                                                                                                        tGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTtct                           >  2:2538431/1‑150 (MQ=255)

                                                                                                                          aaaGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTcc                         <  1:1489925/150‑1 (MQ=255)

                                                                                                                          aaaGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTcc                         >  2:2593827/1‑150 (MQ=255)

                                                                                                                          aaaGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTcc                         >  2:4269165/1‑150 (MQ=255)

                                                                                                                          aaaGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTcc                         <  1:3428857/150‑1 (MQ=255)

                                                                                                                          aaaGGGATGGAGCGTTACGGTTATAGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTcc                         <  2:1947896/150‑1 (MQ=255)

                                                                                                                            aGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTg                       >  2:4817967/1‑150 (MQ=255)

                                                                                                                              ggATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGa                     <  2:4516645/150‑1 (MQ=255)

                                                                                                                                  ggAGCGTTACGGTTATCGTGATGATGCCCTGAATCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGc                 <  2:3874081/150‑1 (MQ=255)

                                                                                                                                   gAGCGTTACGGTTATAGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGcc                <  1:672181/150‑1 (MQ=255)

                                                                                                                                       gTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTATCCTGGAGTGCCgcgc            >  2:1624140/1‑150 (MQ=255)

                                                                                                                                       gTTACGGTTATAGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCgcgc            >  1:2064655/1‑150 (MQ=255)

                                                                                                                                         tACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCAt          <  1:4877634/150‑1 (MQ=255)

                                                                                                                                         tACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCAt          <  1:1560938/150‑1 (MQ=255)

                                                                                                                                         tACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCAt          <  1:2819990/150‑1 (MQ=255)

                                                                                                                                         tACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCAt          >  1:4176734/1‑150 (MQ=255)

                                                                                                                                           cGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCAttt        >  1:3029338/1‑150 (MQ=255)

                                                                                                                                           cGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCAttt        <  2:1852941/150‑1 (MQ=255)

                                                                                                                                           cGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCAttt        <  2:2448866/150‑1 (MQ=255)

                                                                                                                                           cGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCAttt        <  2:862956/150‑1 (MQ=255)

                                                                                                                                            ggTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTg       <  1:369222/150‑1 (MQ=255)

                                                                                                                                               tATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGtat    <  1:3321633/150‑1 (MQ=255)

                                                                                                                                               tATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGtat    <  1:4524311/150‑1 (MQ=255)

TCGATCCTAAAGAGTTCAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATAT  >  CP009273/3224273‑3224567

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