BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	101,803	A→T	6.0%	K208M (AAG→ATG)	prlC →	oligopeptidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	101,803	0	A	T	6.0%	1335.3 / 32.3	447	K208M (AAG→ATG)	prlC	oligopeptidase A
Reads supporting (aligned to +/- strand): ref base A (229/171); new base T (14/17); total (243/204)
Fisher's exact test for biased strand distribution p-value = 1.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 448 total aligned reads displayed.

GCAGAGCAAGCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAGTTCGA  >  NC_002947/101654‑101950

 cAGAGCAAGCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAg                                                                                                                                                    <  1:1984826/150‑1 (MQ=255)

  agagCAAGCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAgg                                                                                                                                                   <  2:267574/150‑1 (MQ=255)

  agagCAAGCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAgg                                                                                                                                                   <  1:986982/150‑1 (MQ=255)

    agCAAGCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGc                                                                                                                                                 <  1:5013079/150‑1 (MQ=255)

    agCAAGCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGc                                                                                                                                                 <  1:1888988/150‑1 (MQ=255)

     gCAAGCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGcc                                                                                                                                                >  1:231795/1‑150 (MQ=255)

       aaGCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTg                                                                                                                                              <  1:7762652/150‑1 (MQ=255)

         gCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCCTTGGCCTggg                                                                                                                                            <  1:8435135/150‑2 (MQ=255)

         gCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGa                                                                                                                                            >  2:607944/1‑150 (MQ=255)

         gCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGa                                                                                                                                            <  1:97453/150‑1 (MQ=255)

         gCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGa                                                                                                                                            >  1:8029339/1‑150 (MQ=255)

         gCTCAGCGAGCTGGACAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGGAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCCTTGGCCTggg                                                                                                                                            <  2:3887811/150‑2 (MQ=255)

          cTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAc                                                                                                                                           <  2:1276445/150‑1 (MQ=255)

          cTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAc                                                                                                                                           <  1:7610888/150‑1 (MQ=255)

           tCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACg                                                                                                                                          >  2:1533908/1‑150 (MQ=255)

           tCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACg                                                                                                                                          >  2:5423661/1‑150 (MQ=255)

             aGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCCTTGGCCTgggctgc                                                                                                                                        <  2:306364/150‑6 (MQ=255)

              gCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCCTTGGCCTgggctgca                                                                                                                                       <  2:6530334/150‑7 (MQ=255)

              gCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCt                                                                                                                                       <  1:6565836/150‑1 (MQ=255)

                  gCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGtt                                                                                                                                   >  2:2808685/1‑150 (MQ=255)

                  gCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGtt                                                                                                                                   <  2:4166410/150‑1 (MQ=255)

                    tGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGa                                                                                                                                 <  1:2874900/150‑1 (MQ=255)

                     gggCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGAt                                                                                                                                >  2:1706348/1‑150 (MQ=255)

                        cAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCAc                                                                                                                             >  1:5234079/1‑150 (MQ=255)

                        cAGCCGTTTCTCCAATCAGCTGCGGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCAc                                                                                                                             <  2:5531503/150‑1 (MQ=255)

                         aGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCAcc                                                                                                                            <  2:4276236/150‑1 (MQ=255)

                         aGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCAcc                                                                                                                            <  2:7935707/150‑1 (MQ=255)

                            cGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCtt                                                                                                                         <  2:3614926/150‑1 (MQ=255)

                            cGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCtt                                                                                                                         <  2:402876/150‑1 (MQ=255)

                             gTTTCTCCAAGCCGCGGCTGGACGCCACCCCGGCCGGGACCACGCAGGTTCCCGCCGCCGGGGGGGGGGCCGGCCGGCCCGAGTCGGCCACGGCGCAGTGGGCCGCTGCAGCCCCGGCCAGGGGCCGGGACGGCTGGTTGATCACCCTTg                                                                                                                        <  1:33980/150‑1 (MQ=255)

                               ttCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGaa                                                                                                                      <  2:1953005/150‑1 (MQ=255)

                               ttCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGaa                                                                                                                      <  2:1410214/150‑1 (MQ=255)

                                  tcCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTc                                                                                                                   >  2:596709/1‑150 (MQ=255)

                                  tcCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTc                                                                                                                   >  1:3772461/1‑150 (MQ=255)

                                  tcCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTc                                                                                                                   >  1:6920908/1‑150 (MQ=255)

                                    cAATCAGCTGCTGGACGCCACCCAGGCCTGGACAAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTccc                                                                                                                 >  2:7724225/1‑150 (MQ=255)

                                     aaTCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTcccc                                                                                                                >  1:3952487/1‑150 (MQ=255)

                                      aTCAGCGGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCa                                                                                                               <  2:5294237/150‑1 (MQ=255)

                                       tCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAg                                                                                                              <  1:7868300/150‑1 (MQ=255)

                                            tgctgGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGctact                                                                                                         >  2:2255160/1‑150 (MQ=255)

                                            tgctgGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGctact                                                                                                         >  2:8326346/1‑150 (MQ=255)

                                            tgctgGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGctact                                                                                                         >  2:5047410/1‑150 (MQ=255)

                                                ggACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTAcgc                                                                                                     >  1:1284421/1‑150 (MQ=255)

                                                  aCGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCgg                                                                                                   <  2:4532515/150‑1 (MQ=255)

                                                  aCGCCACCCAGGCCTGGACAATGATTGTTACGTACTTAGGGTGGGTGTCCGGCTTGACCGATTCGGCCACGGGGCTTTTGGCCGGGGCAGCCGAGGCCATGTGCCGGGTCGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCgg                                                                                                   <  1:5511169/150‑1 (MQ=255)

                                                    gCCACCCAGGCCTGGACCAAGCTTGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGtg                                                                                                 <  1:7015606/150‑1 (MQ=255)

                                                     ccACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCGTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGtga                                                                                                <  2:8101875/150‑1 (MQ=255)

                                                     ccACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGtga                                                                                                >  2:17538/1‑150 (MQ=255)

                                                      cACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTAGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGtgat                                                                                               >  1:6807854/1‑150 (MQ=255)

                                                       aCCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGtgatg                                                                                              >  1:4800980/1‑150 (MQ=255)

                                                           aGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCt                                                                                          >  1:8524683/1‑150 (MQ=255)

                                                             gCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTAc                                                                                        <  1:5605295/150‑1 (MQ=255)

                                                               cTGGACCACGCATGTTACCGACGACGCGGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGc                                                                                      <  1:1824302/150‑1 (MQ=255)

                                                                 ggACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCACGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCa                                                                                    <  1:4159439/150‑1 (MQ=255)

                                                                 ggACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCa                                                                                    <  1:6360604/150‑1 (MQ=255)

                                                                 ggACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCa                                                                                    >  2:2579595/1‑150 (MQ=255)

                                                                   aCCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGc                                                                                  <  2:2073169/150‑1 (MQ=255)

                                                                     cAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGa                                                                                >  2:4942628/1‑150 (MQ=255)

                                                                      aaGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGAc                                                                               >  2:7599528/1‑150 (MQ=255)

                                                                                  gacgaAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTgcgc                                                                   >  2:7457730/1‑150 (MQ=255)

                                                                                     gaAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGaa                                                                <  1:5737632/150‑1 (MQ=255)

                                                                                     gaAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGaa                                                                >  1:1573157/1‑150 (MQ=255)

                                                                                     gaAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGaa                                                                <  1:716560/150‑1 (MQ=255)

                                                                                     gaAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTAGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGaa                                                                <  2:6807854/150‑1 (MQ=255)

                                                                                         cTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGTTGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGc                                                            <  1:4219581/150‑1 (MQ=255)

                                                                                          tGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCt                                                           >  1:8114579/1‑150 (MQ=255)

                                                                                           gcgcTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCta                                                          <  2:7112694/150‑1 (MQ=255)

                                                                                            cgcTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCtat                                                         >  1:7317299/1‑150 (MQ=255)

                                                                                             gcTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCtata                                                        >  1:2535289/1‑150 (MQ=255)

                                                                                                 gccggccTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATAcgcc                                                    >  1:2072589/1‑150 (MQ=255)

                                                                                                 gccggccTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATAcgcc                                                    >  1:5168818/1‑150 (MQ=255)

                                                                                                   cggccTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATAcgccgc                                                  >  2:5846013/1‑150 (MQ=255)

                                                                                                   cggccTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATAcgccgc                                                  <  2:3786947/150‑1 (MQ=255)

                                                                                                   cggccTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATAcgccgc                                                  >  1:6725066/1‑150 (MQ=255)

                                                                                                    ggccTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACgccgcc                                                 <  2:3738139/150‑1 (MQ=255)

                                                                                                        tGACCGCTGCGGCCACGGCGCTGAGGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACt                                             <  1:8289279/150‑1 (MQ=255)

                                                                                                        tGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACt                                             >  1:4504542/1‑150 (MQ=255)

                                                                                                        tGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACt                                             <  2:7227205/150‑1 (MQ=255)

                                                                                                        tGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACt                                             >  2:8584015/1‑150 (MQ=255)

                                                                                                         gACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTg                                            >  1:880423/1‑150 (MQ=255)

                                                                                                             gATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCAcc                                        >  1:2739632/1‑150 (MQ=255)

                                                                                                                tCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCcgc                                     <  2:3385793/150‑1 (MQ=255)

                                                                                                                  ggCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCgcgc                                   >  1:6176214/1‑150 (MQ=255)

                                                                                                                  ggCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCgcgc                                   <  1:8565356/150‑1 (MQ=255)

                                                                                                                        ggCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGa                             <  2:4615714/150‑1 (MQ=255)

                                                                                                                         gcgatgATGGCCGCTGCAGCCCAGGCAAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGAc                            <  1:1533908/145‑1 (MQ=255)

                                                                                                                           gcAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCa                          <  2:1345511/150‑1 (MQ=255)

                                                                                                                             aGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAgg                        <  1:607944/150‑1 (MQ=255)

                                                                                                                                 ggCCGCTGCAGCCCAGGCCACGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGccc                    <  1:4002212/150‑1 (MQ=255)

                                                                                                                                 ggCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGccc                    <  2:8696403/150‑1 (MQ=255)

                                                                                                                                           gtccAGGCCCTTGGCCTGGGCGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGc          >  2:6192082/3‑150 (MQ=255)

                                                                                                                                           gCCCCGGACAAGGGCATGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGTCCTACGCCTGCGACCGCGCCCTGCGCGAAGAGCTATACTCCGCCTACTGCTCCGGCGCCTCCGACCAGGGCCAGAATGCCGGc          <  2:6270237/150‑1 (MQ=255)

                                                                                                                                            tccAGGCCCTTGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGcc         >  1:3196527/2‑150 (MQ=255)

                                                                                                                                              cAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAg       >  1:7496392/1‑150 (MQ=255)

                                                                                                                                              cAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAg       >  1:4670600/1‑150 (MQ=255)

                                                                                                                                              cAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAg       >  1:1014432/1‑150 (MQ=255)

                                                                                                                                                ggCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAGtt     <  2:3173295/150‑1 (MQ=255)

                                                                                                                                                 gCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAGTTc    >  1:1731290/1‑150 (MQ=255)

                                                                                                                                                 gCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAGTTc    >  2:8754534/1‑150 (MQ=255)

                                                                                                                                                   cAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAGTTCGa  <  1:4324970/150‑1 (MQ=255)

GCAGAGCAAGCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAGTTCGA  >  NC_002947/101654‑101950

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