BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,356,985	T→C	5.6%	L64L (CTA→CTG)	PP_2950 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,356,985	0	T	C	5.6%	1030.4 / 29.9	321	L64L (CTA→CTG)	PP_2950	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (165/138); new base C (15/3); total (180/141)
Fisher's exact test for biased strand distribution p-value = 2.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.02e-01

Only 100 of 327 total aligned reads displayed.

AACTCCCAGTTTTCGTTGCCATAGGAGCGAAACCAGTTGCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGC  >  NC_002947/3356836‑3357133

 aCTCCCAGTTTTCGTTGCCATAGGAGCGAAACCAGTTGCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTa                                                                                                                                                     <  1:4222867/150‑1 (MQ=255)

 aCTCCCAGTTTTCGTTGCCATAGGAGCGAAACCAGTTGCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTa                                                                                                                                                     <  1:3432280/150‑1 (MQ=255)

   tCCCAGTTTTCGTTGCCATAGGAGCGAAACCAGTTGCCCGAATCGTCGTGCGACTCGTAGGCGTAACGGACCGCGATGGGGTTGTCGGTGAAGGCGCAGAGCTCTTTGTTCAACCGATAGTCCAGTTCTTTCTTCCTCTTGCGAGTTAgg                                                                                                                                                   >  2:1800280/1‑150 (MQ=255)

   tCCCAGTTTTCGTTGCCATAGGAGCGAAACCAGTTGCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAgg                                                                                                                                                   >  1:1354169/1‑150 (MQ=255)

   tCCCAGTTTTCGTTGCCATAGGAGCGAAACCAGTTGCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAgg                                                                                                                                                   >  1:6832411/1‑150 (MQ=255)

   tCCCAGTTTTCGTTGCCATAGGAGCGAAACCAGTTGCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAgg                                                                                                                                                   >  1:2659993/1‑150 (MQ=255)

   tCCCAGTTTTCGTTGCCATAGGAGCGAAACCAGTTGCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAgg                                                                                                                                                   >  2:2800827/1‑150 (MQ=255)

     ccAGTTTTCGTTGCCATAGGAGCGAAACCAGTTGCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGGGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGaa                                                                                                                                                 >  1:3837772/1‑150 (MQ=255)

     ccAGTTTTCGTTGCCATAGGAGCGAAACCAGTTGCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGaa                                                                                                                                                 >  2:4677061/1‑150 (MQ=255)

     ccAGTTTTCGTTGCCATAGGAGCGAAACCAGTTGCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGaa                                                                                                                                                 <  1:2046865/150‑1 (MQ=255)

        gTTTTCGTTGCCATAGGAGCGAAACCAGTTGCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTc                                                                                                                                              <  1:5401040/150‑1 (MQ=255)

         ttttCGTTGCCATAGGAGCGAAACCAGTTGCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCTGTGAAGGCCCATAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTcc                                                                                                                                             <  2:4011689/150‑1 (MQ=255)

         ttttCGTTGCCATAGGAGCGAAACCAGTTGCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTcc                                                                                                                                             >  2:8012028/1‑150 (MQ=255)

            tagTTTCCATAGTAGCGAAACAATTTGACCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCAACACCTCTTTTATCAACAGATAGTCCAGTTCTTTATTCCACTTGCGAGTTAGGAAGTCCCGa                                                                                                                                          <  2:1714619/148‑1 (MQ=255)

            tCGTTGCCATAGGAGCGAAACCAGTTGCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGa                                                                                                                                          <  2:3752631/150‑1 (MQ=255)

                 gCCATAGGAGCGAAACCAGTTGCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGcttc                                                                                                                                     <  1:7122260/150‑1 (MQ=255)

                  ccATAGGAGCGAAACCAGTTGCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGcttct                                                                                                                                    >  1:2565744/1‑150 (MQ=255)

                   cATAGGAGCGAAACCAGTTGCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGcttctt                                                                                                                                   <  2:1844651/150‑1 (MQ=255)

                   cATAGGAGCGAAACCAGTTGCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGcttctt                                                                                                                                   <  1:4854949/150‑1 (MQ=255)

                     tAGGAGCGAAACCAGTTGCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCa                                                                                                                                 >  1:1155338/1‑150 (MQ=255)

                     tAGGAGCGAAACCAGTTGCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCa                                                                                                                                 >  1:2410153/1‑150 (MQ=255)

                     tAGGAGCGAAACCAGTTGCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCa                                                                                                                                 >  1:1819196/1‑150 (MQ=255)

                     tAGGAGCGAAACCAGTTGCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCa                                                                                                                                 >  2:8114256/1‑150 (MQ=255)

                      aGGAGCGAAACCAGTTGCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCAc                                                                                                                                >  2:8963658/1‑150 (MQ=255)

                              aaCCAGTTGCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGttt                                                                                                                        <  2:7518728/150‑1 (MQ=255)

                                ccAGTTGCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGc                                                                                                                      >  2:8505209/1‑150 (MQ=255)

                                  aGTTGCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCaa                                                                                                                    >  1:4570344/1‑150 (MQ=255)

                                   gTTGCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCaaa                                                                                                                   <  2:9043406/150‑1 (MQ=255)

                                      gCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTc                                                                                                                <  2:6170312/150‑1 (MQ=255)

                                      gCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTc                                                                                                                <  2:3041069/150‑1 (MQ=255)

                                       cccGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCg                                                                                                               >  1:7552385/1‑150 (MQ=255)

                                         cGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGc                                                                                                             >  1:5893172/1‑150 (MQ=255)

                                          gAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCa                                                                                                            >  2:8764398/1‑150 (MQ=255)

                                               gtcgtGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGtt                                                                                                       <  2:7101438/150‑1 (MQ=255)

                                               gtcgtGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGtt                                                                                                       <  2:1478352/150‑1 (MQ=255)

                                               ctcgtGCCACTCGGAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGtt                                                                                                       <  2:8411738/149‑1 (MQ=255)

                                                tcgtGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTa                                                                                                      <  1:5998934/150‑1 (MQ=255)

                                                 cgtGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTAc                                                                                                     <  1:8242526/150‑1 (MQ=255)

                                                 cgtGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTAc                                                                                                     <  1:8356757/150‑1 (MQ=255)

                                                  gtGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACg                                                                                                    <  1:6123952/150‑1 (MQ=255)

                                                      cACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCAt                                                                                                <  1:8883925/150‑1 (MQ=255)

                                                       aCTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCAtt                                                                                               >  1:1906621/1‑150 (MQ=255)

                                                       aCTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCAtt                                                                                               >  1:9079699/1‑150 (MQ=255)

                                                        cTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTg                                                                                              >  1:3071629/1‑150 (MQ=255)

                                                         tCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTgg                                                                                             <  2:3369344/150‑1 (MQ=255)

                                                         tCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTgg                                                                                             <  1:8963699/150‑1 (MQ=255)

                                                         tCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTgg                                                                                             <  1:1565984/150‑1 (MQ=255)

                                                           gTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGgt                                                                                           <  2:3535757/150‑1 (MQ=255)

                                                           gTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGgt                                                                                           <  2:1320565/150‑1 (MQ=255)

                                                                  gaaCGGACCGCGAGGCGGTGGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTGCTTCCACTTGCGAGTTAGTAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCAcc                                                                                    <  2:6912048/149‑1 (MQ=255)

                                                                     cGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGta                                                                                 <  2:1155338/150‑1 (MQ=255)

                                                                      ggACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGtat                                                                                >  2:3539453/1‑150 (MQ=255)

                                                                          cgcgATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGc                                                                            >  1:2821487/1‑150 (MQ=255)

                                                                           gcgATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGcc                                                                           >  2:3884055/1‑150 (MQ=255)

                                                                            cgATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCa                                                                          <  2:4435372/150‑1 (MQ=255)

                                                                            cgATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCa                                                                          <  2:2428689/150‑1 (MQ=255)

                                                                            cgATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCa                                                                          <  1:5281499/150‑1 (MQ=255)

                                                                            cgATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCa                                                                          <  1:6555529/150‑1 (MQ=255)

                                                                            cgATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCa                                                                          >  1:3942210/1‑150 (MQ=255)

                                                                                    ttGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATc                                                                  >  2:71658/1‑150 (MQ=255)

                                                                                    ttGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGGTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCCTTGGGTGTCCACCGTATAGGCCAGGGCAATc                                                                  >  1:5159066/1‑150 (MQ=255)

                                                                                        cGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCATACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCtttt                                                              <  2:2835404/150‑1 (MQ=255)

                                                                                         ggTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTc                                                             <  2:4255422/150‑1 (MQ=255)

                                                                                           tGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCgg                                                           <  1:8414794/150‑1 (MQ=255)

                                                                                           tGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCgg                                                           <  2:5523492/150‑1 (MQ=255)

                                                                                                  ccAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTcgcg                                                    >  1:101337/1‑150 (MQ=255)

                                                                                                     gagCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGt                                                 >  1:2434523/1‑150 (MQ=255)

                                                                                                       gCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTAt                                               >  1:7413657/1‑150 (MQ=255)

                                                                                                           tttGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCa                                           >  1:2888257/1‑150 (MQ=255)

                                                                                                           tttGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCa                                           >  1:5676858/1‑150 (MQ=255)

                                                                                                            ttGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAg                                          >  1:7247732/1‑150 (MQ=255)

                                                                                                                tCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCAt                                      >  2:6299684/1‑150 (MQ=255)

                                                                                                                  aaCCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATcc                                    >  2:8329773/1‑150 (MQ=255)

                                                                                                                  aaCCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATcc                                    >  1:6370372/1‑150 (MQ=255)

                                                                                                                   aCCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCt                                   >  2:6850791/1‑150 (MQ=255)

                                                                                                                   aCCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCATTTCGGGTTCGCGGGTATTCCAGCCATCCt                                   >  2:6152350/1‑150 (MQ=255)

                                                                                                                     cggtactCCAGTTCCTTATTCCACTTGCGGGTCAGGAAGGCTTCGGCTTCGGCACGGTTTGTAACGAATTCAACGCGGTTGCGCCACCGCGTGTCCAAGGTATAGGCGAGAGCGACTTTCGCGGGATCGCGGGTATTCCAGCCATCTTcc                                 <  2:2282083/144‑2 (MQ=255)

                                                                                                                     cGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCa                                 >  1:1774672/1‑150 (MQ=255)

                                                                                                                         aGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGa                             <  2:6587529/150‑1 (MQ=255)

                                                                                                                           tCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGt                           <  1:7800667/150‑1 (MQ=255)

                                                                                                                            ccAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTc                          >  2:5882960/1‑150 (MQ=255)

                                                                                                                                ttctttctTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGAc                      >  1:4389425/1‑150 (MQ=255)

                                                                                                                                   tttcttCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCtt                   <  2:2216847/150‑1 (MQ=255)

                                                                                                                                   tttcttCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCtt                   >  1:2852182/1‑150 (MQ=255)

                                                                                                                                   tttcttCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCtt                   >  2:8635875/1‑150 (MQ=255)

                                                                                                                                   tttcttCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGACAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCtt                   >  2:533153/1‑150 (MQ=255)

                                                                                                                                   cttattcCACTTGCGGGTCAGGAAGGCTTCGGCTTCGGCACGGTTTGTAACGAATTCAACGCGGTTGCGCCACCGCGTGTCCAAGGTATAGGCGAGAGCGACTTTCGCGGGATCGCGGGTATTCCAGCCATCTTCCGCAAGGCGGACCtt                   >  2:1425408/5‑150 (MQ=255)

                                                                                                                                   cttattcCACTTGCGGGTCAGGAAGGCTTCGGCTTCGGCACGGTTTGTAACGAATTCAACGCGGTTGCGCCACCGCGTGTCCAAGGTATAGGCGAGAGCGACTTTCGCGGGATCGCGGGTATTCCAGCCATCTTCCGCAAGGCGGACCtt                   >  2:1619703/5‑150 (MQ=255)

                                                                                                                                   cttattcCACTTGCGGGTCAGGAAGGCTTCGGCTTCGGCACGGTTTGTAACGAATTCAACGCGGTTGCGCCACCGCGTGTCCAAGGTATAGGCGAGAGCGACTTTCGCGGGATCGCGGGTATTCCAGCCATCTTCCGCAAGGCGGACCtt                   >  1:464019/5‑150 (MQ=255)

                                                                                                                                    ttcttcCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGTCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCttt                  <  1:7583285/150‑1 (MQ=255)

                                                                                                                                    ttattcCACTTGCGGGTCAGGAAGGCTTCGGCTTCGGCACGGTTTGTAACGAATTCAACGCGGTTGCGCCACCGCGTGTCCAAGGTATAGGCGAGAGCGACTTTCGCGGGATCGCGGGTATTCCAGCCATCTTCCGCAAGGCGGACCttt                  >  2:7278749/4‑150 (MQ=255)

                                                                                                                                     tattcCACTTGCGGGTCAGGAAGGCTTCGGCTTCGGCACGGTTTGTAACGAATTCAACGCGGTTGCGCCACCGCGTGTCCAAGGTATAGGCGAGAGCGACTTTCGCGGGATCGCGGGTATTCCAGCCATCTTCCGCAAGGCGGACCtttt                 <  2:77904/148‑1 (MQ=255)

                                                                                                                                       ttcCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCa               >  2:9088814/1‑150 (MQ=255)

                                                                                                                                       ttcCACTTGCGAGTTAGGAAGTCCAGAGCTTCTTCACGGTTGGTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGAATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGACATCCTCAGCGAGTCGGACCTTTTCa               >  1:8147628/1‑150 (MQ=255)

                                                                                                                                          cACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATc            >  2:2577482/1‑150 (MQ=255)

                                                                                                                                           aCTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCg           >  2:4880557/1‑150 (MQ=255)

                                                                                                                                               gCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGa       >  2:8825412/1‑150 (MQ=255)

                                                                                                                                                 gAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGAtt     >  2:2354668/1‑150 (MQ=255)

                                                                                                                                                    ttAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTcgc  >  2:5006492/1‑150 (MQ=255)

                                                                                                                                                    ttAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTcgc  <  1:4952041/150‑1 (MQ=255)

AACTCCCAGTTTTCGTTGCCATAGGAGCGAAACCAGTTGCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGC  >  NC_002947/3356836‑3357133

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