BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,737,983	C→A	5.1%	G481C (GGC→TGC) ‡	hutU ←	urocanate hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,737,983	0	C	A	5.1%	519.9 / 11.5	155	G481C (GGC→TGC) ‡	hutU	urocanate hydratase
Reads supporting (aligned to +/- strand): ref base C (67/80); new base A (5/3); total (72/83)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 164 total aligned reads displayed.

GACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCG  >  NC_002947/5737834‑5738131

   gCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGc                                                                                                                                                   >  2:3229312/1‑150 (MQ=255)

    ccAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCg                                                                                                                                                  <  1:946887/150‑1 (MQ=255)

    ccAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGAGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCg                                                                                                                                                  >  1:144064/1‑150 (MQ=255)

     cAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCgg                                                                                                                                                 >  1:824497/1‑150 (MQ=255)

      aGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGgt                                                                                                                                                <  1:3570347/150‑1 (MQ=255)

          ccTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTc                                                                                                                                            <  1:564785/150‑1 (MQ=255)

          ccTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTc                                                                                                                                            <  1:564798/150‑1 (MQ=255)

           cTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTca                                                                                                                                           >  2:4120791/1‑150 (MQ=255)

           cTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTca                                                                                                                                           >  2:2178704/1‑150 (MQ=255)

             gggTCGTTGGTGAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGTCGATGACCTTGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTcagc                                                                                                                                         <  2:2844542/150‑1 (MQ=255)

             gggTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTcagc                                                                                                                                         <  2:2288288/150‑1 (MQ=255)

             gggTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTcagc                                                                                                                                         >  1:3258607/1‑150 (MQ=255)

               gTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTcagcag                                                                                                                                       >  2:3429289/1‑150 (MQ=255)

                tCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAgg                                                                                                                                      <  2:3294865/150‑1 (MQ=255)

                  gTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGc                                                                                                                                    >  1:3774241/1‑150 (MQ=255)

                    tGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGt                                                                                                                                  <  2:1098507/150‑1 (MQ=255)

                       tCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGa                                                                                                                               <  2:3266998/150‑1 (MQ=255)

                        cAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGag                                                                                                                              >  1:4665054/1‑150 (MQ=255)

                         aGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGagc                                                                                                                             >  1:2829432/1‑150 (MQ=255)

                            aCCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGagcagc                                                                                                                          <  2:4382334/150‑1 (MQ=255)

                            aCCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGagcagc                                                                                                                          <  2:5280665/150‑1 (MQ=255)

                                    gATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTc                                                                                                                  >  1:886991/1‑150 (MQ=255)

                                     aTGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCg                                                                                                                 <  2:5150376/150‑1 (MQ=255)

                                       gcgcTGGTCGGCATCTTCCGAACCGTCGCAGCAGATGACCATGCCCGAATGTTGGGAGAAGCCCATGCCGACGCCGCCACCATGATGGAGCGATACCCATGTCGCACCCGATGCGGTGTTGAGCAACGCATTGAGCAGCGGCCAGTCGGa                                                                                                               <  1:3038833/150‑1 (MQ=255)

                                        cgcTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTGGAGCAGCGGCCAGTCGGac                                                                                                              >  2:4556017/1‑149 (MQ=255)

                                         gcTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGaaa                                                                                                             >  2:843330/1‑150 (MQ=255)

                                          cTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAAc                                                                                                            <  1:1201244/150‑1 (MQ=255)

                                           tCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACa                                                                                                           <  1:3696164/150‑1 (MQ=255)

                                           tCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACa                                                                                                           >  1:773737/1‑150 (MQ=255)

                                               cggcCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGt                                                                                                       <  1:4238729/150‑1 (MQ=255)

                                                ggcCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTc                                                                                                      <  1:4214664/150‑1 (MQ=255)

                                                 gcCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTcc                                                                                                     <  2:4819430/150‑1 (MQ=255)

                                                   cTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGa                                                                                                   <  2:2034498/150‑1 (MQ=255)

                                                   cTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGa                                                                                                   <  2:4772313/150‑1 (MQ=255)

                                                   cTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGa                                                                                                   >  1:4482510/1‑150 (MQ=255)

                                                    tCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGTCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGaa                                                                                                  >  2:3706659/1‑150 (MQ=255)

                                                    tCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGaa                                                                                                  >  2:4240856/1‑150 (MQ=255)

                                                       tCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCa                                                                                               <  2:1777024/150‑1 (MQ=255)

                                                        cGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCAt                                                                                              >  1:5239833/1‑150 (MQ=255)

                                                         ggTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATc                                                                                             <  2:5375208/150‑1 (MQ=255)

                                                          gaaCCGTCGCAGCAGATGACCATGCCCGAATGTTGGGAGTAGCCCATGCCGACGCCGCCACCATGATGGAGCGATACCCATGTCGCACCCGATGCGGTGTTGAGCAACGCATTGAGCAGCGGCCAGTCGGAGACTGCGTCCGAGCCGTct                                                                                            >  2:4666671/3‑149 (MQ=255)

                                                          gTACCGTCGCAGACGATGACCATGCCCGAGTGCAGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCa                                                                                            <  2:3407034/150‑1 (MQ=255)

                                                              cGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCa                                                                                        <  1:2084740/150‑1 (MQ=255)

                                                                tCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATg                                                                                      <  2:4230949/150‑1 (MQ=255)

                                                                 cGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCg                                                                                                                                   >  1:3915655/1‑104 (MQ=255)

                                                                 cGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCg                                                                                                                                   <  2:3915655/104‑1 (MQ=255)

                                                                 cGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATgg                                                                                     <  2:3216400/150‑1 (MQ=255)

                                                                      cagATGACCATGCCCGAATGTTGGGAGAAGCCCATGCCGACGCCGCCACCATGATGGAGCGATACCCATGTCGCACCCGATGCGGTGTTGAGCAACGCATTGAGCAGCGGCCAGTCGGAGACTGCGTCCGAGCCGTcttt                                                                                          <  2:2569331/138‑4 (MQ=255)

                                                                        gATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGt                                                                              <  2:3342157/150‑1 (MQ=255)

                                                                             ccATGCCCGAATGTTGGGAGAAGCCCATGCCGACGCCGCCACCATGATGGAGCGATACCCATGTCGCACCCGATGCGGTGTTGAGCAACGCATTGAGCAGCGGCCAGTCGGAGACTGCGTCCGAGCCGTCTTTCATGGATTCGGTTTcac                                                                         >  1:3207414/1‑148 (MQ=255)

                                                                              cATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTcgcg                                                                        <  2:3576568/150‑1 (MQ=255)

                                                                              cATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTcgcg                                                                        <  2:1065557/150‑1 (MQ=255)

                                                                              cATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTcgcg                                                                        <  2:3218485/150‑1 (MQ=255)

                                                                              cATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTcgcg                                                                        >  1:3228503/1‑150 (MQ=255)

                                                                              cATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTcgcg                                                                        >  1:3228499/1‑150 (MQ=255)

                                                                                 gCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGtt                                                                     <  1:1278841/150‑1 (MQ=255)

                                                                                 gCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGtt                                                                     <  1:2193900/150‑1 (MQ=255)

                                                                                   ccGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCg                                                                   >  1:4386620/1‑150 (MQ=255)

                                                                                   ccGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCg                                                                   <  1:2525114/150‑1 (MQ=255)

                                                                                     gAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCggg                                                                 <  1:1161176/150‑1 (MQ=255)

                                                                                      aGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGc                                                                <  1:4029928/150‑1 (MQ=255)

                                                                                           tgcGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGAt                                                           >  1:2082210/1‑150 (MQ=255)

                                                                                           tgcGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGAt                                                           >  1:2868255/1‑150 (MQ=255)

                                                                                             cGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATAc                                                         <  1:2537453/150‑1 (MQ=255)

                                                                                             cGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATAc                                                         <  1:3862453/150‑1 (MQ=255)

                                                                                             cGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATAc                                                         <  1:2537654/150‑1 (MQ=255)

                                                                                              gagaAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGGGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATAcc                                                        <  1:3449227/150‑1 (MQ=255)

                                                                                               agaAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCg                                                       >  1:4722035/1‑150 (MQ=255)

                                                                                               agaAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCg                                                       >  1:2572102/1‑150 (MQ=255)

                                                                                                     ccATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCgg                                                 >  1:2047648/1‑150 (MQ=255)

                                                                                                     ccATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCgg                                                 >  2:1939091/1‑150 (MQ=255)

                                                                                                      cATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGa                                                >  1:2704778/1‑150 (MQ=255)

                                                                                                       aTGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAg                                               >  2:3063647/1‑150 (MQ=255)

                                                                                                           ccACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGa                                           <  2:3989838/150‑1 (MQ=255)

                                                                                                           acACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGTGCAGCGGCCAGTCGGACACAGCGTCCGAACCATCCCGCATGGCCTCGGGTTCGGGGTTCGGGCGGGATTGCGTGCCGGTGTCgg                                           >  2:3962883/2‑149 (MQ=255)

                                                                                                             aCGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAgg                                         >  2:3708379/1‑150 (MQ=255)

                                                                                                              cGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGTGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAggt                                        <  1:412726/150‑1 (MQ=255)

                                                                                                               gCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAggtg                                       >  2:4878646/1‑150 (MQ=255)

                                                                                                               gCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAggtg                                       <  1:1754748/150‑1 (MQ=255)

                                                                                                                   ccGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCa                                   >  2:2448987/1‑150 (MQ=255)

                                                                                                                     gCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGGAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACg                                 >  1:1040255/1‑150 (MQ=255)

                                                                                                                            tGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATc                          <  2:4337675/150‑1 (MQ=255)

                                                                                                                            tGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATc                          <  1:4942476/150‑1 (MQ=255)

                                                                                                                             gCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATca                         <  2:5344781/150‑1 (MQ=255)

                                                                                                                              cAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATcac                        <  2:5118700/150‑1 (MQ=255)

                                                                                                                               aGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTcg                                                                          >  1:1670095/1‑99 (MQ=255)

                                                                                                                               aGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATcacc                       >  2:924061/1‑150 (MQ=255)

                                                                                                                               aGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATcacc                       <  2:2557749/150‑1 (MQ=255)

                                                                                                                                     aCCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGt                 <  1:4976142/150‑1 (MQ=255)

                                                                                                                                       ccAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGc               >  1:1460337/1‑150 (MQ=255)

                                                                                                                                        cAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGcc              <  2:805063/150‑1 (MQ=255)

                                                                                                                                          ggTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGa            <  2:4665054/150‑1 (MQ=255)

                                                                                                                                            tGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACa          <  2:1341183/150‑1 (MQ=255)

                                                                                                                                             ggCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAg         <  2:1409668/150‑1 (MQ=255)

                                                                                                                                              gcgcCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGTAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCa                                   <  2:762044/123‑1 (MQ=255)

                                                                                                                                                gccgccTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTc      >  1:5180529/1‑150 (MQ=255)

                                                                                                                                                gccgccTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTc      >  1:1006310/1‑150 (MQ=255)

                                                                                                                                                gccgccTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTc      >  2:3020100/1‑150 (MQ=255)

                                                                                                                                                    ccTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCTGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTcgccg  <  1:4233215/150‑1 (MQ=255)

                                                                                                                                                    ccTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTcgccg  >  2:2287480/1‑150 (MQ=255)

GACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCG  >  NC_002947/5737834‑5738131

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