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%	1409.5 / NA	359	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 (221/138); total (221/138)
*	NC_008463	3,673,535	0	C	T	95.5%	1551.4 / 31.7	375	intergenic (+84/+14)	fabI/PA14_RS16690	enoyl‑ACP reductase FabI/SurA N‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand): ref base C (11/6); new base T (221/137); total (232/143)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01
*	NC_008463	3,673,536	0	T	C	95.7%	1369.9 / 42.0	375	intergenic (+85/+13)	fabI/PA14_RS16690	enoyl‑ACP reductase FabI/SurA N‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand): ref base T (11/5); new base C (221/138); total (232/143)
Fisher's exact test for biased strand distribution p-value = 6.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 375 total aligned reads displayed.

ATCCTCTACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGNGTAAACGGCCTTTTTCT‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGCTC  >  NC_008463/3673392‑3673680

|||

      tACGTCGATGGCGGCTTCAACACTTCCGCAATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCATTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGtc                                                                                                                                      <  2:5225546/151‑1 (MQ=255)

      tACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCAGGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGtc                                                                                                                                      <  1:1735987/151‑1 (MQ=255)

       aCGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGtct                                                                                                                                     >  1:4783994/1‑151 (MQ=255)

         gTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGCTCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCa                                                                                                                                   <  2:10338705/151‑1 (MQ=255)

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

           cGAGGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGt                                                                                                                                 <  2:9017933/151‑1 (MQ=255)

            gATGGCGGCTTCAACACTACCGCCTTGGGCCCGCTGGACGACGACTGATCACCCCGCTCATTAAAAATGGCCGCTTCGAGCGGCATTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTa                                                                                                                                <  1:5115802/151‑1 (MQ=255)

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

             agggcggcTACAACACTCCCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTAt                                                                                                                               <  1:9506042/149‑1 (MQ=255)

             aTGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAATTGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCTTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTAt                                                                                                                               <  2:3050289/151‑1 (MQ=255)

               ggcggcTACAACACTACCGCCTGGGGCCCGCGGGACGACTACTTAACACCCCGCTCATGAAAATGGGCCGCTTCGCGCGGCATTTTTCAGGGGCGTAATAACGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTAttt                                                                                                                             <  1:7175378/151‑1 (MQ=255)

                gcggcTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAATTGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTtc                                                                                                                            >  1:4821853/1‑151 (MQ=255)

                 cggcTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTAATtct                                                                                                                           >  1:2606256/1‑151 (MQ=255)

                 cggcTTCAACACGACCGCCATGGGCCCGCGGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAAGAAAGTCAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTtct                                                                                                                           <  2:10635206/151‑1 (MQ=255)

                  ggcTTCACCACTACCCCCAGGGGCCCGCTGCACGACGACGGATCACCCCGCACATGAAAATGGGCCGATTCGCGCGGCCTTTTTCATGGGCGTAATAAAGAAACGCCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTtctc                                                                                                                          <  1:3574631/151‑1 (MQ=255)

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

                   gcTTCAACACTACCGCCTTGGGCCCGCTGGACGACGACGGATCACCCCGCTCATGAAAAATGGCCGCTTCTAGCGGCCTTTTTCATGGGCGTAATAAACAAAAGGCCGCGTATGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCg                                                                                                                         <  2:9789932/151‑1 (MQ=255)

                   gcTTCAACACTACCGCCATGGGCCCGCGGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGTGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCg                                                                                                                         <  1:2491370/151‑1 (MQ=255)

                    cTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTTAACACCCCGCACATTAAAAATGGCCGCTTCGAGCGGCATTTTTCTTGGGCGGAATAAACAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGa                                                                                                                        <  1:1009954/151‑1 (MQ=255)

                    cTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGa                                                                                                                        >  1:1925289/1‑151 (MQ=255)

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

                       cAACACTACCGCCATGGGCCCGCTGGCCGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACtt                                                                                                                     <  1:3992470/151‑1 (MQ=255)

                          cacTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCg                                                                                                                  >  2:1350033/1‑151 (MQ=255)

                          cacTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCACATGAAAAATGGCAGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCg                                                                                                                  <  2:5339990/151‑1 (MQ=255)

                           acTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGGATTTCTCGACTTCCGc                                                                                                                 >  2:10573986/1‑151 (MQ=255)

                            cTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGcc                                                                                                                <  1:7706403/151‑1 (MQ=255)

                              aCCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAACGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCtt                                                                                                              <  2:3699154/151‑1 (MQ=255)

                               cccccATGGGCCCGCTGGACGACGACTGAACACCCCGATCATGAAAAATGGCCGCTACGAGCGGCATTTTTAATGGGCGTAATACCCAAAAGGCCGCGTATGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTg                                                                                                             <  1:64136/148‑1 (MQ=255)

                                  ccATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTcc                                                                                                          >  2:4844930/1‑151 (MQ=255)

                                         cccGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACt                                                                                                   >  1:5180526/1‑151 (MQ=255)

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

                                          ccGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCAGGGGCGGAATAAAGCAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTg                                                                                                  <  1:4138682/151‑1 (MQ=255)

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

                                               ggACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAcggc                                                                                             <  1:1416505/151‑1 (MQ=255)

                                                   gacgacTGATCACCCCGCTCATGAAAAAGGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAACGAACAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGaa                                                                                         <  1:3180787/151‑1 (MQ=255)

                                                    acgacTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAg                                                                                        >  1:9161986/1‑151 (MQ=255)

                                                    acgacTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAg                                                                                        >  2:6778324/1‑151 (MQ=255)

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

                                                          gATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGgcggcg                                                                                  >  2:1184460/1‑151 (MQ=255)

                                                            tCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGaa                                                                                >  2:338987/1‑151 (MQ=255)

                                                                   gCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGcc                                                                         >  1:7483820/1‑151 (MQ=255)

                                                                        tcaaaaaTGGCCGCTTCGAGCGGCCTTTTACAGGGGCGTAATAAAGACAAGGCCGCGTCAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTcc                                                                    <  1:8061041/149‑1 (MQ=255)

                                                                         gAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCg                                                                   >  2:7819864/1‑151 (MQ=255)

                                                                           ataaTGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGATTAAAGAAATGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCt                                                                 <  1:6953284/149‑1 (MQ=255)

                                                                             aaTGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGTTCTCAGTATTTCTCGACTTCCGCCTAGTCCTGCAACTGACGGCGGAAGGCGGCGACGTCTGCCTGTCCGCTgc                                                               >  2:5769366/1‑151 (MQ=255)

                                                                             aaTGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTgc                                                               >  1:8953606/1‑151 (MQ=255)

                                                                              aTGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGACTAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTgcg                                                              >  2:5695441/1‑151 (MQ=255)

                                                                                ggCCGCGTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAACGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCACAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGcgcg                                                            <  2:6558212/151‑1 (MQ=255)

                                                                                 gCCGCTTCGCGCGGCCATTTTCATGGGCGGAATAAAGAACAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGa                                                           <  1:1457539/151‑1 (MQ=255)

                                                                                   cGCTTCGAGCGGCATTTTTCATGGGCGGAATAAATAAACGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTAGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAg                                                         <  1:4072746/151‑1 (MQ=255)

                                                                                    gCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTGGGTCTCAGTATTTCTCGACTTCCGCCGTGTCCTGCAACTGACGGCGGAAGGCGGCGTCGTCTGCGTGTCCGCTGCGCGAAGc                                                        >  2:4344036/1‑151 (MQ=255)

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

                                                                                     cTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGcc                                                       >  1:963330/1‑151 (MQ=255)

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

                                                                                        cGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAgg                                                    <  2:3256506/151‑1 (MQ=255)

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

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

                                                                                           gCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAg                                                 <  2:7742334/151‑1 (MQ=255)

                                                                                           gCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAg                                                 >  1:8453625/1‑151 (MQ=255)

                                                                                           gCGGCCTTTTTCATGGGCGGAAGCACGAAACGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCAGCAACTGTCGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAg                                                 <  2:7426261/151‑1 (MQ=255)

                                                                                            cGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGACGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAgc                                                >  2:2603473/1‑151 (MQ=255)

                                                                                               ccTTTGTCATGGGCGGAATAATGAACAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAgcggc                                             <  2:9865038/151‑1 (MQ=255)

                                                                                                cTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCAGCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAgcggcg                                            >  2:166026/1‑151 (MQ=255)

                                                                                                 tttttCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGcggcgg                                           >  2:8066840/1‑151 (MQ=255)

                                                                                                 tttttCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGcggcgg                                           >  2:6194287/1‑151 (MQ=255)

                                                                                                   tttCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTa                                         >  1:5377420/1‑151 (MQ=255)

                                                                                                   tttCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTa                                         >  2:4626659/1‑151 (MQ=255)

                                                                                                    ttCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTAc                                        <  2:2234184/151‑1 (MQ=255)

                                                                                                    ttCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTAc                                        <  1:9752533/151‑1 (MQ=255)

                                                                                                     tcttGGGCGTAATAAAGAAACGGCCGCTTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACa                                       <  1:9011090/148‑1 (MQ=255)

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

                                                                                                        tGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATgg                                    <  2:8781058/151‑1 (MQ=255)

                                                                                                         gggCGGAATAAAGACAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGc                                   <  2:10136333/151‑1 (MQ=255)

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

                                                                                                            cGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCtt                                >  1:1054323/1‑151 (MQ=255)

                                                                                                            cGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCtt                                <  1:5033751/151‑1 (MQ=255)

                                                                                                             ggAATAAAGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGACGCCAGGAAGCGGCGGTACATGGCCTtc                               >  2:199658/1‑151 (MQ=255)

                                                                                                               aaTAAAGCACAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTtctc                             <  2:934642/151‑1 (MQ=255)

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

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

                                                                                                                  aaaGAACAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtc                          <  2:2783813/151‑1 (MQ=255)

                                                                                                                   aaGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtcg                         <  1:8475179/151‑1 (MQ=255)

                                                                                                                   aaGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCGTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCGGCGCGAAGCGCGGCTGCGGCGGTACAGGGCCTTCtcgtcg                         >  2:10030464/1‑151 (MQ=255)

                                                                                                                    aGAAAAGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtcgt                        >  1:3372593/1‑151 (MQ=255)

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

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

                                                                                                                        aaGGCCGCGTAAGCGGCCTTTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGga                    >  1:2348404/1‑151 (MQ=255)

                                                                                                                          ggCCGCGTACGCGGCCTGTTCTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGGCCTGCAACTGACGGCGGAAGGCGGCGCAGTCTGCCTGTCCGATGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGgaga                  >  2:6419001/1‑151 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                                                                                                           ttttCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGCTc  >  1:5962422/1‑151 (MQ=255)

                                                                                                                                           ttttCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGCTc  >  1:1541653/1‑151 (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: `‑`