BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,781,099	A→C	5.1%	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,781,099	0	A	C	5.1%	718.7 / 17.6	230	K152Q (AAG→CAG)	sotB	sugar efflux transporter
Reads supporting (aligned to +/- strand): ref base A (109/107); new base C (7/6); total (116/114)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 234 total aligned reads displayed.

GTTCATTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGATCAAA  >  NC_002947/2780951‑2781248

 ttCATTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCCAg                                                                                                                                                     >  2:3285425/1‑150 (MQ=255)

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

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

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

   cATTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCTGCCGGGCCAGCAGGCCAAGgc                                                                                                                                                   >  2:3377316/1‑150 (MQ=255)

    aTTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCCAgcgc                                                                                                                                                  >  1:1317303/1‑147 (MQ=255)

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

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

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

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

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

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

        tCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCACCGCCGGGCCAGCAGGCCAAGGCGCTgg                                                                                                                                              >  2:2168121/1‑150 (MQ=255)

          aGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGc                                                                                                                                            <  1:2201768/150‑1 (MQ=255)

           gCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGcc                                                                                                                                           >  2:50299/1‑150 (MQ=255)

           gCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGcc                                                                                                                                           >  2:370606/1‑150 (MQ=255)

            ccACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCCGCCGGCCAAGGCGCTGGGccg                                                                                                                                          >  1:1473236/1‑149 (MQ=255)

            ccACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCGGGGCct                                                                                                                                          >  1:734101/1‑150 (MQ=255)

              acgtGCTGTCCTGGCTTTCGCACAGCGTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCAGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCCAGCGCCTTGGCctgc                                                                                                                                        <  1:3369148/147‑1 (MQ=255)

              aCTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCCAGCGCCTTGGCctgc                                                                                                                                        <  1:1074073/150‑1 (MQ=255)

              aCTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCctgc                                                                                                                                        <  1:2899670/150‑1 (MQ=255)

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

                         tGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGt                                                                                                                             >  1:4920120/1‑150 (MQ=255)

                          ggCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGta                                                                                                                            >  2:5000628/1‑150 (MQ=255)

                             tttCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCCGGCCAAGGCGCGGGGCCTGCTGGCCACCGGtactc                                                                                                                         >  1:3219502/1‑149 (MQ=255)

                             tttCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGGCTGCTGGGCACCGGtaata                                                                                                                         >  1:5368858/1‑150 (MQ=255)

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

                                  cAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTAcccag                                                                                                                    >  1:1057065/1‑148 (MQ=255)

                                      gCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCa                                                                                                                >  1:4113537/1‑150 (MQ=255)

                                        ttCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAACGCGCTGGGCCGGCTGGCCACCGGTACTACCCTGGCCATg                                                                                                              >  2:173323/1‑150 (MQ=255)

                                        gtCGCCATGCTGCTGGTCAGCCGCCTCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCCTg                                                                                                              <  2:1395789/149‑1 (MQ=255)

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

                                          cGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGgcg                                                                                                                                                  <  2:3184234/112‑1 (MQ=255)

                                          cGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGgcg                                                                                                                                                  >  1:3184234/1‑112 (MQ=255)

                                             cATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCt                                                                                                         >  2:3628622/1‑150 (MQ=255)

                                             cATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCAGGGCCTGCAGGCCACCGGTACAACCCAGGCCAAGGTcca                                                                                                         >  1:4929482/1‑149 (MQ=255)

                                              aTGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCCCCAGTACAACCCAGGCCAAGGTCCTg                                                                                                        >  2:208780/1‑150 (MQ=255)

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

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

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

                                                       gTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCg                                                                                               >  2:1964840/1‑150 (MQ=255)

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

                                                           gCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTgg                                                                                           <  1:676449/150‑1 (MQ=255)

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

                                                            ccGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCGGCTGGCCACCGGTCCTACCCTGGCCAGGGTCCTGGGCCTCCCGCgggg                                                                                          >  2:878900/1‑150 (MQ=255)

                                                               cATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCg                                                                                       >  2:2715094/1‑150 (MQ=255)

                                                                  cGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCCTCCCCCGGGGCCgggt                                                                                    >  1:3028281/1‑150 (MQ=255)

                                                                      aTTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTgtggtg                                                                                <  1:3428461/150‑1 (MQ=255)

                                                                       ttGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGtggtgg                                                                               <  2:5093473/150‑1 (MQ=255)

                                                                       ttGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCACCCCGCTGGGCCGTGtggtgg                                                                               <  1:3797690/150‑1 (MQ=255)

                                                                           ccTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCAAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGa                                                                           <  2:1664791/150‑1 (MQ=255)

                                                                           ccTGGAGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGa                                                                           <  1:3487083/150‑1 (MQ=255)

                                                                                 gcATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCCAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGGGGTGGGCGAGGCACt                                                                     >  2:4379947/1‑150 (MQ=255)

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

                                                                                    tGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTggg                                                                  <  2:2129668/150‑1 (MQ=255)

                                                                                       ggTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTg                                                               <  1:4665851/150‑1 (MQ=255)

                                                                                       ggTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCAGGCCATGGTCCTGGGCATCCCCCTGGGCCGGGTGGTGGGCGTGGCCCTGGGCTg                                                               >  1:4382162/1‑150 (MQ=255)

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

                                                                                           ttCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGt                                                           <  2:4373294/150‑1 (MQ=255)

                                                                                              tGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGGCCGCAGGCCCAGGCGCGGGGCCGGCTGGCCACCCGTACTACCCAGGCCAAGGTCCTGGGCCAACCCCTGGGCCGTGTGGTGGGCGAGGCCCGGGGCTGGCGGATc                                                        >  1:4766306/1‑150 (MQ=255)

                                                                                              tGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATc                                                        >  2:4596171/1‑150 (MQ=255)

                                                                                              tGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGGGTGGTGGGCGAGGCACTGGGCTGGCGTATc                                                        >  1:4689897/1‑150 (MQ=255)

                                                                                                gTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCCGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCGGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGGATCAc                                                      >  2:2368941/1‑150 (MQ=255)

                                                                                                   cATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCCGCAGGCCAAGGCGCGGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCCTCCCGCGGGGCCGGGTGGTGGGCGAGGCACGGGGCTGGCGTCTCACGtg                                                   >  2:4437682/1‑149 (MQ=255)

                                                                                                   cATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGGGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGtt                                                   >  2:1285512/1‑150 (MQ=255)

                                                                                                    aTCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGGTGGGGGTGCTTGGCACCGGGAATACCCTGGGCATGGTGGTGGGGATACCCCTGGGGCGTGTGGTTGGGGAGGCACTGGGGTGGGGTATCACGttg                                                  >  1:1905267/1‑149 (MQ=255)

                                                                                                     tCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCCGCCGGCCAAGGCGCGGGGCCGGCTGGCCCCCCGGACTACCCAGGCCCAGGTCCGGGGCCTCCCCCGGGGCCCTGGGGTGGGCGAGGCACGGGGCTGGCGGATCACGTTcc                                                 >  1:1994846/1‑150 (MQ=255)

                                                                                                     tCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTcc                                                 >  2:5469299/1‑150 (MQ=255)

                                                                                                     tCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCACCAGGCCAACGCGCTGGGGCGGCTGGCCACCCGTACTACCCGGGCCATGGTCCTGGGCAACCCCCTGGGCCGGGTGGTGGGCGCGGCACGGGGCTGGCGTATCACGGTcc                                                 >  1:79901/1‑150 (MQ=255)

                                                                                                     tCACCGCCTCGCTGGCGGTGCGGGTGGCGCCGCCGGGCCCGCAGGCCAATGCGCTGGGCCTGCTGGCCCCCGGGACTACCCTGGCCAGGGTCCTGGGCATCCCGCGGGGCCGTGGGGTGGGCGAGGCCCGGGGCTGGCGTATCACGTTcc                                                 >  2:3021411/1‑150 (MQ=255)

                                                                                                       aCCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCCCGCCAAGGCGCGGGGCCGGCTGGCCACCCGTACAAACCAGGCCATGGTCCGGGGCATCCCCCGGGGCCGGGTGGTGGGCGAGGCACAGGGCTGGGGTATCACGTTCCTg                                               >  1:3835208/1‑150 (MQ=255)

                                                                                                       aCCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTg                                               <  2:2780753/150‑1 (MQ=255)

                                                                                                           ccTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCGGGGCCTGCTGGCCACCGGTACTACCCTGGCCCAGGTCCTGGGCCTCCCGCTGGGCCGGGTGGTGGGCGAGGCACAGGGCTGGCGTATAACGTTCCAGAGCa                                           >  1:3434670/1‑150 (MQ=255)

                                                                                                                 tGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCg                                     <  1:4478086/150‑1 (MQ=255)

                                                                                                                  ggCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCgg                                    <  2:3918255/150‑1 (MQ=255)

                                                                                                                  ggCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCgg                                    <  2:2346331/150‑1 (MQ=255)

                                                                                                                   gCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCggg                                   <  2:5227533/150‑1 (MQ=255)

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

                                                                                                                              ggCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACGGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTCGc                        >  2:4182971/1‑150 (MQ=255)

                                                                                                                                  ccgccgGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGccacc                    <  1:2938794/150‑1 (MQ=255)

                                                                                                                                  ccgccgGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACGACCCTGGCCAGGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGACTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGccacc                    <  2:2530089/150‑1 (MQ=255)

                                                                                                                                   cgccgGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGccacca                   <  1:2697635/150‑1 (MQ=255)

                                                                                                                                   cgccgGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGccacca                   <  1:5445312/150‑1 (MQ=255)

                                                                                                                                   cgccgGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGccacca                   <  1:426292/150‑1 (MQ=255)

                                                                                                                                     ccgGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATg                 <  1:1613929/150‑1 (MQ=255)

                                                                                                                                     ccgGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATg                 >  2:3018044/1‑150 (MQ=255)

                                                                                                                                      cgGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGc                >  1:659523/1‑150 (MQ=255)

                                                                                                                                          ccAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTg            <  2:4707533/150‑1 (MQ=255)

                                                                                                                                           cagcagGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGAATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGc           <  1:3263322/150‑1 (MQ=255)

                                                                                                                                           cagcagGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACGACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGc           <  1:62731/150‑1 (MQ=255)

                                                                                                                                            agcagGCCCAGCGCCTTGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGcc          >  2:788727/1‑150 (MQ=255)

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

                                                                                                                                            agcagGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGcc          <  1:377934/150‑1 (MQ=255)

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

                                                                                                                                            agcagGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCAGGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGcc          <  2:685586/150‑1 (MQ=255)

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

                                                                                                                                               agGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGa       <  1:329378/150‑1 (MQ=255)

                                                                                                                                               agGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGa       <  2:1221401/150‑1 (MQ=255)

                                                                                                                                                  ccAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGATCa    >  1:1769599/1‑150 (MQ=255)

                                                                                                                                                   cAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGATCaa   >  2:4231556/1‑150 (MQ=255)

GTTCATTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGATCAAA  >  NC_002947/2780951‑2781248

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