BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,169,196	T→A	17.9%	E57V (GAA→GTA)	yqiA ←	acyl CoA esterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,169,196	0	T	A	17.9%	481.2 / 92.8	190	E57V (GAA→GTA)	yqiA	acyl CoA esterase
Reads supporting (aligned to +/- strand): ref base T (78/78); new base A (12/22); total (90/100)
Fisher's exact test for biased strand distribution p-value = 1.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 191 total aligned reads displayed.

TTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAACC  >  CP009273/3169049‑3169341

ttCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGtt                                                                                                                                                   >  1:48086/1‑148 (MQ=255)

 tCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTc                                                                                                                                                  <  1:1687820/148‑1 (MQ=255)

 tCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTc                                                                                                                                                  <  2:2405041/148‑1 (MQ=255)

     aCCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAgg                                                                                                                                              <  1:501208/148‑1 (MQ=255)

     aCCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAgg                                                                                                                                              <  1:2758230/148‑1 (MQ=255)

       cGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGAc                                                                                                                                            >  1:2609791/1‑148 (MQ=255)

        gagaTAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACa                                                                                                                                           >  2:1012210/1‑148 (MQ=255)

          gaTAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAAt                                                                                                                                         <  2:3361234/148‑1 (MQ=255)

           aTAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACGACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATg                                                                                                                                        <  2:963642/148‑1 (MQ=255)

           aTAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATg                                                                                                                                        >  2:2207051/1‑148 (MQ=255)

           aTAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATg                                                                                                                                        >  2:1131153/1‑148 (MQ=255)

           aTAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATg                                                                                                                                        >  2:1305283/1‑148 (MQ=255)

            tAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATgg                                                                                                                                       <  1:3162747/148‑1 (MQ=255)

             aGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGa                                                                                                                                      <  2:2853028/148‑1 (MQ=255)

                ccGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTc                                                                                                                                   >  1:582597/1‑148 (MQ=255)

                      gcagTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCcagcag                                                                                                                             <  1:711443/148‑1 (MQ=255)

                       cagTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCagcagc                                                                                                                            <  2:1494175/148‑1 (MQ=255)

                              aaCGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCcg                                                                                                                     <  1:2469231/148‑1 (MQ=255)

                               aCGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCcgc                                                                                                                    >  2:2932221/1‑148 (MQ=255)

                                  ggCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTc                                                                                                                 >  2:1804055/1‑148 (MQ=255)

                                      cACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGaa                                                                                                             >  2:880670/1‑148 (MQ=255)

                                        ccgccgGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAgg                                                                                                           >  1:1343358/1‑148 (MQ=255)

                                        ccgccgGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAgg                                                                                                           >  2:834286/1‑148 (MQ=255)

                                        ccgccgGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAgg                                                                                                           >  2:1633381/1‑148 (MQ=255)

                                         cgccgGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGa                                                                                                          <  2:770411/148‑1 (MQ=255)

                                           ccgGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATa                                                                                                        >  2:342579/1‑148 (MQ=255)

                                            cgGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATAc                                                                                                       <  2:2713475/148‑1 (MQ=255)

                                                     caccacTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAAc                                                                                              <  1:1756237/148‑1 (MQ=255)

                                                      accacTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACt                                                                                             >  2:1855790/1‑148 (MQ=255)

                                                           tGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGa                                                                                        <  1:2816192/148‑1 (MQ=255)

                                                           tGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGa                                                                                        >  2:2876706/1‑148 (MQ=255)

                                                           tGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGa                                                                                        >  2:447507/1‑148 (MQ=255)

                                                           tGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGa                                                                                        >  2:447573/1‑148 (MQ=255)

                                                              gggCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATg                                                                                     <  1:2199090/148‑1 (MQ=255)

                                                                gcagcaTAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGat                                                                                   <  1:1012210/148‑1 (MQ=255)

                                                                    caTAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGatcatc                                                                               >  2:1060923/1‑148 (MQ=255)

                                                                    caTAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGatcatc                                                                               >  1:1722463/1‑148 (MQ=255)

                                                                     aTAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCAtct                                                                              <  1:186884/148‑1 (MQ=255)

                                                                       aaaaCATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCa                                                                            >  2:2942383/1‑148 (MQ=255)

                                                                       aaaaCATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCa                                                                            <  2:3449134/148‑1 (MQ=255)

                                                                        aaaCATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCaa                                                                           <  2:2995790/148‑1 (MQ=255)

                                                                        aaaCATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCaa                                                                           >  2:1040814/1‑148 (MQ=255)

                                                                         aaCATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAAc                                                                          <  1:636886/148‑1 (MQ=255)

                                                                          aCATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACg                                                                         >  1:2405487/1‑148 (MQ=255)

                                                                           cATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGt                                                                        <  2:2112358/148‑1 (MQ=255)

                                                                           cATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGt                                                                        >  2:3184098/1‑148 (MQ=255)

                                                                           cATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGt                                                                        <  1:2928231/148‑1 (MQ=255)

                                                                                 gACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGt                                                                  <  1:326694/148‑1 (MQ=255)

                                                                                  aCAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTg                                                                 >  1:2962242/1‑148 (MQ=255)

                                                                                     aCCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATg                                                              <  2:2357914/148‑1 (MQ=255)

                                                                                      ccAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGt                                                             >  2:524303/1‑148 (MQ=255)

                                                                                        aGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTc                                                           <  1:2715906/148‑1 (MQ=255)

                                                                                        aGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTc                                                           <  1:2718572/148‑1 (MQ=255)

                                                                                        aGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATNCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTc                                                           <  1:1805652/148‑1 (MQ=255)

                                                                                        aGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATNCCAGCAGCTCGGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTc                                                           <  1:1805667/148‑1 (MQ=255)

                                                                                         ggTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGACTGATCATCTCAACGTCAGGGTGATGTTcc                                                          >  2:1326770/1‑148 (MQ=255)

                                                                                         ggTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTcc                                                          >  2:2507867/1‑148 (MQ=255)

                                                                                          gTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCg                                                         >  1:2430554/1‑148 (MQ=255)

                                                                                            ggCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCgcc                                                       <  2:2568509/148‑1 (MQ=255)

                                                                                              cGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCgccag                                                     >  2:728026/1‑148 (MQ=255)

                                                                                                tAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTTCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCgccagcc                                                   >  1:266802/1‑148 (MQ=255)

                                                                                                 aaTATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCgccagcca                                                  >  2:143381/1‑148 (MQ=255)

                                                                                                 aaTATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCgccagcca                                                  <  2:2641233/148‑1 (MQ=255)

                                                                                                   tatCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCGGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGt                                                <  1:788238/148‑1 (MQ=255)

                                                                                                   tatCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGt                                                >  2:2393460/1‑148 (MQ=255)

                                                                                                    atCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGtt                                               <  1:514593/148‑1 (MQ=255)

                                                                                                     tCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGttt                                              >  2:3079477/1‑148 (MQ=255)

                                                                                                      ccccccAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGtttt                                             <  1:279925/148‑1 (MQ=255)

                                                                                                      ccccccAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACTATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGtttt                                             <  2:1697946/148‑1 (MQ=255)

                                                                                                       cccccAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGttttt                                            >  2:1688177/1‑148 (MQ=255)

                                                                                                       cccccAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGttttt                                            >  2:2805665/1‑148 (MQ=255)

                                                                                                       cccccAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGttttt                                            <  2:1689786/148‑1 (MQ=255)

                                                                                                         cccAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTa                                          >  2:706439/1‑148 (MQ=255)

                                                                                                         cccAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCGGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTa                                          <  1:1378923/148‑1 (MQ=255)

                                                                                                         cccAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTa                                          <  2:1978711/148‑1 (MQ=255)

                                                                                                          ccAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTaa                                         <  2:2609791/148‑1 (MQ=255)

                                                                                                               gACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGc                                    <  1:2335623/148‑1 (MQ=255)

                                                                                                                aCGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCt                                   >  1:339065/1‑148 (MQ=255)

                                                                                                                  gAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCg                                 <  2:3033661/148‑1 (MQ=255)

                                                                                                                    aCCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCt                               >  1:386378/1‑148 (MQ=255)

                                                                                                                        aCAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTg                           <  1:1462265/148‑1 (MQ=255)

                                                                                                                         cAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGc                          >  1:330738/1‑148 (MQ=255)

                                                                                                                          aaTACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCa                         <  1:1071959/148‑1 (MQ=255)

                                                                                                                           aTACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCag                        >  1:2212593/1‑148 (MQ=255)

                                                                                                                            tACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCGGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCaga                       <  1:1010238/148‑1 (MQ=255)

                                                                                                                             aCCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCagag                      <  2:374758/148‑1 (MQ=255)

                                                                                                                              cccAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAgc                     >  2:1029817/1‑148 (MQ=255)

                                                                                                                                 aGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGcgcg                  >  2:2793906/1‑148 (MQ=255)

                                                                                                                                  gCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCgcgg                 >  2:1072469/1‑148 (MQ=255)

                                                                                                                                      aTCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAg             >  2:560535/1‑148 (MQ=255)

                                                                                                                                          ccGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGt         <  2:581491/148‑1 (MQ=255)

                                                                                                                                          cagCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGt         <  2:2309175/146‑1 (MQ=255)

                                                                                                                                           cGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGtt        <  2:2407103/148‑1 (MQ=255)

                                                                                                                                           cGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGtt        <  2:2407158/148‑1 (MQ=255)

                                                                                                                                           cGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGtt        >  1:1636658/1‑148 (MQ=255)

                                                                                                                                            gCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTg       >  2:614736/1‑148 (MQ=255)

                                                                                                                                            gCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTg       <  1:3354217/148‑1 (MQ=255)

                                                                                                                                             ccATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGa      >  1:344107/1‑148 (MQ=255)

                                                                                                                                             ccATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGa      >  1:1116781/1‑148 (MQ=255)

                                                                                                                                              cATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGaa     <  1:630832/148‑1 (MQ=255)

TTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAACC  >  CP009273/3169049‑3169341

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