BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,169,188	C→A	45.8%	G60C (GGT→TGT)	yqiA ←	acyl CoA esterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,169,188	0	C	A	45.8%	131.9 / 222.6	181	G60C (GGT→TGT)	yqiA	acyl CoA esterase
Reads supporting (aligned to +/- strand): ref base C (45/53); new base A (47/36); total (92/89)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 187 total aligned reads displayed.

TCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTT  >  CP009273/3169046‑3169335

tcGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCAt                                                                                                                                                <  2:3287743/148‑1 (MQ=255)

 cGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATg                                                                                                                                               >  2:514538/1‑148 (MQ=255)

 cGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATg                                                                                                                                               <  2:1287196/148‑1 (MQ=255)

 cGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATg                                                                                                                                               >  1:1102567/1‑148 (MQ=255)

  gTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGt                                                                                                                                              >  1:3138213/1‑148 (MQ=255)

  gTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGt                                                                                                                                              <  1:2280493/148‑1 (MQ=255)

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

    tCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTc                                                                                                                                            <  1:1687820/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)

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

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

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

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

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

                  tCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGAtt                                                                                                                              >  1:1849988/1‑148 (MQ=255)

                       cagcagTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCcagc                                                                                                                         >  1:1645008/1‑148 (MQ=255)

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

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

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

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

                            gTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGctc                                                                                                                    <  2:1977984/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)

                                     ggCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTc                                                                                                           >  2:2017813/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)

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

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

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

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

                                                 ggTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATAcgg                                                                                               <  1:543660/148‑1 (MQ=255)

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

                                                          ccacTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTg                                                                                      <  2:1367733/148‑1 (MQ=255)

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

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

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

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

                                                                  ggCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGa                                                                              <  2:2378752/148‑1 (MQ=255)

                                                                   gcagcaTAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGat                                                                             >  2:2678581/1‑148 (MQ=255)

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

                                                                    cagaaTAAAACATTGTGACAACCAGGTTGCGTAAAATCCCAACAGAGACGAACCCACAAAACCCAACGAAACACCGCCATGTACCACTACAATCTATTACATCAGCTCTGAAGCGACTGAAACAAACTGCCGCAACTCAGGaaaaatc                                                                            >  2:1050739/1‑143 (MQ=255)

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

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

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

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

                                                                          aaaaCATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCa                                                                      >  1:592289/1‑148 (MQ=255)

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

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

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

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

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

                                                                                 tgtgACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAg                                                               >  2:1035539/1‑148 (MQ=255)

                                                                                    gACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGt                                                            >  1:1978711/1‑148 (MQ=255)

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

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

                                                                                        acgaGGTGGCGTAATCTCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCGATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCATGGTGATg                                                        <  2:2568577/145‑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)

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

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

                                                                                             gTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCg                                                   >  2:1228014/1‑148 (MQ=255)

                                                                                             gTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCg                                                   >  2:2805681/1‑148 (MQ=255)

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

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

                                                                                                 cGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCgccag                                               <  2:1637045/148‑1 (MQ=255)

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

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

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

                                                                                                      tatCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGt                                          >  2:2393481/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)

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

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

                                                                                                            cccAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCGGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTa                                    <  1:1378923/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)

                                                                                                                       aCCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCt                         <  2:582597/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                >  1:157586/1‑148 (MQ=255)

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

                                                                                                                                    aGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGcgcg            >  2:2793910/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)

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

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

TCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTT  >  CP009273/3169046‑3169335

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