BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	593,437	A→C	5.6%	E391A (GAG→GCG)	PP_0503 →	MFS transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	593,437	0	A	C	5.6%	751.3 / 14.6	242	E391A (GAG→GCG)	PP_0503	MFS transporter
Reads supporting (aligned to +/- strand): ref base A (103/122); new base C (4/13); total (107/135)
Fisher's exact test for biased strand distribution p-value = 8.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 252 total aligned reads displayed.

TTGCTATTCGCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGTGTCGATGTCGACGTACCA  >  NC_002947/593288‑593584

ttGCTATTCGCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGa                                                                                                                                                     >  1:492590/1‑150 (MQ=255)

ttGCTATTCGCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGa                                                                                                                                                     >  1:1688833/1‑150 (MQ=255)

  gCTATTCGCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAgg                                                                                                                                                   <  2:3235021/150‑1 (MQ=255)

  gCTATTCGCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAgg                                                                                                                                                   <  2:10447/150‑1 (MQ=255)

    tATTCGCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGgcc                                                                                                                                                 >  1:121014/1‑150 (MQ=255)

    tATTCGCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGgcc                                                                                                                                                 >  1:1661990/1‑150 (MQ=255)

     aTTCGCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGgccg                                                                                                                                                >  2:1944897/1‑150 (MQ=255)

        cGCTGGGCAGCGGTTTGGGGGCGGTGTTGGCAAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGccc                                                                                                                                             <  2:2119630/150‑1 (MQ=255)

         gCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCATTGCCGCCGCACTGTGCTAATGCATATGCATGCAACGCCCGAGGGCGGCCGAGGCCGCCCt                                                                                                                                            >  2:1558631/1‑150 (MQ=255)

         gCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCTCGGCCGCCCt                                                                                                                                            <  2:4179848/150‑1 (MQ=255)

          cTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTg                                                                                                                                           <  1:3421013/150‑1 (MQ=255)

           tGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGc                                                                                                                                          <  2:121034/150‑1 (MQ=255)

           tGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGc                                                                                                                                          >  2:2043569/1‑150 (MQ=255)

            gggCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCTCGGCCGCCCtcgc                                                                                                                                         <  1:2597137/150‑4 (MQ=255)

            gggCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCa                                                                                                                                         <  1:821961/150‑1 (MQ=255)

               cAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATg                                                                                                                                      >  2:472359/1‑150 (MQ=255)

                  cGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTc                                                                                                                                   <  1:2453674/150‑1 (MQ=255)

                  cGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACCGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTc                                                                                                                                   <  1:1312160/150‑1 (MQ=255)

                   ggTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCa                                                                                                                                  <  2:4298746/150‑1 (MQ=255)

                   ggTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCa                                                                                                                                  <  2:3431294/150‑1 (MQ=255)

                   ggTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCa                                                                                                                                  <  2:2910374/150‑1 (MQ=255)

                   ggTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCa                                                                                                                                  <  2:2534436/150‑1 (MQ=255)

                    gTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAg                                                                                                                                 <  2:2846304/150‑1 (MQ=255)

                            gCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAAc                                                                                                                         <  1:1287900/150‑1 (MQ=255)

                              ggTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACtt                                                                                                                       <  2:4361433/150‑1 (MQ=255)

                              ggTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACtt                                                                                                                       <  2:465752/150‑1 (MQ=255)

                              ggTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACtt                                                                                                                       <  1:1555294/150‑1 (MQ=255)

                               gtgtTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTg                                                                                                                      <  2:3583748/150‑1 (MQ=255)

                                tgtTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGt                                                                                                                     >  2:1792047/1‑150 (MQ=255)

                                 gtTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTa                                                                                                                    <  1:2043569/150‑1 (MQ=255)

                                      cGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCaa                                                                                                               <  2:3795674/150‑1 (MQ=255)

                                      cGAGCTGGGTGTGTGCACACTGGGGCTGGGGTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCaa                                                                                                               <  2:3193555/150‑1 (MQ=255)

                                          cTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAgacc                                                                                                           >  2:2598938/1‑150 (MQ=255)

                                             ggTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCCAgaccgac                                                                                                        >  2:1032522/1‑150 (MQ=255)

                                               tgtTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCa                                                                                                      >  1:4357813/1‑150 (MQ=255)

                                                gtTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCAt                                                                                                     <  2:2031114/150‑1 (MQ=255)

                                                  ttGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGt                                                                                                   <  1:2013352/150‑1 (MQ=255)

                                                     cacaCTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGtata                                                                                                >  2:2358708/1‑150 (MQ=255)

                                                             ggCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACggg                                                                                        >  1:3366739/1‑150 (MQ=255)

                                                              gCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgt                                                                                       >  2:2411700/1‑150 (MQ=255)

                                                               cTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtc                                                                                      <  1:728950/150‑1 (MQ=255)

                                                               cTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtc                                                                                      <  1:343945/150‑1 (MQ=255)

                                                                 gggTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcga                                                                                    >  1:1360344/1‑150 (MQ=255)

                                                                  ggTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgac                                                                                   >  1:3930810/1‑150 (MQ=255)

                                                                     tGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgtc                                                                                <  2:264144/150‑1 (MQ=255)

                                                                     tGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgtc                                                                                <  2:616115/150‑1 (MQ=255)

                                                                       cGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgtcga                                                                              >  2:1305773/1‑150 (MQ=255)

                                                                        ggTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgtcgag                                                                             >  1:1022283/1‑149 (MQ=255)

                                                                            tgtgGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtc                                                                                      <  1:263170/137‑1 (MQ=255)

                                                                            tgtgGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtc                                                                                      >  2:263170/1‑137 (MQ=255)

                                                                             gtgGGCTTGGCATGGGCATCAGTGCCGCCGCAGTGTGCTACTGGCTATGGCTGCACCGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgtcgacgtcga                                                                        <  1:3788480/150‑1 (MQ=255)

                                                                              tgGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgtcgacgtcgac                                                                       <  1:667674/150‑1 (MQ=255)

                                                                              tgGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgtcgacgtcgac                                                                       >  2:1969709/1‑150 (MQ=255)

                                                                               gggCTTGGCATGGGCATCAGTGCCGCCGCACTTTGCTACTGGCTATGGCTGCACCGCGGGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGAcc                                                                      <  1:25546/150‑1 (MQ=255)

                                                                                   ttGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGAccaagg                                                                  >  1:3631745/1‑146 (MQ=255)

                                                                                      gCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCt                                                               <  2:3862322/150‑1 (MQ=255)

                                                                                       cATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCtt                                                              >  1:2201687/1‑150 (MQ=255)

                                                                                         tGGGCATCTGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGc                                                            <  1:4080073/150‑1 (MQ=255)

                                                                                            gCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGt                                                         >  2:3955079/1‑150 (MQ=255)

                                                                                             gcacAGTGCGGCGGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGGCGACCTTGGTCTTGCTGTa                                                        >  1:2901460/4‑150 (MQ=255)

                                                                                             cttCAGTGCCGCCGCAGTGTGCTTGTGGCTATGGCTGCACGGCGCGCGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTa                                                        <  1:4242518/148‑1 (MQ=255)

                                                                                             cATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTAGGTCTTTCTGTa                                                        >  1:3269807/1‑150 (MQ=255)

                                                                                              aTCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAg                                                       <  1:2157934/150‑1 (MQ=255)

                                                                                              aTCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCCGCCGTGGTCTTGGTGTAg                                                       >  1:105840/1‑150 (MQ=255)

                                                                                                cAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGCGGGCGGCCGCGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGcc                                                     <  1:2000264/150‑1 (MQ=255)

                                                                                                  gTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCaa                                                   <  1:968934/150‑1 (MQ=255)

                                                                                                       gccgcACTGTGCTTCTGGCTTTGGCTGCAACGCGCGTGGGGGGCCGAGGCGGCCCTGCAATGTTCAGGTCAGAACTTGTAGCAAAGACGGACCTTGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCg                                              <  2:2540096/150‑1 (MQ=255)

                                                                                                          gcACTGTGCTACTGGCTTTGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCAAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTc                                           <  1:1717701/150‑1 (MQ=255)

                                                                                                               gtgCTAGTGGCTCTGGCTGCCCCGCGCGCGGGCGGCCGCGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGCCGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGcc                                      <  1:3827464/150‑1 (MQ=255)

                                                                                                               gtgCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGcc                                      <  2:1130098/150‑1 (MQ=255)

                                                                                                                 gCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGt                                    <  2:1206634/150‑1 (MQ=255)

                                                                                                                   tACTGGCTTTGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTg                                  <  2:573579/150‑1 (MQ=255)

                                                                                                                    aCTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGa                                 <  2:383807/150‑1 (MQ=255)

                                                                                                                       ggCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGc                              <  2:719143/150‑1 (MQ=255)

                                                                                                                       ggCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGGCGACGGCGACCGTGGTCTTGCGGGAGCCCAGGGCGGGCGGGCCGGTGgcgc                              >  1:105399/1‑150 (MQ=255)

                                                                                                                        gCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCt                             >  2:220636/1‑150 (MQ=255)

                                                                                                                         cTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCACTGCAATGTTCAGGTCAGAACTTTTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCTACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTg                            <  2:861452/150‑1 (MQ=255)

                                                                                                                           aTGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGGCGTCGACCATGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGc                          >  1:356319/1‑150 (MQ=255)

                                                                                                                                tGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGGCGCCGTCGACCAGTGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTgg                     >  1:2734883/1‑150 (MQ=255)

                                                                                                                                   aaCGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGGCGGGCCGTTGACGCTGGCCTTGgtgt                  >  1:1791760/1‑150 (MQ=255)

                                                                                                                                   aaCGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCGTGGTCTTGCTGTAGCCAAGGGCGGGCGGGCCGGTGACGCTGGCCTTGgtgt                  >  1:1732237/1‑150 (MQ=255)

                                                                                                                                    aCGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGtgtc                 >  2:2088694/1‑150 (MQ=255)

                                                                                                                                     cgcgcgcgGGCGGCCGAGGCCGCCCTGAAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGtgtcg                <  1:3206680/150‑1 (MQ=255)

                                                                                                                                     cgcgcgAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGtgtcg                <  2:785789/150‑1 (MQ=255)

                                                                                                                                     cgcgcgAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGtgtcg                <  2:2244942/150‑1 (MQ=255)

                                                                                                                                     cgcgcgAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCATGGTCTTTCTGTAGCCAAGGGCGGTCGGGGCGTTGAAGCTGGCCATGgtttag                >  1:917875/1‑146 (MQ=255)

                                                                                                                                      gcgcgAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGtgtcga               <  1:3955079/150‑1 (MQ=255)

                                                                                                                                      gcgcgAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACCTCCACCGTGGTCTTTCCGTAGCCAACGGCGGTCGGGGCGGGGTAGCGGGCCATGGtgtgga               >  1:2262442/1‑150 (MQ=255)

                                                                                                                                       cgggcgggCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGtgtcgat              <  1:1908837/145‑1 (MQ=255)

                                                                                                                                       cgcgcgGGCGGCCGCGGCCGCCCTGCAATGTTCCGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGtgtcgat              <  1:472241/150‑1 (MQ=255)

                                                                                                                                       cgcgTGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGtgtcgat              <  2:4067528/150‑1 (MQ=255)

                                                                                                                                           aGGGCGGCCTCGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTTCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCAGGCCTTGGtgtcgatgtcg          >  1:1274426/1‑150 (MQ=255)

                                                                                                                                             ggcggcCTCGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGGCGGGCCGTTGACGCTGGCCTTGGTGTCGAGGTCGAc        >  2:1052065/1‑150 (MQ=255)

                                                                                                                                              gcggcCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCATGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGTCGCTGGCCTTGGTGTCGATGTCGACg       >  1:1598585/1‑150 (MQ=255)

                                                                                                                                               cggcCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGTGTCGATGTCGACGt      <  1:1096455/150‑1 (MQ=255)

                                                                                                                                                ggcCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGTGTCGATGTCGACGTa     <  2:1389706/150‑1 (MQ=255)

                                                                                                                                                ggcCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGTGTCGATGTCGACGTa     <  2:3408933/150‑1 (MQ=255)

                                                                                                                                                ggcCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGGTGACGCTGGCCTTGGTGTCGAGGTCGACGTa     >  1:2351912/1‑150 (MQ=255)

                                                                                                                                                 gcCGCGGCCGCCCTGCAATGTTCAGGTCAGAACTGGTAGCCAAGACCGACCATGTAGACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGTGTCGATGTCGACGTAc    <  2:3122359/150‑1 (MQ=255)

                                                                                                                                                   cGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTAGACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGTGTCGATGTCGACGTACCa  <  1:4466577/150‑1 (MQ=255)

TTGCTATTCGCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGTGTCGATGTCGACGTACCA  >  NC_002947/593288‑593584

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