BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,609,484	G→T	30.6%	R22L (CGC→CTC)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,609,484	0	G	T	30.6%	600.2 / 255.4	320	R22L (CGC→CTC)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand): ref base G (126/96); new base T (45/53); total (171/149)
Fisher's exact test for biased strand distribution p-value = 8.87e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.51e-01

Only 100 of 323 total aligned reads displayed.

ACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTC  >  CP009273/4609336‑4609632

 cacCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGc                                                                                                                                                     <  1:2697571/149‑1 (MQ=255)

    cAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACTCttt                                                                                                                                                  <  2:1840636/149‑1 (MQ=255)

    cAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACTCttt                                                                                                                                                  >  2:532690/1‑149 (MQ=255)

    cAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACTCttt                                                                                                                                                  >  1:1500032/1‑149 (MQ=255)

      aCTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTgg                                                                                                                                                >  1:4127656/1‑149 (MQ=255)

      aCTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTgg                                                                                                                                                >  1:1739030/1‑149 (MQ=255)

          tCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACTCTTTGGTGAc                                                                                                                                            <  1:1630571/149‑1 (MQ=255)

           cTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACTCTTTGGTGACg                                                                                                                                           >  1:695819/1‑149 (MQ=255)

            tATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACTCTTTGGTGACGt                                                                                                                                          <  1:4006874/149‑1 (MQ=255)

            tATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGt                                                                                                                                          <  2:1674069/149‑1 (MQ=255)

            tATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGt                                                                                                                                          <  2:4669785/149‑1 (MQ=255)

            tATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGt                                                                                                                                          <  1:1907276/149‑1 (MQ=255)

            tATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGt                                                                                                                                          <  2:4462557/149‑1 (MQ=255)

            tATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGt                                                                                                                                          <  2:710009/149‑1 (MQ=255)

             aTCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACTCTTTGGTGACGTc                                                                                                                                         >  2:3675726/1‑149 (MQ=255)

              tCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCg                                                                                                                                        >  1:993370/1‑149 (MQ=255)

               cGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACTCTTTGGTGACGTCgg                                                                                                                                       <  2:2627374/149‑1 (MQ=255)

               cGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACTCTTTGGTGACGTCgg                                                                                                                                       <  1:538449/149‑1 (MQ=255)

                gCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCggg                                                                                                                                      >  1:1186054/1‑149 (MQ=255)

                  cGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGc                                                                                                                                    >  1:847212/1‑149 (MQ=255)

                    tATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTg                                                                                                                                  >  1:384550/1‑149 (MQ=255)

                      tCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGac                                                                                                                                <  1:925728/149‑1 (MQ=255)

                           gAATAACGGTATACTGTTGTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACTCTTTGGTGACGTCGGGGCTGACACCCt                                                                                                                           <  2:2003553/149‑1 (MQ=255)

                           gAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACTCTTTGGTGACGTCGGGGCTGACACCCt                                                                                                                           <  2:2541429/149‑1 (MQ=255)

                           gAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCt                                                                                                                           >  1:116485/1‑149 (MQ=255)

                           gAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCt                                                                                                                           >  1:3227824/1‑149 (MQ=255)

                             aTAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTgg                                                                                                                         >  1:3417387/1‑149 (MQ=255)

                              tAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTggg                                                                                                                        >  1:3469486/1‑149 (MQ=255)

                               aaCGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGt                                                                                                                       >  2:2999328/1‑149 (MQ=255)

                                 ctgTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCa                                                                                                                     >  2:2574660/3‑149 (MQ=255)

                                  ggTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACTCTTTGGTGACGTCGGGGCTGACACCCTGGGTCat                                                                                                                    <  1:1953966/149‑1 (MQ=255)

                                     ataCTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACTCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATc                                                                                                                 <  1:3692357/149‑1 (MQ=255)

                                     ataCTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACTCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATc                                                                                                                 <  1:2097909/149‑1 (MQ=255)

                                     ataCTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCAGTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATc                                                                                                                 <  2:244198/149‑1 (MQ=255)

                                         tgatctgatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACTCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAg                                                                                                             >  1:1438463/1‑149 (MQ=255)

                                          gatctgatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGa                                                                                                            <  1:4730174/149‑1 (MQ=255)

                                              tgatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCt                                                                                                        >  1:4603847/1‑149 (MQ=255)

                                              tgatcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCt                                                                                                        <  1:1487334/149‑1 (MQ=255)

                                                atcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTtg                                                                                                      <  1:2395669/149‑1 (MQ=255)

                                                 tcATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACTCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTtgt                                                                                                     >  2:1541493/1‑149 (MQ=255)

                                                       aaaTTTGAAGCACTGAGTACGTAGAAAATATGAAACGTTCATTTATTATTGTGCTGTACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACTCTTTGGTGACGTCGGGGCTGACACACTGGGTCATATCGCAGAAGCTTGTGCCaaa                                                                                               <  2:1022939/149‑1 (MQ=255)

                                                       aaaTTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCaaa                                                                                               <  2:2983277/149‑1 (MQ=255)

                                                        aaTTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAg                                                                                              >  1:3520910/1‑149 (MQ=255)

                                                         aTTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCAGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAgg                                                                                             >  2:4574665/1‑149 (MQ=255)

                                                          tttGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGc                                                                                            >  1:1122521/1‑149 (MQ=255)

                                                          tttGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGc                                                                                            <  1:3163831/149‑1 (MQ=255)

                                                          tttGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCGTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGc                                                                                            >  1:780869/1‑149 (MQ=255)

                                                             gAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGaa                                                                                         >  2:1836410/1‑149 (MQ=255)

                                                               aGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACTCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGc                                                                                       <  2:3598610/149‑1 (MQ=255)

                                                               aGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGc                                                                                       >  2:1149611/1‑149 (MQ=255)

                                                               aGCACGGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACTCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGc                                                                                       <  2:2206818/149‑1 (MQ=255)

                                                                gCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCt                                                                                      >  1:722497/1‑149 (MQ=255)

                                                                gCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCt                                                                                      <  2:3192081/149‑1 (MQ=255)

                                                                 cACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTg                                                                                     >  2:804302/1‑149 (MQ=255)

                                                                 cACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTg                                                                                     >  1:3328578/1‑149 (MQ=255)

                                                                    tGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACTCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATa                                                                                  <  1:2540867/149‑1 (MQ=255)

                                                                    tGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATa                                                                                  <  2:1962160/149‑1 (MQ=255)

                                                                     gAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACTCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATaa                                                                                 <  1:2499880/149‑1 (MQ=255)

                                                                        tACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACgg                                                                              <  1:2850251/149‑1 (MQ=255)

                                                                          cGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTc                                                                            <  1:2178831/149‑1 (MQ=255)

                                                                           ggAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCg                                                                           >  2:3356515/1‑149 (MQ=255)

                                                                              gaACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTaa                                                                        >  1:1580392/1‑149 (MQ=255)

                                                                                aCATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAg                                                                      >  1:4387196/1‑149 (MQ=255)

                                                                                  atatGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACTCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGc                                                                    <  2:3490157/149‑1 (MQ=255)

                                                                                    atGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGccc                                                                  <  1:2273932/149‑1 (MQ=255)

                                                                                    atGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGccc                                                                  >  1:2154614/1‑149 (MQ=255)

                                                                                     tGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCg                                                                 >  1:1014877/1‑149 (MQ=255)

                                                                                     tGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCg                                                                 >  1:3800360/1‑149 (MQ=255)

                                                                                      gAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGc                                                                <  2:2347969/149‑1 (MQ=255)

                                                                                      gAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGc                                                                >  2:4366297/1‑149 (MQ=255)

                                                                                      gAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGc                                                                <  2:1781123/149‑1 (MQ=255)

                                                                                         aCGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCa                                                             >  1:450719/1‑149 (MQ=255)

                                                                                              cATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACTCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTg                                                        >  2:2155018/1‑149 (MQ=255)

                                                                                               aTTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGc                                                       >  2:38832/1‑149 (MQ=255)

                                                                                                    ttatGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAAt                                                  >  2:2000954/1‑149 (MQ=255)

                                                                                                        ggTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACTCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGa                                              >  2:3748316/1‑149 (MQ=255)

                                                                                                               gACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCt                                       >  1:4479988/1‑149 (MQ=255)

                                                                                                               gACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCt                                       >  2:731732/1‑149 (MQ=255)

                                                                                                               gACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCt                                       >  1:1690690/1‑149 (MQ=255)

                                                                                                                     ttCGGCATCGGCGCTACAGAAGATGCAGAACTCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCt                                 <  2:1499981/149‑1 (MQ=255)

                                                                                                                     ttCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCt                                 >  1:966804/1‑149 (MQ=255)

                                                                                                                     ttCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCt                                 >  1:636746/1‑149 (MQ=255)

                                                                                                                        ggCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGc                              <  2:1946450/149‑1 (MQ=255)

                                                                                                                        ggCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGAGAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGc                              <  2:1913884/149‑1 (MQ=255)

                                                                                                                         gCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCg                             >  2:2029095/1‑149 (MQ=255)

                                                                                                                           aTCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGaa                           >  1:221994/1‑149 (MQ=255)

                                                                                                                            tCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGaaa                          <  1:1262631/149‑1 (MQ=255)

                                                                                                                            tCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGaaa                          >  2:1457871/1‑149 (MQ=255)

                                                                                                                            tCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGaaa                          <  2:2737546/149‑1 (MQ=255)

                                                                                                                             cGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAg                         <  2:3844605/149‑1 (MQ=255)

                                                                                                                              ggCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGc                        >  2:829409/1‑149 (MQ=255)

                                                                                                                                cgcTACAGAAGATGCAGAACTCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGcac                      <  2:4050740/149‑1 (MQ=255)

                                                                                                                                 gcTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGcaca                     <  1:3500184/149‑1 (MQ=255)

                                                                                                                                  cTACAGAAGATGCAGAACTCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCacac                    <  1:861526/149‑1 (MQ=255)

                                                                                                                                      agaagaTGCAGAACTCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAg                >  2:2888630/1‑149 (MQ=255)

                                                                                                                                       gaagaTGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAgg               >  2:1234022/1‑149 (MQ=255)

                                                                                                                                        aagaTGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGt              <  1:4412492/149‑1 (MQ=255)

                                                                                                                                                gAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCgg      >  2:2813251/1‑149 (MQ=255)

                                                                                                                                                gAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCgg      >  1:4637214/1‑149 (MQ=255)

                                                                                                                                                  aCGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGtt    >  1:4352352/1‑149 (MQ=255)

ACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTC  >  CP009273/4609336‑4609632

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