BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,728,125	T→G	8.6%	Y416S (TAC→TCC)	xylB ←	xylulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,728,125	0	T	G	8.6%	1253.8 / 24.5	444	Y416S (TAC→TCC)	xylB	xylulokinase
Reads supporting (aligned to +/- strand): ref base T (128/266); new base G (46/4); total (174/270)
Fisher's exact test for biased strand distribution p-value = 1.56e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 447 total aligned reads displayed.

GATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTC  >  NC_000913/3727976‑3728273

   aaGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAAt                                                                                                                                                   <  2:622858/150‑1 (MQ=255)

    aGCGGCATAACGCTGCGCATCTGGTGGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATc                                                                                                                                                  <  1:1844300/150‑1 (MQ=255)

     gCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGa                                                                                                                                                <  1:616037/150‑1 (MQ=255)

      cGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGtgcgcgc                                                                                                                                                >  2:1458686/1‑144 (MQ=255)

       ggCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAg                                                                                                                                               <  1:3528929/150‑1 (MQ=255)

        gCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAgc                                                                                                                                              <  1:2707841/150‑1 (MQ=255)

        gCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAgc                                                                                                                                              <  2:1051546/150‑1 (MQ=255)

        gCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAgc                                                                                                                                              <  2:1993992/150‑1 (MQ=255)

         cATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTCCGGTGATCGAgcg                                                                                                                                             >  2:639084/1‑149 (MQ=255)

               gCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGTGCGGTCATCGGggtgcggtc                                                                                                                                       >  2:3672234/1‑146 (MQ=255)

                    gcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTg                                                                                                                                  >  2:91813/1‑150 (MQ=255)

                    gcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTg                                                                                                                                  <  2:1648460/150‑1 (MQ=255)

                         tGGTAGATGCGACTGTTCTAACTGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCTCGATCTGCGCCAGCCTTTCTGCGCCCAGTGCTGGCACCACATCCCCCCCCGTAAGGTACTCGAGCTGCTGACCGGTGttatg                                                                                                                             >  2:838971/1‑149 (MQ=255)

                         tGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCCCCCCCCGTGGGGTGCGCGAGCTGGGGGCCGCTggggtc                                                                                                                             >  2:1406733/1‑145 (MQ=255)

                          ggTAGCTGCGACTGTTCTACCGGTTGTTGCGGCAACAATTCACTGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGGGCTGGCCCCCCATCCCCCCCCGTACGGTAATCGCGCTGCTGACCGCTGATATcc                                                                                                                            <  2:2209871/150‑1 (MQ=255)

                          ggTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATcc                                                                                                                            <  2:1783182/150‑1 (MQ=255)

                           gTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTTCGGGCGTCGGGCTGCTGTCCGCTGAGATcgg                                                                                                                           >  2:2401109/1‑148 (MQ=255)

                           gTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATTTCCg                                                                                                                           >  2:53317/1‑150 (MQ=255)

                              gATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCa                                                                                                                        <  2:349065/150‑1 (MQ=255)

                              gATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCa                                                                                                                        <  1:779585/150‑1 (MQ=255)

                                   gACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATGCCCCCCCGTGGGGGGGTCGAGGTGCGGCCGGCGGGTGTCCGCCAGCGTc                                                                                                                   >  1:247367/1‑150 (MQ=255)

                                    aCTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCt                                                                                                                  <  2:310519/150‑1 (MQ=255)

                                      tGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGa                                                                                                                <  1:2069202/150‑1 (MQ=255)

                                        ttCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGGACGGTACTCGAGCTGCTGACCGCTGCGATCCGCCCGCATCTGGCg                                                                                                              >  2:2629875/1‑150 (MQ=255)

                                           tAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACTTCCCCCGCCGTGCGGGGATCGGGCTGGTGGGCGCGGTTGTGCGCCAGCGGCTGTCGccg                                                                                                           >  1:1709087/1‑149 (MQ=255)

                                           tAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCa                                                                                                           <  2:1556755/150‑1 (MQ=255)

                                           tAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCa                                                                                                           <  2:1218819/150‑1 (MQ=255)

                                              cGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTa                                                                                                        <  1:2791083/150‑1 (MQ=255)

                                                 tAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTc                                                                                                     <  1:3062869/150‑1 (MQ=255)

                                                 tAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTc                                                                                                     <  1:3821670/150‑1 (MQ=255)

                                                   gTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCGCCCCCGCCGGTCGGGAATCGGGCTGCTGACCGCTGGTATCCGCCAGCATCTGGCGCCAGGACTCag                                                                                                   >  1:1935231/1‑149 (MQ=255)

                                                    ttGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGTACGGGGATCGAGCTGCTGACCGCGGATCGCCGCCCGCATCGGACGCCAGTGCTcgcg                                                                                                  >  2:2802733/1‑147 (MQ=255)

                                                    ttGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTTATCGAGCTGCTGACCGCTGATATCCGCCAGCCTCTGACGCCTGTACTCGCt                                                                                                  >  2:1881500/1‑150 (MQ=255)

                                                    ttGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACt                                                                                                  <  1:1373681/150‑1 (MQ=255)

                                                        ggCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACGTCCGCCCCCGTACGGGGGGCGAGCTGCTGGGCGCGGAGGTCCGCGGGCGGCTGGCGCCAGGCCTCAcgccgg                                                                                              >  1:2977746/1‑145 (MQ=255)

                                                          caacaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACgcgc                                                                                            <  1:4057112/150‑1 (MQ=255)

                                                           aacaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAGTCGAGCTGCTGGCCGCTGATATCCGCCCGCATCGGTCGCCGGGGCTCGCTAGGcgcc                                                                                           >  2:773833/1‑150 (MQ=255)

                                                              aaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgcccca                                                                                        <  2:2345178/150‑2 (MQ=255)

                                                                 tCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccccgccc                                                                                     <  1:91813/150‑1 (MQ=255)

                                                                    aTGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCa                                                                                  <  1:3216335/150‑1 (MQ=255)

                                                                      gAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCCCATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAAt                                                                                <  2:3368163/150‑1 (MQ=255)

                                                                       aGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATc                                                                               <  2:343334/150‑1 (MQ=255)

                                                                         cGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCGCCCCCCCCGTGCGGGAGGCGAGCGGCGGGCCGCGGGGAGCCGCCCGCGTGGGGCGCCAGGGCTCGCGGCGCGCCCCGcccccgcgccg                                                                             >  1:3912016/1‑144 (MQ=255)

                                                                         cGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCaa                                                                             <  1:2374205/150‑1 (MQ=255)

                                                                            tgtctcTGGCTTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTTGTGGTGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCAAGTACTCACTACGCGCCCCGCCCCCACTCAACGt                                                                          <  2:2837325/148‑1 (MQ=255)

                                                                              tctcTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTaa                                                                        <  1:2130195/150‑1 (MQ=255)

                                                                                tctGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGGGCCCCGCCCCCAATCAACGTAaca                                                                      <  2:2561052/150‑1 (MQ=255)

                                                                                     aTTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTg                                                                 <  2:2994851/150‑1 (MQ=255)

                                                                                      ttCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGc                                                                <  2:3994024/150‑1 (MQ=255)

                                                                                          ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTGCTCTCTACGCGCCCCGCCGCCTATCAACGTAACACGTTGCGGtt                                                            >  2:3224527/1‑150 (MQ=255)

                                                                                          ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGtt                                                            <  1:1463424/150‑1 (MQ=255)

                                                                                          ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGtt                                                            <  2:263836/150‑1 (MQ=255)

                                                                                          ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGtt                                                            <  1:1255176/150‑1 (MQ=255)

                                                                                          ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGttt                                                           <  2:1920763/150‑1 (MQ=255)

                                                                                           cgcgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACGGCCCCCCCCGGGCGGGAATGGAGCTGCTGACCGCTGGTGTCCGCCAGCATCTGACGCGGGTACTCGCTGCGCGCCCCGCCCCCAGTCAACGGGACGCTGTGCggggg                                                           >  1:3616214/1‑147 (MQ=255)

                                                                                           cgcgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGttt                                                           <  2:4275185/150‑1 (MQ=255)

                                                                                             cgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTaa                                                         <  2:89487/150‑1 (MQ=255)

                                                                                              gATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAAt                                                        <  1:1805416/150‑1 (MQ=255)

                                                                                               aTCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATa                                                       <  1:2882833/150‑1 (MQ=255)

                                                                                                  tGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCg                                                    <  1:1082466/150‑1 (MQ=255)

                                                                                                  tGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCg                                                    <  1:2256335/150‑1 (MQ=255)

                                                                                                   gcgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCGGTGAGCGTGCTGCTGGCCGCTGGTATCCGCCGGCGTCTGACGCCAGTTCTCACTACGCGCCCCGCCCCCGGTCCACGGAACACTTTGCGGGGTATTACCGc                                                   >  2:207672/1‑150 (MQ=255)

                                                                                                    cgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCTTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCa                                                  <  2:246032/150‑1 (MQ=255)

                                                                                                    cgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCa                                                  <  1:3763548/150‑1 (MQ=255)

                                                                                                         gCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCAt                                             <  2:443901/150‑1 (MQ=255)

                                                                                                            ttGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCCATCAACGTAACACTGTGCGGGTTAATACCGCAGGCATGcc                                          >  2:2929883/1‑149 (MQ=255)

                                                                                                            ttGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCa                                          <  2:3045256/150‑1 (MQ=255)

                                                                                                             tgctgcGCCCAGTGCTGGCCCCACATCCACCCCCGTACGGTAATCGAGCTGCTGACCTCTGATATCCGCCTGCATCTGACGCCAGTACTCAGTACGCGCCACTCCCCCTATCAACGTTTCACTTTGCGGTTTTATACCGCAGGCATGCAc                                         <  2:1711170/150‑1 (MQ=255)

                                                                                                               ctgcGCCCAGTGCTGGCCCCGCAGCCCCCCCCGGGCGGGCGTCGAGCTGCTGACCGCTGGTCTCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTACCACTTTGCGGTTTAATACCGCAGGCATGCACGa                                       >  2:3814106/1‑150 (MQ=255)

                                                                                                               ctgcGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGa                                       <  1:2612662/150‑1 (MQ=255)

                                                                                                               ctgcGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGa                                       >  2:1855470/1‑150 (MQ=255)

                                                                                                                tgcGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTCCTCGAGCTGCTGGCCGCTGATATCCGCCAGCATCTGACGCCAGTGCTCACTACGCGCCCCGCCCCCCCTCCACGTAACACTTTGCGGTTTAATACCGCAGGCGTGCACGAc                                      >  2:2288312/1‑150 (MQ=255)

                                                                                                                 gcgcCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACa                                     <  1:3111121/150‑1 (MQ=255)

                                                                                                                   gcCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATc                                   <  1:1424624/150‑1 (MQ=255)

                                                                                                                     ccAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCa                                 <  1:2488226/150‑1 (MQ=255)

                                                                                                                     ccAGTGCTGGCCCCACATCCCCCCCCGGACGGTCATCGAGCTGCTGACCGCGGGTATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCCAGCAACGGAACACGTTGCGGTGTAAGACCGCCGGCCTGCACGACATccc                                 >  1:1178110/1‑149 (MQ=255)

                                                                                                                     ccAGTGCTGGCCCCACATCCCCCCCCGGACGGTCATCGAGCTGCGGACCGCGGATATCCGCCAGCATCTGGCGCCCGGACTCACTACGCGCCCCGCCCCCCATCAACGTAACACTTTGCGGTTTTATACCGCAGGCATGCACGACATCCa                                 >  1:684971/1‑150 (MQ=255)

                                                                                                                      cAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCGGACCGCGGATATCCGCCAGCATCTGACGCCCGTACGCACTCCGCGCCCCGCCCCCAATCAACGGAACACTTTGCGGTGTAATACCGCAGGCATGCACGACATCCAt                                >  1:1519877/1‑150 (MQ=255)

                                                                                                                       aGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGGGTAATACCGCAGGCATGCACGACCTCCATg                               >  1:1925967/1‑150 (MQ=255)

                                                                                                                            tGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATc                          >  2:3623011/1‑150 (MQ=255)

                                                                                                                             ggCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTCCTCACTACGCGCCCCGCCCCCAATCAACGTACCACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCt                         >  1:4150490/1‑150 (MQ=255)

                                                                                                                             ggCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCt                         <  1:1406845/150‑1 (MQ=255)

                                                                                                                              gCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCCACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTg                        >  2:3557489/1‑150 (MQ=255)

                                                                                                                              gCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTg                        >  1:1503622/1‑150 (MQ=255)

                                                                                                                                cccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGcc                      >  1:2033259/1‑150 (MQ=255)

                                                                                                                                cccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGcc                      <  1:571352/150‑1 (MQ=255)

                                                                                                                                 ccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCAACTGccc                     >  1:3344023/1‑149 (MQ=255)

                                                                                                                                   acaTCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAgc                   <  2:4238523/150‑1 (MQ=255)

                                                                                                                                   acaTCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTATGCGGTTTAATACCGCAGGCATGCACGACAACCATCCCATCTGCCAgc                   >  1:4258125/1‑150 (MQ=255)

                                                                                                                                         ccccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATaa             >  2:4109840/1‑150 (MQ=255)

                                                                                                                                          cccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAAc            >  2:390908/1‑150 (MQ=255)

                                                                                                                                             ccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCa         >  1:111213/1‑150 (MQ=255)

                                                                                                                                             ccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACCTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCa         >  1:1001041/1‑150 (MQ=255)

                                                                                                                                              cGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCAc        <  2:4054056/150‑1 (MQ=255)

                                                                                                                                                   ggTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCtt   >  1:1677582/1‑150 (MQ=255)

                                                                                                                                                   ggTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCtt   >  1:1601785/1‑150 (MQ=255)

                                                                                                                                                   ggTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCtt   >  2:3716146/1‑150 (MQ=255)

                                                                                                                                                   ggTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCtt   >  2:1277215/1‑150 (MQ=255)

                                                                                                                                                    gTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTc  <  1:3719458/150‑1 (MQ=255)

                                                                                                                                                    gTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTc  >  1:2335837/1‑150 (MQ=255)

GATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTC  >  NC_000913/3727976‑3728273

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