BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,778,072	A→C	5.4%	K152Q (AAG→CAG)	sotB →	sugar efflux transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,778,072	0	A	C	5.4%	569.7 / 13.4	181	K152Q (AAG→CAG)	sotB	sugar efflux transporter
Reads supporting (aligned to +/- strand): ref base A (81/89); new base C (7/4); total (88/93)
Fisher's exact test for biased strand distribution p-value = 3.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 184 total aligned reads displayed.

GTTCATTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGATCAAA  >  NC_002947/2777924‑2778221

 ttCATTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAg                                                                                                                                                     <  1:2830878/150‑1 (MQ=255)

  tCATTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAgg                                                                                                                                                    <  2:2651292/150‑1 (MQ=255)

  tCATTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAgg                                                                                                                                                    >  1:1856222/1‑150 (MQ=255)

    aTTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGgcg                                                                                                                                                  >  1:2237868/1‑150 (MQ=255)

    aTTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGgcg                                                                                                                                                  >  2:2360504/1‑150 (MQ=255)

       gTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTg                                                                                                                                               >  1:4399529/1‑150 (MQ=255)

       gTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTg                                                                                                                                               <  2:4519051/150‑1 (MQ=255)

       gTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTg                                                                                                                                               <  1:3185131/150‑1 (MQ=255)

        tCAGGCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTgg                                                                                                                                              <  1:2602567/150‑1 (MQ=255)

        tCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCCGGCCAAGGCGCTgg                                                                                                                                              >  1:1690803/1‑150 (MQ=255)

        tCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTgg                                                                                                                                              <  2:2546391/150‑1 (MQ=255)

          aGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGc                                                                                                                                            >  2:3665571/1‑150 (MQ=255)

          aGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGc                                                                                                                                            >  1:2852690/1‑150 (MQ=255)

          aGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGc                                                                                                                                            >  1:4481380/1‑150 (MQ=255)

          aGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCACGCCAAGGCGCTGGGc                                                                                                                                            >  1:3869582/1‑150 (MQ=255)

             cACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCctg                                                                                                                                         <  1:2055425/150‑1 (MQ=255)

             cACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCctg                                                                                                                                         <  1:3648562/150‑1 (MQ=255)

             cACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCACGCCCAGGCGCTGGGccgg                                                                                                                                         >  1:100113/1‑148 (MQ=255)

              aCTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCctgc                                                                                                                                        >  1:2231646/1‑150 (MQ=255)

               cTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCctgct                                                                                                                                       <  2:1056816/150‑1 (MQ=255)

               cTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCctgct                                                                                                                                       <  2:2768538/150‑1 (MQ=255)

               cTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCctgct                                                                                                                                       <  2:4052158/150‑1 (MQ=255)

                ttGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCCAGCGCCTTGGCctgctg                                                                                                                                      <  2:354128/150‑1 (MQ=255)

                     gTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCAc                                                                                                                                 <  1:4036250/150‑1 (MQ=255)

                      tCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCAcc                                                                                                                                >  1:4275844/1‑150 (MQ=255)

                       ccTGGCTTTCGCATAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCTCTGGGCCTGCTGGCCACCg                                                                                                                               <  2:866079/150‑1 (MQ=255)

                         tGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGGCAAGGCGCTGGGCCTGCTGGCCACCGGt                                                                                                                             >  1:838159/1‑150 (MQ=255)

                         tGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCCCCGGt                                                                                                                             >  1:2379468/1‑150 (MQ=255)

                           gCTTTCGCAAAGCGTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGtac                                                                                                                           <  2:834992/150‑1 (MQ=255)

                                 gCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCt                                                                                                                     >  2:4432595/1‑150 (MQ=255)

                                  cAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTg                                                                                                                    >  2:3703550/1‑150 (MQ=255)

                                  cAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTg                                                                                                                    <  2:1493732/150‑1 (MQ=255)

                                   aaaGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTgg                                                                                                                   <  2:871392/150‑1 (MQ=255)

                                   aaaGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTgg                                                                                                                   >  2:1825922/1‑150 (MQ=255)

                                    aaGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGc                                                                                                                  <  1:4510179/150‑1 (MQ=255)

                                         tCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCTGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATgg                                                                                                             <  1:2394400/150‑1 (MQ=255)

                                         tCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATgg                                                                                                             <  1:2986588/150‑1 (MQ=255)

                                            ccATGCTGCTGGTCAGCCGCATCGGCAGTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCGGCTGGCCACCGGTACTACCCTGGCCATGGTcc                                                                                                          >  2:4208353/1‑150 (MQ=255)

                                             cATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCAGGGTCCt                                                                                                         <  2:3179830/150‑1 (MQ=255)

                                               tgctgcTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTgg                                                                                                       <  1:182530/150‑1 (MQ=255)

                                               tgctgcTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTgg                                                                                                       <  1:2719106/150‑1 (MQ=255)

                                                gctgctGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTggg                                                                                                      <  2:4227720/150‑1 (MQ=255)

                                                  tgctgGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCa                                                                                                    >  2:2455377/1‑150 (MQ=255)

                                                  tgctgGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCGGGGCCTGCTGGCCACCGGTACTACCCTGGCCCTGGTCCTGGGCa                                                                                                    >  2:1579317/1‑150 (MQ=255)

                                                    ctgGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCCGGCCAAGGCGCAGGGCCTGCTGGCCACCCGTACTACCCAGGCCAAGGTCCTGGGCCTc                                                                                                  >  1:2670701/1‑150 (MQ=255)

                                                    ctgGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATc                                                                                                  <  1:4292794/150‑1 (MQ=255)

                                                      ggTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAACGCGCTGGGCCTGCTGGCCACCGGTACTACCCAGGCCATGGTCCTGGGCCTccc                                                                                                >  2:956979/1‑150 (MQ=255)

                                                         cAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCt                                                                                             <  2:3807627/150‑1 (MQ=255)

                                                          aGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTg                                                                                            <  1:1896559/150‑1 (MQ=255)

                                                            ccGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTggg                                                                                          <  1:1642762/150‑1 (MQ=255)

                                                                aTCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCAGGCCATGGTCCAGGGCATCCCGCTGGGCCgg                                                                                      >  1:4303862/1‑149 (MQ=255)

                                                                  cGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCgtgt                                                                                    <  1:1279766/150‑1 (MQ=255)

                                                                  cGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCgtgt                                                                                    <  2:2111760/150‑1 (MQ=255)

                                                                   ggCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCACCAGGCCCACGCGCGGGGCCTGGGGGCCACCGGTACTACCCTGGCCATGGGCCTGGGCATCCCCCGGGGCCGtgtg                                                                                   >  1:3771349/1‑150 (MQ=255)

                                                                    gCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCGGGGCATCCCGCTGGGCCGTgggg                                                                                  >  2:1211080/1‑150 (MQ=255)

                                                                     cATTGCCCTGGCGCATGCGGTGTTCTGGTCCCTCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTgtggt                                                                                 <  1:3398930/150‑1 (MQ=255)

                                                                           ccTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGa                                                                           <  1:2198015/150‑1 (MQ=255)

                                                                             tGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAgg                                                                         <  2:2441721/150‑1 (MQ=255)

                                                                                 gcATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACt                                                                     >  2:591703/1‑150 (MQ=255)

                                                                                  cATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTg                                                                    <  2:3309364/150‑1 (MQ=255)

                                                                                  cATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTg                                                                    >  2:3039296/1‑150 (MQ=255)

                                                                                         tgtTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTgggcgggc                                                             >  1:697454/1‑150 (MQ=255)

                                                                                         tgtTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGc                                                             <  2:4276099/150‑1 (MQ=255)

                                                                                         cgtgcgggTCCAGCCCCGCCCCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCAGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGc                                                             <  1:254105/144‑1 (MQ=255)

                                                                                              tGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGGCCCCGGTACTACCCTGGCCAAGGGCCTGGGCATCCCGCGGGGCCGTGTGGGGGGCGGGGGCCTGGGCTGGCGTATc                                                        >  1:3901843/1‑150 (MQ=255)

                                                                                              tGGTCCATCACCGCCTCGCTGGCGGTGCGGGTGGCGCCGCCGCGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATc                                                        <  1:3048758/150‑1 (MQ=255)

                                                                                               ggTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCCCCCGGCCCAGGCGCTGGGGCGGCTGGCCCCCGGTACTACCCTGGCCCAGGGCCGGGGCAACCCGCGGGGCCGGGTGGTGGGCGAGGCACAGGGCGGGCGGATCa                                                       >  1:4237077/1‑150 (MQ=255)

                                                                                               ggTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCa                                                       <  1:2273487/150‑1 (MQ=255)

                                                                                                  ccATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGt                                                    >  2:2745564/1‑150 (MQ=255)

                                                                                                   cATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGtt                                                   <  2:3109352/150‑1 (MQ=255)

                                                                                                   cATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCACGCCAAGGCGCTGGGCCGGCTTGCCCCCCGTACTACCCTGGCCATGGTCCAGGGCCTACCCCGGGGCCGGGGGGTGGGCGAGGCACTGGGCTGGCGGATAACGtt                                                   >  1:3786497/1‑150 (MQ=255)

                                                                                                   acgcACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCACCCCGCCAAGGCGCTGGGCCGGCTGGCCCCCCGTACTACCCGGGCCCTGGTCCGGGGCCTCCCGCGGGGCCGGGTGGTGGGCGAGGGACTGGGCTGGCGTCTCACGtt                                                   >  1:101164/4‑150 (MQ=255)

                                                                                                    aTCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTc                                                  <  2:3867000/150‑1 (MQ=255)

                                                                                                     tCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGTGGCACTGGGCTGGCGTATCACGTTcc                                                 >  2:3255852/1‑150 (MQ=255)

                                                                                                     tCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGttac                                                 >  2:923051/1‑148 (MQ=255)

                                                                                                        ccGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGa                                              <  1:1161869/150‑1 (MQ=255)

                                                                                                        ccGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGa                                              <  2:3052128/150‑1 (MQ=255)

                                                                                                        ccGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGa                                              <  2:168011/150‑1 (MQ=255)

                                                                                                          gCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCGGGGCCGGCTGGCCACCGGTACTACCCAGGCCCTGGTCCTGGGCATCCCGCTGGGGCGGGGGGTGGGCGAGGCCCTGGGCTGGCGTATCACGGTCCAGAGc                                            >  1:562040/1‑150 (MQ=255)

                                                                                                             tCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGGCCTGGGCATCCCGCTGGGCCGTGTGGGGGGCGTGGCACTGGGCTGGCGGATCACGTTCCTGAGCAtg                                         >  1:4111199/1‑149 (MQ=255)

                                                                                                                     ggTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCGGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGt                                 >  1:1446646/1‑150 (MQ=255)

                                                                                                                        gCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAgc                              <  2:1360597/150‑1 (MQ=255)

                                                                                                                         cGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCAGGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAgcg                             <  1:1153971/150‑1 (MQ=255)

                                                                                                                            gtgGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTa                          <  1:1706712/150‑1 (MQ=255)

                                                                                                                            gtgGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTa                          <  2:1925887/150‑1 (MQ=255)

                                                                                                                              ggCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGc                        <  2:4387243/150‑1 (MQ=255)

                                                                                                                               gcgcCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCGGCTGGCCAACCGTACTACCCTGGCCATGGTCCTGGGCCTCCCGCAGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATAACGTTCCTGAGCATTGCCGGGGTCGCGCTAGcc                       >  2:1740629/1‑150 (MQ=255)

                                                                                                                                    gccgGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCAt                  <  1:3703550/150‑1 (MQ=255)

                                                                                                                                    gccgGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCAt                  <  2:280254/150‑1 (MQ=255)

                                                                                                                                         gCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCttt             <  1:964609/150‑1 (MQ=255)

                                                                                                                                          ccAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTCGCGCTAGCCACCATGCTTTg            >  1:3895699/1‑150 (MQ=255)

                                                                                                                                           cagcagGCCCAGCGCCTTGGCCTGCTGGCCCCCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGc           >  2:587779/1‑150 (MQ=255)

                                                                                                                                           cagcagGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGc           <  1:2878145/150‑1 (MQ=255)

                                                                                                                                            agcagGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTCGCCACCATGCTTTGcc          >  2:314163/1‑150 (MQ=255)

                                                                                                                                            agcagGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGcc          >  2:3684020/1‑150 (MQ=255)

                                                                                                                                            agcagGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGcc          >  1:843724/1‑150 (MQ=255)

                                                                                                                                             gcagGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCt         <  2:4112217/150‑1 (MQ=255)

                                                                                                                                              cagGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTg        >  1:3532626/1‑150 (MQ=255)

                                                                                                                                                    aaGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGATCaaa  <  1:1945638/150‑1 (MQ=255)

                                                                                                                                                    aaGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGATCaaa  >  2:1314630/1‑150 (MQ=255)

GTTCATTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGATCAAA  >  NC_002947/2777924‑2778221

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