BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_008463	3,673,534	0	T	C	99.7%	1413.4 / ‑4.1	361	intergenic (+83/+15)	fabI/PA14_RS16690	enoyl‑ACP reductase FabI/SurA N‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (206/154); total (207/154)
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.91e-01
*	NC_008463	3,673,535	0	C	T	97.6%	1531.1 / 12.9	369	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/2); new base T (206/154); total (213/156)
Fisher's exact test for biased strand distribution p-value = 3.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01
*	NC_008463	3,673,536	0	T	C	97.6%	1408.4 / 19.6	369	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/2); new base C (206/154); total (213/156)
Fisher's exact test for biased strand distribution p-value = 3.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

Only 100 of 369 total aligned reads displayed.

TCCTCTACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGNGTAAACGGCCTTTT‑TCT‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGCT  >  NC_008463/3673393‑3673679

|||

tCCTCTACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCAtc                                                                                                                                           <  2:456458/150‑1 (MQ=255)

tCCTCTACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCAtc                                                                                                                                           <  2:10603579/150‑1 (MQ=255)

  ctctCCGTCGATGGCGGATACAACACTACCGCCATGGGCCCGCTGGACTACTACTGAACACCCCGAACATTAAAATTGGCCGCTCCTAGCGGCCTTTTTCATGGGCGTAATAACGAAAAGCCCGCGTAAGCGGCCTTTT‑CTC‑TGCAtctc                                                                                                                                         <  1:5554928/150‑1 (MQ=255)

     tACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGt                                                                                                                                      >  2:1352614/1‑150 (MQ=255)

     tACGTCGAGGGCGGCTTCAACACAACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAACGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGt                                                                                                                                      <  1:8413964/150‑1 (MQ=255)

        gTCGAGGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGAACACCCCGCTCATGAAAATGGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAACTAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGtctc                                                                                                                                   <  1:3333304/150‑1 (MQ=255)

          cGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAg                                                                                                                                 <  1:4285888/150‑1 (MQ=255)

              ggcggcTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTAtt                                                                                                                             >  1:3886753/1‑150 (MQ=255)

               gcggcTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAATTGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTAttt                                                                                                                            <  1:4348973/150‑1 (MQ=255)

                  gcTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGACTAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTtctc                                                                                                                         >  1:220135/1‑150 (MQ=255)

                  gcTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTCTTtctc                                                                                                                         >  2:7750813/1‑150 (MQ=255)

                      cAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACt                                                                                                                     >  2:4243911/1‑150 (MQ=255)

                        acacTCCCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATTAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATACAGAACAGCCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACtac                                                                                                                   <  1:1613720/150‑3 (MQ=255)

                         cacTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTcc                                                                                                                  <  2:5700248/150‑1 (MQ=255)

                          accaCCGCCAGGGGCCCGCTGGACGACGACTGGTCACCCCGCACATTAAAAATGGCCGCTTCGAGCGGCCTTTTTCAGGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCg                                                                                                                 <  1:9180591/147‑1 (MQ=255)

                            tACCGCCTTGGGCCCGCTGGACGACGACTGAACACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAACGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGcc                                                                                                               <  1:967097/150‑1 (MQ=255)

                            tACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGcc                                                                                                               >  1:6315812/1‑150 (MQ=255)

                              ccGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCtt                                                                                                             >  1:7756733/1‑150 (MQ=255)

                               cGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCTTGGGCGGAAAAAACAAACGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTg                                                                                                            <  1:4539970/150‑1 (MQ=255)

                                gCCATGGGCCCGCTGGACGACGACTGATCACCCCGCCCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAACAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGt                                                                                                           <  1:4388719/150‑1 (MQ=255)

                                 ccATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAACGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTc                                                                                                          <  2:7896296/150‑1 (MQ=255)

                                 ccATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTc                                                                                                          >  2:926727/1‑150 (MQ=255)

                                   aTGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGGATTTCTCGACTTCCGCCTTGTCCt                                                                                                        >  1:5870165/1‑150 (MQ=255)

                                    tGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTg                                                                                                       >  2:9797222/1‑150 (MQ=255)

                                      ggCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCAGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAATGGCCGCGGAAGCGGCCTGTT‑CTC‑TGCATCTCGGTCTCAGTAGTTCTCGTCTTCCGCCTGGTCCTGCa                                                                                                     >  2:9138792/1‑150 (MQ=255)

                                          cGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTg                                                                                                 >  2:5309624/1‑150 (MQ=255)

                                           gCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGa                                                                                                >  1:2108927/1‑150 (MQ=255)

                                           gCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCATCGAGCGGCCTTTTTCTTGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGa                                                                                                <  1:9838470/150‑1 (MQ=255)

                                           gCTGGACGACGACTGATCACCCCGCTCATGAAAAAGGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAACGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGa                                                                                                <  1:9212528/150‑1 (MQ=255)

                                            cTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAc                                                                                               >  1:1267071/1‑150 (MQ=255)

                                             tGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAcg                                                                                              <  2:4708011/150‑1 (MQ=255)

                                             tGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAACAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCACCTGAcg                                                                                              <  1:10078690/150‑1 (MQ=255)

                                              ggACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGTCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAcgg                                                                                             >  2:7298944/1‑150 (MQ=255)

                                              ggACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCACCTGAcgg                                                                                             >  2:10258141/1‑150 (MQ=255)

                                                acgacgacTGATCACCCCGCTCATGAAAATTGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCGCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGAcggcg                                                                                           >  2:5930396/1‑150 (MQ=255)

                                                  gacgacTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGa                                                                                         <  1:10493622/150‑1 (MQ=255)

                                                   acgacTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGaa                                                                                        >  1:2662065/1‑150 (MQ=255)

                                                      acTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAATGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAggc                                                                                     >  2:239/1‑150 (MQ=255)

                                                      acTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAggc                                                                                     >  2:8040685/1‑150 (MQ=255)

                                                      acTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAggc                                                                                     >  1:8738363/1‑150 (MQ=255)

                                                         gATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAggcggc                                                                                  >  2:8348281/1‑150 (MQ=255)

                                                           tCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGACAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGa                                                                                <  2:7500733/150‑1 (MQ=255)

                                                                ccGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATATAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGGCGTGCACCTGTCGGCGGAAGGCGGCGCAGGCt                                                                           >  1:7424360/1‑150 (MQ=255)

                                                                  gCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGACGTCTGc                                                                         >  2:10655303/1‑150 (MQ=255)

                                                                  gCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGc                                                                         <  2:177814/150‑1 (MQ=255)

                                                                       tGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTc                                                                    >  1:10674867/1‑150 (MQ=255)

                                                                        gAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTcc                                                                   >  2:860687/1‑150 (MQ=255)

                                                                           aaaTGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGACAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCt                                                                <  1:9388017/150‑1 (MQ=255)

                                                                           aaaTGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCt                                                                >  2:581569/1‑150 (MQ=255)

                                                                            aaTGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATACAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTg                                                               <  2:307877/150‑1 (MQ=255)

                                                                                gCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGcgcg                                                           >  2:9504431/1‑150 (MQ=255)

                                                                                gCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGcgcg                                                           >  1:7079034/1‑150 (MQ=255)

                                                                                gCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGcgcg                                                           <  2:2211374/150‑1 (MQ=255)

                                                                                     ttCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGcc                                                      >  1:10318954/1‑150 (MQ=255)

                                                                                     ttCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGcc                                                      >  1:5407151/1‑150 (MQ=255)

                                                                                        gAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTNTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAgg                                                   >  1:5667916/1‑150 (MQ=255)

                                                                                         aGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGa                                                  >  2:6438532/1‑150 (MQ=255)

                                                                                         aGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGa                                                  <  2:368827/150‑1 (MQ=255)

                                                                                             gCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAgcg                                              >  2:3965522/1‑150 (MQ=255)

                                                                                               cTTTTTCATGGGCGGAATAAAGAAATGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCGCAGTAGTTCTCGACTTCCGCCTTGTCCTGCCACTGACGGCGGAAGGCGGCGAAGTCTGCCGGTCCGCTGCGCGAAGCCCGGAAgcggc                                            >  2:7026787/1‑150 (MQ=255)

                                                                                               cTTTTTCATGGGCGGAATAAAGAAACGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAgcggc                                            <  2:6033895/150‑1 (MQ=255)

                                                                                                 ttttCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGcggcgg                                          >  1:3183896/1‑150 (MQ=255)

                                                                                                   ttCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTa                                        <  1:8911537/150‑1 (MQ=255)

                                                                                                    tCATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGACGCGGCGGTAc                                       >  1:6977859/1‑150 (MQ=255)

                                                                                                     cATGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACa                                      >  2:1925815/1‑150 (MQ=255)

                                                                                                      aTGGGCGGAATCAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACAt                                     <  2:1168111/150‑1 (MQ=255)

                                                                                                       tGGGCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATg                                    <  2:6947773/150‑1 (MQ=255)

                                                                                                        gggCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATgg                                   >  1:8157563/1‑150 (MQ=255)

                                                                                                         ggCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGc                                  >  2:3529900/1‑150 (MQ=255)

                                                                                                          gCGGAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGcc                                 <  1:1741194/150‑1 (MQ=255)

                                                                                                             gAATAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTtc                              >  1:4973335/1‑150 (MQ=255)

                                                                                                              aaTACACAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCACGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTtct                             <  1:2514620/150‑1 (MQ=255)

                                                                                                              aaTAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTtct                             <  1:3807659/150‑1 (MQ=255)

                                                                                                               aTAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTtctc                            >  1:5773554/1‑150 (MQ=255)

                                                                                                               aTAAAGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTtctc                            >  1:2850599/1‑150 (MQ=255)

                                                                                                               aTAAACAAACGGCCGCGTACGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTtctc                            <  1:1076912/150‑1 (MQ=255)

                                                                                                                  aaGAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtc                         >  1:6626990/1‑150 (MQ=255)

                                                                                                                    gAAAAGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtcgt                       >  1:6107604/1‑150 (MQ=255)

                                                                                                                     aaaaGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtcgtc                      <  2:8815459/150‑1 (MQ=255)

                                                                                                                     aaaaGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtcgtc                      >  1:902120/1‑150 (MQ=255)

                                                                                                                     aaaaGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCACGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCtcgtcgtc                      <  2:5651511/150‑1 (MQ=255)

                                                                                                                       aaGGCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCgg                    >  1:5104832/1‑150 (MQ=255)

                                                                                                                         ggCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGgag                  <  2:955800/150‑1 (MQ=255)

                                                                                                                          gCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGgagg                 >  1:6497612/1‑149 (MQ=255)

                                                                                                                          gCCGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGgaga                 >  1:523384/1‑150 (MQ=255)

                                                                                                                           ccGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGAt                >  2:8985946/1‑150 (MQ=255)

                                                                                                                           ccGCGTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTATCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGAt                >  2:1590968/1‑150 (MQ=255)

                                                                                                                              cgTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGc             >  1:4085236/2‑150 (MQ=255)

                                                                                                                               gTAAGCGGCCTTTTCTCT‑T‑CATCTCGGTCTCAGTATTTAACGACTTCCGCCTTGACCTGCAACTGACTGCGGAAGTCCGCGAAGTCTGCCAGACCGCTGCGCGAAGCCAGGAAGAGGCGGTACATGTCCTTCTCGTCGTCGGAGATAGcc            >  1:1887376/1‑150 (MQ=255)

                                                                                                                               gTAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcc            >  1:3048391/1‑150 (MQ=255)

                                                                                                                                tAAGCGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcct           >  1:3418398/1‑150 (MQ=255)

                                                                                                                                  agcGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcctcc         >  2:3171937/3‑150 (MQ=255)

                                                                                                                                  agcGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcctcc         >  2:2873729/3‑150 (MQ=255)

                                                                                                                                   gcGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcctcct        >  1:7649376/2‑150 (MQ=255)

                                                                                                                                   gcGGCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATGACTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCGGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGcctcct        <  1:4605221/149‑1 (MQ=255)

                                                                                                                                     ggCCTTTT‑CTC‑TGCAGCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCTCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTAGTCGGAGATAGCAtcctcc      <  2:7153839/150‑1 (MQ=255)

                                                                                                                                      gCCTTTT‑CTC‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCg     <  2:2528970/150‑1 (MQ=255)

                                                                                                                                         cttt‑tCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGc   >  2:5141977/2‑150 (MQ=255)

                                                                                                                                         cttt‑tCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGc   <  2:3759272/149‑1 (MQ=255)

                                                                                                                                         cttt‑tCTCTGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGc   >  2:9681629/2‑150 (MQ=255)

|||

TCCTCTACGTCGATGGCGGCTTCAACACTACCGCCATGGGCCCGCTGGACGACGACTGATCACCCCGCTCATGAAAAATGGCCGCTTCGAGCGGCCTTTTTCATGGGCGGAATAAAGAAAAGGCCGNGTAAACGGCCTTTT‑TCT‑TGCATCTCGGTCTCAGTATTTCTCGACTTCCGCCTTGTCCTGCAACTGACGGCGGAAGGCGGCGAAGTCTGCCTGTCCGCTGCGCGAAGCCAGGAAGCGGCGGTACATGGCCTTCTCGTCGTCGGAGATAGCCTCCTCCGGCT  >  NC_008463/3673393‑3673679

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