BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP009974	5,513,866	G→A	A266V (GCA→GTA)	RPPX_24600 ←	integrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009974	5,513,866	0	G	A	99.5%	633.5 / NA	187	A266V (GCA→GTA)	RPPX_24600	integrase
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (84/102); minor base C (0/1); total (84/103)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 187 total aligned reads displayed.

GTCCTGGTTCGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTGCAGCCTGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGAAAAACTCAGA  >  CP009974/5513718‑5514008

 tCCTGGTTCGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAg                                                                                                                                             <  2:2967006/151‑1 (MQ=255)

  ccTGGTTCGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGc                                                                                                                                            >  2:813050/1‑151 (MQ=255)

  ccTGGTTCGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGc                                                                                                                                            >  1:4400658/1‑151 (MQ=255)

    tGGTTCGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGccc                                                                                                                                          <  1:1120664/151‑2 (MQ=255)

      gTTCGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGa                                                                                                                                        >  2:3899777/1‑151 (MQ=255)

           agaCTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCt                                                                                                                                   <  1:3726856/151‑1 (MQ=255)

           agaCTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCt                                                                                                                                   <  1:3529004/151‑1 (MQ=255)

             aCTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAGTTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGt                                                                                                                                 <  2:5090597/151‑1 (MQ=255)

             aCTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGt                                                                                                                                 >  1:5332714/1‑151 (MQ=255)

             aCTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGt                                                                                                                                 >  2:2567546/1‑151 (MQ=255)

               ttgtGGCTCGGGCAGGAGTGGCGGGCGGGTATCAAGACGATGGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTAt                                                                                                                               <  2:4429469/148‑1 (MQ=37)

               ttATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTAt                                                                                                                               <  1:1131818/151‑1 (MQ=255)

                tATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTAtt                                                                                                                              <  1:2423656/151‑1 (MQ=255)

                   ggcggaggcaggAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGGTCGCCCTCGCGCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCca                                                                                                                           <  1:2096223/146‑1 (MQ=37)

                   ggCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCca                                                                                                                           >  2:2174131/1‑151 (MQ=255)

                      tCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAg                                                                                                                        >  2:1520548/1‑151 (MQ=255)

                      tCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAg                                                                                                                        <  2:3618218/151‑1 (MQ=255)

                        aggcaggAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAggg                                                                                                                      >  2:229309/1‑151 (MQ=255)

                        aggcaggAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAggg                                                                                                                      >  2:425891/1‑151 (MQ=255)

                          gcaggGGTCGCAGGCGGGTATCAATACGATTGAATGGACCGGCGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCt                                                                                                                    <  1:4684753/151‑1 (MQ=25)

                           caggAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTc                                                                                                                   >  2:4628180/1‑151 (MQ=255)

                            aggAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCa                                                                                                                  >  2:1737856/1‑151 (MQ=255)

                            aggAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCCCa                                                                                                                  >  2:3000520/1‑151 (MQ=255)

                             ggAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAg                                                                                                                 <  1:1553634/151‑1 (MQ=255)

                              gAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAgg                                                                                                                >  1:4831041/1‑151 (MQ=255)

                                 tGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCa                                                                                                             <  1:2040182/151‑1 (MQ=255)

                                  ggCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCAt                                                                                                            >  1:5124876/1‑151 (MQ=255)

                                         gggTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGtct                                                                                                     <  1:1505555/151‑1 (MQ=255)

                                           gTATCAATACGATTGAATGACCCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCa                                                                                                   <  2:3376589/151‑1 (MQ=255)

                                           gTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCa                                                                                                   <  2:545228/151‑1 (MQ=255)

                                           gTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCa                                                                                                   >  1:2670475/1‑151 (MQ=255)

                                              tCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCt                                                                                                <  2:3042659/151‑1 (MQ=255)

                                               cAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCtt                                                                                               >  2:1235540/1‑151 (MQ=255)

                                                   aCGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATg                                                                                           >  1:4907860/1‑151 (MQ=255)

                                                     gATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCg                                                                                         <  1:376894/151‑1 (MQ=255)

                                                        tgaatgaaCCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATc                                                                                      <  1:1736989/151‑1 (MQ=255)

                                                        tgaatgaaCCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATc                                                                                      <  2:1339487/151‑1 (MQ=255)

                                                          aatgaaCCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCAt                                                                                    <  1:1234849/151‑1 (MQ=255)

                                                          aatgaaCCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCCCAGGCCATTGAGTCTCAGCTTGATGCGATCAt                                                                                    >  1:1383100/1‑151 (MQ=255)

                                                            tgaaCCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCACTGAGTCCCAGCTTGATGCGATCATgg                                                                                  >  2:2253181/1‑151 (MQ=255)

                                                              aaCCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGtt                                                                                <  1:3533903/151‑1 (MQ=255)

                                                              aaCCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGtt                                                                                <  1:676515/151‑1 (MQ=255)

                                                                   gggCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTgtagt                                                                           <  2:2376002/151‑1 (MQ=255)

                                                                    ggCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTgtagta                                                                          <  2:2190176/151‑1 (MQ=255)

                                                                     gccgccTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGtagtag                                                                         <  2:4140085/151‑1 (MQ=255)

                                                                      ccgccTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTagtagt                                                                        <  1:1182498/151‑1 (MQ=255)

                                                                        gccTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGa                                                                      >  2:2878327/1‑151 (MQ=255)

                                                                        gccTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGa                                                                      >  1:5333362/1‑151 (MQ=255)

                                                                        gccTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGa                                                                      <  1:4493847/151‑1 (MQ=255)

                                                                         ccTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGaa                                                                     <  1:2245132/151‑1 (MQ=255)

                                                                         ccTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGaa                                                                     >  1:3986673/1‑151 (MQ=255)

                                                                          cTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAAt                                                                    >  2:1010898/1‑151 (MQ=255)

                                                                            tGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGt                                                                  <  2:2672990/151‑1 (MQ=255)

                                                                            tGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGt                                                                  <  1:190771/151‑1 (MQ=255)

                                                                            tGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGt                                                                  <  1:3334249/151‑1 (MQ=255)

                                                                              cGGCCGTTCGCGGGCTCGCCCTCGGTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTAc                                                                <  2:5009721/151‑1 (MQ=255)

                                                                                gCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGTCAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTc                                                              <  1:3193404/151‑1 (MQ=255)

                                                                                gCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGGTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTc                                                              <  2:5168596/151‑1 (MQ=255)

                                                                                 cccTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGa                                                                                        <  1:904886/124‑1 (MQ=255)

                                                                                 cccTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGa                                                                                        >  2:904886/1‑124 (MQ=255)

                                                                                    ttCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATcc                                                          <  1:1068871/151‑1 (MQ=255)

                                                                                       ccaggCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGa                                                       <  1:301487/148‑1 (MQ=255)

                                                                                        cgGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGAc                                                      <  1:1308554/151‑1 (MQ=255)

                                                                                        cgGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGAc                                                      <  1:1307773/151‑1 (MQ=255)

                                                                                         gggCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGAcc                                                     >  1:2053412/1‑151 (MQ=255)

                                                                                          ggCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGAccc                                                    >  1:3858881/1‑151 (MQ=255)

                                                                                           gCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCg                                                   <  1:662717/151‑1 (MQ=255)

                                                                                             tCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGct                                                 >  1:832037/1‑151 (MQ=255)

                                                                                             tCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGct                                                 <  2:3986673/151‑1 (MQ=255)

                                                                                              cGCCCTCGCTCCCACAATTTGGCCTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGctc                                                <  1:537437/151‑1 (MQ=255)

                                                                                                 ccTCGCTCCCATAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCa                                             <  1:4968509/151‑1 (MQ=255)

                                                                                                 ccTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCa                                             <  1:529797/151‑1 (MQ=255)

                                                                                                   tCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGt                                                                  <  2:1219889/128‑1 (MQ=255)

                                                                                                   tCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGt                                                                  >  1:1219889/1‑128 (MQ=255)

                                                                                                     gCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCtc                                         >  2:3362726/1‑151 (MQ=255)

                                                                                                     ccTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGTCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCTGCtc                                         <  1:2749432/150‑1 (MQ=255)

                                                                                                      cTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCtca                                        >  1:2009331/1‑151 (MQ=255)

                                                                                                      cTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCtca                                        >  1:4544400/1‑151 (MQ=255)

                                                                                                       tCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCtcat                                       <  1:2866/151‑1 (MQ=255)

                                                                                                         ccACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCtcatca                                     >  1:5331482/1‑151 (MQ=255)

                                                                                                          cacaATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCaa                                    >  2:1545146/1‑151 (MQ=255)

                                                                                                              aTTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCCCAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCACCCAGACCCCCTCTCAGCTCACCAATGcc                                >  2:1061291/1‑150 (MQ=11)

                                                                                                                      tGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAAc                        <  2:3166609/151‑1 (MQ=12)

                                                                                                                      tGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAAc                        <  2:4116542/151‑1 (MQ=9)

                                                                                                                      tGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAAc                        >  2:4711543/1‑151 (MQ=12)

                                                                                                                       gCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACg                       >  1:3630212/1‑151 (MQ=12)

                                                                                                                       gCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACg                       <  2:3289899/151‑1 (MQ=14)

                                                                                                                         cccccAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCt                     <  1:4601781/151‑1 (MQ=12)

                                                                                                                             cAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAg                 <  1:425891/151‑1 (MQ=17)

                                                                                                                              aaaaGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAgg                >  2:511505/1‑151 (MQ=17)

                                                                                                                              aaaaGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAgg                >  2:4486554/1‑151 (MQ=12)

                                                                                                                                aaGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTa              >  1:2915044/1‑151 (MQ=12)

                                                                                                                                aaGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTa              <  2:2086299/151‑1 (MQ=14)

                                                                                                                                   ttGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGaa           <  1:1060976/151‑1 (MQ=12)

                                                                                                                                    tGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGaaa          <  1:1430395/151‑1 (MQ=17)

                                                                                                                                     ggAAAGTACTAAGCTACAGCCAGATCCCTGTGGTGCACAGGGCTCAGGCCATGGAGTCCCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTAATCCAGTCCCGCTCTCAGCTCATCAATGCTCTCGAAACGGCTCAGGTAGaaaa         <  2:4392461/151‑1 (MQ=11)

                                                                                                                                       aCAGATCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGAAAAAc       <  2:498374/151‑1 (MQ=9)

                                                                                                                                        cAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGAAAAACt      <  2:1450865/151‑1 (MQ=14)

                                                                                                                                         aGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCACTGCTCTCGAAACGCTTCAGGTAGAAAAACTc     >  2:3551180/1‑151 (MQ=9)

                                                                                                                                            tCTAAGCTCCAGCCGGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGAAAAACTCAGa  <  2:4052886/151‑1 (MQ=14)

GTCCTGGTTCGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTGCAGCCTGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGAAAAACTCAGA  >  CP009974/5513718‑5514008

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