BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,388,209	A→T	8.0%	L16Q (CTG→CAG)	PP_1204 ←	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,388,209	0	A	T	8.0%	664.6 / 30.9	281	L16Q (CTG→CAG)	PP_1204	lipoprotein
Reads supporting (aligned to +/- strand): ref base A (104/128); new base T (8/17); total (123/158)
Fisher's exact test for biased strand distribution p-value = 2.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

Only 100 of 295 total aligned reads displayed.

GCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCTTGTACTCCA  >  NC_002947/1388060‑1388358

gcACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGcaca                                                                                                                                                       <  2:5760284/150‑1 (MQ=255)

gcACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGcaca                                                                                                                                                       <  1:3769493/150‑1 (MQ=255)

   cacGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTg                                                                                                                                                    >  2:5139551/1‑150 (MQ=255)

       tCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTgcggc                                                                                                                                                >  2:4989372/1‑150 (MQ=255)

         gCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCAAAGTGCGgcgc                                                                                                                                              >  1:4790997/1‑150 (MQ=255)

          cTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACTGTGCGGCGcc                                                                                                                                             <  2:7138836/150‑1 (MQ=255)

          cTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACTGTGCGGCGcc                                                                                                                                             <  2:2758500/150‑1 (MQ=255)

            gCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGcccc                                                                                                                                           <  1:5819362/150‑1 (MQ=255)

               gAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCaa                                                                                                                                        <  2:1337056/150‑1 (MQ=255)

               gAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCaa                                                                                                                                        >  1:6358763/1‑150 (MQ=255)

                aaGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCGGGGCGCCGCACTGTGCGGCGccctgctc                                                                                                                                       <  1:2615124/150‑6 (MQ=255)

                 aGCTCCATCACCTTGGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAgg                                                                                                                                      <  1:7151907/150‑1 (MQ=255)

                 aGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCCAgg                                                                                                                                      >  2:2653311/1‑150 (MQ=255)

                 aGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAgg                                                                                                                                      <  2:4639778/150‑1 (MQ=255)

                      cATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTaaa                                                                                                                                 <  2:6330143/150‑1 (MQ=255)

                       aTCACCTTCGGCTCGCCGACATCGTTATGCAGGGTATGCTCTAGCAGTTTGCGCTCTTCGCGCTACTCATGGTCGCTGCTATGTGACATCTGGTTGGCGCAGCCGCTGTCTTGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAg                                                                                                                                <  2:1115631/150‑1 (MQ=255)

                            cTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACt                                                                                                                           >  2:2960906/1‑150 (MQ=255)

                                ggCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGc                                                                                                                       <  2:6685294/150‑1 (MQ=255)

                                 gCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCTTGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCt                                                                                                                      <  2:1455598/150‑1 (MQ=255)

                                 gCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCCTGGTCGGTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCt                                                                                                                      <  1:4813542/150‑1 (MQ=255)

                                 gCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCCCACTGAGGCGCCCCCCAGGGTAAAGGTACTTAGCt                                                                                                                      >  1:1243397/1‑150 (MQ=255)

                                  cTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCtt                                                                                                                     >  1:613382/1‑150 (MQ=255)

                                   tCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGGGGCGCCCCCCCGGGGAACGGTACTTCGCtgg                                                                                                                    >  1:2048032/1‑148 (MQ=255)

                                   tCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCCAGGGTAAAGGTACTTCGCGTg                                                                                                                    >  1:4049839/1‑150 (MQ=255)

                                   tCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTg                                                                                                                    >  2:6402508/1‑150 (MQ=255)

                                       cGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACa                                                                                                                >  2:5406642/1‑150 (MQ=255)

                                        gaattCGATCTGCAGGGTTTGCTCGAGCAGTTGGCGCTCGACGCGCTCCTCATGGTCGCGGCGCTGGGACATCTGGTTGGCGCCGCCGGTGAGTAGCAGGGCGCCGCACCTGGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACAt                                                                                                               <  1:6454210/146‑1 (MQ=255)

                                        gACATCGATCTGCTGGTTATGCTCGAGCTGTTTGCGCGCGACGGGCGCCTCTTGGTCGCTGCGCTGGGACTTCTGGTTGGCGCAGCCGCTGACGAGCTGGGCGCCGAACAGTGGGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACAt                                                                                                               <  1:3452715/150‑1 (MQ=255)

                                        gACATCGATCTGCCGGGTCTGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACTTCTGGGTGGCGCAGCCGCTGCCGAGCTGGGCGCCGCACCGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACAt                                                                                                               <  2:653434/150‑1 (MQ=255)

                                          aaTCGATCTGCTGGGTATGCTCGAGCATTTTGCGCTCGACGGGCTCCTCTTGGTCGCTGGGGTGGGACATGTGGTTGGCGCAGCCGCTTACGAGCTGGGCGCCGCACTGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGa                                                                                                             <  1:4555137/149‑1 (MQ=255)

                                           aTCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGAt                                                                                                            <  1:3006357/150‑1 (MQ=255)

                                              gATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCAc                                                                                                         >  2:4308029/1‑150 (MQ=255)

                                                tcggCAGGGGATGCTCGAGCCGTGGGCCCGGGCCGCGCTCCTCAGGGTCGCTGCGCTGGGACATCTGGTGGGGGCTGCGGCTGAGCAGACGGGCGCCGCACATTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTc                                                                                                       <  1:3732665/147‑1 (MQ=255)

                                                  tGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGccc                                                                                                                                            >  2:219888/1‑111 (MQ=255)

                                                  tGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGt                                                                                                     <  2:5520187/150‑1 (MQ=255)

                                                  tGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGt                                                                                                     >  1:2278682/1‑150 (MQ=255)

                                                      gggTATGCTCGAGCAGTTTGCGCTCGCCGCGCTCCTCTTGGTCGCTGCGCGGGGACATCTGGTTGGCGCAGCCGCTGACGAGCTGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCAt                                                                                                 <  1:4451835/150‑1 (MQ=255)

                                                       ggTATGCTCGAGCAGGTGGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTGGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCAtt                                                                                                <  1:6617458/150‑1 (MQ=255)

                                                        gTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGGCGCACAGTGCGGCGCCCCCAAGGGTTAAGGTACTTAGGTTGTTGATGATGAATCGTGGATTg                                                                                               >  1:2349996/1‑150 (MQ=255)

                                                         tATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTgg                                                                                              <  2:4756383/150‑1 (MQ=255)

                                                              tCGCGCTGTTGGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCGGGTGGGCGCAGCCGCTGACGAGCAGGGCGCCGCACCGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATc                                                                                         <  1:1031707/150‑1 (MQ=255)

                                                              tCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACCGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACCTGATCACTCGTGCATTGGTTAta                                                                                         >  1:3613033/1‑149 (MQ=255)

                                                               cgtgctgTTTGCGCTCGCCGCGCTCCTCTTGGTGGCTGCGCGGGGACATGTGGTTGGCGCAGCCGCTGACGAGCTGGGCGCCGCACTGTGGGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCa                                                                                        <  1:151925/144‑1 (MQ=255)

                                                                gAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACATCGCTTGTACATGATCACTCGTGCATTGGTTATCAg                                                                                       >  1:4476521/1‑150 (MQ=255)

                                                                  gCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAgcg                                                                                     <  1:2449492/150‑1 (MQ=255)

                                                                  gCAGTTTGCGCTCGACGCGCTCCTCAGGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAgcg                                                                                     <  1:4239539/150‑1 (MQ=255)

                                                                   cAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAgcgg                                                                                    >  2:5951692/1‑150 (MQ=255)

                                                                   cAGTTTGCGCTCGACGCGCGCCTCATGGTCGCTGCGGTGGGACTTCTGGTTGGCGCAGCCGCTGACGAGCTGGGCGCCGCACCGGGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAgcgg                                                                                    <  1:21499/150‑1 (MQ=255)

                                                                      tttGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGa                                                                                 >  2:5428078/1‑150 (MQ=255)

                                                                      tttGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCGTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGa                                                                                 >  2:6999784/1‑150 (MQ=255)

                                                                       ttGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCagcagccag                                                                                <  2:6189208/150‑4 (MQ=255)

                                                                            cTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCCCAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGGGCATTGGTTATCAGCGGCGAACGCgg                                                                           >  2:6428485/1‑150 (MQ=255)

                                                                              cGACGCGCTCCTCATGGTCGCTGCGCTGGGACTTCTGGTTGGCGCAGCCGGTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGc                                                                         <  1:3703265/150‑1 (MQ=255)

                                                                              cGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGc                                                                         <  1:4466603/150‑1 (MQ=255)

                                                                               gACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGGGCGGCGCCCCCCAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGcc                                                                        >  2:5874522/1‑150 (MQ=255)

                                                                                 cgagcTCCTCATGGTCGCTGCGCTGGCACATCTGGTTTGCGCAGCCGCTGACGAGCAGGGCGCCGCACCGAGCGGCGCCCCCAAGGGTAAAGGTACTTCGCGTGTACATGATCACTCGTCCAGTGGTTATCAGCGGCGAAAGCGGGACtc                                                                      >  2:7344851/4‑149 (MQ=255)

                                                                                  gcgcTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCAGTGGTTATCAGCGGCGTACGCGGGCCTGc                                                                     >  2:303729/1‑150 (MQ=255)

                                                                                   cgcTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCa                                                                    <  2:4945467/150‑1 (MQ=255)

                                                                                   cgcTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCa                                                                    <  1:6767103/150‑1 (MQ=255)

                                                                                   cgcTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCa                                                                    <  2:613382/150‑1 (MQ=255)

                                                                                    gcTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACACTGCGGCGCCCCCCAGGGTAAAGGTACGTCGCGTGTACAGGATCACGCGTGCATTGGTTATAAGCGGCGAACGCGGGCCTGCAg                                                                   >  1:2230993/1‑150 (MQ=255)

                                                                                      tcctcTTGGGCGCTGCGCGGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGa                                                                 <  1:921529/150‑1 (MQ=255)

                                                                                           aTGGTCGCTGCGCTGGGACACCTGGTTGGCGCAGCCGCTGCCGAGCTGGGCGCCGCACCGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTc                                                            <  1:4662481/150‑1 (MQ=255)

                                                                                               tCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCa                                                        >  2:3685779/1‑150 (MQ=255)

                                                                                                     cgctgGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTAGCTTGTACATGATGACTCGTGGATTGGTTATCAGCGGCGAAAGCGGGCCTGCAGGAAGGTCAGCACTTCAt                                                  >  2:777472/1‑150 (MQ=255)

                                                                                                        tgGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGa                                               <  2:2670360/150‑1 (MQ=255)

                                                                                                        tgGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGa                                               <  2:2443399/150‑1 (MQ=255)

                                                                                                          ggACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGc                                             >  2:3855803/1‑150 (MQ=255)

                                                                                                            aCATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCgg                                           <  2:6833652/150‑1 (MQ=255)

                                                                                                               tCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACACGTGCATTTGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCaa                                        >  1:2929342/1‑150 (MQ=255)

                                                                                                                 tGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACg                                      >  2:3974999/1‑150 (MQ=255)

                                                                                                                 tGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACg                                      >  1:5554332/1‑150 (MQ=255)

                                                                                                                 tGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACGTCATTGAGCGGCAACg                                      >  2:1822960/1‑150 (MQ=255)

                                                                                                                      ggCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTgg                                 <  2:3697221/150‑1 (MQ=255)

                                                                                                                      ggCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTgg                                 <  1:5246863/150‑1 (MQ=255)

                                                                                                                      ggCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGATTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTgg                                 <  1:7123497/150‑1 (MQ=255)

                                                                                                                        cgcAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATAAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGc                               >  2:3181209/1‑150 (MQ=255)

                                                                                                                           aGCCGCTGACGTGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATTATCACTCTTGCATTGGTTATCAGCGGCGAACTCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTc                            <  1:1161566/150‑1 (MQ=255)

                                                                                                                           aGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTc                            <  1:7013818/150‑1 (MQ=255)

                                                                                                                           aGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTc                            <  1:4665776/150‑1 (MQ=255)

                                                                                                                           aGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTc                            <  2:3558651/150‑1 (MQ=255)

                                                                                                                             ccGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCt                          <  2:3459187/150‑1 (MQ=255)

                                                                                                                                    ccttgaGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTcggt                   >  2:2138567/6‑150 (MQ=255)

                                                                                                                                     gAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATAACACGTGCATTGGTTATCAGCGGCGAACCCGGGCCAGCAGGAAGGTCAGCACTTCATTGAGCGGCAACCCCAGGGCGTCCTGCTcggtg                  >  1:4956397/1‑150 (MQ=255)

                                                                                                                                      ttgggggCGCCGCACTGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTcggtgc                 >  1:1736256/5‑150 (MQ=255)

                                                                                                                                       gCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTcggtgcg                >  2:1696231/1‑150 (MQ=255)

                                                                                                                                           ggCGCCGCACTGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCggtgcggtgc            <  1:2653311/150‑1 (MQ=255)

                                                                                                                                           ggCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCAt                                                                                                 <  2:7434355/65‑1 (MQ=255)

                                                                                                                                           ggCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCAt                                                                                                 >  1:7434355/1‑65 (MQ=255)

                                                                                                                                           ggCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCggtgcggtgc            >  1:39042/1‑150 (MQ=255)

                                                                                                                                             cgccgcACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCtt          <  2:1981351/150‑1 (MQ=255)

                                                                                                                                             cgccgcACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCtt          <  2:4128930/150‑1 (MQ=255)

                                                                                                                                             cgccgaACTGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCtt          >  2:2507039/1‑150 (MQ=255)

                                                                                                                                              gccgcACTGTGCGGCGCCCTCAAGGGTAAAGGTACTTCGCTTGTACAGGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCAGTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCTTg         >  2:6360716/1‑150 (MQ=255)

                                                                                                                                              gccgcACCGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCTTg         <  1:1876486/150‑1 (MQ=255)

                                                                                                                                               ccgcACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCTTGt        >  1:3347565/1‑150 (MQ=255)

                                                                                                                                                  cacaGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCTTGTACt     >  1:248865/1‑150 (MQ=255)

                                                                                                                                                   acaGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCTTGTACTc    <  2:3376812/150‑1 (MQ=255)

                                                                                                                                                    caGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCTTGTACTcc   <  1:4000940/150‑1 (MQ=255)

                                                                                                                                                     aGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCTTGTACTCCa  <  1:1694885/150‑1 (MQ=255)

GCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCTTGTACTCCA  >  NC_002947/1388060‑1388358

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