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%	1627.0 / NA	421	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 (222/199); total (222/199)
*	NC_008463	3,673,535	0	C	T	97.2%	1752.3 / 19.7	433	intergenic (+84/+14)	fabI/PA14_RS16690	enoyl‑ACP reductase FabI/SurA N‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand): ref base C (10/2); new base T (221/199); total (232/201)
Fisher's exact test for biased strand distribution p-value = 4.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01
*	NC_008463	3,673,536	0	T	C	97.2%	1602.6 / 27.9	433	intergenic (+85/+13)	fabI/PA14_RS16690	enoyl‑ACP reductase FabI/SurA N‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand): ref base T (10/2); new base C (222/199); total (232/201)
Fisher's exact test for biased strand distribution p-value = 4.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 433 total aligned reads displayed.

ATCCTCTACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGNGTAAACGGCCTTTTTCT‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGCTC  >  NC_008463/3673392‑3673680

|||

ctccgatACGTCGATGGCGGCTTCAACACTCCCGCCATGGGCCCGCTGGACGCCGACTGATCACCCCGCTCAGGAAAATTGGCCGCTTCGAGCGGCCTTTTTCTTGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCAtc                                                                                                                                            <  1:9488264/145‑1 (MQ=255)

 tCCTCTACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGTATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCAtct                                                                                                                                           >  1:4896467/1‑151 (MQ=255)

   cactACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCACATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCAGGGGCGGAATAAAGACAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCg                                                                                                                                         <  2:6468077/149‑1 (MQ=255)

    tctcCGTCGATGGCGGCTACAACACTACCCCCATGGGCCCGCTGGACGACGCCTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAACGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCgg                                                                                                                                        <  1:6350300/151‑1 (MQ=255)

        cGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCGCGGtctc                                                                                                                                    >  1:4466948/1‑151 (MQ=255)

        cGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCAGCTCGGtctc                                                                                                                                    >  2:10017264/1‑151 (MQ=255)

        cGTCGATGGCGGCTACAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAACGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGtctc                                                                                                                                    <  2:10210455/151‑1 (MQ=255)

         gTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCAAATAGGTCTCa                                                                                                                                   >  1:3219331/1‑151 (MQ=255)

            gATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAATGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTa                                                                                                                                >  2:457103/1‑151 (MQ=255)

            gATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTa                                                                                                                                >  1:7104611/1‑151 (MQ=255)

            gATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCGCGGTCTCAGTa                                                                                                                                >  1:1930747/1‑151 (MQ=255)

             aTGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTAt                                                                                                                               <  2:9667871/151‑1 (MQ=255)

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

                 cggcTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTtct                                                                                                                           <  2:5538786/151‑1 (MQ=255)

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

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

                       cAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACtt                                                                                                                     <  2:1949737/151‑1 (MQ=255)

                        aaCACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTc                                                                                                                    >  1:10028260/1‑151 (MQ=255)

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

                           ccaaCCGCCCTGGGCCCGCTGGACGACGCCTGAACACCCCCCTCATCAAAAATGGCCGCTTCGAGCGGCATTTTTCATGGGCGTAAAAAACAAAAGGCCGCTTACGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGc                                                                                                                 <  1:1018378/148‑1 (MQ=255)

                           aatACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCCTGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGCATAAAGAAAAGGCCGAGTAAGCGGCCTTTTCTC‑TGCATCACGGGCTCAGTATTTCTCGACTTCAGc                                                                                                                 >  2:5194473/3‑151 (MQ=255)

                             tACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCATTTTTCATGGGCGGAATAAAGAAACGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCt                                                                                                               <  1:4094521/151‑1 (MQ=255)

                              aCCGCCATGGGCCCGCTGGACGACGACTTATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCtt                                                                                                              <  1:5134899/151‑1 (MQ=255)

                                   cATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAATAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCt                                                                                                         <  2:3034409/151‑1 (MQ=255)

                                   cATGGGCCCGCTGGACGACGACTGATCACCCCGCTCAGGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGCAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCt                                                                                                         <  1:2064443/151‑1 (MQ=255)

                                    aTGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTg                                                                                                        >  2:1167198/1‑151 (MQ=255)

                                     tGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAATAAAAGGCCGCGTAAGAGGCCTTTTCTC‑TGCAACACGGGATCAGTATTTCTCGACTTCCGCCTTGTCCTGc                                                                                                       >  1:790912/1‑151 (MQ=255)

                                         cccGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACt                                                                                                   <  2:2720284/151‑1 (MQ=255)

                                           cGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGa                                                                                                 >  2:787778/1‑151 (MQ=255)

                                           agCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTGTCATGGGCGGAATAAAGAAAAGGCCGCGTGAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTAGTTCTCGACTTCCGCCTTGTCCTGCTACTGa                                                                                                 >  2:9861687/2‑151 (MQ=255)

                                            gCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAc                                                                                                <  2:1558790/151‑1 (MQ=255)

                                             cTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAcg                                                                                               >  2:10595140/1‑151 (MQ=255)

                                               ggACTACGACTGATCACCCCGCTAATCAAAAATGGCCGCTTCGAGCGGCCTTTTTAATGGGCGTAATAAATAAACGGCCGCTTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAcggc                                                                                             <  1:9267165/151‑1 (MQ=255)

                                               ggACGACGACTGATCACCCCGCACATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAcggc                                                                                             <  1:2353620/151‑1 (MQ=255)

                                                  cgacgacTGATCACCCCGCTCCTGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAAGCCAGAACAGGCCGCGTACGCGGCCATTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGa                                                                                          <  2:8068639/151‑1 (MQ=255)

                                                  cgacgacTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGa                                                                                          >  2:628856/1‑151 (MQ=255)

                                                   gactacTGTTCACCCCGAACTTGAAAATTGGCCTCTACGCGCTGCCTTTTTCTTGGGGGGAATAACGAACCGGCCGCGTCAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGaa                                                                                         <  1:3093323/151‑1 (MQ=255)

                                                   gacgacTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGaa                                                                                         <  2:10760645/151‑1 (MQ=255)

                                                     cgacTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAgg                                                                                       >  1:2583382/1‑151 (MQ=255)

                                                          gATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATACAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGgcggcg                                                                                  <  1:9088341/151‑1 (MQ=255)

                                                          gATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACGTCCGCCTTGTCCTGCAACTGACGGCGGAAGgcggcg                                                                                  <  1:1087567/151‑1 (MQ=255)

                                                           aTCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTGCAGGGGACGCATAAAGACAAGGCCGCGTAAGCGGCCGTGTCTC‑AGCATCTCGGACTCAGTATTTCTCGCATTCCGCCTTGTCCAGCAACTGACGGCGGAAGGCGGCGa                                                                                 >  2:7162757/1‑151 (MQ=255)

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

                                                                cccGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCt                                                                            >  1:223344/1‑151 (MQ=255)

                                                                    cTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGACGGCGGCGAAGTCTGCCt                                                                        >  2:8417613/1‑151 (MQ=255)

                                                                    cTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCt                                                                        >  1:5848177/1‑151 (MQ=255)

                                                                      cATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTGTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGt                                                                      >  1:6067725/1‑151 (MQ=255)

                                                                       aTGAAAATTGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTc                                                                     >  1:5545636/1‑151 (MQ=255)

                                                                        tGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTcc                                                                    >  2:2904453/1‑151 (MQ=255)

                                                                           aaaaTGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCt                                                                 >  1:1333589/1‑151 (MQ=255)

                                                                               tGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTgcgc                                                             <  1:4322142/151‑1 (MQ=255)

                                                                                ggCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGcgcg                                                            >  1:2199533/1‑151 (MQ=255)

                                                                                ggCCGCTTCGAGCGGCCTATTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGcgcg                                                            >  2:4184937/1‑151 (MQ=255)

                                                                                    gCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGc                                                        >  1:3643023/1‑151 (MQ=255)

                                                                                       tCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAg                                                     <  1:729597/151‑1 (MQ=255)

                                                                                          aGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGaa                                                  <  1:1199808/151‑1 (MQ=255)

                                                                                            cGGCCTTTTTCATGGGCGGAATAACGAAATGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAgc                                                <  1:65414/151‑1 (MQ=255)

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

                                                                                                cTTTTTCAGGGGCGGAATCACGACAAGGCCGCGTAAGCGGCCATTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAgcggcg                                            <  2:6669618/151‑1 (MQ=255)

                                                                                                 tttttCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGcggcgg                                           <  1:1798794/151‑1 (MQ=255)

                                                                                                     tCATGGGCGGAATAAAGAACAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACa                                       <  2:1792315/151‑1 (MQ=255)

                                                                                                       aTGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATg                                     >  1:3512488/1‑151 (MQ=255)

                                                                                                       aTGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATg                                     <  1:8656560/151‑1 (MQ=255)

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

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

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

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

                                                                                                           gCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCt                                 <  2:4916295/151‑1 (MQ=255)

                                                                                                            cgcaaTAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCtt                                >  2:7544065/4‑151 (MQ=255)

                                                                                                              gAATAAAGAAACGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTtct                              <  1:8511360/151‑1 (MQ=255)

                                                                                                                aTAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcg                            <  2:7912212/151‑1 (MQ=255)

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

                                                                                                                   aaGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACATCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGACTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtcg                         >  2:8539091/1‑151 (MQ=255)

                                                                                                                    aGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGACGCGGCGGTACATGGCCTTCtcgtcgt                        >  2:5758936/1‑151 (MQ=255)

                                                                                                                    aGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGACGCGGCGGTACATGGCCTTCtcgtcgt                        >  2:9168749/1‑151 (MQ=255)

                                                                                                                     gAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtcgtc                       >  1:2753144/1‑151 (MQ=255)

                                                                                                                     caaaaGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGACTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtcgtc                       <  2:3665379/150‑1 (MQ=255)

                                                                                                                      aaaaGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtcgtcg                      <  1:7711224/151‑1 (MQ=255)

                                                                                                                       aaaGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCgg                     <  1:9033678/151‑1 (MQ=255)

                                                                                                                       aaaGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCgg                     >  2:2413428/1‑151 (MQ=255)

                                                                                                                         aGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGgag                   >  1:11069204/1‑151 (MQ=255)

                                                                                                                          ggCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGgaga                  >  1:2864949/1‑151 (MQ=255)

                                                                                                                          ggCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGgaga                  >  2:4369019/1‑151 (MQ=255)

                                                                                                                               cgTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcc             >  2:2144486/2‑151 (MQ=255)

                                                                                                                                 tAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcctc           >  1:8402620/1‑151 (MQ=255)

                                                                                                                                  aagcGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcctcc          >  2:5078532/4‑151 (MQ=255)

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

                                                                                                                                  aagcGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcctcc          >  2:6060202/4‑151 (MQ=255)

                                                                                                                                  aagcGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcctcc          >  2:706926/4‑151 (MQ=255)

                                                                                                                                  aagcGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcctcc          >  2:10998397/4‑151 (MQ=255)

                                                                                                                                  aagcGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcctca          >  2:5897655/4‑150 (MQ=255)

                                                                                                                                  aagcGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAACCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGCAGATAGcctcc          >  2:10125958/4‑151 (MQ=255)

                                                                                                                                   cgcGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcctcct         <  2:3172987/149‑1 (MQ=255)

                                                                                                                                     cGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCtcctcc       >  2:8311463/1‑151 (MQ=255)

                                                                                                                                      tgCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCg      >  2:1427473/2‑151 (MQ=255)

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

                                                                                                                                      ggCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCg      >  1:3378440/1‑151 (MQ=255)

                                                                                                                                      ggCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCg      >  1:4974382/1‑151 (MQ=255)

                                                                                                                                          cttttCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGCt   >  1:4652836/2‑151 (MQ=255)

                                                                                                                                           ttttCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGCTc  <  2:3843784/151‑1 (MQ=255)

|||

ATCCTCTACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGNGTAAACGGCCTTTTTCT‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGCTC  >  NC_008463/3673392‑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: `‑`