BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP009974	5,597,121	C→T	T614T (ACC→ACT)	RPPX_24915 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009974	5,597,121	0	C	T	99.7%	1009.3 / NA	360	T614T (ACC→ACT)	RPPX_24915	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (160/199); minor base G (1/0); total (161/199)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 372 total aligned reads displayed.

CTTTAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACCACTGCGCGTGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTGATGCACCGGAATACACCCTG  >  CP009974/5597013‑5597268

cTTTAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTcc                                                                                                           <  2:3837670/151‑1 (MQ=2)

  ttAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCcgg                                                                                                         >  1:976669/1‑151 (MQ=9)

  ttAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCcgg                                                                                                         <  1:2006807/151‑1 (MQ=9)

    aGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCcggcg                                                                                                       >  2:3528790/1‑151 (MQ=9)

     gCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCcggcgg                                                                                                      <  1:1165313/151‑1 (MQ=12)

     gCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCcggcgg                                                                                                      <  1:1489350/151‑1 (MQ=12)

         ggaggTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCggtggg                                                                                                  >  1:3281892/1‑150 (MQ=9)

          gaggTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCggtggtg                                                                                                 <  2:553388/151‑1 (MQ=14)

           aggTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCggtggtgg                                                                                                <  2:82606/151‑1 (MQ=12)

           aggTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCggtggtgg                                                                                                >  1:3772782/1‑151 (MQ=12)

            ggTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCggtggtggt                                                                                               >  1:5019347/1‑151 (MQ=9)

             gTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCggtggtggtg                                                                                              >  1:3123340/1‑151 (MQ=12)

                cAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGtt                                                                                           >  1:717647/1‑151 (MQ=12)

                cAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGtt                                                                                           <  2:2962627/151‑1 (MQ=12)

                cAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGtt                                                                                           >  1:2962627/1‑151 (MQ=17)

                cAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGGGGGGGTGGtt                                                                                           >  2:1364108/1‑151 (MQ=255)

                 aGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGttt                                                                                          <  1:4622037/151‑1 (MQ=12)

                 aGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGttt                                                                                          >  1:2326234/1‑151 (MQ=17)

                  gCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTa                                                                                         <  1:1750644/151‑1 (MQ=12)

                  gCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTa                                                                                         >  2:3544373/1‑151 (MQ=12)

                   cGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTaa                                                                                        <  1:5071947/151‑1 (MQ=12)

                    ggTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCACGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTaaa                                                                                       >  1:2589601/1‑151 (MQ=9)

                     gTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAAc                                                                                      <  1:2513210/151‑1 (MQ=12)

                     gTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAAc                                                                                      >  1:3847999/1‑151 (MQ=9)

                      tATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACg                                                                                     <  2:1715581/151‑1 (MQ=12)

                       aTTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGc                                                                                    <  1:2732238/151‑1 (MQ=12)

                         ttgcggggCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTAGCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAg                                                                                  <  2:5220610/144‑1 (MQ=11)

                           gAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAggg                                                                                >  1:645541/1‑151 (MQ=255)

                           gAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAggg                                                                                <  2:2107215/151‑1 (MQ=19)

                            aaGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGa                                                                               <  2:3292827/151‑1 (MQ=19)

                            aaGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGa                                                                               <  1:4124291/151‑1 (MQ=19)

                            aaGAGCGGCAGGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGa                                                                               <  2:2674996/151‑1 (MQ=255)

                             agagCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGAt                                                                              <  2:87900/151‑1 (MQ=255)

                               agCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGTGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAGCCCAACTTCCTTGGGTGGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGAtgt                                                                            <  2:4033870/151‑1 (MQ=255)

                                  ggCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGAtgctgt                                                                         <  1:5171079/151‑1 (MQ=16)

                                  ggCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGAtgctgt                                                                         <  1:4857697/151‑1 (MQ=255)

                                   gCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTtgtg                                                                        >  1:4939500/1‑151 (MQ=255)

                                     aTGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGt                                                                      >  1:536545/1‑151 (MQ=255)

                                          acaGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAAc                                                                 >  1:1441182/1‑151 (MQ=255)

                                           caGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACa                                                                >  2:2322236/1‑151 (MQ=9)

                                           caGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGCGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACa                                                                <  2:3737922/151‑1 (MQ=255)

                                            aGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACAt                                                               >  2:1538631/1‑151 (MQ=9)

                                             ggTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATc                                                              >  1:3216646/1‑151 (MQ=255)

                                               ttGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCtt                                                            <  2:1233205/151‑1 (MQ=255)

                                                tGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTa                                                           >  2:3937794/1‑151 (MQ=255)

                                                tGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTa                                                           >  2:4798109/1‑151 (MQ=255)

                                                 gAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTAt                                                          >  2:4714138/1‑151 (MQ=255)

                                                  aaGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATc                                                         <  1:1267151/151‑1 (MQ=255)

                                                     agaTGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCt                                                      <  1:1680073/151‑1 (MQ=255)

                                                       aTGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCa                                                    <  2:3069608/151‑1 (MQ=255)

                                                        tGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCaa                                                   >  1:3050937/1‑151 (MQ=255)

                                                          aaaGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCaaca                                                 <  1:425642/151‑3 (MQ=255)

                                                          aaaGCGTTGCGCATCTCGGGTAGCCGCACGGCGTCTATCGAGCTGACCACTACTGCGCGCGTGAAAACCAACTTCCTTGGGTTGCTCACCTCCTGCGGGGGTGGTGGGTTAAACGCAGGGATGATGTGGTTAAACCTGTTATCGCTCAAGa                                                 >  2:4476117/1‑151 (MQ=11)

                                                           aaGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGAc                                                <  1:4012825/151‑1 (MQ=255)

                                                            aGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACg                                               <  1:1039589/151‑1 (MQ=255)

                                                             gCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACgg                                              <  1:1575380/151‑1 (MQ=255)

                                                             gCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACgg                                              <  1:4316694/151‑1 (MQ=255)

                                                               gTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGcc                                            >  1:4798607/1‑151 (MQ=11)

                                                               gTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGcc                                            >  1:4364268/1‑151 (MQ=11)

                                                                ttGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCt                                           >  1:4515/1‑151 (MQ=255)

                                                                ttGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGGTCAACATCTTATAGCTCAACACGGCCt                                           >  2:4589985/1‑151 (MQ=11)

                                                                 tGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTa                                          >  2:445453/1‑151 (MQ=11)

                                                                     cATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTacaac                                      <  2:607126/151‑1 (MQ=11)

                                                                        ctcGGGTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGc                                   >  1:3369216/1‑151 (MQ=255)

                                                                            ggTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGc                               <  2:159013/151‑1 (MQ=255)

                                                                             gTAGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCa                              <  1:831346/151‑1 (MQ=255)

                                                                               aGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCaga                            >  2:423475/1‑151 (MQ=11)

                                                                               aGCCGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGGGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCaga                            >  2:2477103/1‑151 (MQ=255)

                                                                                  cGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCagaaga                         <  2:1171626/151‑1 (MQ=255)

                                                                                  cGCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCagaaga                         >  1:2781437/1‑151 (MQ=14)

                                                                                   gCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAagag                        <  2:4324170/151‑1 (MQ=255)

                                                                                   gCACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAagag                        >  2:3190462/1‑151 (MQ=11)

                                                                                    cACGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGt                       >  2:2756394/1‑151 (MQ=255)

                                                                                     aCGGTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAagagcg                      <  2:4906369/151‑3 (MQ=11)

                                                                                       ggTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAagagcgac                    >  2:2159144/1‑147 (MQ=11)

                                                                                       ggTGTCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAagagcgac                    >  1:1770307/1‑147 (MQ=11)

                                                                                        gtgtCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAAGAGCGacg                   >  1:3028783/1‑149 (MQ=11)

                                                                                         tgtCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTGATGc                  <  2:4182914/151‑1 (MQ=255)

                                                                                         tgtCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAAGAGCGACGc                  >  2:4864932/1‑151 (MQ=11)

                                                                                         ggtCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAAGAGCGACGc                  <  2:4104649/150‑1 (MQ=11)

                                                                                           tCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTGATGCAc                >  1:4210738/1‑151 (MQ=255)

                                                                                           tCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAAGAGCGACGCAc                >  1:2907308/1‑151 (MQ=11)

                                                                                           tCTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAAGAGCGACGCAc                >  1:2863537/1‑151 (MQ=11)

                                                                                            cTATCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTGATGCAcc               <  2:1088268/151‑1 (MQ=255)

                                                                                              aTCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAAGAGCGACGCACCgg             >  2:3462451/1‑151 (MQ=11)

                                                                                              aTCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAAGAGCGACGCACCgg             >  1:2438831/1‑151 (MQ=11)

                                                                                              aTCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATACTGTGGTTCAACACGTTATCGCTCAACACGGCCTACAACAGCCTGCAGAAGAGCGACGCACCgg             >  2:5022819/1‑151 (MQ=11)

                                                                                               tCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAAGAGCGACGCACCGGa            <  2:2781437/151‑1 (MQ=11)

                                                                                               tCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAAGAGCGACGCACCGGa            <  2:4770708/151‑1 (MQ=11)

                                                                                               tCGAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAAGAGCGACGCACCGGa            >  1:3706597/1‑151 (MQ=11)

                                                                                                 gAGCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAAGAGCGACGCACCGGAAt          <  1:4703554/151‑1 (MQ=11)

                                                                                                   gCTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTGATGCACCGGAATac        <  1:498223/151‑1 (MQ=255)

                                                                                                    cTGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTGATGCACCGGAATaca       <  1:3631520/151‑1 (MQ=255)

                                                                                                     tGACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAAGAGCGACGCACCGGAATacac      >  2:1812322/1‑151 (MQ=11)

                                                                                                      gACCACTACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAAGAGCGACGCACCGGACTACAcc     >  2:5147696/1‑151 (MQ=11)

                                                                                                       accactACTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTGATGCACCGGAATACAccc    >  1:1352848/1‑151 (MQ=255)

                                                                                                         cactacTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCATGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACATCCTGCAGAAGATCGACGCACCGGAATACACCCTg  >  1:889002/1‑151 (MQ=11)

                                                                                                         cactacTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTGATGCACCGGAATACACCCTg  <  1:1890326/151‑1 (MQ=255)

                                                                                                         cactacTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAAGAGCGACGCACCGGAATACACCCTg  <  2:2382183/151‑1 (MQ=11)

                                                                                                         cactacTGCGCGCGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGCTGTGGTTCAACATCTTATCGCTCAACACGGCCTACAACAGCCTGCAGAAGAGCGACGCACCGGAATACACCCTg  <  2:2829550/151‑1 (MQ=11)

CTTTAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCACAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACGGTGTCTATCGAGCTGACCACCACTGCGCGTGTGAAACCCAACTTCCTTGGGTTGCTCACCTCCGGCGGTGGTGGTGGTTTAAACGCAGGGATGTTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGTGATGCACCGGAATACACCCTG  >  CP009974/5597013‑5597268

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