BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,929,196	C→A	6.8%	T85K (ACG→AAG)	fucI →	L‑fucose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,929,196	0	C	A	6.8%	1194.4 / 45.0	323	T85K (ACG→AAG)	fucI	L‑fucose isomerase
Reads supporting (aligned to +/- strand): ref base C (165/136); new base A (14/8); total (179/144)
Fisher's exact test for biased strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 329 total aligned reads displayed.

AGCTACGGCCGCACTGCTGACCGAGAAACTGCGCCATGCCTGCGGAGCTGCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAG  >  CP009273/2929047‑2929345

   tACGGCCGCACTGCTGACCGAGAAACTGCGCCATGCCTGCGGAGCTGCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGt                                                                                                                                                    >  1:715766/1‑150 (MQ=255)

   tACGGCCGCACTGCTGACCGAGAAACTGCGCCATGCCTGCGGAGCTGCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGt                                                                                                                                                    >  2:1112170/1‑150 (MQ=255)

   tACGGCCGCACTGCTGACCGAGAAACTGCGCCATGCCTGCGGAGCTGCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGt                                                                                                                                                    <  2:287477/150‑1 (MQ=255)

      ggCCGCACTGCTGACCGAGAAACTGCGCCATGCCTGCGGAGCTGCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAAc                                                                                                                                                 >  2:661922/1‑150 (MQ=255)

      ggCCGCACTGCTGACCGAGAAACTGCGCCATGCCTGCGGAGCTGCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAAc                                                                                                                                                 <  1:3652408/150‑1 (MQ=255)

      ggCCGCACTGCTGACCGAGAAACTGCGCCATGCCTGCGGAGCTGCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAAc                                                                                                                                                 >  1:294314/1‑150 (MQ=255)

           cACTGCTGACCGAGAAACTGCGCCATGCCTGCGGAGCTGCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCtt                                                                                                                                            <  2:485186/150‑1 (MQ=255)

                ctgACCGAGAAACTGCGCCATGCCTGCGGAGCTGCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTgg                                                                                                                                       <  1:3570813/150‑1 (MQ=255)

                ctgACCGAGAAACTGCGCCATGCCTGCGGAGCTGCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTgg                                                                                                                                       >  2:1894874/1‑150 (MQ=255)

                 tgACCGAGAAACTGCGCCATGCCTGCGGAGCTGCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGt                                                                                                                                      <  2:3318242/150‑1 (MQ=255)

                 tgACCGAGAAACTGCGCCATGCCTGCGGAGCTGCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGt                                                                                                                                      >  1:4620861/1‑150 (MQ=255)

                     cGAGAAACTGCGCCATGCCTGCGGAGCTGCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTa                                                                                                                                  <  1:1854530/150‑1 (MQ=255)

                     cGAGAAACTGCGCCATGCCTGCGGAGCTGCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTa                                                                                                                                  <  2:1215717/150‑1 (MQ=255)

                      gagaAACTGCGCCATGCCTGCGGAGCTGCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTAt                                                                                                                                 <  2:3822843/150‑1 (MQ=255)

                         aaaCTGCGCCATGCCTGCGGAGCTGCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGc                                                                                                                              >  2:263728/1‑150 (MQ=255)

                          aaCTGCGCCATGCCTGCGGAGCTGCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCa                                                                                                                             >  2:3791263/1‑150 (MQ=255)

                          aaCTGCGCCATGCCTGCGGAGCTGCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTAAGGTACCGCCTTGCTGGTGCTATGGCa                                                                                                                             >  2:3369574/1‑150 (MQ=255)

                            cTGCGCCATGCCTGCGGAGCTGCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGt                                                                                                                           >  1:4899061/1‑150 (MQ=255)

                            cTGCGCCATGCCTGCGGAGCTGCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGt                                                                                                                           <  1:4339278/150‑1 (MQ=255)

                              gcgcCATGCCTGCGGAGCTGCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGa                                                                                                                         <  2:2028661/150‑1 (MQ=255)

                              gcgcCATGCCTGCGGAGCTGCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGa                                                                                                                         <  2:5227469/150‑1 (MQ=255)

                              gcgcCATGCCTGCGGAGCTGCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGa                                                                                                                         <  1:1058291/150‑1 (MQ=255)

                               cgcCATGCCTGCGGAGCTGCCGTCGAGTGTGTCATTTCCGTTACCTGTATCGCGGGTATGGATGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAAGGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGaa                                                                                                                        >  2:823149/1‑150 (MQ=255)

                                 ccATGCCTGCGGAGCTGCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATTTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAAc                                                                                                                      >  2:2619842/1‑150 (MQ=255)

                                 ccATGCCTGCGGAGCTGCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAAc                                                                                                                      <  1:2070118/150‑1 (MQ=255)

                                     gCCTGCGGAGCTGCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATc                                                                                                                  <  2:143649/150‑1 (MQ=255)

                                     gCCTGCGGAGCTGCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATc                                                                                                                  <  2:2679518/150‑1 (MQ=255)

                                       cTGCGGAGCTGCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGa                                                                                                                >  2:3697589/1‑150 (MQ=255)

                                         gCGGAGCTGCCGTCGAGTGTGTCATGTCCGATCCCGGTATCGCGGGAATGGCTGAAGCCGCTGCTTGCGAAGAACAATTCAGCCGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCAACGACa                                                                                                              <  1:4238963/150‑1 (MQ=255)

                                             aGCTGCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGa                                                                                                          >  2:2400929/1‑150 (MQ=255)

                                                tGCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCAGCTTGCGAAGAAAAATTCAGCAGACAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATcc                                                                                                       >  2:1415402/1‑150 (MQ=255)

                                                 gCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCa                                                                                                      >  1:4519033/1‑150 (MQ=255)

                                                  ccGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCaa                                                                                                     <  1:2873411/150‑1 (MQ=255)

                                                  ccGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCaa                                                                                                     <  2:3960101/150‑1 (MQ=255)

                                                      cGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCAATACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAAcccg                                                                                                 >  2:4400499/1‑150 (MQ=255)

                                                         gtgtgtCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAAcccgccc                                                                                              <  2:3412675/150‑1 (MQ=255)

                                                         gtgtgtCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAAcccgccc                                                                                              >  1:3450813/1‑150 (MQ=255)

                                                         gtgtgtCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTAAGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAAcccgccc                                                                                              <  2:2125775/150‑1 (MQ=255)

                                                          tgtgtCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAAcccgcccg                                                                                             <  1:609968/150‑1 (MQ=255)

                                                                aTTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGcc                                                                                       >  1:4852931/1‑150 (MQ=255)

                                                                aTTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGcc                                                                                       <  1:1305255/150‑1 (MQ=255)

                                                                aTTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACAAGGATCCAACCCGCCCGAAGGcc                                                                                       >  2:4944448/1‑150 (MQ=255)

                                                                aTTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTAAGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGcc                                                                                       >  2:402892/1‑150 (MQ=255)

                                                                  ttCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCAt                                                                                     <  1:1889314/150‑1 (MQ=255)

                                                                    ccGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCAttt                                                                                   >  2:2168871/1‑150 (MQ=255)

                                                                    ccGATACCAGTATCGCGTCTAACGCAGAAGCAGCCGCGAGCGAAGAAAAATTCAGCAGTCAGAATGTACGCCTCACCAATACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAAACCGCCCGAAGGCCAttt                                                                                   <  2:3185323/150‑1 (MQ=255)

                                                                    ccGAGACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCAttt                                                                                   <  2:4727309/150‑1 (MQ=255)

                                                                     cGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAACCCATCGACATGGATCCAACACGCCCGAAGGCCATTTg                                                                                  >  2:3703207/1‑150 (MQ=255)

                                                                     cGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTg                                                                                  >  1:2395775/1‑150 (MQ=255)

                                                                       aTACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTggg                                                                                >  1:3992353/1‑150 (MQ=255)

                                                                       aTACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTggg                                                                                >  1:3978921/1‑150 (MQ=255)

                                                                         aCCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGc                                                                              <  1:1171238/150‑1 (MQ=255)

                                                                          ccTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCt                                                                             <  2:1644577/150‑1 (MQ=255)

                                                                           cTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCtt                                                                            >  1:5094477/1‑150 (MQ=255)

                                                                             gTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTa                                                                          <  2:1973443/150‑1 (MQ=255)

                                                                             gTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTa                                                                          >  1:4037311/1‑150 (MQ=255)

                                                                             gTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTa                                                                          >  1:4181081/1‑150 (MQ=255)

                                                                              tATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTaa                                                                         <  1:3091923/150‑1 (MQ=255)

                                                                               aTCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAAc                                                                        <  2:85900/150‑1 (MQ=255)

                                                                               aTCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCAGTTGGGGCTTTAAc                                                                        >  2:1357019/1‑150 (MQ=255)

                                                                                   cgGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCa                                                                    <  2:2093507/150‑1 (MQ=255)

                                                                                    gggTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCAc                                                                   >  1:1956294/1‑150 (MQ=255)

                                                                                     ggTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACt                                                                  <  1:327408/150‑1 (MQ=255)

                                                                                          ggCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACg                                                             <  1:4411753/150‑1 (MQ=255)

                                                                                          ggCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACg                                                             <  2:4608625/150‑1 (MQ=255)

                                                                                            cTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGcc                                                           >  1:4387162/1‑150 (MQ=255)

                                                                                              gAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGcccc                                                         <  1:211015/150‑1 (MQ=255)

                                                                                              gAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGcccc                                                         <  1:1472076/150‑1 (MQ=255)

                                                                                              gAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGcccc                                                         >  2:3103040/1‑150 (MQ=255)

                                                                                                        ctTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGttt                                               >  2:3446228/1‑150 (MQ=255)

                                                                                                        ctTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTCGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGttt                                               >  1:2196934/1‑150 (MQ=255)

                                                                                                            cGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCt                                           <  2:3571063/150‑1 (MQ=255)

                                                                                                             gaagaaAAATTCAGCAGTCAGAATGTATGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTg                                          <  2:2596066/150‑1 (MQ=255)

                                                                                                             gaagaaAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTg                                          >  2:4542715/1‑150 (MQ=255)

                                                                                                              aagaaAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTgg                                         <  1:3259809/150‑1 (MQ=255)

                                                                                                              aagaaAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTgg                                         >  2:1043820/1‑150 (MQ=255)

                                                                                                              aagaaAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTgg                                         <  1:1826643/150‑1 (MQ=255)

                                                                                                                gaaAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCa                                       >  1:1641403/1‑150 (MQ=255)

                                                                                                                gaaAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCGTGCTGGTGCGATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCa                                       <  2:1776760/150‑1 (MQ=255)

                                                                                                                 aaaaaTTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCAg                                      >  2:2374381/1‑150 (MQ=255)

                                                                                                                   aaaTTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCg                                    <  1:4971945/150‑1 (MQ=255)

                                                                                                                        cagcagTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGctct                               >  2:812529/1‑150 (MQ=255)

                                                                                                                           cagTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGc                            >  1:4751856/1‑150 (MQ=255)

                                                                                                                           cagTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGc                            <  2:3083478/150‑1 (MQ=255)

                                                                                                                                aGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTc                       >  1:100678/1‑150 (MQ=255)

                                                                                                                                aGAATGTAGGCCTCACCATTAAGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTc                       >  1:3870047/1‑150 (MQ=255)

                                                                                                                                 taaTGTAGGACTCACAATTACGGTTACTCCTTGCTGGTGATATGGCAGTGAATCCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTca                      <  2:2772825/149‑1 (MQ=255)

                                                                                                                                 gAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTca                      >  2:4974938/1‑150 (MQ=255)

                                                                                                                                      tAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGcc                 >  2:4662885/1‑150 (MQ=255)

                                                                                                                                      tAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGcc                 >  1:2420527/1‑150 (MQ=255)

                                                                                                                                       aGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCa                >  1:855090/1‑150 (MQ=255)

                                                                                                                                       aGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCa                <  1:5183531/150‑1 (MQ=255)

                                                                                                                                           cTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGaaa            >  1:985378/1‑150 (MQ=255)

                                                                                                                                             cACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAgg          >  1:4435591/1‑150 (MQ=255)

                                                                                                                                             cACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAgg          <  2:2201545/150‑1 (MQ=255)

                                                                                                                                              aCCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCTACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGTAAGGc         <  1:1230492/150‑1 (MQ=255)

                                                                                                                                                cATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCAt       <  1:3791263/150‑1 (MQ=255)

                                                                                                                                                cATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCAt       <  1:4733568/150‑1 (MQ=255)

                                                                                                                                                 aTTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATc      >  1:3167805/1‑150 (MQ=255)

                                                                                                                                                    aCGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCa   <  2:2912902/150‑1 (MQ=255)

AGCTACGGCCGCACTGCTGACCGAGAAACTGCGCCATGCCTGCGGAGCTGCCGTCGAGTGTGTCATTTCCGATACCTGTATCGCGGGTATGGCTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAG  >  CP009273/2929047‑2929345

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