BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,805,420	C→A	17.1%	A212E (GCG→GAG)	emrA →	multidrug efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,805,420	0	C	A	17.1%	874.2 / 114.9	280	A212E (GCG→GAG)	emrA	multidrug efflux system
Reads supporting (aligned to +/- strand): ref base C (128/104); new base A (24/24); total (152/128)
Fisher's exact test for biased strand distribution p-value = 5.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 310 total aligned reads displayed.

CGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACG  >  CP009273/2805271‑2805568

 gCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGgcg                                                                                                                                                     >  1:2110645/1‑150 (MQ=255)

  cccGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGgcgc                                                                                                                                                    >  1:1146847/1‑150 (MQ=255)

   ccGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCt                                                                                                                                                   >  1:454344/1‑150 (MQ=255)

   ccGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCt                                                                                                                                                   >  1:3016024/1‑150 (MQ=255)

    cgcgACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGAGTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTg                                                                                                                                                  >  2:3729489/1‑150 (MQ=255)

        aCGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGc                                                                                                                                              >  1:2238197/1‑150 (MQ=255)

             gTCACCAGTGCCCAGGCGCCACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACt                                                                                                                                         <  1:4311525/150‑1 (MQ=255)

             gTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGAGCTGGAGCGTACt                                                                                                                                         >  1:3847433/1‑150 (MQ=255)

              tCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTc                                                                                                                                        <  2:3448848/150‑1 (MQ=255)

               cACCAGTGCCCATGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCg                                                                                                                                       <  2:1018882/150‑1 (MQ=255)

                aCCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGt                                                                                                                                      <  1:1596549/150‑1 (MQ=255)

                aCCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGt                                                                                                                                      >  2:2482598/1‑150 (MQ=255)

                 ccAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGta                                                                                                                                     <  1:474742/150‑1 (MQ=255)

                 ccAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGta                                                                                                                                     <  1:761030/150‑1 (MQ=255)

                   aGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGtatt                                                                                                                                   >  2:1880282/1‑150 (MQ=255)

                       cccAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCa                                                                                                                               <  2:585134/150‑1 (MQ=255)

                        ccAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAg                                                                                                                              <  2:2583746/150‑1 (MQ=255)

                        ccAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAg                                                                                                                              >  2:2710264/1‑150 (MQ=255)

                         cAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGAGCTGGAGCGTACTCGTATTATCAGt                                                                                                                             >  1:1894698/1‑150 (MQ=255)

                               cAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATg                                                                                                                       <  2:1421225/150‑1 (MQ=255)

                                aaCTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGa                                                                                                                      >  2:3362374/1‑150 (MQ=255)

                                  cTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGAcc                                                                                                                    >  1:1606639/1‑150 (MQ=255)

                                    ggACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTCTTATCAGTCCGATGACCgg                                                                                                                  >  2:4431131/1‑150 (MQ=255)

                                    ggACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCgg                                                                                                                  >  2:734358/1‑150 (MQ=255)

                                     gACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGt                                                                                                                 >  1:534930/1‑150 (MQ=255)

                                         tCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATg                                                                                                             >  2:4544027/1‑150 (MQ=255)

                                          cgcgATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGt                                                                                                            >  1:1948460/1‑150 (MQ=255)

                                          cgcgATTCAACAATAAAATGCCAATCAGGCGATGATTCTGTGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTATCGCATGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGTCCGGTTATGt                                                                                                            <  2:41777/150‑1 (MQ=255)

                                             gATCCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGtctc                                                                                                         <  1:3097000/150‑1 (MQ=255)

                                              aTTCAACAATACACTGCCTATCAGGCGATGATGCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTATCGCCTGGCTGGCGCTGGATCGTACTAGTATTATGAGTCCGATGACCGGTTATGTCTcc                                                                                                        <  2:1377976/150‑1 (MQ=255)

                                                 caacaaTACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCccgc                                                                                                     >  1:3783584/1‑150 (MQ=255)

                                                  aacaaTACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCccgcc                                                                                                    >  1:2929146/1‑150 (MQ=255)

                                                         caatGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGt                                                                                             >  2:2873802/1‑150 (MQ=255)

                                                          aatGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTa                                                                                            >  1:3267485/1‑150 (MQ=255)

                                                            tGCCAATCAGGCGATGATTCTGGGGACTACACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGAGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACa                                                                                          <  1:246062/150‑1 (MQ=255)

                                                              ccAATCAGGCGATGATTCTGTGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGc                                                                                        >  1:2888249/1‑150 (MQ=255)

                                                                aaTCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCt                                                                                      <  1:2955719/150‑1 (MQ=255)

                                                                aaTCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCt                                                                                      <  1:2101041/150‑1 (MQ=255)

                                                                   cAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTggg                                                                                   >  2:4115148/1‑150 (MQ=255)

                                                                       cGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGgcgc                                                                               >  2:921518/1‑150 (MQ=255)

                                                                        gatgatTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCa                                                                              >  1:4195159/1‑150 (MQ=255)

                                                                         atgatTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGAGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCaa                                                                             >  1:3288018/1‑150 (MQ=255)

                                                                           gatTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAAt                                                                           <  2:2131592/150‑1 (MQ=255)

                                                                            atTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAAtt                                                                          >  2:662357/1‑150 (MQ=255)

                                                                              tCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAg                                                                        >  1:995615/1‑150 (MQ=255)

                                                                              tCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGAGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAg                                                                        >  2:550891/1‑150 (MQ=255)

                                                                                tGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGcc                                                                      >  1:2715024/1‑150 (MQ=255)

                                                                                tGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGcc                                                                      <  1:3137959/150‑1 (MQ=255)

                                                                                  gggACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCa                                                                    >  2:819129/1‑150 (MQ=255)

                                                                                   ggACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGAGCTGGAGCGTACTCTTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCaa                                                                   <  2:2465361/150‑1 (MQ=255)

                                                                                    gACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAac                                                                  <  2:3048862/150‑1 (MQ=255)

                                                                                        aaaCTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAacgacg                                                              >  2:2736012/1‑150 (MQ=255)

                                                                                          aCTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGAcgcc                                                            <  2:4578946/150‑1 (MQ=255)

                                                                                            tGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGAGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGAcgccgc                                                          <  2:1043449/150‑1 (MQ=255)

                                                                                             ggAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCt                                                         >  2:4284605/1‑150 (MQ=255)

                                                                                                aGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGAt                                                      <  1:1456576/150‑1 (MQ=255)

                                                                                                aGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGAt                                                      >  1:970002/1‑150 (MQ=255)

                                                                                                 gATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATg                                                     >  1:1267105/1‑150 (MQ=255)

                                                                                                 gATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGAGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATg                                                     <  2:2671809/150‑1 (MQ=255)

                                                                                                  aTCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATgg                                                    >  1:1738621/1‑150 (MQ=255)

                                                                                                   tCAGCCCGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGc                                                   >  2:1235189/1‑150 (MQ=255)

                                                                                                   tCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGc                                                   <  2:1798707/150‑1 (MQ=255)

                                                                                                     agccagccGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCgg                                                 >  2:2311799/1‑150 (MQ=255)

                                                                                                      gccagccGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGt                                                <  2:2121212/150‑1 (MQ=255)

                                                                                                       ccagccGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTc                                               <  1:3370544/150‑1 (MQ=255)

                                                                                                       ccagccGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTc                                               <  2:756965/150‑1 (MQ=255)

                                                                                                       ccagccGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGAGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTc                                               <  1:1816950/150‑1 (MQ=255)

                                                                                                        cagccGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCg                                              >  1:2797444/1‑150 (MQ=255)

                                                                                                        cagccGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCg                                              <  2:3841092/150‑1 (MQ=255)

                                                                                                          gccGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGAGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGtt                                            >  2:3008190/1‑150 (MQ=255)

                                                                                                          gccGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGAGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGtt                                            >  2:4388794/1‑150 (MQ=255)

                                                                                                            cGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTcc                                          >  2:4518923/1‑150 (MQ=255)

                                                                                                             gTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCa                                         >  2:1079171/1‑150 (MQ=255)

                                                                                                                 aaCAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGcca                                     >  1:2140248/1‑150 (MQ=255)

                                                                                                                   cAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGccacc                                   <  2:2927439/150‑1 (MQ=255)

                                                                                                                      gCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAat                                <  1:588732/150‑1 (MQ=255)

                                                                                                                           cACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATAtgtg                           <  2:3676153/150‑1 (MQ=255)

                                                                                                                             ccGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATgtggg                         <  1:4602666/150‑1 (MQ=255)

                                                                                                                              cGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATgtgggt                        <  2:2816950/150‑1 (MQ=255)

                                                                                                                                  gTACGTAACGCCTGGCTGGAGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGAt                    >  2:392438/1‑150 (MQ=255)

                                                                                                                                     cGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGcc                 >  2:4317099/1‑150 (MQ=255)

                                                                                                                                     cGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGcc                 <  2:3355453/150‑1 (MQ=255)

                                                                                                                                        aaCGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAAc              >  2:3376898/1‑150 (MQ=255)

                                                                                                                                        aaCGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAAc              >  1:3556635/1‑150 (MQ=255)

                                                                                                                                         aCGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACt             <  2:1590800/150‑1 (MQ=255)

                                                                                                                                          cGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACtt            >  1:1879650/1‑150 (MQ=255)

                                                                                                                                           gCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACttt           >  1:203300/1‑150 (MQ=255)

                                                                                                                                           gCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACttt           <  1:3649573/150‑1 (MQ=255)

                                                                                                                                            ccTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTa          >  2:4635272/1‑150 (MQ=255)

                                                                                                                                            ccTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTa          >  1:1378887/1‑150 (MQ=255)

                                                                                                                                             ctggctggCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTaa         <  1:3724731/150‑1 (MQ=255)

                                                                                                                                              tggctggcGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTaaa        >  1:1009777/1‑150 (MQ=255)

                                                                                                                                              tggctggcGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTaaa        >  2:4270477/1‑150 (MQ=255)

                                                                                                                                              tggctggcGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTaaa        >  2:1484061/1‑150 (MQ=255)

                                                                                                                                              tggctggcGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTaaa        >  1:40654/1‑150 (MQ=255)

                                                                                                                                               ggctggaGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAag       <  2:3991024/150‑1 (MQ=255)

                                                                                                                                                gctggcGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAaga      >  2:1138086/1‑150 (MQ=255)

                                                                                                                                                gctggcGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAaga      >  1:4257229/1‑150 (MQ=255)

                                                                                                                                                  tggagctggagcGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAgaga    >  1:2579043/1‑150 (MQ=255)

                                                                                                                                                    gcgcTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACg  <  2:3582279/150‑1 (MQ=255)

CGCCCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACG  >  CP009273/2805271‑2805568

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