BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,753,293	G→T	48.9%	V85V (GTC→GTA)	sufE ←	sulfur acceptor protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,753,293	0	G	T	48.9%	128.5 / inf	257	V85V (GTC→GTA)	sufE	sulfur acceptor protein
Reads supporting (aligned to +/- strand): ref base G (64/67); new base T (64/62); total (128/129)
Fisher's exact test for biased strand distribution p-value = 8.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

Only 100 of 277 total aligned reads displayed.

GCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTTGTGGACTTCTGTCTTCGTCGCGTAA  >  CP009273/1753144‑1753439

  ggCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAtac                                                                                                                                                  >  1:4159002/1‑150 (MQ=255)

  ggCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAgac                                                                                                                                                  <  2:2040252/150‑1 (MQ=255)

   gCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAtacg                                                                                                                                                 >  2:1554932/1‑150 (MQ=255)

   gCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAtacg                                                                                                                                                 >  1:2578871/1‑150 (MQ=255)

   gCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAtacg                                                                                                                                                 <  1:2266638/150‑1 (MQ=255)

   gCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAgacg                                                                                                                                                 <  2:11549/150‑1 (MQ=255)

    cTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAtacga                                                                                                                                                >  1:2098090/1‑150 (MQ=255)

    cTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAgacga                                                                                                                                                >  2:3932290/1‑150 (MQ=255)

     tttGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAtacgac                                                                                                                                               <  1:876781/150‑1 (MQ=255)

         gcgcGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCa                                                                                                                                           >  2:4594506/1‑150 (MQ=255)

            cgAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATa                                                                                                                                        >  2:1642686/1‑150 (MQ=255)

            cgAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATa                                                                                                                                        >  2:2023462/1‑150 (MQ=255)

             gAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATaa                                                                                                                                       >  2:177259/1‑150 (MQ=255)

              aaTTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAg                                                                                                                                      <  1:2798423/150‑1 (MQ=255)

                ttGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGcc                                                                                                                                    <  2:3067052/150‑1 (MQ=255)

                    gcgAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCttt                                                                                                                                >  1:3607401/1‑150 (MQ=255)

                     cgAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTc                                                                                                                               >  2:445403/1‑150 (MQ=255)

                      gAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGTTGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTca                                                                                                                              <  2:503150/150‑1 (MQ=255)

                        atcatcGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTcaca                                                                                                                            <  2:447955/150‑1 (MQ=255)

                          catcGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAAt                                                                                                                          >  2:2296815/1‑150 (MQ=255)

                           atcGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATc                                                                                                                         <  1:1399520/150‑1 (MQ=255)

                              gCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATcgcc                                                                                                                      <  2:621987/150‑1 (MQ=255)

                                  ccAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCAt                                                                                                                  >  2:1081980/1‑150 (MQ=255)

                                       ccTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTg                                                                                                             >  2:4184532/1‑150 (MQ=255)

                                       ccTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTg                                                                                                             >  2:3085918/1‑150 (MQ=255)

                                        cTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGt                                                                                                            >  1:1649260/1‑150 (MQ=255)

                                        cTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGt                                                                                                            >  2:735850/1‑150 (MQ=255)

                                        cTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGt                                                                                                            >  1:335575/1‑150 (MQ=255)

                                         ttGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTc                                                                                                           <  1:3412950/150‑1 (MQ=255)

                                          tgtgAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCg                                                                                                          >  1:4365546/1‑150 (MQ=255)

                                           gtgAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGc                                                                                                         >  2:1819165/1‑150 (MQ=255)

                                                 gagaTGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTa                                                                                                   <  2:3358002/150‑1 (MQ=255)

                                                 gagaTGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTa                                                                                                   <  2:2005080/150‑1 (MQ=255)

                                                    aTGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAAtt                                                                                                <  2:3247525/150‑1 (MQ=255)

                                                     tGGGGTGAGTTGTTGGGTGAGTGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTc                                                                                               <  1:1264533/150‑1 (MQ=255)

                                                      ggggTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCa                                                                                              <  1:2384610/150‑1 (MQ=255)

                                                      ggggTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCa                                                                                              <  1:2300846/150‑1 (MQ=255)

                                                      ggggTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCa                                                                                              <  1:4503407/150‑1 (MQ=255)

                                                         gTGAGATGTTGGGTGAGCGCCATTTTTTCANACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCaata                                                                                           >  1:4281777/1‑150 (MQ=255)

                                                          tGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCaataa                                                                                          <  1:2566503/150‑1 (MQ=255)

                                                            agaTGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAAtt                                                                                        >  2:3596523/1‑150 (MQ=255)

                                                               tGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTccc                                                                                     <  1:1700973/150‑1 (MQ=255)

                                                                gTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCt                                                                                    >  2:1518564/1‑150 (MQ=255)

                                                                gTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCt                                                                                    >  1:2768839/1‑150 (MQ=255)

                                                                  tGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTgg                                                                                  >  1:451671/1‑150 (MQ=255)

                                                                    ggTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGc                                                                                >  2:1213678/1‑150 (MQ=255)

                                                                    ggTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGc                                                                                <  2:2216521/150‑1 (MQ=255)

                                                                     gTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCa                                                                               <  2:3990252/150‑1 (MQ=255)

                                                                      tGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCAt                                                                              >  1:1466670/1‑150 (MQ=255)

                                                                      tGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCAt                                                                              <  1:3556533/150‑1 (MQ=255)

                                                                           gcCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTgg                                                                         <  1:1019313/150‑1 (MQ=255)

                                                                              aTTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGgcgc                                                                      >  1:1589339/1‑150 (MQ=255)

                                                                               ttttttCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCa                                                                     <  1:784765/150‑1 (MQ=255)

                                                                                 ttttCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATg                                                                   <  1:2981807/150‑1 (MQ=255)

                                                                                    tCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACa                                                                <  2:2145027/150‑1 (MQ=255)

                                                                                      aaaCCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAAt                                                              <  1:1845801/150‑1 (MQ=255)

                                                                                        aCCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATcc                                                            <  1:4447891/150‑1 (MQ=255)

                                                                                         ccACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCca                                                           <  1:4103484/150‑1 (MQ=255)

                                                                                            cGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCacac                                                        <  2:3436319/150‑1 (MQ=255)

                                                                                             ggACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACAcc                                                       >  1:3955082/1‑150 (MQ=255)

                                                                                             ggACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACAcc                                                       >  1:3342049/1‑150 (MQ=255)

                                                                                               aCGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTg                                                     <  2:3207620/150‑1 (MQ=255)

                                                                                               aCGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTg                                                     <  2:4176778/150‑1 (MQ=255)

                                                                                                  cacaTCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGAct                                                  >  1:135421/1‑150 (MQ=255)

                                                                                                   acaTCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGActc                                                 >  2:3000066/1‑150 (MQ=255)

                                                                                                   acaTCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGActc                                                 >  2:2369489/1‑150 (MQ=255)

                                                                                                   acaTCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGActc                                                 <  2:1125700/150‑1 (MQ=255)

                                                                                                    caTCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGActct                                                <  2:1160126/150‑1 (MQ=255)

                                                                                                    caTCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGActct                                                <  2:3172394/150‑1 (MQ=255)

                                                                                                    caTCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGActct                                                >  1:1691196/1‑150 (MQ=255)

                                                                                                       cGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGAc                                             >  1:3978391/1‑150 (MQ=255)

                                                                                                       cGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGAc                                             >  1:101037/1‑150 (MQ=255)

                                                                                                        gAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACa                                            >  1:2583353/1‑150 (MQ=255)

                                                                                                          aaTTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGc                                          >  1:522126/1‑150 (MQ=255)

                                                                                                          aaTTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGc                                          >  1:2595729/1‑150 (MQ=255)

                                                                                                            ttGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGccc                                        <  1:3197827/150‑1 (MQ=255)

                                                                                                             tGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCt                                       <  1:1929768/150‑1 (MQ=255)

                                                                                                                  atatCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATg                                  >  1:1080907/1‑150 (MQ=255)

                                                                                                                      ccTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTAt                              >  2:3505010/1‑150 (MQ=255)

                                                                                                                       cTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTAtt                             <  2:2852019/150‑1 (MQ=255)

                                                                                                                       cTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAAAAATTCCCTGGGCATTCTGGCGCATGACAAACCACACCTGACTCTGACAGCCCTGAATGCTAtt                             >  1:1908592/1‑150 (MQ=255)

                                                                                                                         gcggcgTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTAtttt                           <  2:4084838/150‑1 (MQ=255)

                                                                                                                          cggcgTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTtg                          >  2:833611/1‑150 (MQ=255)

                                                                                                                           ggcgTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTtgt                         <  1:923811/150‑1 (MQ=255)

                                                                                                                              gTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTTGTGGa                      >  1:2388960/1‑150 (MQ=255)

                                                                                                                                 aTCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTTGTGGACtt                   <  1:58791/150‑1 (MQ=255)

                                                                                                                                 aTCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTTGTGGACtt                   >  1:2701868/1‑150 (MQ=255)

                                                                                                                                  tCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTTGTGGACTTc                  >  1:2824686/1‑150 (MQ=255)

                                                                                                                                   cTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTTGTGGACTTCt                 <  1:3855006/150‑1 (MQ=255)

                                                                                                                                    tGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTTGTGGACTTCTg                >  2:1804495/1‑150 (MQ=255)

                                                                                                                                    tGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTTGTGGACTTCTg                <  1:4065070/150‑1 (MQ=255)

                                                                                                                                    tGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTTGTGGACTTCTg                <  1:4017204/150‑1 (MQ=255)

                                                                                                                                      aTCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTTGTGGACTTCTGTc              >  2:1518943/1‑150 (MQ=255)

                                                                                                                                       tCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTTGTGGACTTCTGTCt             <  2:3787919/150‑1 (MQ=255)

                                                                                                                                        cGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTTGTGGACTTCTGTCtt            >  2:2688073/1‑150 (MQ=255)

                                                                                                                                         gTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTTGTGGACTTCTGTCTtc           >  1:4550589/1‑150 (MQ=255)

                                                                                                                                            gagaATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTTGTGGACTTCTGTCTtcgtc        <  1:3784681/150‑1 (MQ=255)

                                                                                                                                             agaATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTTGTGGACTTCTGTCTtcgtcg       >  1:880730/1‑150 (MQ=255)

                                                                                                                                             agaATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTTGTGGACTTCTGTCTtcgtcg       >  1:2450489/1‑150 (MQ=255)

                                                                                                                                               aaTAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTTGTGGACTTCTGTCTTCGTcgcg     >  2:194579/1‑150 (MQ=255)

GCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTTGTGGACTTCTGTCTTCGTCGCGTAA  >  CP009273/1753144‑1753439

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