BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP009974	5,513,872	T→C	Q264R (CAG→CGG)	RPPX_24600 ←	integrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009974	5,513,872	0	T	C	99.6%	820.2 / ‑5.4	261	Q264R (CAG→CGG)	RPPX_24600	integrase
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (134/126); total (134/127)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 261 total aligned reads displayed.

TTCGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTGCAGCCTGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGAAAAACTCAGA  >  CP009974/5513725‑5514008

  cGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGGCTGCCCCCCAAAAGTTGGACACTTCTATGCTACAGCCCGAGcc                                                                                                                                     >  1:4420187/1‑151 (MQ=38)

   gagaCTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGccc                                                                                                                                    >  2:3600713/1‑151 (MQ=255)

   gagaCTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGccc                                                                                                                                    <  1:2147822/151‑1 (MQ=255)

   gagaCTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGccc                                                                                                                                    >  1:5540840/1‑151 (MQ=255)

      aCTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGt                                                                                                                                 <  1:898469/151‑1 (MQ=255)

       cTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTa                                                                                                                                <  2:1108840/151‑1 (MQ=255)

          aTGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTc                                                                                                                             <  1:2528139/151‑1 (MQ=255)

          aTGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGCTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTc                                                                                                                             >  2:4816034/1‑151 (MQ=255)

            ggCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTAtt                                                                                                                              <  1:1242711/148‑1 (MQ=255)

            ggCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTAtt                                                                                                                              >  2:1242711/1‑148 (MQ=255)

            ggCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCCCGTCCTCGCTCCCACAAGTTGGACTGCCCCCCAAAAGTTTGACAGTTCTAAGTTACAGCCCGAGCCCTGTATTAca                                                                                                                           >  1:2227570/1‑151 (MQ=255)

             gCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCcac                                                                                                                          >  2:3128078/1‑151 (MQ=255)

             gCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCcac                                                                                                                          >  1:2338514/1‑151 (MQ=255)

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

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

                caggcaggAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAgg                                                                                                                       <  1:2916964/151‑1 (MQ=255)

                caggcaggAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAgg                                                                                                                       >  1:568435/1‑151 (MQ=255)

                   gcaggAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCGCTCGCTCCCACAACTTGGACTGCCCCCCAAAAGTTGGACGGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCt                                                                                                                    >  1:593884/1‑151 (MQ=255)

                   gcaggAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCCCCCACAATTTGGACTGCCCCCCAAAAGTCGCACAGTTCTAAGCTACAGCCCGAGCCCTCTATTCCACAGGGcc                                                                                                                    >  2:508457/1‑150 (MQ=255)

                     aggAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCa                                                                                                                  <  1:3773546/151‑1 (MQ=255)

                        aGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGACCCCTGTATTCCACAGGGCTCAGGc                                                                                                               >  1:1414057/1‑151 (MQ=255)

                         gtcgCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGcc                                                                                                              <  2:4570232/148‑1 (MQ=255)

                          tGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCa                                                                                                             >  2:2274765/1‑151 (MQ=255)

                          tGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCa                                                                                                             >  1:1061727/1‑151 (MQ=255)

                           ggCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCAt                                                                                                            <  2:6338969/151‑1 (MQ=255)

                               ggCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAg                                                                                                        >  1:3325748/1‑151 (MQ=255)

                               ggCGGGTATCAATACGATGGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAg                                                                                                        <  2:1722575/151‑1 (MQ=255)

                                 cGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGtc                                                                                                      <  2:2443449/151‑1 (MQ=255)

                                  gggTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGtct                                                                                                     >  1:5330069/1‑151 (MQ=255)

                                    gTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGtctcc                                                                                                   >  2:4708099/1‑150 (MQ=255)

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

                                     tATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAg                                                                                                  >  1:5765132/1‑151 (MQ=255)

                                     tATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAg                                                                                                  >  1:3404321/1‑151 (MQ=255)

                                      aTCAATACGATTGAATGAACCGGGGCGGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGc                                                                                                 <  2:4527238/151‑1 (MQ=255)

                                      aTCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGc                                                                                                 >  2:1544422/1‑151 (MQ=255)

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

                                        cAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCtt                                                                                               <  1:258511/151‑1 (MQ=255)

                                         aaTACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTg                                                                                              <  1:3474352/151‑1 (MQ=255)

                                         aaTACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTg                                                                                              >  1:2901240/1‑151 (MQ=255)

                                           tACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGAt                                                                                            >  2:6072465/1‑151 (MQ=255)

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

                                            aCGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATg                                                                                           <  1:3047957/151‑1 (MQ=255)

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

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

                                                 tgaatgaaCCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTAGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATc                                                                                      <  1:4699213/151‑1 (MQ=255)

                                                  gaatgaaCCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCa                                                                                     <  1:1890399/151‑1 (MQ=255)

                                                   aatgaaCCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCAt                                                                                    >  1:4656671/1‑151 (MQ=255)

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

                                                         ccGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTgt                                                                              <  1:1676191/151‑1 (MQ=255)

                                                          cGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCCTGGTTgga                                                                             >  2:662042/1‑149 (MQ=255)

                                                           ggggCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTgtag                                                                            <  1:4156841/151‑1 (MQ=255)

                                                            gggCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTgtagt                                                                           >  1:5469225/1‑151 (MQ=255)

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

                                                             ggCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTACTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTgtagta                                                                          <  1:448781/151‑1 (MQ=255)

                                                                cgccTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTg                                                                       <  2:5983836/151‑1 (MQ=255)

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

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

                                                                    ttGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATg                                                                   <  2:486512/151‑1 (MQ=255)

                                                                     tGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGt                                                                  >  1:4348333/1‑151 (MQ=255)

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

                                                                       cGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTAc                                                                <  2:2302092/151‑1 (MQ=255)

                                                                        ggCCCGTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACt                                                               <  2:3181047/151‑1 (MQ=255)

                                                                          cccTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTCGTAGTAGTGAATGTACTCa                                                             >  2:4977424/1‑151 (MQ=255)

                                                                           ccTTCGCGGGCTCGCCATCGCTCCTACAATTTGGCATGCCACCCAAAAGTTGGACAGTTCTAAGCTACAGACCGAGCACTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCAt                                                            <  2:577659/151‑1 (MQ=17)

                                                                               cgcgGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAg                                                        <  2:4096761/151‑1 (MQ=255)

                                                                                gcgGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGa                                                       <  1:6287052/151‑1 (MQ=255)

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

                                                                                    gCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCg                                                   >  1:2219342/1‑151 (MQ=255)

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

                                                                                      tCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTCGTGAATTTACTCATCTCGTCCCGct                                                 >  1:4900857/1‑151 (MQ=14)

                                                                                        gCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGctct                                               <  1:205509/151‑1 (MQ=255)

                                                                                         cccTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCtctc                                              >  1:161464/1‑151 (MQ=255)

                                                                                          ccTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGGATGTACTCATCCTGACCAGTtctcc                                             >  1:4095766/1‑150 (MQ=14)

                                                                                          ccTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATCGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCa                                             <  1:4332424/151‑1 (MQ=255)

                                                                                            tCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGc                                           <  2:2822106/151‑1 (MQ=255)

                                                                                             cGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCt                                          <  1:564035/151‑1 (MQ=255)

                                                                                             cGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCAGCTCTCAGCt                                          >  2:5253267/1‑151 (MQ=255)

                                                                                                 cccACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCtcatc                                      >  1:746834/1‑151 (MQ=255)

                                                                                                   cacaATTTGGACTGCCCCCCCAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCaa                                    >  2:2698995/1‑151 (MQ=255)

                                                                                                      aaTTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGc                                 >  1:6194375/1‑151 (MQ=255)

                                                                                                       aTTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGct                                <  2:2902816/151‑1 (MQ=16)

                                                                                                       aTTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGct                                >  2:5914997/1‑151 (MQ=19)

                                                                                                       aTTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGct                                >  2:2507974/1‑151 (MQ=19)

                                                                                                          tGGACTGCGACCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTACTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCtctc                             <  2:4082874/151‑1 (MQ=255)

                                                                                                          tGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCtctc                             <  2:2190596/151‑1 (MQ=12)

                                                                                                           ggACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCg                            <  2:4421413/151‑1 (MQ=9)

                                                                                                               tGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAAc                        >  2:6510408/1‑151 (MQ=9)

                                                                                                               ggCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAAc                        <  2:5580209/150‑1 (MQ=9)

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

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

                                                                                                                   cccaAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACCGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCtt                    <  1:3253312/151‑1 (MQ=255)

                                                                                                                     ccAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCa                  >  2:5346463/1‑151 (MQ=12)

                                                                                                                      cAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAg                 >  1:6459782/1‑151 (MQ=12)

                                                                                                                        aaaGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGt               <  2:5897979/151‑1 (MQ=14)

                                                                                                                        aaaGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGt               <  2:2022967/151‑1 (MQ=17)

                                                                                                                          aGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGCCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAg             <  2:1889268/151‑1 (MQ=9)

                                                                                                                              ggACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAACCGCTTCACGTAGaaaa         >  1:5164470/1‑151 (MQ=9)

                                                                                                                               gACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGaaaaa        <  1:4550661/151‑1 (MQ=14)

                                                                                                                                    ttCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGAAAAACTCAg   >  1:4994059/1‑151 (MQ=9)

                                                                                                                                     tCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGAAAAACTCAGa  <  1:6043909/151‑1 (MQ=9)

TTCGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTGCAGCCTGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGAAAAACTCAGA  >  CP009974/5513725‑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: `‑`