BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	549,805	C→G	5.2%	S159S (TCC→TCG)	tufB →	elongation factor Tu‑B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	549,805	0	C	G	5.2%	561.5 / 15.1	172	S159S (TCC→TCG)	tufB	elongation factor Tu‑B
Reads supporting (aligned to +/- strand): ref base C (83/80); new base G (3/6); total (86/86)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 174 total aligned reads displayed.

GATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGCCTACATCCCTGAGCCAGTTCGTG  >  NC_002947/549656‑549953

 aTGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCAGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCa                                                                                                                                                     <  1:591917/150‑1 (MQ=255)

  tGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCAc                                                                                                                                                    <  2:385246/150‑1 (MQ=39)

  tGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCAc                                                                                                                                                    <  2:350656/150‑1 (MQ=39)

       cGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACg                                                                                                                                               >  1:914205/1‑150 (MQ=39)

         aTGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGAc                                                                                                                                             <  2:2850856/150‑1 (MQ=39)

              gCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTccc                                                                                                                                        >  2:1557656/1‑150 (MQ=38)

              gCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTccc                                                                                                                                        <  2:2006377/150‑1 (MQ=38)

                    ccGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGcga                                                                                                                                  >  1:585500/1‑150 (MQ=37)

                       tGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAGCTCGTAGAACTCGAAGTTCGCGAACTGCTTTCGAAGTACGACTTCCCAGGcgacga                                                                                                                               <  1:3175683/150‑1 (MQ=2)

                          gCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGacac                                                                                                                            >  2:4007078/1‑150 (MQ=37)

                           cacaTCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACt                                                                                                                           >  1:432520/1‑150 (MQ=37)

                              aTCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCg                                                                                                                        >  1:4951601/1‑150 (MQ=37)

                              aTCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCg                                                                                                                        >  1:2340943/1‑150 (MQ=37)

                               tCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGa                                                                                                                       >  2:3089821/1‑150 (MQ=37)

                               tCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGa                                                                                                                       >  2:1035901/1‑150 (MQ=37)

                                      gTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGatcatcat                                                                                                                <  1:2672153/150‑1 (MQ=34)

                                      cgctcGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAACTCGTAGAACTCGAAGTTCGCGAACTGCTTTCGAAGTACGACTTCCCAGGCGACGAAGTTCCAatcatcaa                                                                                                                <  2:3112262/146‑2 (MQ=2)

                                      cgctcGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAACTCGTAGAACTCGAAGTTCGCGAACTGCTTTCGAAGTACGACTTCCCAGGCGACGAAGTTCCAatcatcaa                                                                                                                <  1:4039689/146‑2 (MQ=2)

                                         ccGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCgg                                                                                                             >  2:4554927/1‑150 (MQ=34)

                                            tCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTc                                                                                                          <  2:1487984/150‑1 (MQ=34)

                                                gTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCt                                                                                                      >  2:518861/1‑150 (MQ=34)

                                                      gTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATg                                                                                                <  1:5053137/150‑1 (MQ=34)

                                                       ttCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCATATgg                                                                                               >  1:2078410/1‑150 (MQ=21)

                                                           gTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCt                                                                                           <  1:3960075/150‑1 (MQ=37)

                                                           gTACATCGTGGTCGTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCt                                                                                           <  2:2180445/150‑1 (MQ=37)

                                                            tACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGATCCGCTCGTATGGCCCTg                                                                                          >  1:2793184/1‑150 (MQ=37)

                                                              cATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGa                                                                                        <  1:1097917/150‑1 (MQ=37)

                                                               aTCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGaa                                                                                       <  1:722588/150‑1 (MQ=37)

                                                                tCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAg                                                                                      <  2:2142348/150‑1 (MQ=37)

                                                                    ggTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCaa                                                                                  <  2:1188918/150‑1 (MQ=37)

                                                                    ggTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCaa                                                                                  <  2:1188997/150‑1 (MQ=37)

                                                                      tCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCCTCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAg                                                                                >  1:4487655/1‑150 (MQ=37)

                                                                      tCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAg                                                                                >  1:885793/1‑150 (MQ=37)

                                                                       cTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAga                                                                               <  1:2562115/150‑1 (MQ=37)

                                                                       cTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAga                                                                               <  2:430204/150‑1 (MQ=37)

                                                                       cTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAga                                                                               <  2:3158865/150‑1 (MQ=37)

                                                                       cTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAga                                                                               <  2:2532607/150‑1 (MQ=37)

                                                                       cTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAga                                                                               <  1:3607060/150‑1 (MQ=37)

                                                                       cTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAga                                                                               >  2:1825209/1‑150 (MQ=37)

                                                                        ttCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGa                                                                                                                       >  1:1241140/1‑109 (MQ=35)

                                                                        ttCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAgac                                                                              >  2:4398249/1‑150 (MQ=37)

                                                                          ccTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAgacga                                                                            <  2:2359316/150‑1 (MQ=37)

                                                                          ccTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAgacga                                                                            >  2:5429590/1‑150 (MQ=37)

                                                                           cTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAgacgac                                                                           <  2:1337406/150‑1 (MQ=37)

                                                                            tGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACa                                                                          >  2:1467009/1‑150 (MQ=37)

                                                                              aaCAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAAc                                                                        >  1:4101659/1‑150 (MQ=37)

                                                                                cAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGa                                                                      >  1:4601648/1‑150 (MQ=37)

                                                                                 aaGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTAc                                                                                                                                                >  1:2509401/1‑75 (MQ=39)

                                                                                   ggCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAAt                                                                   <  1:4802643/150‑1 (MQ=37)

                                                                                     ccGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATgg                                                                 <  1:1206205/150‑1 (MQ=35)

                                                                                      cGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATggg                                                                >  1:1632570/1‑150 (MQ=35)

                                                                                      cGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATggg                                                                <  1:1363742/150‑1 (MQ=35)

                                                                                      cGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATggg                                                                <  2:3577833/150‑1 (MQ=35)

                                                                                      cGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATggg                                                                >  2:3855933/1‑150 (MQ=35)

                                                                                        aCCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAg                                                                                      <  1:3421160/126‑1 (MQ=35)

                                                                                        aCCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAg                                                                                      >  2:3421201/1‑126 (MQ=35)

                                                                                                  cgacgCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGt                                                    <  2:3380190/150‑1 (MQ=34)

                                                                                                  cgacgCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGt                                                    <  1:2084572/150‑1 (MQ=34)

                                                                                                   gacgCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTTTCCACCTACGACTTCCCAGGCGACGACACTc                                                                                                                          <  1:4633362/79‑1 (MQ=21)

                                                                                                   gacgCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGatcatc                                                                                                                  <  2:3683676/87‑1 (MQ=33)

                                                                                                   gacgCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGatcatc                                                                                                                  >  1:3683676/1‑87 (MQ=33)

                                                                                                   gacgCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATgg                                                                 <  1:2991824/136‑1 (MQ=34)

                                                                                                   gacgCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATgg                                                                 >  2:2991824/1‑136 (MQ=34)

                                                                                                      gCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCaag                                                >  2:4394230/1‑150 (MQ=37)

                                                                                                        tGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCTGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCaagaa                                              >  2:1323727/1‑150 (MQ=37)

                                                                                                        tGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCaagaa                                              >  1:1122967/1‑150 (MQ=37)

                                                                                                            ctgctgGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATgg                                                                 <  2:882663/127‑1 (MQ=35)

                                                                                                              gctgGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGt                                        >  2:2070401/1‑150 (MQ=37)

                                                                                                               ctgGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATggg                                                                >  2:2280471/1‑125 (MQ=35)

                                                                                                               ctgGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATggg                                                                <  1:2280471/125‑1 (MQ=35)

                                                                                                                tgGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCaagaa                                              >  2:4968286/1‑142 (MQ=37)

                                                                                                                tgGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCaagaa                                              >  1:1013607/1‑142 (MQ=37)

                                                                                                                tgGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAg                                      >  2:794988/1‑150 (MQ=37)

                                                                                                                tgGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAg                                      >  1:3710944/1‑150 (MQ=37)

                                                                                                                 ggAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGa                                     >  1:712814/1‑150 (MQ=37)

                                                                                                                   aaCTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGaaa                                   <  2:2755705/150‑1 (MQ=35)

                                                                                                                        gTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTg                              <  2:652031/150‑1 (MQ=34)

                                                                                                                         tCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTgg                             >  1:4278227/1‑150 (MQ=34)

                                                                                                                         tCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTgg                             >  1:1472407/1‑150 (MQ=34)

                                                                                                                         tCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTgg                             >  1:794577/1‑150 (MQ=34)

                                                                                                                          cGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGa                            <  1:2070401/150‑1 (MQ=34)

                                                                                                                           gagaTGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGAt                           >  2:2564757/1‑150 (MQ=35)

                                                                                                                           gagaTGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGAt                           >  1:4742711/1‑150 (MQ=35)

                                                                                                                             gaTGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGc                         >  1:4103344/1‑150 (MQ=37)

                                                                                                                             gaTGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGc                         <  1:2143582/150‑1 (MQ=37)

                                                                                                                              aTGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGt                                                    >  1:851410/1‑122 (MQ=32)

                                                                                                                                 gAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGCCTAc                     >  2:2672960/1‑150 (MQ=40)

                                                                                                                                   aGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGCCTACAt                   <  2:4907852/150‑1 (MQ=255)

                                                                                                                                     ttCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGCCTACATcc                 >  1:2575956/1‑150 (MQ=255)

                                                                                                                                      tCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGCCTACATccc                <  2:4986312/150‑1 (MQ=255)

                                                                                                                                       cgcgACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGCCTACATCCCt               <  2:1666801/150‑1 (MQ=255)

                                                                                                                                       cgcgACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGCCTACATCCCt               <  1:146997/150‑1 (MQ=255)

                                                                                                                                       cgcgACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGCCTACATCCCt               <  2:1934870/150‑1 (MQ=255)

                                                                                                                                             ctgctgTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGCCTACATCCCTGAGCCa         <  1:1282437/150‑1 (MQ=255)

                                                                                                                                              tgctgTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTgg                                         <  1:2200616/117‑1 (MQ=35)

                                                                                                                                                 tgTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGCCTACATCCCTGAGCCAGTTc     <  1:1510989/150‑1 (MQ=255)

                                                                                                                                                   tCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGCCTACATCCCTGAGCCAGTTCGt   <  1:3459174/150‑1 (MQ=255)

                                                                                                                                                   tCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGCCTACATCCCTGAGCCAGTTCGt   <  1:5003752/150‑1 (MQ=255)

                                                                                                                                                   tCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGCCTACATCCCTGAGCCAGTTCGt   <  2:2685745/150‑1 (MQ=255)

                                                                                                                                                    ccACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGCCTACATCCCTGAGCCAGTTCGTg  <  2:747200/150‑1 (MQ=255)

GATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGCCTACATCCCTGAGCCAGTTCGTG  >  NC_002947/549656‑549953

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