BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,504,963	A→T	6.0%	Q98L (CAG→CTG)	speB →	agmatinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,504,963	0	A	T	6.0%	599.8 / 17.3	184	Q98L (CAG→CTG)	speB	agmatinase
Reads supporting (aligned to +/- strand): ref base A (78/95); new base T (7/4); total (85/99)
Fisher's exact test for biased strand distribution p-value = 3.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

Only 100 of 187 total aligned reads displayed.

CCGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTGGCGCTGATCCATTTC  >  NC_002947/2504814‑2505111

ccGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCa                                                                                                                                                      <  1:1693792/150‑1 (MQ=255)

ccGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCCGTCCAGCAGGCCTGCGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCCACTAAGGCGACTGACCCTTCGACAGCGGCACGCCGCa                                                                                                                                                      >  2:3063627/1‑150 (MQ=255)

 cGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAg                                                                                                                                                     >  1:3008907/1‑150 (MQ=255)

  gCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACTGCGGCGTGCCGCTGt                                                                                                                                                    <  2:4796950/150‑1 (MQ=255)

   ccccGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCCGGCCTGGGCCCGCCACTGGCCGTGGGCGGTCGACCCGATCGACCACCTGGCGGTCATCGACTCTGGCGACTGTCCCGTCGACAGCGGCACGCCGCAGTc                                                                                                                                                   >  2:3331423/1‑150 (MQ=255)

   ccccGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTc                                                                                                                                                   <  2:706309/150‑1 (MQ=255)

    cccGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACTGCGGCGTGCCGCTGTCg                                                                                                                                                  <  2:779817/150‑1 (MQ=255)

    cccGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCg                                                                                                                                                  >  2:4312744/1‑150 (MQ=255)

      cGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACTGCGGCGTGCCGCTGTCgaa                                                                                                                                                <  2:2365619/150‑3 (MQ=255)

        gcgcgcgcTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTAc                                                                                                                                              <  1:4010111/150‑1 (MQ=255)

         cgcgcgcTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGG‑AACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCg                                                                                                                                            <  2:2610875/150‑1 (MQ=255)

            gcgcTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGa                                                                                                                                          <  1:4320738/150‑1 (MQ=255)

              gcTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACa                                                                                                                                        <  1:1524042/150‑1 (MQ=255)

                    gCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTTGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCg                                                                                                                                  >  2:4024247/1‑150 (MQ=255)

                        gcgGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGc                                                                                                                              <  2:446903/150‑1 (MQ=255)

                          gggCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATc                                                                                                                                   >  1:1108933/1‑143 (MQ=255)

                          gggCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATc                                                                                                                                   <  2:1108933/143‑1 (MQ=255)

                          gggCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGccc                                                                                                                            <  1:2678906/150‑1 (MQ=255)

                              gATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGc                                                                                                                        <  1:3416877/150‑1 (MQ=255)

                               aTCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGcc                                                                                                                       <  2:4756886/150‑1 (MQ=255)

                                tCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCg                                                                                                                      >  2:1446530/1‑150 (MQ=255)

                                 ccGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGa                                                                                                                     <  1:2384465/150‑1 (MQ=255)

                                   gTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGc                                                                                                                   >  1:2485260/1‑150 (MQ=255)

                                     gctgcCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGcac                                                                                                                 >  2:3203589/1‑150 (MQ=255)

                                     gctgcCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGcac                                                                                                                 >  1:526525/1‑150 (MQ=255)

                                     gctgcCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGcac                                                                                                                 >  2:3952852/1‑150 (MQ=255)

                                         ccTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGc                                                                                                                   >  2:409598/1‑144 (MQ=255)

                                         ccTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGc                                                                                                                   <  1:409598/144‑1 (MQ=255)

                                         ccTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTc                                                                                                             <  2:1240087/150‑1 (MQ=255)

                                         ccTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTc                                                                                                             >  1:5208632/1‑150 (MQ=255)

                                         ccTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTc                                                                                                             <  2:894962/150‑1 (MQ=255)

                                          cTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCt                                                                                                            <  1:2881521/150‑1 (MQ=255)

                                          cTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCt                                                                                                            <  1:4813024/150‑1 (MQ=255)

                                           tCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGGTCGACAACCTGTCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTg                                                                                                           >  1:2210255/1‑150 (MQ=255)

                                            cGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGc                                                                                                          <  1:4574776/150‑1 (MQ=255)

                                             ggTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCa                                                                                                         <  2:2532380/150‑1 (MQ=255)

                                             ggTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCa                                                                                                         <  2:352693/150‑1 (MQ=255)

                                             ggTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCa                                                                                                         <  2:3564413/150‑1 (MQ=255)

                                               tCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAgg                                                                                                       >  1:369995/1‑150 (MQ=255)

                                                ccAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGc                                                                                                      <  1:3663844/150‑1 (MQ=255)

                                                 cagcagGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCg                                                                                                     <  2:3072235/150‑1 (MQ=255)

                                                 cagcagGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCg                                                                                                     <  2:5005762/150‑1 (MQ=255)

                                                  agcagGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCgg                                                                                                    <  1:1256357/150‑1 (MQ=255)

                                                  agcagGCCTGGGCCCGCCACGGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCgg                                                                                                    <  2:665378/150‑1 (MQ=255)

                                                       gCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTgc                                                                                               <  1:1318594/150‑1 (MQ=255)

                                                            ggCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCAt                                                                                          <  2:2001779/150‑1 (MQ=255)

                                                              cccGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGc                                                                                        >  1:784012/1‑150 (MQ=255)

                                                                   cACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGAcc                                                                                   >  2:3665450/1‑150 (MQ=255)

                                                                    aCTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGAccc                                                                                  <  2:5399406/150‑1 (MQ=255)

                                                                      tctccGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCTTCGTCTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCtt                                                                                <  2:2534944/147‑1 (MQ=255)

                                                                      tGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCtt                                                                                >  1:595201/1‑150 (MQ=255)

                                                                      tGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCtt                                                                                >  1:595302/1‑150 (MQ=255)

                                                                       ggCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTtg                                                                               <  2:1862999/150‑1 (MQ=255)

                                                                       ggCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTtg                                                                               <  1:5162411/150‑1 (MQ=255)

                                                                       ggCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTtg                                                                               <  2:4159327/150‑1 (MQ=255)

                                                                         ccGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTtggt                                                                             <  2:3148644/150‑1 (MQ=255)

                                                                         ccGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTtggt                                                                             <  2:5208632/150‑1 (MQ=255)

                                                                          cGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTtggtg                                                                            >  1:1385358/1‑150 (MQ=255)

                                                                          cGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCAACGACTATGGCGACTGTCCCAACGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCAAGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTtgaga                                                                            >  1:2534944/1‑147 (MQ=255)

                                                                            tGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGc                                                                          <  2:1669196/150‑1 (MQ=255)

                                                                            tGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGc                                                                          >  2:1941145/1‑150 (MQ=255)

                                                                            tGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGc                                                                          >  2:1097733/1‑150 (MQ=255)

                                                                            tGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGc                                                                          <  1:1719424/150‑1 (MQ=255)

                                                                            tGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGGACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGc                                                                          <  2:2039043/150‑1 (MQ=255)

                                                                              ggCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGa                                                                        >  2:2737972/1‑150 (MQ=255)

                                                                              ggCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGCCGACTCTCCCTTCTACAGAGGCACGCCGCAGTCGGTACCGGACAGCATCGAGCCCCAACCCGAGCACATTCTGCAGGCGGGATGCGCCAAGCTGACCCTTGGTGGCGa                                                                        >  2:3262058/1‑150 (MQ=255)

                                                                                 gttcgacccgttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCa                                                                     >  2:854872/1‑150 (MQ=255)

                                                                                 gttcgacccgttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCa                                                                     <  1:2737972/150‑1 (MQ=255)

                                                                                 gttcgacccgttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTGGGTGGCGACCa                                                                     <  2:266275/150‑1 (MQ=255)

                                                                                 gttcgacccgttcgaccACCTGGCAGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCa                                                                     <  1:4225615/150‑1 (MQ=255)

                                                                                   tcgacccgttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACt                                                                   <  2:1486810/150‑1 (MQ=255)

                                                                                       cccgttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTtcat                                                               <  1:1259569/150‑1 (MQ=255)

                                                                                       cccgttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTtcat                                                               <  2:4949950/150‑1 (MQ=255)

                                                                                          gttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAg                                                            >  2:4398648/1‑150 (MQ=255)

                                                                                                         ggTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCaa                                             <  2:2822553/150‑1 (MQ=255)

                                                                                                          gTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAg                                            <  1:1977200/150‑1 (MQ=255)

                                                                                                              tCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTc                                        >  1:2678855/1‑150 (MQ=255)

                                                                                                               cGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCa                                       >  1:5264468/1‑150 (MQ=255)

                                                                                                                gACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCAc                                      <  1:5398268/150‑1 (MQ=255)

                                                                                                                      ggCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCt                                                                                                            <  1:3831724/74‑1 (MQ=255)

                                                                                                                      ggCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCt                                                                                                            >  2:3831724/1‑74 (MQ=255)

                                                                                                                        cGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCccgccg                              >  1:325301/1‑150 (MQ=255)

                                                                                                                          aCTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCgccgcc                            <  1:5334672/150‑1 (MQ=255)

                                                                                                                            tGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCAt                          >  2:4231340/1‑150 (MQ=255)

                                                                                                                              tCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCa                                                                     >  1:4846648/1‑105 (MQ=255)

                                                                                                                               ggtaccGACTGCGGCGTGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGc                       >  1:4467177/6‑150 (MQ=255)

                                                                                                                                ccTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGcc                      >  1:684507/1‑150 (MQ=255)

                                                                                                                                ccTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGcc                      <  1:3570830/150‑1 (MQ=255)

                                                                                                                                   ccGACTGCGGCGTGCCGCTGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGc                   >  2:1043225/2‑150 (MQ=255)

                                                                                                                                     gACTGCGGCGTGCCGCTGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGa                                                                                     >  1:2164976/1‑82 (MQ=255)

                                                                                                                                     gACTGCGGCGTGCCGCTGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGa                                                                                     <  2:2164976/82‑1 (MQ=255)

                                                                                                                                     gACTGCGGCGTGCCGCTGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTg                 >  1:3035355/1‑150 (MQ=255)

                                                                                                                                     gACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTg                 >  1:1780358/1‑150 (MQ=255)

                                                                                                                                      aCAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTgg                >  2:3938864/1‑150 (MQ=255)

                                                                                                                                       cAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTGgc               <  1:2947295/150‑1 (MQ=255)

                                                                                                                                         gCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTGgcgc             <  1:854872/150‑1 (MQ=255)

                                                                                                                                         gCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTGgcgc             <  1:854824/150‑1 (MQ=255)

                                                                                                                                             cACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTGGCGCTGAt         >  2:3492511/1‑150 (MQ=255)

                                                                                                                                             cACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTGGCGCTGAt         >  2:3461160/1‑150 (MQ=255)

                                                                                                                                                    cAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTGGCGCTGATCCTTTTc  <  2:5147932/150‑1 (MQ=255)

CCGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTGGCGCTGATCCATTTC  >  NC_002947/2504814‑2505111

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