BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,241,983	C→A	5.2%	D284Y (GAT→TAT)	treA ←	periplasmic trehalase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,241,983	0	C	A	5.2%	942.9 / 16.5	248	D284Y (GAT→TAT)	treA	periplasmic trehalase
Reads supporting (aligned to +/- strand): ref base C (110/125); new base A (7/6); total (117/131)
Fisher's exact test for biased strand distribution p-value = 7.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 253 total aligned reads displayed.

ACGATGCTGGTGGTGCGTAAGGTATTTAACTGCTGCGGGTTGTCCATCCAGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGCAGGTTTTCAACACCGTCCATCCAGTAAGCATATTCTTTTTGC  >  CP009273/1241835‑1242131

aCGATGCTGGTGGTGCGTAAGGTATTTAACTGCTGCGGGTTGTCCATCCAGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATc                                                                                                                                                      <  2:499512/149‑1 (MQ=255)

    tGCTGGTGGTGCGTAAGGTATTTAACTGCTGCGGGTTGTCCATCCAGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTcc                                                                                                                                                  >  1:679176/1‑149 (MQ=255)

     gCTGGTGGTGCGTAAGGTATTTAACTGCTGCGGGTTGTCCATCCAGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCa                                                                                                                                                 <  2:3794385/149‑1 (MQ=255)

      cTGGTGGTGCGTAAGGTATTTAACTGCTGCGGGTTGTCCATCCAGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCAc                                                                                                                                                <  1:3651192/149‑1 (MQ=255)

       tggtggTGCGTAAGGTATTTAACTGCTGCGGGTTGTCCATCCAGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCAcc                                                                                                                                               >  1:1484219/1‑149 (MQ=255)

              gCGTAAGGTATTTAACTGCTGCGGGTTGTCCATCCAGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGAct                                                                                                                                        >  2:1762066/1‑149 (MQ=255)

              gCGTAAGGTATTTAACTGCTGCGGGTTGTCCATCCAGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGAct                                                                                                                                        >  1:2475109/1‑149 (MQ=255)

                 tAAGGTATTTAACTGCTGCGGGTTGTCCATCCAGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTg                                                                                                                                     >  1:3497733/1‑149 (MQ=255)

                   aGGTATTTAACTGCTGCGGGTTGTCCATCCAGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGt                                                                                                                                   >  2:3238625/1‑149 (MQ=255)

                     gTATTTAACTGCTGCGGGTTGACCATCCAGCGCGAGCTGAACTCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTACATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCg                                                                                                                                 <  2:3585615/149‑1 (MQ=255)

                      tATTTAACTGCTGCGGGTTGTCCATCCAGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGt                                                                                                                                >  1:446995/1‑149 (MQ=255)

                       aTTTAACTGCTGCGGGTTGTCCATCCAGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTg                                                                                                                               >  1:106870/1‑149 (MQ=255)

                       aTTTAACTGCTGCGGGTTGTCCATCCAGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTg                                                                                                                               >  1:2702264/1‑149 (MQ=255)

                        tttAACTGCTGCGGGTTGTCCATCCAGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTgg                                                                                                                              <  1:3227432/149‑1 (MQ=255)

                           aaCTGCTGCGGGTTGTCCATCCAGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGt                                                                                                                           <  1:2644076/149‑1 (MQ=255)

                            aCTGCTGCGGGTTGTCCATCCAGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTa                                                                                                                          <  1:2606177/149‑1 (MQ=255)

                             ctgctgCGGGTTGTCCATCCAGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTAt                                                                                                                         <  2:2705872/149‑1 (MQ=255)

                             ctgctgCGGGTTGTCCATCCAGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTAt                                                                                                                         >  2:3260388/1‑149 (MQ=255)

                             ctgctgCGGGTTGTCCATCCAGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTAt                                                                                                                         >  2:2988712/1‑149 (MQ=255)

                             ctgctgCGGGTTGTCCATCCAGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGGGGCAGGTCGAGTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTAt                                                                                                                         >  1:572672/1‑149 (MQ=255)

                                 tgcGGGTTGTCCATCCAGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATcgcg                                                                                                                     >  2:2563240/1‑149 (MQ=255)

                                 tgcGGGTTGTCCATCCAGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATcgcg                                                                                                                     >  2:1890820/1‑149 (MQ=255)

                                     ggTTGTCCATCCAGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCg                                                                                                                 >  1:64513/1‑149 (MQ=255)

                                       ttGTCCATCCAGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTc                                                                                                               <  1:3384122/149‑1 (MQ=255)

                                        tGTCCATCCAGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTcc                                                                                                              <  1:3403487/149‑1 (MQ=255)

                                          tccatccaGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTACATTTCAGTGGCAGGTCGCGTCGGAGTGCTGTTGGCGGTGGCAAGATCTTCCACCCATGACTCTGGTCGTGGGGTAGCGCGATCGTCCCa                                                                                                            >  2:1223485/1‑149 (MQ=255)

                                          tccatccaGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCa                                                                                                            <  1:2923388/149‑1 (MQ=255)

                                            catccaGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGt                                                                                                          <  1:2182566/149‑1 (MQ=255)

                                              tccaGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAg                                                                                                        >  1:2922726/1‑149 (MQ=255)

                                               ccaGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATATTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGc                                                                                                       <  1:1332276/149‑1 (MQ=255)

                                                caGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTAGCGTATCGCGATCGTCCCAGTAGCg                                                                                                      <  1:1645023/149‑1 (MQ=255)

                                                 aGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCgg                                                                                                     >  2:2387213/1‑149 (MQ=255)

                                                  gcgcGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGt                                                                                                    >  2:1263703/1‑149 (MQ=255)

                                                  gcgcGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGt                                                                                                    <  1:3823861/149‑1 (MQ=255)

                                                  gcgcGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGt                                                                                                    >  1:2197740/1‑149 (MQ=255)

                                                    gcgAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTg                                                                                                  <  1:1826311/149‑1 (MQ=255)

                                                    gcgAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTg                                                                                                  <  2:3674055/149‑1 (MQ=255)

                                                      gAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTgag                                                                                                >  1:2269243/1‑149 (MQ=255)

                                                        gCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGaa                                                                                              <  1:3323812/149‑1 (MQ=255)

                                                          tGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAgg                                                                                            >  1:315881/1‑149 (MQ=255)

                                                                cccAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCa                                                                                      <  2:1684055/149‑1 (MQ=255)

                                                                 ccagccagACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCAt                                                                                     >  1:1642036/1‑149 (MQ=255)

                                                                 ccagccagACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCAt                                                                                     <  2:3743317/149‑1 (MQ=255)

                                                                  cagccagACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATc                                                                                    <  2:930341/149‑1 (MQ=255)

                                                                  cagccagACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATc                                                                                    <  1:111274/149‑1 (MQ=255)

                                                                    gccagACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCt                                                                                  >  2:3432383/1‑149 (MQ=255)

                                                                        gACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAg                                                                              >  1:266231/1‑149 (MQ=255)

                                                                         aCGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGt                                                                             <  1:1116407/149‑1 (MQ=255)

                                                                          cGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGtt                                                                            <  1:2046415/149‑1 (MQ=255)

                                                                           gCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAATGGTACCATCCTGAAGttt                                                                           >  2:2157114/1‑149 (MQ=255)

                                                                           gCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCACAGTAGCGGTTGAGAAGGGTACCATCCTGAAGttt                                                                           >  2:76556/1‑149 (MQ=255)

                                                                            cAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTg                                                                          >  2:758124/1‑149 (MQ=255)

                                                                            cAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTg                                                                          >  2:1382700/1‑149 (MQ=255)

                                                                             aGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGa                                                                         >  1:1313540/1‑149 (MQ=255)

                                                                              gCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGCTCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTGTGAc                                                                        >  2:3056022/1‑149 (MQ=255)

                                                                              gCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGAc                                                                        <  2:110602/149‑1 (MQ=255)

                                                                                     agCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGt                                                                 <  2:2868803/149‑1 (MQ=255)

                                                                                     agCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGt                                                                 <  1:362533/149‑1 (MQ=255)

                                                                                      gcgcAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAAcgcggt                                                                >  1:3089336/1‑147 (MQ=255)

                                                                                        gcAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGtttt                                                              >  2:1484764/1‑149 (MQ=255)

                                                                                         cAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGttttt                                                             <  2:1829552/149‑1 (MQ=255)

                                                                                          aGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTc                                                            >  2:899979/1‑149 (MQ=255)

                                                                                            gTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCt                                                          <  2:3145817/149‑1 (MQ=255)

                                                                                             tCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCtg                                                         <  1:1180306/149‑1 (MQ=255)

                                                                                              cgcgGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATAGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCtgt                                                        <  1:399745/149‑1 (MQ=255)

                                                                                               gcgGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTCTCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCtgtt                                                       >  2:2481191/1‑149 (MQ=255)

                                                                                                  gTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTc                                                    >  1:3959304/1‑149 (MQ=255)

                                                                                                  gTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTc                                                    >  1:3499480/1‑149 (MQ=255)

                                                                                                  gTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTc                                                    >  2:2010065/1‑149 (MQ=255)

                                                                                                      aTTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGGCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGc                                                <  2:628947/149‑1 (MQ=255)

                                                                                                          cAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGc                                            <  1:2815329/149‑1 (MQ=255)

                                                                                                          cAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGc                                            <  1:3773877/149‑1 (MQ=255)

                                                                                                            gTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGCAg                                          >  2:836737/1‑149 (MQ=255)

                                                                                                             tGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGCAgg                                         <  2:3456897/149‑1 (MQ=255)

                                                                                                               gCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATATTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGCAGGtt                                       <  2:244341/149‑1 (MQ=255)

                                                                                                                  ggTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGCAGGTTTTc                                    <  2:480117/149‑1 (MQ=255)

                                                                                                                  ggTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGCAGGTTTTc                                    >  2:2680603/1‑149 (MQ=255)

                                                                                                                      gATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGCAGGTTTTCAaca                                <  1:3144641/149‑1 (MQ=255)

                                                                                                                      gATTCGGATTGCTTTTGGCGGTGGCAATATATTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGCAGGTTTTCAaca                                >  2:1865059/1‑149 (MQ=255)

                                                                                                                       aTTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGCAGGTTTTCAacac                               >  2:2538797/1‑149 (MQ=255)

                                                                                                                       aTTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGCAGGTTTTCAacac                               >  1:2117080/1‑149 (MQ=255)

                                                                                                                       aTTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGCAGGTTTTCAacac                               >  2:2511118/1‑149 (MQ=255)

                                                                                                                          cGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGCAGGTTTTCAACACCGt                            <  1:2989047/149‑1 (MQ=255)

                                                                                                                            gATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGCAGGTTTTCAACACCGtcc                          <  1:1279866/149‑1 (MQ=255)

                                                                                                                             aTTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGCAGGTTTTCAACACCGtcca                         <  1:937055/149‑1 (MQ=255)

                                                                                                                               tGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGCAGGTTTTCAACACCGtccatc                       <  2:1253392/149‑1 (MQ=255)

                                                                                                                                        cGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGCAGGTTTTCAACACCGTCCATCCAGTAAGCa              <  2:1042855/149‑1 (MQ=255)

                                                                                                                                        cGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGCAGGTTTTCAACACCGTCCATCCAGTAAGCa              <  1:1141618/149‑1 (MQ=255)

                                                                                                                                          gTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGCAGGTTTTCAACACCGTCCATCCAGTCAGCata            >  1:843299/1‑149 (MQ=255)

                                                                                                                                          gTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGCAGGTTTTCAACACCGTCCATCCAGTAAGCata            >  2:1619787/1‑149 (MQ=255)

                                                                                                                                           tGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGCAGGTTTTCAACACCGTCCATCCAGTAAGCatat           >  2:1948680/1‑149 (MQ=255)

                                                                                                                                            ggCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGCAGGTTTTCAACACCGTCCATCCAGTAAGCATAtt          <  1:1032393/149‑1 (MQ=255)

                                                                                                                                              cAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGCAGGTTTTCAACACCGTCCATCCAGTAAGCATATTCt        <  1:3178347/149‑1 (MQ=255)

                                                                                                                                              cAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGCAGGTTTTCAACACCGTCCATCCAGTAAGCATATTCt        <  2:3497733/149‑1 (MQ=255)

                                                                                                                                              cAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGCAGGTTTTCAACACCGTCCATCCAGTAAGCATATTCt        <  2:315881/149‑1 (MQ=255)

                                                                                                                                               aaTATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGCAGGTTTTCAACACCGTCCATCCAGTAAGCATATTCtt       >  2:3327430/1‑149 (MQ=255)

                                                                                                                                               aaTATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGCAGGTTTTCAACACCGTCCATCCAGTAAGCATATTCtt       <  1:773408/149‑1 (MQ=255)

                                                                                                                                                atatCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGCAGGTTTTCAACACCGTCCATCCAGTAAGCATATTCttt      >  2:1497469/1‑149 (MQ=255)

                                                                                                                                                atatCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGCAGGTTTTCAACACCGTCCATCCAGTAAGCATATTCttt      <  1:1762066/149‑1 (MQ=255)

                                                                                                                                                 tatCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGCAGGTTTTCAACACCGTCCATCCAGTAAGCATATTCtttt     <  2:2856610/149‑1 (MQ=255)

ACGATGCTGGTGGTGCGTAAGGTATTTAACTGCTGCGGGTTGTCCATCCAGCGCGAGCTGAAATCCCAGCCAGACGCAGCGGCAGAGCGCAGGTCGCGGTAAATTTCAGTGGCAGGTCGATTCGGATTGCTTTTGGCGGTGGCAATATCTTCCACCCATGACTCTGGTCGTGGCGTATCGCGATCGTCCCAGTAGCGGTTGAGAAGGGTACCATCCTGAAGTTTGACAACGCGTTTTTCCTGTTGTCCGGCTTGCAGGTTTTCAACACCGTCCATCCAGTAAGCATATTCTTTTTGC  >  CP009273/1241835‑1242131

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