BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,971,716	T→C	8.2%	K519R (AAG→AGG)	PP_5236 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,971,716	0	T	C	8.2%	884.1 / 45.6	291	K519R (AAG→AGG)	PP_5236	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (123/144); new base C (6/18); total (129/162)
Fisher's exact test for biased strand distribution p-value = 5.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

Only 100 of 304 total aligned reads displayed.

GGCGTAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAGGTTCAGTTCGGCGCAG  >  NC_002947/5971567‑5971865

ggCGTAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACt                                                                                                                                                       <  2:4492100/150‑1 (MQ=255)

 gCGTAACAATATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACtt                                                                                                                                                      <  1:7267833/150‑1 (MQ=255)

 gCGTAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACtt                                                                                                                                                      <  2:6382203/150‑1 (MQ=255)

 gCGTAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACtt                                                                                                                                                      >  2:2818890/1‑150 (MQ=255)

   gTAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTgg                                                                                                                                                    <  2:6528811/150‑1 (MQ=255)

     aaCAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcc                                                                                                                                                  >  1:3400242/1‑150 (MQ=255)

      aCAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcct                                                                                                                                                 >  2:3106561/1‑150 (MQ=255)

      aCAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAAGTACAACCTGgcct                                                                                                                                                 <  1:4663859/150‑1 (MQ=255)

      aCAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAAGTACAACCTGgcct                                                                                                                                                 <  2:5109096/150‑1 (MQ=255)

      aCAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAAGTACAACCTGgcct                                                                                                                                                 <  1:7104107/150‑1 (MQ=255)

      aCAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGAGCAGGCAGGCCAAGTACAACCTGgcct                                                                                                                                                 <  2:3576741/150‑1 (MQ=255)

        aaaaTCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcctgc                                                                                                                                               <  2:2968960/150‑1 (MQ=255)

        aaaaTCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcctgc                                                                                                                                               <  1:6446136/150‑1 (MQ=255)

        aaaaTCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcctgc                                                                                                                                               <  1:2813195/150‑1 (MQ=255)

         aaaTCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcctgcc                                                                                                                                              <  1:102451/150‑1 (MQ=255)

         aaaTCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGAGCAGGCAGGCCAAGTACAACCTGgcctgcc                                                                                                                                              <  1:1726672/150‑1 (MQ=255)

          aaTCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcctgcct                                                                                                                                             <  1:4194780/150‑1 (MQ=255)

           aTCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGcctgcctg                                                                                                                                            <  2:1529789/150‑1 (MQ=255)

               ccAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTGGTACTGGGCCTGCCtgctgc                                                                                                                                        >  2:6192346/1‑150 (MQ=255)

                 aGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGTTTGTACTTGGCCTGCCTGCTGccc                                                                                                                                      <  1:3264717/150‑1 (MQ=255)

                 aGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGccc                                                                                                                                      <  1:2549617/150‑1 (MQ=255)

                        aCCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCTGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCa                                                                                                                               <  1:1209496/150‑1 (MQ=255)

                         ccGCGCACACGCTTTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCATGTTTGTACTTGGCCTGCCTGCTGCCCGCCTCCAc                                                                                                                              <  1:880972/150‑1 (MQ=255)

                          cgcgCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCAcc                                                                                                                             <  2:1768194/150‑1 (MQ=255)

                               acacGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACg                                                                                                                        >  2:208622/1‑150 (MQ=255)

                                cacGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGc                                                                                                                       <  2:2196605/150‑1 (MQ=255)

                                    cATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCa                                                                                                                   <  2:6600105/150‑1 (MQ=255)

                                     aTTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGCCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTgg                                                                                                                                                    <  1:4024860/116‑1 (MQ=255)

                                     aTTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGCCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTgg                                                                                                                                                    >  2:4024860/1‑116 (MQ=255)

                                        ggTCTGTCCGTTGCCCAGCGCGTCGACACGCCGGCGCCTGACCGCCAGCTTGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCa                                                                                                               <  1:6947206/150‑1 (MQ=255)

                                        ggTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCa                                                                                                               <  2:1636114/150‑1 (MQ=255)

                                        ggTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCa                                                                                                               <  1:4945810/150‑1 (MQ=255)

                                          tCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGt                                                                                                             >  1:5645661/1‑150 (MQ=255)

                                          tCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGt                                                                                                             >  2:6909738/1‑150 (MQ=255)

                                           cTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTa                                                                                                            <  1:8968211/150‑1 (MQ=255)

                                           cTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTa                                                                                                            <  2:1318749/150‑1 (MQ=255)

                                            tGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTAc                                                                                                           >  1:9151024/1‑150 (MQ=255)

                                              cccGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTAccc                                                                                                         >  1:2702574/1‑150 (MQ=255)

                                                cGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCgg                                                                                                       >  1:1968730/1‑150 (MQ=255)

                                                cGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCgg                                                                                                       <  1:7376072/150‑1 (MQ=255)

                                                  ttGCCCAGCGCGTCGCCAAGCTGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCt                                                                                                     <  2:2388060/150‑1 (MQ=255)

                                                      ccAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCt                                                                                                 <  1:3577524/150‑1 (MQ=255)

                                                       cAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTc                                                                                                >  2:7591737/1‑150 (MQ=255)

                                                       cAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTc                                                                                                <  2:259407/150‑1 (MQ=255)

                                                       cAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTc                                                                                                >  1:2472636/1‑150 (MQ=255)

                                                          cgcgTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTc                                                                                             <  1:5105345/150‑1 (MQ=255)

                                                              tCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATcc                                                                                         <  2:8463935/150‑1 (MQ=255)

                                                                  cAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGt                                                                                     >  1:3934150/1‑150 (MQ=255)

                                                                   aaGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTa                                                                                    >  1:3671645/1‑150 (MQ=255)

                                                                   aaGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTa                                                                                    >  1:2838698/1‑150 (MQ=255)

                                                                    aGCAGGCGCCTGACCGCCAGCTTGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCTTGCGCAGGCTGGCCTGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTaa                                                                                   <  1:2268847/150‑1 (MQ=255)

                                                                     gCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAAt                                                                                  <  2:7141633/150‑1 (MQ=255)

                                                                     gCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAAt                                                                                  >  2:8916747/1‑150 (MQ=255)

                                                                      cAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATg                                                                                 >  2:5365185/1‑150 (MQ=255)

                                                                      cAGGCGCCTGACCGCCAGCATGGCACTGCGCACTTTCAGATCCTCAACGCCGTCCATGCGCAGGCAGGCCAGTTTGTACTTGCCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATg                                                                                 <  1:8443113/150‑1 (MQ=255)

                                                                         gcgcCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCCGGCCCGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCCCCGCACGCTGCAGCCAGTACCCGGCTGCCGCGTCATCCAGGGAATgcgc                                                                              >  2:8390032/1‑150 (MQ=255)

                                                                            ccTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCAGCGCCGTCCATGCGCAGGCAGGCCAGGTTTTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCgg                                                                           <  1:5343250/150‑1 (MQ=255)

                                                                             cTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCggg                                                                          <  1:6651970/150‑1 (MQ=255)

                                                                             cTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCggg                                                                          >  2:1999947/1‑150 (MQ=255)

                                                                             cTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCggg                                                                          <  2:2347771/150‑1 (MQ=255)

                                                                             cTGACCGCCAGCAGGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCggg                                                                          <  1:321484/150‑1 (MQ=255)

                                                                              tGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCgggg                                                                         <  1:1486908/150‑1 (MQ=255)

                                                                                   gCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTg                                                                    <  1:6260982/150‑1 (MQ=255)

                                                                                    ccAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTgg                                                                   >  1:3423880/1‑150 (MQ=255)

                                                                                    ccAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTgg                                                                   >  1:5031754/1‑150 (MQ=255)

                                                                                    ccAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGCAATGCGCCGGGGTGTTgg                                                                   >  1:6066702/1‑150 (MQ=255)

                                                                                     cAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCACACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGc                                                                  <  1:6716696/150‑1 (MQ=255)

                                                                                        cATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAg                                                               <  1:1341484/150‑1 (MQ=255)

                                                                                         aTGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAgg                                                              <  1:7972792/150‑1 (MQ=255)

                                                                                           ggCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGgcc                                                            <  2:6045578/150‑1 (MQ=255)

                                                                                           ggCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGACAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGgcc                                                            >  1:7835065/1‑150 (MQ=255)

                                                                                             cACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGgccgc                                                          >  2:7847570/1‑150 (MQ=255)

                                                                                              aCTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGgccgcc                                                         >  1:687244/1‑150 (MQ=255)

                                                                                                    cTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTg                                                   <  2:7960264/150‑1 (MQ=255)

                                                                                                     tCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGc                                                  <  1:4521603/150‑1 (MQ=255)

                                                                                                      cGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCa                                                 <  2:7229318/150‑1 (MQ=255)

                                                                                                        ttCAGATCCTCATGGCCGTCCATGCGCTGGCAGGCCAGGTTGTTCTTGGCCTGCCTGCTGCCCGCCTCCACCGCACTCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGc                                               <  1:2342409/150‑1 (MQ=255)

                                                                                                        ttCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGc                                               >  1:1756917/1‑150 (MQ=255)

                                                                                                           aGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCa                                            <  1:1172351/150‑1 (MQ=255)

                                                                                                           aGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCa                                            <  1:4329904/150‑1 (MQ=255)

                                                                                                               ccTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACa                                        >  1:8653476/1‑150 (MQ=255)

                                                                                                                cTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACaa                                       >  1:2126228/1‑150 (MQ=255)

                                                                                                                cTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACaa                                       <  1:6001843/150‑1 (MQ=255)

                                                                                                                 tCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAg                                      >  2:4422883/1‑150 (MQ=255)

                                                                                                                 tCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAg                                      >  2:5598338/1‑150 (MQ=255)

                                                                                                                     cGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCa                                  >  1:8336596/1‑150 (MQ=255)

                                                                                                                     cGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCa                                  <  2:8338344/150‑1 (MQ=255)

                                                                                                                       ccGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAgg                                >  1:6354793/1‑150 (MQ=255)

                                                                                                                          tCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGgca                             <  1:5921427/150‑1 (MQ=255)

                                                                                                                                  caggcaggCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCgg                     <  1:303492/150‑1 (MQ=255)

                                                                                                                                   aggcaggcCAAGTACAACCTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGc                    >  1:6520856/1‑150 (MQ=255)

                                                                                                                                     gcaggcCAAGTACAACCTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAg                  >  2:8636042/1‑150 (MQ=255)

                                                                                                                                      caggcCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCCCCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAgg                 <  2:4887623/150‑1 (MQ=255)

                                                                                                                                      caggcCAAGTACAACCTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAgg                 <  1:1666103/150‑1 (MQ=255)

                                                                                                                                       aggcCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCCGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAGGt                <  2:737438/150‑1 (MQ=255)

                                                                                                                                              gTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAGGTTCAGTTc         >  2:3799334/1‑150 (MQ=255)

                                                                                                                                                tGTACTTGGCCGGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAGGTTCAGTTCgg       <  2:3648545/150‑1 (MQ=255)

                                                                                                                                                     ttGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAGGTTCAGTTCGGCGCAg  >  1:4872421/1‑150 (MQ=255)

                                                                                                                                                     ttGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAGGTTCAGTTCGGCGCAg  >  2:4616750/1‑150 (MQ=255)

                                                                                                                                                     ttGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAGGTTCAGTTCGGCGCAg  >  1:8116931/1‑150 (MQ=255)

GGCGTAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAGGTTCAGTTCGGCGCAG  >  NC_002947/5971567‑5971865

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