BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_008463	3,673,534	3 bp→CTC	intergenic (+83/+13)	fabI → / ← PA14_RS16690	enoyl‑ACP reductase FabI/SurA N‑terminal domain‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_008463	3,673,534	0	T	C	99.8%	1638.3 / ‑4.2	422	intergenic (+83/+15)	fabI/PA14_RS16690	enoyl‑ACP reductase FabI/SurA N‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (229/191); total (230/192)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01
*	NC_008463	3,673,535	0	C	T	96.5%	1758.4 / 24.4	436	intergenic (+84/+14)	fabI/PA14_RS16690	enoyl‑ACP reductase FabI/SurA N‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand): ref base C (12/3); new base T (229/190); total (242/194)
Fisher's exact test for biased strand distribution p-value = 6.47e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	NC_008463	3,673,536	0	T	C	96.6%	1620.7 / 35.7	436	intergenic (+85/+13)	fabI/PA14_RS16690	enoyl‑ACP reductase FabI/SurA N‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand): ref base T (12/3); new base C (230/191); total (242/194)
Fisher's exact test for biased strand distribution p-value = 6.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 436 total aligned reads displayed.

CCTCTACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGNGTAAACGGCCTTTTTCT‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGCT  >  NC_008463/3673394‑3673679

|||

ccTCTACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCAtct                                                                                                                                          >  2:8773162/1‑150 (MQ=255)

     aCGTCGATGGCGGCTTCAACCCTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGtc                                                                                                                                     <  2:11108141/150‑1 (MQ=255)

      cGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCATTTTTCATGGGCGGAATAAAGAAAAGGCCGCTTAAGCGGCCTTTTCTC‑TGCATCTCGGtct                                                                                                                                    <  1:10875234/150‑1 (MQ=255)

       gTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAATTGGCCGCTTCTAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGtctc                                                                                                                                   <  1:5956770/150‑1 (MQ=255)

        tCGATGGCGGCTTCAACACTCCCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCa                                                                                                                                  <  2:5249105/150‑1 (MQ=255)

        tCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCa                                                                                                                                  >  1:5147722/1‑150 (MQ=255)

        tCGATGGCGGCTTCAACACAACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCa                                                                                                                                  >  2:163634/1‑150 (MQ=255)

        tCGAGGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCCTGAAAAATGGCCGCTTCGAGCTGCCTTTTTCATGGGCGGACTAAAGACAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCa                                                                                                                                  <  2:2938099/150‑1 (MQ=255)

          gATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGt                                                                                                                                <  1:9538655/150‑1 (MQ=255)

          gATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTACTTGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGt                                                                                                                                >  1:7119743/1‑150 (MQ=255)

           aTGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTa                                                                                                                               <  1:7198532/150‑1 (MQ=255)

            tGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTAt                                                                                                                              >  1:7899913/1‑150 (MQ=255)

                 gcTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTtctc                                                                                                                         >  2:1883652/1‑150 (MQ=255)

                 gcTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTtctc                                                                                                                         >  1:6339803/1‑150 (MQ=255)

                    tCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGAc                                                                                                                      <  1:1616976/150‑1 (MQ=255)

                    tCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGAc                                                                                                                      <  1:5531471/150‑1 (MQ=255)

                    tCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGAc                                                                                                                      >  2:8259574/1‑150 (MQ=255)

                      aaCACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCTTGGGCGTAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACtt                                                                                                                    <  1:6610581/150‑1 (MQ=255)

                       acacTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCACATTAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAACCAAAAGGCCGCGTACGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTc                                                                                                                   <  1:10339667/150‑1 (MQ=255)

                           tACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGcc                                                                                                               <  1:5862710/150‑1 (MQ=255)

                              cGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTg                                                                                                            >  2:2033835/1‑150 (MQ=255)

                                   tGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTg                                                                                                       >  2:4199539/1‑150 (MQ=255)

                                     ggCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCa                                                                                                     <  2:4663760/150‑1 (MQ=255)

                                       cccGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAAc                                                                                                   <  2:9189175/150‑1 (MQ=255)

                                        ccGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACt                                                                                                  >  1:2337704/1‑150 (MQ=255)

                                         cGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTg                                                                                                 >  1:8356279/1‑150 (MQ=255)

                                           cTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAc                                                                                               <  2:6039492/150‑1 (MQ=255)

                                            tGGACGACGACTGATCACCCCGCTCATGAAATATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCATTGTCCTGCAACTGAcg                                                                                              >  2:8891663/1‑150 (MQ=255)

                                              gacgacgacTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAcggc                                                                                            >  2:7424227/1‑150 (MQ=255)

                                              gacgacgacTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAcggc                                                                                            <  1:6017667/150‑1 (MQ=255)

                                                  acgacTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCCGCAACTGACGGCGGaa                                                                                        <  2:7884095/150‑1 (MQ=255)

                                                     acTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAggc                                                                                     <  1:7722440/150‑1 (MQ=255)

                                                     acTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGAGGACTAAAGAAAAGTCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCTTCTTCCGCCTTGTCCTGCACCTGACGGCGGAAggc                                                                                     >  2:988164/1‑150 (MQ=255)

                                                        gATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAggcggc                                                                                  <  2:229820/150‑1 (MQ=255)

                                                        gATCAAACCGCTCAAGAAAATTGGCCGCTTCGAGCGGCCTATATCATCTGCGGAATAAAGAAAAGGCCGCGTAAGAGGCCTATTCTC‑TGCATCTCGGTCTCAGTATTACTCGACTTCCGACTGGTCATGCAAAAGAGGGAGGAAggcggc                                                                                  >  2:10830050/1‑150 (MQ=255)

                                                             ccccGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGt                                                                             <  2:686820/150‑1 (MQ=255)

                                                             ccccGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTCTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGt                                                                             >  1:3777995/1‑150 (MQ=255)

                                                              cccGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAACAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTc                                                                            <  1:6740239/150‑1 (MQ=255)

                                                              cccGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTc                                                                            >  1:8277161/1‑150 (MQ=255)

                                                                 gCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGc                                                                         >  1:424121/1‑150 (MQ=255)

                                                                  cTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAACGCAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGcc                                                                        <  1:3400883/150‑1 (MQ=255)

                                                                  cTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTACGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTCCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGcc                                                                        <  1:1511448/150‑1 (MQ=255)

                                                                  cTCATGAAAAAGGGCCGCTACGAGCGGCCTATTTCATGGGCGGAATACCGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCGCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGcc                                                                        <  2:7165760/150‑1 (MQ=255)

                                                                   tCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAACGAAAAGGCCGCGTACGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCt                                                                       <  1:1279058/150‑1 (MQ=255)

                                                                   tCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCt                                                                       >  1:7443731/1‑150 (MQ=255)

                                                                      tGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTc                                                                    <  2:3922471/150‑1 (MQ=255)

                                                                      tGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTc                                                                    >  2:6277434/1‑150 (MQ=255)

                                                                          aaaTGGCCGCTTCGAGCGGCATTTTTCTTGGGCGGAATAAAGACAAGGCCGCTTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCt                                                                <  1:3535333/150‑1 (MQ=255)

                                                                           aaTGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATCAAGAAAAGGCCGCGGAAGCGGCCTTTTCTC‑TGCATCGCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGACGTCTGCCTGTCCGCTg                                                               <  2:6310849/150‑1 (MQ=255)

                                                                            aTGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTgc                                                              >  1:8153746/1‑150 (MQ=255)

                                                                                   cTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGc                                                       >  2:10109738/1‑150 (MQ=255)

                                                                                    ttCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGcc                                                      <  2:9313802/150‑1 (MQ=255)

                                                                                    ttCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGcc                                                      >  2:10425911/1‑150 (MQ=255)

                                                                                     tCGAGCGGCCTTTTTCTTGGGCGTAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCa                                                     <  1:3152336/150‑1 (MQ=255)

                                                                                       gAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAgg                                                   >  2:4889164/1‑150 (MQ=255)

                                                                                        aGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGa                                                  >  2:2636529/1‑150 (MQ=255)

                                                                                        aGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGa                                                  <  1:3686909/150‑1 (MQ=255)

                                                                                        aGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGa                                                  <  2:1541381/150‑1 (MQ=255)

                                                                                         gCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGaa                                                 <  1:8598773/150‑1 (MQ=255)

                                                                                           ggCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAgc                                               <  1:6316587/150‑1 (MQ=255)

                                                                                           ggCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAgc                                               >  1:3093717/1‑150 (MQ=255)

                                                                                            gCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAgcg                                              >  1:4719979/1‑150 (MQ=255)

                                                                                               tttttCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAgcggcg                                           >  1:4648807/1‑150 (MQ=255)

                                                                                                ttttCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGcggcgg                                          >  1:6854914/1‑150 (MQ=255)

                                                                                                     aTGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACAt                                     <  2:6715398/150‑1 (MQ=255)

                                                                                                      tGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATg                                    >  1:5388719/1‑150 (MQ=255)

                                                                                                        ggCGGAATCAAGAAACGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGc                                  <  2:9452980/150‑1 (MQ=255)

                                                                                                          cgcaaTAAAGAACAGGCCGCTTAAGCGGCATTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCt                                <  1:9132467/147‑1 (MQ=255)

                                                                                                             aaTAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTtct                             >  1:6466235/1‑150 (MQ=255)

                                                                                                              aTAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGACGCGGCGGTACATGGCCTtctc                            >  2:10300191/1‑150 (MQ=255)

                                                                                                              aTAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTtctc                            >  2:3969831/1‑150 (MQ=255)

                                                                                                              aTAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCGtctc                            >  2:965060/1‑150 (MQ=255)

                                                                                                               tAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcg                           >  2:5684026/1‑150 (MQ=255)

                                                                                                                acaGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgt                          <  1:10729643/148‑1 (MQ=255)

                                                                                                                aaccaaacGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgt                          <  1:4890623/146‑1 (MQ=255)

                                                                                                                aaaGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgt                          >  1:285641/1‑150 (MQ=255)

                                                                                                                 aaGAACAGGCCGCGTAAGCGGCCTTTACTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGACTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtc                         <  1:2113076/150‑1 (MQ=255)

                                                                                                                   gAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtcgt                       >  2:8076841/1‑150 (MQ=255)

                                                                                                                    aaaaGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtcgtc                      >  2:1724163/1‑150 (MQ=255)

                                                                                                                     aaaGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtcgtcg                     >  1:10562350/1‑150 (MQ=255)

                                                                                                                     aaaGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtcgtcg                     <  2:10962483/150‑1 (MQ=255)

                                                                                                                     aaaGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCAAGGAAGCGGCGGTACATGGCCTTCtcgtcgtcg                     >  2:3321157/1‑150 (MQ=255)

                                                                                                                     aaaGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTAGTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtcgtcg                     <  1:9788070/150‑1 (MQ=255)

                                                                                                                       aGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGga                   >  2:169855/1‑150 (MQ=255)

                                                                                                                         gCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGgaga                 >  2:3097575/1‑150 (MQ=255)

                                                                                                                            gCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAg              >  1:8834800/1‑150 (MQ=255)

                                                                                                                              gTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcc            <  2:7702295/150‑1 (MQ=255)

                                                                                                                              gTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcc            <  1:899249/150‑1 (MQ=255)

                                                                                                                                 agcGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcctcc         <  2:10655201/148‑1 (MQ=255)

                                                                                                                                  gcGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcctcct        <  1:5735061/149‑1 (MQ=255)

                                                                                                                                  gcGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcctcct        <  2:3007054/149‑1 (MQ=255)

                                                                                                                                  gcGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcctcct        >  1:11188609/2‑150 (MQ=255)

                                                                                                                                  gcGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcctcct        >  1:4036117/2‑150 (MQ=255)

                                                                                                                                    ggCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCtcctcc      >  2:4409602/1‑150 (MQ=255)

                                                                                                                                     gCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCg     <  2:36010/150‑1 (MQ=255)

                                                                                                                                     gCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCg     >  1:1935836/1‑150 (MQ=255)

                                                                                                                                     gCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCg     >  1:9258699/1‑150 (MQ=255)

                                                                                                                                     gCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCg     >  2:9772060/1‑150 (MQ=255)

                                                                                                                                         ttttCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGCt  >  1:10132757/1‑150 (MQ=255)

                                                                                                                                         ttttCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGCt  >  1:9641206/1‑150 (MQ=255)

|||

CCTCTACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGNGTAAACGGCCTTTTTCT‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGCT  >  NC_008463/3673394‑3673679

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`