BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,728,109	A→C	6.9%	D421E (GAT→GAG)	xylB ←	xylulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,728,109	0	A	C	6.9%	722.2 / 2.4	263	D421E (GAT→GAG)	xylB	xylulokinase
Reads supporting (aligned to +/- strand): ref base A (74/161); new base C (21/4); total (97/166)
Fisher's exact test for biased strand distribution p-value = 4.27e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 272 total aligned reads displayed.

GTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCG  >  NC_000913/3727960‑3728257

gTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCcaca                                                                                                                                                       >  1:4223069/1‑150 (MQ=255)

   tctcGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACCT‑cc                                                                                                                                                   >  1:268074/1‑150 (MQ=255)

   tctcGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑cc                                                                                                                                                   >  1:1501342/1‑150 (MQ=255)

    ctcGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTACCGGTAGTTGCGGCAACAAGTCAATGAGCGCGTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑ccc                                                                                                                                                  <  2:834740/150‑1 (MQ=255)

    ctcGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACCT‑ccc                                                                                                                                                  >  2:2818514/1‑150 (MQ=255)

     tcgtcgTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCC‑GTGCTGGCCCCCCAG‑cccc                                                                                                                                                 >  2:3459945/1‑150 (MQ=255)

     tcgtcgTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑cccc                                                                                                                                                 >  2:1510173/1‑150 (MQ=255)

      cgtcgtGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑ccccc                                                                                                                                                <  2:2595682/150‑1 (MQ=255)

        tcgtGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑ccccccc                                                                                                                                              <  1:4144708/150‑1 (MQ=255)

         cgtGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑cccccccc                                                                                                                                             <  2:2449948/150‑1 (MQ=255)

          gtGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CGCGCCCCg                                                                                                                                            >  1:1776078/1‑150 (MQ=255)

             gCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCCCCG‑CCCCCCCCggcc                                                                                                                                         >  2:1168966/1‑147 (MQ=255)

              cTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACg                                                                                                                                        <  1:271111/150‑1 (MQ=255)

                gATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGt                                                                                                                                      <  2:3901224/150‑1 (MQ=255)

                      cGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGa                                                                                                                                <  1:1797641/150‑1 (MQ=255)

                      cGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGa                                                                                                                                <  1:1563132/150‑1 (MQ=255)

                      cGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGGACGGGACTCGa                                                                                                                                >  2:3309519/1‑150 (MQ=255)

                       ggCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAg                                                                                                                               <  1:3979894/150‑1 (MQ=255)

                        gCATAACTCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTATTCGAgc                                                                                                                              <  1:3525557/150‑1 (MQ=255)

                        gCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAgc                                                                                                                              <  1:1670158/150‑1 (MQ=255)

                        gCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAgc                                                                                                                              <  1:2452178/150‑1 (MQ=255)

                         cATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAgct                                                                                                                             <  2:3795683/150‑1 (MQ=255)

                         cATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAgct                                                                                                                             >  2:3615107/1‑150 (MQ=255)

                          aTAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAgctg                                                                                                                            <  2:3360947/150‑1 (MQ=255)

                          aTAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAG‑CCCCCCCCGTCCGGTACTCGAgggg                                                                                                                            >  2:1424780/1‑147 (MQ=255)

                           tAACGCTGCGCATCTGGTAGATGCGACTGTCCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑‑CCCCCCCGTACGGTAATCGAgctgct                                                                                                                          >  1:3135471/1‑150 (MQ=255)

                             aCGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑ACCCCCCCGTACGGTAATCGAGctgctg                                                                                                                         <  2:3347565/150‑1 (MQ=255)

                             aCGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATGTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGctgctg                                                                                                                         <  2:2631439/150‑1 (MQ=255)

                              cGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGa                                                                                                                        <  1:2956963/150‑1 (MQ=255)

                               gCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACCT‑CCCCCCCCGTACGGGGATCGAGctgctgtc                                                                                                                       >  2:493775/1‑148 (MQ=255)

                                 tGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCg                                                                                                                     <  1:533995/150‑1 (MQ=255)

                                  tcgcATCTACCAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATCAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTCCGCGGGCTGCTGGCCGc                                                                                                                    >  2:4288634/2‑150 (MQ=255)

                                  gcgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCg                                                                                                                     <  1:2041980/150‑1 (MQ=255)

                                   cgcATCTGGTAGAGGCGACTGTTCTCACGGTAGGTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAG‑CCCCCCCCGTACGGTACTCGCGCTGCTGACCGCt                                                                                                                   <  2:617028/150‑1 (MQ=255)

                                         tGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGGGCGGGCGTCGGGCTGCTGACCGCTGAGATc                                                                                                             >  2:3745304/1‑150 (MQ=255)

                                          ggTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATcc                                                                                                            <  1:289179/150‑1 (MQ=255)

                                          ggTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATcc                                                                                                            <  2:4332538/150‑1 (MQ=255)

                                          ggTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATcc                                                                                                            <  2:539721/150‑1 (MQ=255)

                                          ggTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATcc                                                                                                            <  1:110239/150‑1 (MQ=255)

                                           gTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCg                                                                                                           <  2:971156/150‑1 (MQ=255)

                                           gTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCg                                                                                                           <  1:1452614/150‑1 (MQ=255)

                                              gATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCa                                                                                                        <  1:1425574/150‑1 (MQ=255)

                                              gATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCa                                                                                                        <  2:3396056/150‑1 (MQ=255)

                                              gATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCa                                                                                                        <  2:1445269/150‑1 (MQ=255)

                                              gATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCa                                                                                                        <  1:1266028/150‑1 (MQ=255)

                                               ttgggccTGTTCGTACGGTTGGTGGGGCAACAATTCAATGCGCGCTTTCCCGTGTTTCGGCGGGGTCTTGGGCAGGCTTTCTGGGGCCAGGGGCGGGCCCAACTC‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCa                                                                                                        <  2:101471/144‑1 (MQ=255)

                                               aTGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAg                                                                                                       <  1:1882627/150‑1 (MQ=255)

                                               aTGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAg                                                                                                       <  1:3067775/150‑1 (MQ=255)

                                                    aCTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCt                                                                                                  <  2:2843260/150‑1 (MQ=255)

                                                    aCTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCt                                                                                                  <  1:3604305/150‑1 (MQ=255)

                                                    aCTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAG‑CCCCCCCCGTGCGGGGGTCGTGCTGCTGGCCGCGGGTTTCCGCCGGCTTCt                                                                                                  >  2:1275107/1‑150 (MQ=255)

                                                         tCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGc                                                                                             <  1:2503673/150‑1 (MQ=255)

                                                          cTAACGTTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCTCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCCCAG‑CCCCCCCCGGACGGTAAGCGTGCGGCTGAGCGGGGATATCCGCCAGGAACTGCCGca                                                                                            >  1:94504/1‑149 (MQ=255)

                                                            aaCGGTAGTTGCGGCAACAATTCACTGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAg                                                                                          <  2:4139361/150‑1 (MQ=255)

                                                              cGGTAGTTGCGGCAACAATTCAATGAGCGATGTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTa                                                                                        <  1:1741122/150‑1 (MQ=255)

                                                               ggTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTGCGGTACGCGAGCTGCTGAGCGCGGAGGGCCGCCCGGCGCTGACGGCAGtcc                                                                                       >  2:1707656/1‑148 (MQ=255)

                                                                        ggCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGGACGGTAATCGAGCTGCTGACCGCTGCTATCCGCCAGCAGCTGACGCCCGTACGCACTAcgc                                                                              >  2:3407975/1‑150 (MQ=255)

                                                                          caacaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACgcgc                                                                            <  1:1343391/150‑1 (MQ=255)

                                                                          caacaaTTAAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCAAGCATCTGACGCCATTACTCACTACgcgc                                                                            <  1:11563/150‑1 (MQ=255)

                                                                             caaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACCT‑CCCCCCCCGGACGGTCATCGGGCTGCTGGCCGCTGCTCTCCGCCCGCGGGGGCCGCCAGTCCTCAGTTCGcgcccg                                                                         >  2:1955650/1‑149 (MQ=255)

                                                                              aaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTGCTCGGGCTGCTGGCCGCTGATGTCCGCCCGCATCTGAGGCCGGTACTCACGGCGcgccccg                                                                        >  2:291070/1‑150 (MQ=255)

                                                                               aTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccccgc                                                                       <  2:1500793/150‑1 (MQ=255)

                                                                                ttCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACCT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccccgcc                                                                      <  2:2040286/150‑1 (MQ=255)

                                                                                    aTGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACCG‑CCCCCCCCGTGCGGTGGTGGGGGTGCTGGCCGGGGGTGTCCGCCAGCGTCTGACGCCAGTGCTCACTAGGCGGCCCGgccccc                                                                  >  1:2361607/1‑149 (MQ=255)

                                                                                         cGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGTCCGCGGATATCCGCCAGCATCTGTCGCCAGTCCTCACTACGCGCCCCGCCCCCAAGCac                                                             >  2:2325740/1‑149 (MQ=255)

                                                                                         cGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCaa                                                             <  1:3654074/150‑1 (MQ=255)

                                                                                         cGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGGAATCGAGCTGCTGACCGCTGATATCCGCCAGCCTCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATcgc                                                             >  2:4168196/1‑148 (MQ=255)

                                                                                         cGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGGAATCGAGCTGCTGACCGCGGATATCCGCCCGCATCTGACGCCAGTGCTCACTACGCGCCCCGcccccccgcaa                                                             >  2:1513239/1‑144 (MQ=255)

                                                                                                    gATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTACTCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGGAACACttt                                                  >  2:4215502/1‑150 (MQ=255)

                                                                                                         gccgcGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTGCGGTGATGGGGCTGCTGGCCGGTGATATCCGCCCGCAGCTGGCGCCGGTACTCACTACGCGCCCCGCCCCCCAGCAACGGCACACGTTGCggg                                             >  2:3000910/1‑149 (MQ=255)

                                                                                                             cgATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTaa                                         <  2:3960861/150‑1 (MQ=255)

                                                                                                                   gcgcCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGc                                   <  1:2180865/150‑1 (MQ=255)

                                                                                                                     gcCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAg                                 <  1:2857125/150‑1 (MQ=255)

                                                                                                                     gcCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAg                                 <  1:215338/150‑1 (MQ=255)

                                                                                                                     gcCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCACGTACGGTTATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCACTCAACGTAACACTTTGCGGTATTAAACCGCaa                                 <  2:825262/150‑2 (MQ=255)

                                                                                                                         gCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCAt                             <  1:65874/150‑1 (MQ=255)

                                                                                                                         gCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCAt                             <  1:3863600/150‑1 (MQ=255)

                                                                                                                               ctgcGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGa                       <  2:3670470/150‑1 (MQ=255)

                                                                                                                                  cgcCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTCCGGTAATCGAGCTGCTGGCCGCGGATAGCCGCCCGCCTCGGCCGCCAGTACGCACTACGCGCCCCGCCCCCCCTCACCGGAACACGTTGCGGGTGAATACCGCCGGCCTGCACGACag                    >  2:1857842/1‑149 (MQ=255)

                                                                                                                                  cgcCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACAt                    <  1:69217/150‑1 (MQ=255)

                                                                                                                                  cgcCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGAGATCCGCCAGCGTCTGACGCCAGTACGCACGACGCGCCCCGCCCCCCATCAACGTACCACGGTGCGGGGTATTGCGGCAGGCCTGCACGagcg                    >  2:2368525/1‑147 (MQ=255)

                                                                                                                                  cgcCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGGAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACGCACTACGCGCCCCGCCCCCACGCAACGTAACACTTTGCGGTTTAATACCGCAGGCAGGCACGACAt                    >  2:4058225/1‑150 (MQ=255)

                                                                                                                                   gcCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTACTCGAGCTGCTGACCGCTGATATCCGCCCGCGTCTGACGCCCGTCCTCACTACGCGCCCCGCCCCCAAGCAACGTAACACGGTGCGGTTTAATACCGCTGGCATGCGCGACGTc                   >  2:926776/1‑150 (MQ=255)

                                                                                                                                     ccA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCa                 <  2:2942167/150‑1 (MQ=255)

                                                                                                                                       a‑GTGCTGGCCCCACAT‑CCCCCCCCGTTCGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCCTCTGACGCCAGTTCTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGGTTAATACCGCAGGCCTGCACGACATCCATg               >  2:2548848/1‑150 (MQ=255)

                                                                                                                                       a‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCCATCAACGTAACGCGTTGCGGGTGAATACGGCAGGCATGCACGACATCCTTg               >  2:771224/1‑150 (MQ=255)

                                                                                                                                       a‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATg               <  1:1301083/150‑1 (MQ=255)

                                                                                                                                         gTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGc              <  2:3636794/150‑1 (MQ=255)

                                                                                                                                              ggCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCAAGCACGACAACCATGCCATCt         >  1:2247425/1‑150 (MQ=255)

                                                                                                                                               gCCCCACCT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACATTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTg        >  1:4215018/1‑150 (MQ=255)

                                                                                                                                               gCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCAACTGACGCCCGTACTCACTCCGCGCCCCGCCCCCCATCAACGGGACACTTTGCGGTGTAATACCGCAGGCATGCACGACATCCATGCCATCTg        >  1:3385112/1‑150 (MQ=255)

                                                                                                                                                ccccACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCTGGCATGCACGACATCCATGCCATCTGc       <  1:4287004/150‑1 (MQ=255)

                                                                                                                                                 cccACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCAGATATCCGCCAGCATCTGACGCCAGTACACACTACGCGCCCCGCCCCCAAACAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCAAGCCATCTGcc      >  1:3809976/1‑150 (MQ=255)

                                                                                                                                                  ccACCT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCACTCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCa     <  1:2201550/150‑1 (MQ=255)

                                                                                                                                                  ccACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCa     >  2:4015373/1‑150 (MQ=255)

                                                                                                                                                  ccACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCa     <  1:2318287/150‑1 (MQ=255)

                                                                                                                                                   cacaT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAg    >  1:4073370/1‑150 (MQ=255)

                                                                                                                                                    acaT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAgc   <  2:3007491/150‑1 (MQ=255)

                                                                                                                                                     caT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCCGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAgcg  <  2:3108275/150‑1 (MQ=255)

                                                                                                                                                     caT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAgcg  <  2:11396/150‑1 (MQ=255)

GTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCA‑GTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCG  >  NC_000913/3727960‑3728257

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`