BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,753,293	G→T	13.1%	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	13.1%	818.8 / 70.9	249	V85V (GTC→GTA)	sufE	sulfur acceptor protein
Reads supporting (aligned to +/- strand): ref base G (94/122); new base T (16/17); total (110/139)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.80e-01

Only 100 of 272 total aligned reads displayed.

GCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTTGTGGACTTCTGTCTTCGTCGCGTAATT  >  CP009273/1753144‑1753441

 cGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACAGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAga                                                                                                                                                     <  1:2497545/150‑1 (MQ=255)

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

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

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

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

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

     tttGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAgacgac                                                                                                                                                 >  1:3813244/1‑150 (MQ=255)

     tttGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAgacgac                                                                                                                                                 >  1:2822372/1‑150 (MQ=255)

     tttGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAAAATCCTGCGGCGTCATCTGATCGTAGAGAATAAAgacgac                                                                                                                                                 >  1:3498216/1‑150 (MQ=255)

      ttGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGAcc                                                                                                                                                >  1:3946079/1‑150 (MQ=255)

           gcgAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAAt                                                                                                                                           <  1:177225/150‑1 (MQ=255)

              aaTTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAg                                                                                                                                        >  1:2213738/1‑150 (MQ=255)

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

                  gcgcGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCt                                                                                                                                    >  2:1193601/1‑150 (MQ=255)

                   cgcgAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCtt                                                                                                                                   >  2:934344/1‑150 (MQ=255)

                     cgAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTc                                                                                                                                 >  1:3625490/1‑150 (MQ=255)

                       aaTCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTcac                                                                                                                               >  2:1265154/1‑150 (MQ=255)

                       aaTCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTcac                                                                                                                               >  2:194643/1‑150 (MQ=255)

                       aaTCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTcac                                                                                                                               >  2:2510440/1‑150 (MQ=255)

                          catcGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAAt                                                                                                                            >  2:776272/1‑150 (MQ=255)

                            tcGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATcg                                                                                                                          >  2:1015630/1‑150 (MQ=255)

                             cGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGACTAAAGACGACCGCAATAAGCCCTTTCACAATcgc                                                                                                                         >  2:3539507/1‑150 (MQ=255)

                             cGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATcgc                                                                                                                         >  1:553034/1‑150 (MQ=255)

                              gCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATcgcc                                                                                                                        >  1:4339112/1‑150 (MQ=255)

                              gCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATcgcc                                                                                                                        <  1:3668711/150‑1 (MQ=255)

                              gCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATcgcc                                                                                                                        <  1:4085393/150‑1 (MQ=255)

                              gCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATcgcc                                                                                                                        <  2:3780218/150‑1 (MQ=255)

                                ttCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATcgccgc                                                                                                                      <  1:3722604/150‑1 (MQ=255)

                                 tCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCa                                                                                                                     <  1:2564480/150‑1 (MQ=255)

                                  ccAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCAt                                                                                                                    <  2:4244164/150‑1 (MQ=255)

                                    aGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCg                                                                                                                  <  1:3816496/150‑1 (MQ=255)

                                    aGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCg                                                                                                                  <  1:3594178/150‑1 (MQ=255)

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

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

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

                                            tgAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGcc                                                                                                          >  2:1336500/1‑150 (MQ=255)

                                            tgAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGcc                                                                                                          >  2:313071/1‑150 (MQ=255)

                                             gAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGccc                                                                                                         >  2:4263790/1‑150 (MQ=255)

                                                cGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGt                                                                                                      <  1:4084550/150‑1 (MQ=255)

                                                  agaTGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAACTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTaa                                                                                                    >  2:1076367/1‑150 (MQ=255)

                                                     tGGGGTGTGCTGTTGGGTGTGCGCCATTTTTTCAAACCCCGGACGCACATCGAAATTGACAATATCCTGGGGGGTCATCTGATCGTATAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTc                                                                                                 <  1:424996/150‑1 (MQ=255)

                                                      ggggTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCCCGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCGGTTCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCa                                                                                                <  2:2005124/150‑1 (MQ=255)

                                                        ggTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCaat                                                                                              <  2:956139/150‑1 (MQ=255)

                                                        ggTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCaat                                                                                              <  1:2899282/150‑1 (MQ=255)

                                                          tGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGTACGCACATCGAAATTGACAATATCCTGCGGAGTCATCTGATCTTAGATAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCaataa                                                                                            <  1:2741313/150‑1 (MQ=255)

                                                            agaTGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAAtt                                                                                          >  1:949439/1‑150 (MQ=255)

                                                             gaTGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTc                                                                                         <  2:1478692/150‑1 (MQ=255)

                                                               tGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTccc                                                                                       >  1:2546171/1‑150 (MQ=255)

                                                                 ttGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTg                                                                                     >  2:3052928/1‑150 (MQ=255)

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

                                                                       gAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCAtt                                                                               <  1:3139733/150‑1 (MQ=255)

                                                                             cATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGgcg                                                                         >  2:3447974/1‑150 (MQ=255)

                                                                              aTTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGgcgc                                                                        >  2:4205045/1‑150 (MQ=255)

                                                                                 ttttCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATg                                                                     >  2:1457566/1‑150 (MQ=255)

                                                                                   ttCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGAc                                                                   >  1:637145/1‑150 (MQ=255)

                                                                                   ttCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGAc                                                                   >  2:371385/1‑150 (MQ=255)

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

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

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

                                                                                          cACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCcac                                                            >  1:889218/1‑150 (MQ=255)

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

                                                                                             ggACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACAcc                                                         >  1:3254548/1‑150 (MQ=255)

                                                                                             ggACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACAcc                                                         >  2:4125968/1‑150 (MQ=255)

                                                                                             ggACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACAcc                                                         >  2:1442344/1‑150 (MQ=255)

                                                                                             ggACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACAcc                                                         <  2:708119/150‑1 (MQ=255)

                                                                                               aCGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTg                                                       <  1:3425781/150‑1 (MQ=255)

                                                                                                cGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGa                                                      >  1:3165443/1‑150 (MQ=255)

                                                                                                 gCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGAc                                                     >  2:4050630/1‑150 (MQ=255)

                                                                                                 gCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGAc                                                     <  2:310284/150‑1 (MQ=255)

                                                                                                 gCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGAc                                                     <  2:4246568/150‑1 (MQ=255)

                                                                                                 gCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGAc                                                     >  2:1785080/1‑150 (MQ=255)

                                                                                                  cacaTCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGAct                                                    >  2:2705453/1‑150 (MQ=255)

                                                                                                  cacaTCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGAct                                                    <  2:2822328/150‑1 (MQ=255)

                                                                                                  cacaTCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGAct                                                    <  2:1755732/150‑1 (MQ=255)

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

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

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

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

                                                                                                       cGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGAc                                               <  2:2664273/150‑1 (MQ=255)

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

                                                                                                        gAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACa                                              <  1:3053433/150‑1 (MQ=255)

                                                                                                           aTTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGcc                                           <  1:1015630/150‑1 (MQ=255)

                                                                                                             tGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCt                                         >  2:3857800/1‑150 (MQ=255)

                                                                                                             tGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCCTTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCt                                         >  1:4147233/1‑150 (MQ=255)

                                                                                                              gACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTg                                        >  2:406535/1‑150 (MQ=255)

                                                                                                                 aaTATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAAt                                     <  2:2640563/150‑1 (MQ=255)

                                                                                                                  atatCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATg                                    >  1:3292313/1‑150 (MQ=255)

                                                                                                                  atatCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATg                                    <  2:4339112/150‑1 (MQ=255)

                                                                                                                         gcggcgTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTAtttt                             <  2:2185029/150‑1 (MQ=255)

                                                                                                                             cgTCATCTGATCGTAGAGAATAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTTGTgg                         <  2:461251/150‑1 (MQ=255)

                                                                                                                               tCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTTGTGGAc                       >  1:1904192/1‑150 (MQ=255)

                                                                                                                                  tCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTTGTGGACTTc                    <  1:83039/150‑1 (MQ=255)

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

                                                                                                                                       tCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTTGTGGACTTCTGTCt               <  1:3157572/150‑1 (MQ=255)

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

                                                                                                                                          tAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTTGTGGACTTCTGTCTtcg            <  2:685821/150‑1 (MQ=255)

                                                                                                                                             agaATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTTGTGGACTTCTGTCTtcgtcg         <  2:1607793/150‑1 (MQ=255)

                                                                                                                                                aTAAATACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTTGTGGACTTCTGTCTTCGTCGCGt      <  1:4502343/150‑1 (MQ=255)

                                                                                                                                                 tAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTTGTGGACTTCTGTCTTCGTCGCGTa     >  1:1370596/1‑150 (MQ=255)

                                                                                                                                                   aaGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTTGTGGACTTCTGTCTTCGTCGCGTAAt   >  1:3848122/1‑150 (MQ=255)

GCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAATTCCCTGGGCATTCTGGCGCATGACAATCCACACCTGACTCTGACAGCCCTGAATGCTATTTTGTGGACTTCTGTCTTCGTCGCGTAATT  >  CP009273/1753144‑1753441

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