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%	1637.2 / NA	418	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 (227/191); total (227/191)
*	NC_008463	3,673,535	0	C	T	97.0%	1777.4 / 22.5	431	intergenic (+84/+14)	fabI/PA14_RS16690	enoyl‑ACP reductase FabI/SurA N‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand): ref base C (7/6); new base T (227/189); total (234/197)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01
*	NC_008463	3,673,536	0	T	C	97.0%	1619.6 / 32.0	431	intergenic (+85/+13)	fabI/PA14_RS16690	enoyl‑ACP reductase FabI/SurA N‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand): ref base T (7/6); new base C (227/191); total (234/197)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

Only 100 of 431 total aligned reads displayed.

ATCCTCTACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGNGTAAACGGCCTTTTTCT‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGCTC  >  NC_008463/3673392‑3673680

|||

aTCCTCTACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCTTGGGCGTAATAAATAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCAtc                                                                                                                                            <  1:8860778/151‑1 (MQ=255)

   ctctACGTCGATGGCGGCTTCAACACTACCGCCACGGGCCCGCTTGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAACGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCg                                                                                                                                         <  2:8000095/151‑1 (MQ=255)

    tctACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCgg                                                                                                                                        <  1:8623846/151‑1 (MQ=255)

     ctACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGt                                                                                                                                       <  1:1697196/151‑1 (MQ=255)

      tACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCCTGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGtc                                                                                                                                      >  1:9065650/1‑151 (MQ=255)

       aCGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGACAAGGCCGCGTAAGCGGCCATTTCTC‑TGCATCTCGGtct                                                                                                                                     <  2:6840819/151‑1 (MQ=255)

       aCGTCGATGGCGGCTTAAACACTACCGCCATGGGCCCGCTGGACGACGACTTATCACCCCGCTCATGAAAATTGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAACAAAAGGCCGCGTACGCGGCCTTTTCTC‑TGCATCTCGGtct                                                                                                                                     <  1:1060884/151‑1 (MQ=255)

        cGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAAAAAACAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGtctc                                                                                                                                    <  1:10289055/151‑1 (MQ=255)

        cGTCGAGGGCGGCTTCAACCCAACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAACGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGtctc                                                                                                                                    <  2:196104/151‑1 (MQ=255)

         gTCGATGGCGGCTACACCACTCCCGCCATGGGCCCGCTGGACGACGACGGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAAAAAAGAAACGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCa                                                                                                                                   <  2:7821510/151‑1 (MQ=255)

          tCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGtctctg                                                                                                                                  >  2:4388041/1‑149 (MQ=255)

          tCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAg                                                                                                                                  <  2:6509639/151‑1 (MQ=255)

           cGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGt                                                                                                                                 >  2:2811627/1‑151 (MQ=255)

            gATGGCGGCTTCAACCCGACCGCCATGGGCCCGCTGGACGACGACATATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCATTTCTC‑TGCATCTCGGTCTCAGTa                                                                                                                                <  1:3502037/151‑1 (MQ=255)

              tGGCGGCTTCAACCATACCGCCATGGGCCCTCTGGACGACGACTGATCACCCCTCTCAGGAAAAATGGCCGCTTCCCGCGGCCTTTTTAATGGTCGTAATAAAGAAACGGCCGCTAAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTAtt                                                                                                                              <  1:93509/151‑1 (MQ=255)

                gcggcTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAACGAAATGGCCGCGTACGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTtc                                                                                                                            <  1:4904034/151‑1 (MQ=255)

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

                   gcTACAACACTACCGCCTTGGGCCCGCTGGACGACGACTGATCACCCCCCTAATGATAAATGGCGGCTTCTTGCGGCATTTTACATGGGCTGAATAACCAAACGCCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCg                                                                                                                         <  1:1521027/151‑1 (MQ=255)

                     ttCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGAc                                                                                                                       <  1:7173204/151‑1 (MQ=255)

                      tCAACCCTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAACAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACt                                                                                                                      <  2:9649208/151‑1 (MQ=255)

                       cAACACTACCTCCAAGGGCCCGCTGGTCGACGACAGATCACCCCGCTCATGAAAAATGGCCGCTACGAGCGGCCATTTTCATGGGCGGAATAAAGAAAAGGCCGCGAAAGCGGCCTTTTCAC‑TGCATCTCGGTCTCAGTATTTCTCGACtt                                                                                                                     <  2:10643148/151‑1 (MQ=255)

                           acaaCCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGATCATCAAAAAGGGCCGCTTCGCGCGGCATTTTTCATGGGCGGAATAAACAAAAGGCCGCTTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGc                                                                                                                 <  1:5368704/148‑1 (MQ=255)

                               ccGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTACGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTg                                                                                                             <  1:9573223/151‑1 (MQ=255)

                                 gCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTc                                                                                                           >  2:6007937/1‑151 (MQ=255)

                                   cATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCt                                                                                                         >  1:10003991/1‑151 (MQ=255)

                                   cATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCt                                                                                                         >  1:1847372/1‑151 (MQ=255)

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

                                       ggCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCaa                                                                                                     >  1:3860226/1‑151 (MQ=255)

                                          ccGCTGGACGACGACTGATCCCCCCGCTCATCAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTg                                                                                                  <  2:6348306/151‑1 (MQ=255)

                                          ccGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTg                                                                                                  >  1:9297523/1‑151 (MQ=255)

                                            gCTGGACGACGACTGATCACCCCGCCCCTGACAAATGGCCGCGTCGAGCGGCCAATTACATGGGCGCAATAAAGCACAGGCCGCGTACGCGGCCTTTTCTC‑TGCATCTCGGTCTCACTATTTCTCGACTTCCGCCTGGGCCTGCAACTGAc                                                                                                <  2:5889989/151‑1 (MQ=255)

                                              tGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAcgg                                                                                              >  2:2705383/1‑151 (MQ=255)

                                               ggACGACGACTTATCACCCCGATCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGTAATAAACAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAcggc                                                                                             <  1:3189886/151‑1 (MQ=255)

                                               ggACGACGACTGATCACCCCGCTCATGAAAATTGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAcggc                                                                                             >  2:309112/1‑151 (MQ=255)

                                                 acgacgacTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAcggcgg                                                                                           <  1:8740847/151‑1 (MQ=255)

                                                      gacTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATCAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAggc                                                                                      <  2:9577856/151‑1 (MQ=255)

                                                         tGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTgacggcggacggcgtc                                                                                   >  2:8892751/1‑149 (MQ=255)

                                                             cACCCCGCACATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGTAATAAAGAAAAGGCCGCGAAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAg                                                                               <  1:2649263/151‑1 (MQ=255)

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

                                                                cccGCTCCTGAAAAATGGCCGCTTCGAGCGGCCTTTGTCATGGGCGGAATACAGACACGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACGTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCt                                                                            <  2:5812223/151‑1 (MQ=255)

                                                                       aTGACAAATGGCCGCTACGAGCGGCCTATTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTgcc                                                                     <  1:5073140/151‑3 (MQ=255)

                                                                           aaaaTGGCCGCTTCGAGCGGCCTTTTTCATGGGCGTAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCt                                                                 <  2:3944399/151‑1 (MQ=255)

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

                                                                               tGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTgcgc                                                             >  1:8444703/1‑151 (MQ=255)

                                                                                ggCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAACGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGcgcg                                                            <  1:9811433/151‑1 (MQ=255)

                                                                                  ccGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGACGTCTGCCTGTCCGCGGCGCGaa                                                          >  2:7515606/1‑151 (MQ=255)

                                                                                  ccGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGaa                                                          >  1:6383521/1‑151 (MQ=255)

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

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

                                                                                    gCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTACTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGc                                                        >  1:1273053/1‑151 (MQ=255)

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

                                                                                       tCGAGCGGCCTTTTTCATGGGCGGACTAAACACAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAg                                                     <  2:10027309/151‑1 (MQ=255)

                                                                                       tCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAg                                                     >  2:6625698/1‑151 (MQ=255)

                                                                                       tCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAg                                                     >  1:10437898/1‑151 (MQ=255)

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

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

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

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

                                                                                                cTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAgcggcg                                            >  1:10335540/1‑151 (MQ=255)

                                                                                                cTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAgcggcg                                            >  1:5870730/1‑151 (MQ=255)

                                                                                                     tCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACa                                       >  2:6156314/1‑151 (MQ=255)

                                                                                                      caggggCGGAATCCAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACAt                                      <  1:2901846/148‑1 (MQ=255)

                                                                                                        tGGGCGGAATAACGAAACGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATgg                                    <  1:4554761/151‑1 (MQ=255)

                                                                                                           gCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCt                                 >  2:2266583/1‑151 (MQ=255)

                                                                                                             ggAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTtc                               >  2:8420887/1‑151 (MQ=255)

                                                                                                             ggAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTtc                               >  2:5036828/1‑151 (MQ=255)

                                                                                                             ggAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTtc                               >  2:4594522/1‑151 (MQ=255)

                                                                                                              gAATAAAGACACGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTtct                              <  2:1754801/151‑1 (MQ=255)

                                                                                                              gAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTtct                              >  1:5984830/1‑151 (MQ=255)

                                                                                                              gAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTtct                              <  2:8615689/151‑1 (MQ=255)

                                                                                                               aaTCAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTtctc                             <  1:7956589/151‑1 (MQ=255)

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

                                                                                                                aTAAAGACAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcg                            <  2:1160205/151‑1 (MQ=255)

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

                                                                                                                 tAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTCTTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGGCTGCCTGTCCGCTGCGCGAAGCCAGGACGCGGCGGTACATGGCCTTCtcgt                           >  2:525739/1‑151 (MQ=255)

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

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

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

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

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

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

                                                                                                                         aGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGgag                   <  2:1378331/151‑1 (MQ=255)

                                                                                                                           gCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGACGCGGCGGTACATGGCCTTCTCGTCGTCGGAGAt                 >  2:9318593/1‑151 (MQ=255)

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

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

                                                                                                                              gCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGc              <  1:7254438/151‑1 (MQ=255)

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

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

                                                                                                                                gTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcct            <  2:10286609/151‑1 (MQ=255)

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

                                                                                                                                 tAAGCGGCCTTTTCTC‑TGAATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcctc           <  1:4238728/151‑1 (MQ=255)

                                                                                                                                   agcGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcctcct         >  2:4833913/3‑151 (MQ=255)

                                                                                                                                    gcGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCctcctc        >  2:9173694/2‑151 (MQ=255)

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

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

                                                                                                                                       gCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCACGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCgg     <  1:5637341/151‑1 (MQ=255)

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

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

                                                                                                                                          cttttCTCTGCATCTCGGTCTAAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGCt   <  1:9172345/150‑1 (MQ=255)

                                                                                                                                           ttttCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGCTc  <  1:4067599/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: `‑`