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	100.0%	1896.3 / NA	485	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/0); new base C (238/247); total (238/247)
*	NC_008463	3,673,535	0	C	T	95.8%	2047.2 / 41.6	506	intergenic (+84/+14)	fabI/PA14_RS16690	enoyl‑ACP reductase FabI/SurA N‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand): ref base C (14/7); new base T (237/245); total (252/254)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01
*	NC_008463	3,673,536	0	T	C	95.9%	1851.5 / 55.8	506	intergenic (+85/+13)	fabI/PA14_RS16690	enoyl‑ACP reductase FabI/SurA N‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand): ref base T (14/7); new base C (238/247); total (252/254)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 506 total aligned reads displayed.

TCCTCTACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGNGTAAACGGCCTTTTTCT‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGCTC  >  NC_008463/3673393‑3673680

|||

tCCTCTACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAATAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCAtct                                                                                                                                           >  2:3437553/1‑151 (MQ=255)

tCCTCTACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTATTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCAtct                                                                                                                                           >  1:1921308/1‑151 (MQ=255)

 ccactACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCAtctc                                                                                                                                          <  2:2453630/148‑1 (MQ=255)

 ccTCTACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAACAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCAtctc                                                                                                                                          <  2:1919312/151‑1 (MQ=255)

 ccTCTACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTACGCGGCCTTTTCTC‑TGCAtctc                                                                                                                                          <  1:7186915/151‑1 (MQ=255)

   tctAGGTCGATGGCGGCTACAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATCAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCgg                                                                                                                                        <  1:8057351/151‑1 (MQ=255)

   tctACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCgg                                                                                                                                        >  1:1514105/1‑151 (MQ=255)

    ctACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGt                                                                                                                                       >  1:9269678/1‑151 (MQ=255)

     tACGTCGATGGCGGCTACAACACTACCGCCATGGGCCCGCTGGACGACGACTGAACACCCCGCACATTAAAATTGGCCGCTGCGCGCGGCATTTTTCTTGGGCGTAATAAACAAACGCCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGtc                                                                                                                                      <  1:2095073/151‑1 (MQ=255)

       cGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCAGGACTACGACTTATCACCCCCCTCATGAAAATTTGCCGCTTCGAGCGGCCTTTTTCAGGGGCGGAATAAAGAAATGGCAGCGTAAGCGGCCTTGTCTC‑TGCATCTCGGtctc                                                                                                                                    >  1:5408786/1‑151 (MQ=255)

        gTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCa                                                                                                                                   >  2:11296489/1‑151 (MQ=255)

           gATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTa                                                                                                                                <  1:11099673/151‑1 (MQ=255)

            aTGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGTTCACCCCGCTCATGAAATATGGCCGCTTCGAGCTGCCTTTTTCATGGGCGGAATATAGAAAAGGCCGCGTAAGCGGCCTGTTCGC‑TGGATCTCGGTCTCAGTAt                                                                                                                               >  2:10597345/1‑151 (MQ=255)

              ggcggcTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTAttt                                                                                                                             <  2:6565100/151‑1 (MQ=255)

               gcggcTTCAACCCTACCGCCAGGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTACGAGCGGCATTTTTTATGGGCGTAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTtc                                                                                                                            <  1:2141155/151‑1 (MQ=255)

               gcggcTTCAACACTACCGCCAGGGGCCCGCTGTACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTACGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTtc                                                                                                                            <  1:9682500/151‑1 (MQ=255)

                 ggcTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTtctc                                                                                                                          >  1:2449026/1‑151 (MQ=255)

                 ggcTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTtctc                                                                                                                          <  2:1553153/151‑1 (MQ=255)

                    ttCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGAc                                                                                                                       >  1:3836922/1‑151 (MQ=255)

                    ttCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGAc                                                                                                                       >  2:5419058/1‑151 (MQ=255)

                    tacAACACTCCCGCCTTGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGTAATAACGAAACGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGAc                                                                                                                       <  1:7958326/149‑1 (MQ=255)

                        acacTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTcc                                                                                                                   <  2:5802528/151‑1 (MQ=255)

                        acacTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTcc                                                                                                                   <  2:5318649/151‑1 (MQ=255)

                         cacTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTACCg                                                                                                                  <  2:5875252/151‑1 (MQ=255)

                            tACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCACATGAAAAATGGCCGCTTCGAGCGGCCTTTTACATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCt                                                                                                               <  1:6706234/151‑1 (MQ=255)

                               cGCCATGGGCCCGCTGGACGACGAATGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGTAATAAACAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGt                                                                                                            <  1:4442573/151‑1 (MQ=255)

                                  cATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAACGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCt                                                                                                         <  1:1764194/151‑1 (MQ=255)

                                  cATGGGCCCGCGGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTACGAGCGGCCTTTTTCATGGGCGGAATAACGAACCGGCCGCGTACGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCt                                                                                                         <  1:8245022/151‑1 (MQ=255)

                                     gggCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCa                                                                                                      >  2:3030832/1‑151 (MQ=255)

                                     gggCCCGCTGGACGACGACTGATCACCCCGCACATGAAAAATGGCCGCTTCGAGCGGCATTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCa                                                                                                      <  1:10737284/151‑1 (MQ=255)

                                          cGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGACTAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGa                                                                                                 >  2:7733115/1‑151 (MQ=255)

                                           gCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAc                                                                                                >  1:9381736/1‑151 (MQ=255)

                                           gCTGGACGACGACTGATCACCCCGCTCATCAAAAATGGCCGATTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAACAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAc                                                                                                <  1:5072744/151‑1 (MQ=255)

                                                     gacTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGACggc                                                                                      >  1:6895173/1‑151 (MQ=255)

                                                      acTGATCACCCCGCTAATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTAGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAggcg                                                                                     <  2:5053817/151‑1 (MQ=255)

                                                         taTCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAACGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGgcggcg                                                                                  <  1:2419981/150‑1 (MQ=255)

                                                          aTCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGa                                                                                 >  2:4564045/1‑151 (MQ=255)

                                                            cACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAg                                                                               >  2:10059093/1‑151 (MQ=255)

                                                             aCCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGt                                                                              >  1:76034/1‑151 (MQ=255)

                                                                   cTCATGACAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTTAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCt                                                                        <  1:11499555/151‑1 (MQ=255)

                                                                    tCAAGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATACCGAACAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTg                                                                       <  1:1707894/151‑1 (MQ=255)

                                                                     cATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAACGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGt                                                                      <  2:3136098/151‑1 (MQ=255)

                                                                      atcaaaaaTGCCCCCCTCGAGCGGCCTATTTCATGGGCGGAATAAACAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTc                                                                     <  1:1333289/148‑1 (MQ=255)

                                                                       tGAACAATGGCCGCTTCGAGCGGCCTTTGTCATGGGCGGAATAACGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTcc                                                                    <  2:6339470/151‑1 (MQ=255)

                                                                          aaaaTGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCt                                                                 <  1:6474430/151‑1 (MQ=255)

                                                                           aaaTGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATACAGCAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTg                                                                <  1:6256529/151‑1 (MQ=255)

                                                                           aaaTGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTg                                                                <  2:1793116/151‑1 (MQ=255)

                                                                             aTGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGACGGCGGCGAAGTCTGCCTGTCCGCTgcg                                                              >  2:6121547/1‑151 (MQ=255)

                                                                             aTGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTgcg                                                              <  1:10078848/151‑1 (MQ=255)

                                                                              tGGCCGCTTCGAGCGGCCATTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTgcgc                                                             <  1:10579487/151‑1 (MQ=255)

                                                                               ggCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATCAAGCAACGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGcgcg                                                            <  2:930026/151‑1 (MQ=255)

                                                                               ggCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTTTTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGcgcg                                                            >  2:7277699/1‑151 (MQ=255)

                                                                               ggCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGcgcg                                                            <  2:4046208/151‑1 (MQ=255)

                                                                               ggCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCGTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGcgcg                                                            >  1:2668809/1‑151 (MQ=255)

                                                                                 ccGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAACCGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGaa                                                          <  2:3721774/151‑1 (MQ=255)

                                                                                  cGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAg                                                         >  2:9492172/1‑151 (MQ=255)

                                                                                    cgtCGAGCGGCCTATTTCATGGGCGGAATAAAGAACAGGCCGCGTACGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGcc                                                       <  1:8540429/149‑1 (MQ=255)

                                                                                     ttCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCa                                                      >  1:9776068/1‑151 (MQ=255)

                                                                                       cGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAgg                                                    >  1:3366443/1‑151 (MQ=255)

                                                                                        gAGCGGCCTTTTTCATGGGCGGAATAAAGACAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGa                                                   <  2:10153095/151‑1 (MQ=255)

                                                                                        gAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGACGCCAGGa                                                   >  2:9689759/1‑151 (MQ=255)

                                                                                         aGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGaa                                                  >  1:9564606/1‑151 (MQ=255)

                                                                                            ggCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAgcg                                               >  2:4230337/1‑151 (MQ=255)

                                                                                             gCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAgcgg                                              <  2:11125191/151‑1 (MQ=255)

                                                                                             gCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAgcgg                                              >  2:3518478/1‑151 (MQ=255)

                                                                                             gCCTTGTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCAGTCCGCTGCGCGAAGCCAGGAAgcgg                                              >  2:7095910/1‑151 (MQ=255)

                                                                                              ccTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAgcggc                                             >  1:11169090/1‑151 (MQ=255)

                                                                                                tttttCATGGGCGCAATAAAGAAAAGGCCGCGTAAGCGCCCAATTCTC‑AGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGcggcgg                                           <  1:4198332/151‑1 (MQ=255)

                                                                                                     cATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACAt                                      >  2:10975070/1‑151 (MQ=255)

                                                                                                       tGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATgg                                    >  2:3553882/1‑151 (MQ=255)

                                                                                                        gggCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGc                                   >  2:8541676/1‑151 (MQ=255)

                                                                                                         ggCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGcc                                  >  1:77877/1‑151 (MQ=255)

                                                                                                              aaTAACGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTtctc                             <  1:10860813/151‑1 (MQ=255)

                                                                                                              aaTAAAGAAACGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTtctc                             <  1:1562531/151‑1 (MQ=255)

                                                                                                               aTAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcg                            >  2:2723110/1‑151 (MQ=255)

                                                                                                               aTAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGACTCAGGATTTCTCGACTTCCGCCTTGGCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcg                            <  1:3554667/151‑1 (MQ=255)

                                                                                                                tAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgt                           >  2:10347265/1‑151 (MQ=255)

                                                                                                                tAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgt                           >  2:11580776/1‑151 (MQ=255)

                                                                                                                 acaGAAAAGCCCGCGTAAGCGGCCTTTTCAC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtc                          <  1:8861611/149‑1 (MQ=255)

                                                                                                                 aaaGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtc                          <  1:908040/151‑1 (MQ=255)

                                                                                                                    gAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtcgtc                       <  1:570971/151‑1 (MQ=255)

                                                                                                                     aaacGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtcgtcg                      <  2:7586040/151‑1 (MQ=255)

                                                                                                                     aaaaGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtcgtcg                      >  2:4880970/1‑151 (MQ=255)

                                                                                                                         ggCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGgaga                  <  1:3041144/151‑1 (MQ=255)

                                                                                                                         ggCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCACGTCGTCGgaga                  >  2:7802505/1‑151 (MQ=255)

                                                                                                                          gCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGAt                 >  1:4797875/1‑151 (MQ=255)

                                                                                                                          gCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGACTCAGAATTTCTCGACTTCCGCCTTGTCCTGCCACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGAt                 <  1:4934516/151‑1 (MQ=255)

                                                                                                                           ccGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATa                >  1:5893025/1‑151 (MQ=255)

                                                                                                                           ccGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATa                >  2:8965297/1‑151 (MQ=255)

                                                                                                                            cgcgTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAg               <  2:6794424/151‑1 (MQ=255)

                                                                                                                            cgcgTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAg               >  2:8439764/1‑151 (MQ=255)

                                                                                                                             gCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGc              >  1:10203342/1‑151 (MQ=255)

                                                                                                                               gTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcct            >  2:1301967/1‑151 (MQ=255)

                                                                                                                               gTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcct            >  2:8898373/1‑151 (MQ=255)

                                                                                                                                 aagcGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcctcc          >  1:3960722/4‑151 (MQ=255)

                                                                                                                                    cGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCtcctcc       <  1:1002721/151‑1 (MQ=255)

                                                                                                                                     ggCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCg      <  2:3926686/151‑1 (MQ=255)

                                                                                                                                      gCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCgg     >  2:8088748/1‑151 (MQ=255)

                                                                                                                                         cttttCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGCt   <  1:10305761/150‑1 (MQ=255)

                                                                                                                                          ttttCTCTGCATCTCGGTCACAGTATTTCTCGACTTACGCCTTGTCCAGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGACTCCTCCGGCTc  >  1:6577935/1‑151 (MQ=255)

|||

TCCTCTACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGNGTAAACGGCCTTTTTCT‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGCTC  >  NC_008463/3673393‑3673680

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