BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,097,262	C→T	44.6%	F297F (TTC→TTT)	mutY →	adenine DNA glycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,097,262	0	C	T	44.6%	170.7 / inf	332	F297F (TTC→TTT)	mutY	adenine DNA glycosylase
Reads supporting (aligned to +/- strand): ref base C (77/107); new base T (73/75); total (150/182)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 340 total aligned reads displayed.

AGCGTCCGCCGAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACT  >  CP009273/3097114‑3097409

aGCGTCCGCCGAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTcc                                                                                                                                                    >  1:481132/1‑150 (MQ=255)

aGCGTCCGCCGAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTcc                                                                                                                                                    <  1:1558734/150‑1 (MQ=255)

    tCCGCCGAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACtt                                                                                                                                                <  1:739558/150‑1 (MQ=255)

     ccgccgAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTa                                                                                                                                               <  2:4753547/150‑1 (MQ=255)

      cgccgAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAg                                                                                                                                              >  2:2336696/1‑150 (MQ=255)

      cgccgAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAg                                                                                                                                              <  2:45421/150‑1 (MQ=255)

        ccgAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGat                                                                                                                                            <  2:3777855/150‑1 (MQ=255)

          gAGCGGATTGTGGGGCGGTTTGTACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGatat                                                                                                                                          <  2:144427/150‑1 (MQ=255)

           aGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATAtt                                                                                                                                         >  1:2836856/1‑150 (MQ=255)

            gCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATtg                                                                                                                                        <  2:561906/150‑1 (MQ=255)

              ggATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATtgtg                                                                                                                                      >  2:1568150/1‑150 (MQ=255)

                  tgtgGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTa                                                                                                                                  >  2:2179050/1‑150 (MQ=255)

                   gtgGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTAt                                                                                                                                 <  2:857391/150‑1 (MQ=255)

                    tgGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTAtg                                                                                                                                >  2:3651700/1‑150 (MQ=255)

                    tgGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTAtg                                                                                                                                >  1:3155515/1‑150 (MQ=255)

                    tgGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTAtg                                                                                                                                >  1:1554060/1‑150 (MQ=255)

                     ggggCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTAtgt                                                                                                                               >  2:630896/1‑150 (MQ=255)

                      gggCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTAtgtg                                                                                                                              >  1:2628257/1‑150 (MQ=255)

                        gCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGc                                                                                                                            <  1:4659975/150‑1 (MQ=255)

                          ggTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCtt                                                                                                                          >  2:4704401/1‑150 (MQ=255)

                            tttATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTcc                                                                                                                        <  2:2828317/150‑1 (MQ=255)

                            tttATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTcc                                                                                                                        >  1:2318010/1‑150 (MQ=255)

                              tataCTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCg                                                                                                                      <  2:1801304/150‑1 (MQ=255)

                               ataCTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCgt                                                                                                                     <  2:3433232/150‑1 (MQ=255)

                               ataCTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCgt                                                                                                                     <  2:4257474/150‑1 (MQ=255)

                                 aCTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCgtgt                                                                                                                   <  1:3980374/150‑1 (MQ=255)

                                 aCTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCgtgt                                                                                                                   <  1:4030358/150‑1 (MQ=255)

                                   tGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTgtcg                                                                                                                 <  1:2368437/150‑1 (MQ=255)

                                     tttCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTgtcgtc                                                                                                               >  2:3136459/1‑150 (MQ=255)

                                            cAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCAcc                                                                                                        <  2:1795542/150‑1 (MQ=255)

                                             aGTTTGCCGACGAAGAAAGTTTGCGGCAGGGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCg                                                                                                       <  2:394573/150‑1 (MQ=255)

                                                ttGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCt                                                                                                    <  1:3521761/150‑1 (MQ=255)

                                                 tGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTg                                                                                                   >  2:3193638/1‑150 (MQ=255)

                                                  gCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGc                                                                                                  <  2:2618372/150‑1 (MQ=255)

                                                   ccGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCa                                                                                                 >  1:3602015/1‑150 (MQ=255)

                                                      acgaAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATgg                                                                                              <  1:1641713/150‑1 (MQ=255)

                                                       cgaAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGa                                                                                             <  1:2454287/150‑1 (MQ=255)

                                                        gaagaaAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGAt                                                                                            >  1:2982397/1‑150 (MQ=255)

                                                        gaagaaAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGAt                                                                                            <  1:962154/150‑1 (MQ=255)

                                                        gaagaaAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGAt                                                                                            >  1:4340484/1‑150 (MQ=255)

                                                          agaaAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGa                                                                                          >  1:2221643/1‑150 (MQ=255)

                                                             aaGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAgg                                                                                       >  1:1383490/1‑150 (MQ=255)

                                                               gTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCa                                                                                     >  1:2169774/1‑150 (MQ=255)

                                                                tttGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCaa                                                                                    <  1:1596381/150‑1 (MQ=255)

                                                                 ttGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAAt                                                                                   <  2:2513252/150‑1 (MQ=255)

                                                                 ttGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAAt                                                                                   <  2:4207554/150‑1 (MQ=255)

                                                                 ttGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAAt                                                                                   <  2:3454663/150‑1 (MQ=255)

                                                                 ttGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAAt                                                                                   <  2:2734770/150‑1 (MQ=255)

                                                                       cAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGctc                                                                             >  2:2221826/1‑150 (MQ=255)

                                                                       cAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGctc                                                                             <  1:4524971/150‑1 (MQ=255)

                                                                          tggctggcGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTgg                                                                          >  1:2044353/1‑150 (MQ=255)

                                                                             ctggcGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGtat                                                                       >  2:769959/1‑150 (MQ=255)

                                                                              tggcGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGtata                                                                      >  1:1083888/1‑150 (MQ=255)

                                                                              tggcGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGtata                                                                      >  1:3490326/1‑150 (MQ=255)

                                                                              tggcGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGtata                                                                      <  2:846924/150‑1 (MQ=255)

                                                                                gcgcAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAAc                                                                    <  2:4590726/150‑1 (MQ=255)

                                                                                  gcAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTTAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACtt                                                                  >  1:3332102/1‑150 (MQ=255)

                                                                                  gcAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACtt                                                                  >  2:30181/1‑150 (MQ=255)

                                                                                    aaCGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAg                                                                >  1:3873244/1‑150 (MQ=255)

                                                                                       ggCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAgcgc                                                             >  2:1405960/1‑150 (MQ=255)

                                                                                       ggCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAgcgc                                                             >  2:648331/1‑150 (MQ=255)

                                                                                        gCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCa                                                            <  1:899627/150‑1 (MQ=255)

                                                                                         cAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCaa                                                           >  2:1519170/1‑150 (MQ=255)

                                                                                           gATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAAcc                                                         >  1:81421/1‑150 (MQ=255)

                                                                                            aTTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAAccg                                                        >  2:2017425/1‑150 (MQ=255)

                                                                                             ttGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAAccgc                                                       >  1:3704122/1‑150 (MQ=255)

                                                                                               gctgcCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAAccgccg                                                     <  2:3832619/150‑1 (MQ=255)

                                                                                                  gcCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCa                                                  <  1:4118353/150‑1 (MQ=255)

                                                                                                       tAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTgg                                             >  1:1226943/1‑150 (MQ=255)

                                                                                                          ccTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCt                                          >  2:3072064/1‑150 (MQ=255)

                                                                                                           cTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTa                                         <  1:4566160/150‑1 (MQ=255)

                                                                                                            tGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAg                                        <  1:2336695/150‑1 (MQ=255)

                                                                                                            tGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAg                                        >  2:4655131/1‑150 (MQ=255)

                                                                                                              aCGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCg                                      <  1:4459883/150‑1 (MQ=255)

                                                                                                              aCGCAACTGACCGCGTTTCGGCATACCTTCAGCCAATTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACAGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCg                                      >  2:4067138/1‑150 (MQ=255)

                                                                                                               cGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCgg                                     <  1:630896/150‑1 (MQ=255)

                                                                                                                  aaCTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTc                                  <  1:3828706/150‑1 (MQ=255)

                                                                                                                  aaCTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTc                                  <  1:2664818/150‑1 (MQ=255)

                                                                                                                  aaCTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTc                                  <  1:1869662/150‑1 (MQ=255)

                                                                                                                   aCTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTcc                                 >  1:1014154/1‑150 (MQ=255)

                                                                                                                   aCTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTcc                                 <  2:428526/150‑1 (MQ=255)

                                                                                                                     tGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCg                               >  2:1041469/1‑150 (MQ=255)

                                                                                                                     tGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCg                               >  2:1722351/1‑150 (MQ=255)

                                                                                                                     tGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCg                               <  1:3240978/150‑1 (MQ=255)

                                                                                                                      gACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGt                              >  1:1135572/1‑150 (MQ=255)

                                                                                                                        ccGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTgg                            >  2:1947016/1‑150 (MQ=255)

                                                                                                                           cgTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGc                         <  1:586603/150‑1 (MQ=255)

                                                                                                                            gTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCg                        >  2:1249891/1‑150 (MQ=255)

                                                                                                                              ttCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGtt                      <  2:3873244/150‑1 (MQ=255)

                                                                                                                                 ggCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGt                   <  1:2372825/150‑1 (MQ=255)

                                                                                                                                  gCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGtt                  <  2:4670577/150‑1 (MQ=255)

                                                                                                                                   cATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTa                 <  1:1186548/150‑1 (MQ=255)

                                                                                                                                   cATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTa                 >  2:2464210/1‑150 (MQ=255)

                                                                                                                                     tACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTAca               >  2:3349456/1‑150 (MQ=255)

                                                                                                                                       ccTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTAcagc             >  1:2100729/1‑150 (MQ=255)

                                                                                                                                          tCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGt          <  2:2939419/150‑1 (MQ=255)

                                                                                                                                           cAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGtt         >  1:440566/1‑150 (MQ=255)

                                                                                                                                           cAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGtt         <  1:2152870/150‑1 (MQ=255)

                                                                                                                                                 tttCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCAc   >  1:3140402/1‑150 (MQ=255)

                                                                                                                                                  ttCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACt  <  2:3571156/150‑1 (MQ=255)

AGCGTCCGCCGAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACT  >  CP009273/3097114‑3097409

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