BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP009974	1,796,148	C→T	Y243Y (TAC→TAT)	RPPX_08005 →	integrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009974	1,796,148	0	C	T	99.2%	639.9 / ‑4.4	267	Y243Y (TAC→TAT)	RPPX_08005	integrase
Reads supporting (aligned to +/- strand): ref base C (0/2); new base T (132/133); total (132/135)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 274 total aligned reads displayed.

TAAAGCAGAGCATGTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTACTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAGCGCGG  >  CP009974/1796000‑1796293

    gCAGAGCATGTCGCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGACACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAAc                                                                                                                                             <  2:6903037/151‑1 (MQ=12)

     cAGAGCATGTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGCAAGCTTCTTCGGCACGCCCCAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATGGATGAGCTGAGCGCGGGTCTGGATGAGTACATTCACTATTACAAcc                                                                                                                                            >  1:3249891/1‑151 (MQ=11)

     cAGAGCATGTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAAcc                                                                                                                                            <  1:3851147/151‑1 (MQ=17)

      agagCATGTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCa                                                                                                                                           >  1:4086481/1‑151 (MQ=17)

      agagCATGTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCa                                                                                                                                           >  1:4005801/1‑151 (MQ=17)

        agCATGTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATg                                                                                                                                         >  1:1336067/1‑151 (MQ=17)

            tGTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCg                                                                                                                                     >  2:3523590/1‑151 (MQ=17)

             gTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGc                                                                                                                                    <  2:3552628/151‑1 (MQ=17)

              tctcGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCa                                                                                                                                   >  1:2565984/1‑151 (MQ=14)

                tcGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATc                                                                                                                                 <  2:109614/151‑1 (MQ=17)

                  gAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTTAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTGTGGAGGAGGACAGTCACTATTACGAGCATGGTCTCaggaa                                                                                                                               >  2:328700/1‑147 (MQ=11)

                    aaaGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACATCCATGATCGCATCAAGc                                                                                                                             <  1:4309694/151‑1 (MQ=17)

                     aaGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCt                                                                                                                            >  1:4203285/1‑151 (MQ=17)

                      aGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTg                                                                                                                           <  1:5605356/151‑1 (MQ=17)

                       ggCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTga                                                                                                                          <  2:4956183/151‑1 (MQ=17)

                         cAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTgaga                                                                                                                        >  1:6583712/1‑151 (MQ=17)

                          aaTTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGAc                                                                                                                       >  2:5701135/1‑151 (MQ=17)

                          aaTTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGAc                                                                                                                       >  1:4448217/1‑151 (MQ=17)

                           aTTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACt                                                                                                                      <  1:475341/151‑1 (MQ=17)

                            ttGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTc                                                                                                                     >  2:937251/1‑151 (MQ=17)

                            ttGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTc                                                                                                                     >  1:1652465/1‑151 (MQ=17)

                              gCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCaa                                                                                                                   >  1:5016595/1‑151 (MQ=14)

                               ccTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAAt                                                                                                                  >  2:2737439/1‑151 (MQ=17)

                                  ggACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGc                                                                                                               <  1:5337237/151‑1 (MQ=17)

                                   gACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGcc                                                                                                              >  1:3526124/1‑151 (MQ=17)

                                    aCAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCt                                                                                                             <  2:732226/151‑1 (MQ=17)

                                    aCAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCt                                                                                                             >  1:1781711/1‑151 (MQ=17)

                                        tGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGc                                                                                                         >  1:4187923/1‑151 (MQ=17)

                                        tGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGc                                                                                                         <  2:2226591/151‑1 (MQ=17)

                                         gCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGcc                                                                                                        <  1:6570614/151‑1 (MQ=14)

                                            gCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTTCTTGTTTTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCtg                                                                                                     >  2:62028/1‑151 (MQ=11)

                                            gCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCtg                                                                                                     <  2:1931690/151‑1 (MQ=17)

                                               aTGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTgg                                                                                                  <  2:2994418/151‑1 (MQ=17)

                                               aTGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTgg                                                                                                  >  1:2260292/1‑151 (MQ=17)

                                                tGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGa                                                                                                 >  1:1902872/1‑151 (MQ=17)

                                                 ggAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGaa                                                                                                >  1:5545001/1‑151 (MQ=17)

                                                 ggAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGaa                                                                                                >  2:4422456/1‑151 (MQ=14)

                                                  gAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAAt                                                                                               <  1:6197023/151‑1 (MQ=17)

                                                   aaaGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATa                                                                                              >  2:4664593/1‑151 (MQ=17)

                                                     aGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCCAGCTGAGACTCAATTGCCTGAGCCCTGTGGAAGACa                                                                                            >  2:3650741/1‑151 (MQ=9)

                                                       cttcttCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAgg                                                                                          <  2:5378720/151‑1 (MQ=17)

                                                       cttcttCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAgg                                                                                          <  1:749448/151‑1 (MQ=14)

                                                          cttcGGCACGCTCAAATCTGAGTTTTTCTACCTTAAGTGTTTCGAGTGCATTGATGATCTGAGAGCGGTTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCt                                                                                       >  1:1030031/1‑151 (MQ=255)

                                                          cttcGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGa                                                                                                 <  2:2410917/141‑1 (MQ=17)

                                                          cttcGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCt                                                                                       >  1:3628572/1‑151 (MQ=17)

                                                           ttcGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTc                                                                                      <  1:2038665/151‑1 (MQ=17)

                                                              ggCCCGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGACTACAGGGCTGAgg                                                                                   >  2:3222134/1‑151 (MQ=255)

                                                                cACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCt                                                                                 <  1:4878332/151‑1 (MQ=17)

                                                                 aCGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTg                                                                                >  1:5395330/1‑151 (MQ=14)

                                                                   gCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCa                                                                              >  1:3780004/1‑151 (MQ=17)

                                                                     tCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGc                                                                            >  2:4734738/1‑151 (MQ=17)

                                                                      cAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCt                                                                           <  1:3706921/151‑1 (MQ=17)

                                                                      cAAATCTGAGTTTTTCTACCTGAAGTGCTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCACTGGCCTGAGCCCTGTCGAATACAGGGCTCCGGCTGCAGCt                                                                           >  2:4970123/1‑151 (MQ=11)

                                                                        aaTCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCACGCTGCCGCTTa                                                                         >  2:5234765/1‑151 (MQ=9)

                                                                         aTCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTaa                                                                        <  2:662942/151‑1 (MQ=14)

                                                                           cTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTaaaa                                                                      >  1:4552368/1‑151 (MQ=17)

                                                                            tGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAAc                                                                     <  2:5651861/151‑1 (MQ=14)

                                                                             gAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACt                                                                    <  2:1062976/151‑1 (MQ=17)

                                                                             aaGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACt                                                                    <  1:1612839/150‑1 (MQ=14)

                                                                              aGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTg                                                                   <  1:5215512/151‑1 (MQ=17)

                                                                                tttttCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTc                                                                 <  2:1807689/151‑1 (MQ=17)

                                                                                 ttttCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTcc                                                                <  2:4985780/151‑1 (MQ=17)

                                                                                 ttttCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTcc                                                                <  1:4964457/151‑1 (MQ=17)

                                                                                   ttCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCaa                                                              >  2:2120711/1‑151 (MQ=14)

                                                                                    tCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCCCAGGCTGCAGCTTAAAACTGTCCAAc                                                             >  2:810349/1‑151 (MQ=17)

                                                                                     cTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACt                                                            <  1:5740319/151‑1 (MQ=14)

                                                                                        ccTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACtttt                                                         >  1:3322128/1‑151 (MQ=17)

                                                                                        ccTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACtttt                                                         >  2:4254579/1‑151 (MQ=17)

                                                                                         cTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTg                                                        >  2:373941/1‑151 (MQ=17)

                                                                                               tgtTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGc                                                  <  2:1556594/151‑1 (MQ=11)

                                                                                               tgtTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGc                                                  <  1:988155/151‑1 (MQ=11)

                                                                                                       gagCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTTAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATa                                          <  2:5999591/151‑1 (MQ=11)

                                                                                                        agCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATAc                                         <  2:3412012/151‑1 (MQ=14)

                                                                                                          cATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCt                                       >  1:6743483/1‑151 (MQ=17)

                                                                                                           aTTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCtg                                      >  1:6425575/1‑151 (MQ=17)

                                                                                                           aTTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCtg                                      <  1:6297028/151‑1 (MQ=14)

                                                                                                             tgatgaGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCtgtg                                    >  2:5799414/1‑151 (MQ=17)

                                                                                                              gatgaGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGCGACTCAATGGCCTGCGCCCTGTGGAATACAGAGCTCAAGCCGCCGCTTCAAACTGCTCAACTTTTGGGGGGCAGTCCATAcccggtg                                   >  2:7126937/1‑146 (MQ=11)

                                                                                                               atgaGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTggg                                  >  2:2320901/1‑151 (MQ=14)

                                                                                                                   gCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCg                              <  1:309914/151‑1 (MQ=14)

                                                                                                                    cTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCgg                             >  1:5896207/1‑151 (MQ=17)

                                                                                                                     tGAGAGCGGGTCTGGATGAGAACATTCACTATTACAACCCTGATCGCATCAAGCTGAGACTCAATGGCCTGATCCCTGGGGAATACGGGGCGCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCggg                            <  2:6528361/151‑1 (MQ=11)

                                                                                                                      gagaGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGc                           <  2:3091816/151‑1 (MQ=14)

                                                                                                                        gagCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCaa                         >  1:3446198/1‑151 (MQ=18)

                                                                                                                         agCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAg                        >  1:7042978/1‑151 (MQ=21)

                                                                                                                         agCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAg                        <  1:527458/151‑1 (MQ=21)

                                                                                                                          gCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGc                       >  2:7008878/1‑151 (MQ=25)

                                                                                                                          gCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGc                       <  2:18427/151‑1 (MQ=25)

                                                                                                                            gggTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGccc                     <  1:1059208/151‑1 (MQ=37)

                                                                                                                            gggTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGccc                     <  1:2136553/151‑1 (MQ=37)

                                                                                                                             ggTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGAGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGTGCAAGTCCATACCCGTGGGAGCGGGCAAGCAcg                    >  2:4144385/1‑151 (MQ=11)

                                                                                                                              gTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCcgc                   >  2:5417691/1‑151 (MQ=25)

                                                                                                                               tCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCTGCTTTAAACTGTCCAACTTTTGGGGTGCAGTTCATACCCGGGGGAGCGGGCAAGCCcgcg                  >  1:2400737/1‑151 (MQ=11)

                                                                                                                               tCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCcgcg                  <  1:6022083/151‑1 (MQ=255)

                                                                                                                                      gAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGcc           <  1:2788570/151‑1 (MQ=38)

                                                                                                                                        gTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGc         <  2:1335140/151‑1 (MQ=255)

                                                                                                                                          aCATTCACTATTAGAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAGCTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAg       <  2:3500101/151‑1 (MQ=255)

                                                                                                                                           cATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGTGCCGCAgc      >  2:104994/1‑151 (MQ=255)

                                                                                                                                           cATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAgc      >  2:4010281/1‑151 (MQ=255)

                                                                                                                                               cACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAGCGCgg  >  1:2946243/1‑151 (MQ=255)

TAAAGCAGAGCATGTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTACTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAGCGCGG  >  CP009974/1796000‑1796293

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