BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,728,109	A→C	8.1%	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	8.1%	397.4 / 2.6	167	D421E (GAT→GAG)	xylB	xylulokinase
Reads supporting (aligned to +/- strand): ref base A (54/89); new base C (15/2); total (75/92)
Fisher's exact test for biased strand distribution p-value = 1.06e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 171 total aligned reads displayed.

TTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGC  >  NC_000913/3727961‑3728256

tttCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAt                                                                                                                                                    <  2:1508895/150‑1 (MQ=255)

  tctcGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCTCACATcc                                                                                                                                                  >  2:1253614/1‑150 (MQ=255)

  tctcGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATcc                                                                                                                                                  <  1:1085247/150‑1 (MQ=255)

    tcgtcgTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTcgcc                                                                                                                                                >  2:1329961/1‑150 (MQ=255)

    tcgtcgTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCGcccc                                                                                                                                                >  2:1329985/1‑150 (MQ=255)

      gtcgtGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACGTcccccc                                                                                                                                              >  2:3445281/1‑150 (MQ=255)

       tcgtGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACGTcccccgc                                                                                                                                             >  1:3199735/1‑148 (MQ=255)

       tcgtGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccc                                                                                                                                             >  2:2356588/1‑150 (MQ=255)

       tcgtGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATccccccc                                                                                                                                             >  2:3529988/1‑150 (MQ=255)

        cgtGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATcccccccc                                                                                                                                            <  1:374001/150‑1 (MQ=255)

         gtGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCg                                                                                                                                           >  2:1468462/1‑150 (MQ=255)

         gtGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCg                                                                                                                                           <  1:1840541/150‑1 (MQ=255)

         gtGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGATGCGCCCAGTGCTGGCCCCACGAGCCCCCCCg                                                                                                                                           >  1:965780/1‑150 (MQ=255)

           ggCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCggg                                                                                                                                         >  2:291641/1‑148 (MQ=255)

            gCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGGCCCCCCCgggc                                                                                                                                        >  2:1311058/1‑147 (MQ=255)

              tGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTCCgg                                                                                                                                      >  2:2065354/1‑150 (MQ=255)

              tGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCCCCCCCCGTCCgg                                                                                                                                      >  2:2065379/1‑150 (MQ=255)

               gATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAAATCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCccgccggcgggt                                                                                                                                     >  1:3033813/1‑144 (MQ=255)

                   aGCGGCATAACGCTGCGCATCTGGTTGATTCGTCTGTTCTTAGTGTAGTTGCGGCTACAATTGTTTTAGGGTTTTCTCTGGTTTGGCGTCGTTCTGCGCGTGGGTTGCTGCGCCCAGTGCTGGCCCCACTTGCCCCCCCGTTAGGTTATc                                                                                                                                 <  2:3896873/150‑1 (MQ=255)

                   aGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATc                                                                                                                                 <  2:1943264/150‑1 (MQ=255)

                   aGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATc                                                                                                                                 <  2:2908488/150‑1 (MQ=255)

                     cGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCCCGTCGGCCCCCGTGCgggccgcgc                                                                                                                               >  1:1404765/1‑143 (MQ=255)

                     cGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGa                                                                                                                               <  1:725272/150‑1 (MQ=255)

                     cGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGa                                                                                                                               <  1:838634/150‑1 (MQ=255)

                     cGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCGCCCGCCGTGCggggagcga                                                                                                                               >  1:1459172/1‑143 (MQ=255)

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

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

                        cATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCGCCCCCCCCGGACGGTGATCgcgcg                                                                                                                            >  1:1807897/1‑149 (MQ=255)

                          tAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAgctgc                                                                                                                          <  1:3868964/150‑1 (MQ=255)

                                tGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCGCCCCCCGGACGGTGATCGGgctggggggcg                                                                                                                    >  1:940728/1‑143 (MQ=255)

                                tGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACgg                                                                                                                                      >  2:2784400/1‑132 (MQ=255)

                                   gcATCTGGTAGATTCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTg                                                                                                                 <  2:1411616/150‑1 (MQ=255)

                                   gcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCCCCCCCCGGCGGGTGGTCGGGCTGCTGGCCGcgg                                                                                                                 >  2:1116670/1‑148 (MQ=255)

                                      tCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGata                                                                                                              <  1:1351042/150‑1 (MQ=255)

                                        tGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCCCCGCCCCCCCCGGGCGGGCGTCGGGCTGCTGACCGCGGATGTc                                                                                                            >  2:1211960/1‑150 (MQ=255)

                                         ggTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGGACGCGGGCTGCTGACCGCGGATAGcc                                                                                                           >  2:2010153/1‑150 (MQ=255)

                                          gTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGAGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCGGTAATCGAGCTGCTGACCGCGGATATGCg                                                                                                          >  2:1781058/1‑150 (MQ=255)

                                            aGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGTACGGTCCTCGGGCTGCTGGCCGCTGATTTCCGcc                                                                                                        >  2:2182749/1‑150 (MQ=255)

                                                 cGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCAt                                                                                                   <  1:1424240/150‑1 (MQ=255)

                                                  gACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCCCCTCCCCCCCCGTGCGGGGAGCGCGGTGGTGGCCGCGGATGTCCGCCGGCagg                                                                                                  >  2:3564176/1‑148 (MQ=255)

                                                  gACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATc                                                                                                  <  2:489826/150‑1 (MQ=255)

                                                  gACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGGCATCGAGCTGCTGACCGCTGATATCCGCCCGCATc                                                                                                  >  2:3688707/1‑150 (MQ=255)

                                                  gACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCt                                                                                                 >  1:3797314/1‑150 (MQ=255)

                                                     tGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTACTCGAGCGGCTGACCGCTGATATCCGCCCGCCTCTGa                                                                                               >  2:1563929/1‑150 (MQ=255)

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

                                                      gTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCGCCCCCGTGCGGGGCGCGGGCTGCTGGCCGCGGGTAGCCGCCGGCATCTggc                                                                                              >  1:419213/1‑148 (MQ=255)

                                                      gTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGGACGGTAGTCGAGCTGCTGTCCGCGGAGATGCGCCAGCGGCTGAc                                                                                              >  1:2780521/1‑150 (MQ=255)

                                                                  gTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCAc                                                                                  <  2:871321/150‑1 (MQ=255)

                                                                  gTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCAc                                                                                  <  1:2639979/150‑1 (MQ=255)

                                                                  gTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCAc                                                                                  <  1:3529988/150‑1 (MQ=255)

                                                                  gTTGCGGCAACAATTCAATGAGCGAGTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCAc                                                                                  <  1:2673332/150‑1 (MQ=255)

                                                                   ttGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCGGTAATCGGGCTGCTGGCCGCTGGTATCCGCCGGCGGCTGGCGCCGGTAGTcgcg                                                                                 >  2:425113/1‑147 (MQ=255)

                                                                      cGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCGGGACTCGGGCTGCTGACCGCTGTGGTCCGCCAGCATCTGACGCCCGTGCTCAcgggg                                                                              >  2:874892/1‑146 (MQ=255)

                                                                          aacaaTTCAATGAGAGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCGGGAAGCGAGCTGCTGACCGCTGAGAGGCGCCGGCCTCTGACGCCCGGCCTCACTACGcgcc                                                                          >  2:3041620/1‑150 (MQ=255)

                                                                               ttCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGGTTCTGACGCCAGTACTCACTACGcgccccgcc                                                                     <  1:3874021/150‑1 (MQ=255)

                                                                               ttCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccccgcc                                                                     <  1:2268024/150‑1 (MQ=255)

                                                                               ttCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccccgcc                                                                     <  1:3588563/150‑1 (MQ=255)

                                                                                  aaTGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTCATCGAGCTGCTGACCGCTGATATCCGCCAGCAGCTGACGCCAGTACTCACTACGCGCCCCGccccc                                                                  >  2:2079602/1‑150 (MQ=255)

                                                                                  aaTGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCACCCCGccccc                                                                  <  2:1901029/150‑1 (MQ=255)

                                                                                      aGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCCCATCCCCCCCCGGACGGGAATCGAGCTGCTGACCGCTGCTATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATc                                                              >  2:3557131/1‑150 (MQ=255)

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

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

                                                                                         gttttCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAAc                                                           <  1:1147550/148‑1 (MQ=255)

                                                                                         gATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAAc                                                           <  1:1147632/150‑1 (MQ=255)

                                                                                         gATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAAc                                                           <  1:3699035/150‑1 (MQ=255)

                                                                                         gATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAAc                                                           <  1:794501/150‑1 (MQ=255)

                                                                                         gATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAAc                                                           <  1:3928918/150‑1 (MQ=255)

                                                                                               tctGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAaca                                                     <  1:892082/150‑1 (MQ=255)

                                                                                                   gATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGAGATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCACTCAACGTAACACgtt                                                 >  2:2152518/1‑150 (MQ=255)

                                                                                                   gATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTGATCGAGCTGCTGACCGCTGATGTCCGCCAGCGGCTGACGCCAGTGCGCCCTGCGCGCCCCGCCCCCCAGCACCGTCACACttt                                                 >  2:3748428/1‑150 (MQ=255)

                                                                                                       cgccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCgg                                             <  1:3482183/150‑1 (MQ=255)

                                                                                                         ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGTGCGGTCATCGGGCTGCTGGCCGCTGCTATCCGCCAGCCTCTGGCGCCCGTACGCGCGCCGCGCCCCGCCCCCCCGCGACGGACCACGTTGCgggg                                           >  2:2737814/1‑148 (MQ=255)

                                                                                                         ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTGCGCGAGCTCCTGCCCGCGGAGATCCGCCGGCATCTGGCGCCCGTACGCGCGTCGCGCCCCGCCGCCCATCAGCGGAACACGTGGCGGtg                                           >  2:2737969/1‑149 (MQ=255)

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

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

                                                                                                                cTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATAcc                                    <  1:186019/150‑1 (MQ=255)

                                                                                                                 tGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGGGCTGCTGACCGCGGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAGTCAACGTAACACTTTGCGGTGTAATACCg                                   >  2:1257602/1‑150 (MQ=255)

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

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

                                                                                                                 tGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGGAATCGGGCTGCTGACCGCTGATATCCGCCCGCGTCTGACGCCAGTACTCACTACGCGCCCCGCCCCCACTCTGCGTAACACTTTGCGGTTTAGTACCg                                   >  2:2251626/1‑150 (MQ=255)

                                                                                                                   cgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGGACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTCCGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCa                                 <  2:2228799/150‑1 (MQ=255)

                                                                                                                     ccAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTCCGGGAATCGAGCTGCTGGCCGCTGATATCCGCCCGCATCTGGCGCCGGTACTCACTACGCGCCCCGCCCCCCATCAACGTAACCCTTTGCGGTTTAATACCGCCgg                               >  2:493976/1‑150 (MQ=255)

                                                                                                                     ccAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCGGCTATCCGCCAGCCTCTGACGCCCGTGCTCACTGCGCGCCCCGCCCCCCCGCAACGGCCCACTGTGCGGTGTAATGCCGCAgg                               >  2:2607226/1‑150 (MQ=255)

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

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

                                                                                                                           ttGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGGCCGCTGCTATCCGCCAGCGTCTGACGCCAGTACGCACTGCGCGCCCCGCCCCCCCGCCCCGTAACACTTTGCGGGGTAATACCGCGGGCAGGcc                         >  2:3806044/1‑149 (MQ=255)

                                                                                                                             gctgcGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACg                       <  1:2244593/150‑1 (MQ=255)

                                                                                                                             gctgcGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACg                       <  1:3751209/150‑1 (MQ=255)

                                                                                                                               tgcGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCCGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGAc                     <  2:905288/150‑1 (MQ=255)

                                                                                                                                   cccAGTGCTGGCCCCACCTCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGGTATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATcc                 >  1:879577/1‑150 (MQ=255)

                                                                                                                                     ctgTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTTCGCGCCCCGCCCCCAATCAACGTAACTCTTTGCGGTTTAATACCGCAGGCTTGCACGACATCCAt               <  2:3152064/148‑1 (MQ=255)

                                                                                                                                     cAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCAt               <  1:1329961/150‑1 (MQ=255)

                                                                                                                                     cAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCAt               <  2:3423364/150‑1 (MQ=255)

                                                                                                                                     cAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCAt               <  2:954059/150‑1 (MQ=255)

                                                                                                                                      aGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATg              <  1:2093410/150‑1 (MQ=255)

                                                                                                                                      aGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATg              <  1:1329028/150‑1 (MQ=255)

                                                                                                                                             gCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTg       <  1:1644198/150‑1 (MQ=255)

                                                                                                                                             gCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTg       <  1:3184590/150‑1 (MQ=255)

                                                                                                                                               cccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCACCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCAACTGcc     >  2:886904/1‑150 (MQ=255)

                                                                                                                                                  acaTCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAgc  >  2:3528597/1‑150 (MQ=255)

TTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGC  >  NC_000913/3727961‑3728256

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