BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,778,072	A→C	6.9%	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	6.9%	463.7 / 14.2	159	K152Q (AAG→CAG)	sotB	sugar efflux transporter
Reads supporting (aligned to +/- strand): ref base A (81/66); new base C (9/3); total (90/69)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 161 total aligned reads displayed.

TTCATTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGATCAA  >  NC_002947/2777925‑2778220

ttCATTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAg                                                                                                                                                     >  2:3473863/1‑150 (MQ=255)

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

ttCATTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAg                                                                                                                                                     >  1:2979051/1‑150 (MQ=255)

    ttGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGgcgc                                                                                                                                                 >  1:675725/1‑150 (MQ=255)

      gTCAGCCACTTGCTGTCCTGGCTTTCGCATAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCtt                                                                                                                                               <  1:263486/150‑2 (MQ=255)

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

       tCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTgg                                                                                                                                              >  2:3485491/1‑150 (MQ=255)

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

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

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

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

              cTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCCAGCGCCTTGGCctgct                                                                                                                                       <  2:757099/150‑1 (MQ=255)

              aatGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCctgct                                                                                                                                       >  2:1834736/3‑150 (MQ=255)

               ttGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCctgctg                                                                                                                                      >  1:3732226/1‑150 (MQ=255)

                 gCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCCTGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGc                                                                                                                                    <  2:1107124/150‑1 (MQ=255)

                 gCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCa                                                                                                                                                       >  2:3625265/1‑131 (MQ=255)

                 gCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCa                                                                                                                                                       <  1:3625265/131‑1 (MQ=255)

                 gCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCCAGCGCCTTGGCCTGCTGGc                                                                                                                                    <  1:4488095/150‑1 (MQ=255)

                   tGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCa                                                                                                                                  >  2:1201635/1‑150 (MQ=255)

                       cTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCCGCAGGCCAAGGCGCTGGGCCGGCGGGCCCCCgg                                                                                                                              >  1:4617893/1‑150 (MQ=255)

                       cTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCgg                                                                                                                              >  1:2445674/1‑150 (MQ=255)

                        tGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCCGCAGGCCAAGGCGCGGGGCCGGCTGGCCCCCGGt                                                                                                                             >  1:2006843/1‑150 (MQ=255)

                         ggCTTTCGCAAAGCTTCGCCATGCTGCTGGTAAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCCAGGCGCTGGGCCGGCTGGCCACCGGta                                                                                                                            >  1:3310955/1‑150 (MQ=255)

                          gCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCCGCCCGCCAACGCGCTGGGCCGGCTGGCCACCCGtac                                                                                                                           >  1:1526902/1‑150 (MQ=255)

                          gCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGtac                                                                                                                           >  2:3581340/1‑150 (MQ=255)

                            tttCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGtacta                                                                                                                         >  1:4685871/1‑150 (MQ=255)

                              tCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTAcc                                                                                                                       <  1:4372412/150‑1 (MQ=255)

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

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

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

                                 cAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCCGCCGGCCCAGGCGCTGGGCCGGCTGGCCCCCGGTACTAcccag                                                                                                                    >  1:1847526/1‑148 (MQ=255)

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

                                  aaaGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCAgg                                                                                                                   >  1:2211294/1‑150 (MQ=255)

                                   aaGCGTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑c                                                                                                                 <  2:3179934/150‑1 (MQ=255)

                                      cTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCAt                                                                                                              >  2:4450913/1‑150 (MQ=255)

                                       ttCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCCTg                                                                                                             >  2:4551261/1‑150 (MQ=255)

                                       ttCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATg                                                                                                             >  1:5412737/1‑150 (MQ=255)

                                            gaTGCTGCTGGTCAGCCGCATCGGCGTTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCctgctg                                                                                                                                      >  2:3903118/2‑121 (MQ=255)

                                            gaTGCTGCTGGTCAGCCGCATCGGCCTTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCctgctg                                                                                                                                      <  1:3903118/120‑1 (MQ=255)

                                              tgctgcTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTgg                                                                                                      <  1:4042118/150‑1 (MQ=255)

                                              tgctgcTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCGGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTgg                                                                                                      <  1:4791887/150‑1 (MQ=255)

                                                 tgctgGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCa                                                                                                   >  2:3587754/1‑150 (MQ=255)

                                                 tgctgGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCACCAGGCCAAAGCGCTGGGGCAGCTGGCCACACGTACTACACAGG‑CCAAGGTCCATGGCa                                                                                                   >  1:2342563/1‑150 (MQ=255)

                                                  gctgGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCCAGGCGCTGGGCCGGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCAt                                                                                                  >  2:3809313/1‑150 (MQ=255)

                                                   ctgGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCACGCCAAGGCGCGGGGCCGGCTGGCCACCCGTACAACCCAGG‑CCATGGTCCTGGGCATc                                                                                                 >  1:1862559/1‑150 (MQ=255)

                                                    tgGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATcc                                                                                                <  2:2211294/150‑1 (MQ=255)

                                                        cAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCt                                                                                            <  1:1772764/150‑1 (MQ=255)

                                                        cAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCt                                                                                            <  1:1772623/150‑1 (MQ=255)

                                                           ccGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTggg                                                                                         >  2:1151725/1‑150 (MQ=255)

                                                           ccGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGGGTGGCGCCGCCGGGCCCGCAGGCCCCGGCGCGGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCGGGGCCGCCCGCTggg                                                                                         >  2:1164433/1‑150 (MQ=255)

                                                             gCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCGGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCGGGGCCTCCCGCTGGGcc                                                                                       >  1:3914468/1‑150 (MQ=255)

                                                              cATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTggg                                                                                                                                             <  1:1294929/96‑1 (MQ=255)

                                                              cATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTggg                                                                                                                                             >  2:1294929/1‑96 (MQ=255)

                                                               ataggCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCgt                                                                                     <  2:2050867/147‑1 (MQ=255)

                                                                tCGGCATTTCCCTGGAGCATGCGGTGTTCTGGTCCATCAACGACTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCAAGCAGGCCAAGGCGCTGGGCCTGCTGGCACCCCGGACAACCCAGG‑CCCAGGTCCTGGGCAAACCGCTGGGCCggg                                                                                    >  2:1612778/1‑148 (MQ=255)

                                                                 cGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCCAGCGCCTTGGCctgctg                                                                                                                                      >  1:2555554/1‑100 (MQ=38)

                                                                 cGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCCAGCGCCTTGGCctgctg                                                                                                                                      <  2:2555554/100‑1 (MQ=38)

                                                                  ggCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGtgtg                                                                                  >  2:3945285/1‑150 (MQ=255)

                                                                    cATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTgtggt                                                                                <  2:5319426/150‑1 (MQ=255)

                                                                     aTTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTgtggtg                                                                               <  1:3129980/150‑1 (MQ=255)

                                                                        gCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGc                                                                            <  2:2672494/150‑1 (MQ=255)

                                                                          ccTGGCGCATGCGTTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGa                                                                          <  2:1882688/150‑1 (MQ=255)

                                                                          ccTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGa                                                                          <  2:2338467/150‑1 (MQ=255)

                                                                          ccTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGa                                                                          <  1:5084318/150‑1 (MQ=255)

                                                                          ccTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGa                                                                          >  2:3154313/1‑150 (MQ=255)

                                                                           cTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAg                                                                         <  2:675725/150‑1 (MQ=255)

                                                                            tGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAgg                                                                        <  2:4495751/150‑1 (MQ=255)

                                                                                gcATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCAGGGCCGGCTGGCCCCCGGTAAAAACCAGGCCCATGG‑CCAGGGCATCCCCCAGGGCCCGGTGGGGGGCGGGGCAcg                                                                    >  1:2274394/1‑149 (MQ=255)

                                                                                   tGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTggg                                                                 >  2:4225458/1‑150 (MQ=255)

                                                                                   tGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTggg                                                                 >  1:5305206/1‑150 (MQ=255)

                                                                                    gCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCGGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGc                                                                >  2:4934621/1‑150 (MQ=255)

                                                                                    gCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAAGCGCGGGGCCTGCTGGCCACCGGTACTACCCAGG‑CCCTGGTACTGGGCATCCCCCAGGGCCGGGTGGTGGGCGAGGCCCTGGGc                                                                >  1:1597360/1‑150 (MQ=255)

                                                                                        tgtTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGc                                                            <  1:2734976/150‑1 (MQ=255)

                                                                                          ttCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCGGCTGGCCACCCGTACAACCCAGG‑CCAAGGTCCTGGGCAGACCCCAGGGCCGGGTGGTGGGCGAGGCACAGGGCTGGCGt                                                          >  1:1452210/1‑150 (MQ=255)

                                                                                           tCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTa                                                         <  1:3954329/150‑1 (MQ=255)

                                                                                             tGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTgg                                                                                                      <  2:4493954/103‑1 (MQ=255)

                                                                                             tGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCAGGCCCCCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGGGGTGGGCGAGGCACGGGGCTGGCGTATc                                                       >  1:2251383/1‑150 (MQ=255)

                                                                                                  cATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCGGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCCTCCCGCTGGGCCGTGTGGTGGGCGAGGCCCGGGGCTGGCGTATCACGtt                                                  >  1:3087558/1‑150 (MQ=255)

                                                                                                   aTCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTc                                                 <  1:838625/150‑1 (MQ=255)

                                                                                                   aTCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTc                                                 <  1:2880239/150‑1 (MQ=255)

                                                                                                   aTCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTc                                                 <  1:2880363/150‑1 (MQ=255)

                                                                                                   aTCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTc                                                 >  2:3564755/1‑150 (MQ=255)

                                                                                                    tCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTcc                                                >  1:3634278/1‑150 (MQ=255)

                                                                                                       ccGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGa                                             <  2:2348485/150‑1 (MQ=255)

                                                                                                       ccGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGa                                             <  2:2222922/150‑1 (MQ=255)

                                                                                                        cGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAg                                            <  1:1451094/150‑1 (MQ=255)

                                                                                                           cTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCAGG‑CCCTGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATAACGTTACTGAGCAt                                         >  2:1471362/1‑150 (MQ=255)

                                                                                                                tGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCg                                    >  1:4999502/1‑150 (MQ=255)

                                                                                                                tGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCAGG‑CCATGGTCCTGGGCATACCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCAGTGCCg                                    >  2:5298494/1‑150 (MQ=255)

                                                                                                                  gCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGc                                                                                        <  2:2556609/96‑1 (MQ=255)

                                                                                                                         gtgtGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAgcgc                           <  2:1552386/150‑1 (MQ=255)

                                                                                                                           gtgGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTa                         <  1:868465/150‑1 (MQ=255)

                                                                                                                             ggCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGc                       <  1:5058477/150‑1 (MQ=255)

                                                                                                                             ggCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGc                       <  2:2857851/150‑1 (MQ=255)

                                                                                                                                gccgccGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGccac                    <  2:4617893/150‑1 (MQ=255)

                                                                                                                                gccgccGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGccac                    <  1:2643724/150‑1 (MQ=255)

                                                                                                                                          cagcagGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGc          <  2:5071201/150‑1 (MQ=255)

                                                                                                                                           agcagGCCCAGCGCCTTGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGcc         >  1:861855/1‑150 (MQ=255)

                                                                                                                                            gcagGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCt        <  1:3471056/150‑1 (MQ=255)

                                                                                                                                                gCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGATc    >  1:4745409/1‑150 (MQ=255)

TTCATTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGG‑CCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGATCAA  >  NC_002947/2777925‑2778220

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