BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,508,146	T→G	5.0%	P33P (CCT→CCG)	oprC →	copper receptor OprC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,508,146	0	T	G	5.0%	1477.8 / 45.1	419	P33P (CCT→CCG)	oprC	copper receptor OprC
Reads supporting (aligned to +/- strand): ref base T (202/196); new base G (16/5); total (218/201)
Fisher's exact test for biased strand distribution p-value = 2.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.19e-01

Only 100 of 425 total aligned reads displayed.

ACGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCGACTACCTGAAGAC  >  NC_002947/5507997‑5508295

aCGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCt                                                                                                                                                       <  1:2522109/150‑1 (MQ=255)

  gAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTgg                                                                                                                                                     <  2:3734546/150‑1 (MQ=255)

  gAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTgg                                                                                                                                                     >  1:2753073/1‑150 (MQ=255)

   aaGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGc                                                                                                                                                    <  1:5296851/150‑1 (MQ=255)

    aGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGcc                                                                                                                                                   <  2:2473187/150‑1 (MQ=255)

     gCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCa                                                                                                                                                  <  2:7431408/150‑1 (MQ=255)

      ccAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCAt                                                                                                                                                 <  1:5646268/150‑1 (MQ=255)

      ccAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCAt                                                                                                                                                 <  2:1548746/150‑1 (MQ=255)

       cAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATg                                                                                                                                                >  2:7077115/1‑150 (MQ=255)

        aggaggCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGa                                                                                                                                               >  2:6442095/1‑150 (MQ=255)

        aggaggCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGa                                                                                                                                               >  2:2693346/1‑150 (MQ=255)

         ggaggCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGaa                                                                                                                                              >  1:4521728/1‑150 (MQ=255)

         ggaggCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGACCGCCGGCCCTGGCCATGaa                                                                                                                                              >  1:2898770/1‑150 (MQ=255)

             gCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGAcc                                                                                                                                          <  2:1112147/150‑1 (MQ=255)

             gCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGAcc                                                                                                                                          <  1:1296063/150‑1 (MQ=255)

              cTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCa                                                                                                                                         <  2:6289312/150‑1 (MQ=255)

              cTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCa                                                                                                                                         >  2:5735751/1‑150 (MQ=255)

               ttCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTCCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCAt                                                                                                                                        <  1:5103879/150‑1 (MQ=255)

                 cgcgCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCAtgt                                                                                                                                      <  1:5805988/150‑1 (MQ=255)

                 cgcgCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCAtgt                                                                                                                                      <  2:4730676/150‑1 (MQ=255)

                 cgcgCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCAtgt                                                                                                                                      <  2:5828119/150‑1 (MQ=255)

                   cgcTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGc                                                                                                                                    >  1:4127405/1‑150 (MQ=255)

                       aaCTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCacga                                                                                                                                >  2:2859936/1‑150 (MQ=255)

                          tCATGACTGATCGAGTGGAATCATTATGTCCGGCTGCACGGGTGTTTTTGGCTTGTCCCTGCGTGGGACCATCGCCGCGCTGTGGGGCTCGCTGCTCGCGCGCATGGGGCTGGGCGCCGGCCCTGGCCATGTAGTGCATGTGCACGACGc                                                                                                                             >  2:5991001/1‑150 (MQ=255)

                          tCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGc                                                                                                                             >  1:2938207/1‑150 (MQ=255)

                           cATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCa                                                                                                                            <  2:606228/150‑1 (MQ=255)

                              gACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCAccc                                                                                                                         <  1:3618119/150‑1 (MQ=255)

                               aCTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCg                                                                                                                        <  1:45905/150‑1 (MQ=255)

                                 tGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAg                                                                                                                      >  1:1912109/1‑150 (MQ=255)

                                   aTCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGtt                                                                                                                    <  2:2029132/150‑1 (MQ=255)

                                   aTCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGtt                                                                                                                    <  2:2584637/150‑1 (MQ=255)

                                    tCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTg                                                                                                                   >  1:1798438/1‑150 (MQ=255)

                                     cGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGa                                                                                                                  <  2:1588432/150‑1 (MQ=255)

                                       aCTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGc                                                                                                                <  2:5083127/150‑1 (MQ=255)

                                       aCTGGAATCATGATGTCCGGCGGCACCCCTGTTTGTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAGGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGc                                                                                                                <  2:11777/150‑1 (MQ=255)

                                        cTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGATCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGcc                                                                                                               <  1:3197185/150‑1 (MQ=255)

                                        cTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGcc                                                                                                               >  1:7331278/1‑150 (MQ=255)

                                         tGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGccc                                                                                                              <  1:4975753/150‑1 (MQ=255)

                                           gAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGa                                                                                                            <  2:1427778/150‑1 (MQ=255)

                                           gAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGa                                                                                                            <  1:5512189/150‑1 (MQ=255)

                                                aTTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTg                                                                                                       <  1:7421875/150‑1 (MQ=255)

                                                aTTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTg                                                                                                       <  1:1353907/150‑1 (MQ=255)

                                                 ttATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGa                                                                                                      <  1:3068027/150‑1 (MQ=255)

                                                   aTGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATc                                                                                                    >  1:4321629/1‑150 (MQ=255)

                                                     gTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCAc                                                                                                  >  1:1000367/1‑150 (MQ=255)

                                                       ccGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCAccg                                                                                                >  1:3964890/1‑150 (MQ=255)

                                                         ggCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCAccgcc                                                                                              <  2:2804858/150‑1 (MQ=255)

                                                         ggCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCAGGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCAccgcc                                                                                              <  2:4204638/150‑1 (MQ=255)

                                                           cTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGt                                                                                            >  1:7173965/1‑150 (MQ=255)

                                                             gCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTgg                                                                                          <  2:4321629/150‑1 (MQ=255)

                                                                   cTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGa                                                                                    <  2:1945104/150‑1 (MQ=255)

                                                                   cTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGa                                                                                    <  1:3299791/150‑1 (MQ=255)

                                                                   cTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGa                                                                                    <  1:477409/150‑1 (MQ=255)

                                                                     gTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGt                                                                                  <  1:4885629/150‑1 (MQ=255)

                                                                      tttgggCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGtt                                                                                 <  2:167338/145‑1 (MQ=255)

                                                                         ttGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTcgc                                                                              >  1:6224937/1‑150 (MQ=255)

                                                                         ttGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTcgc                                                                              <  1:2243147/150‑1 (MQ=255)

                                                                           gCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTcgccg                                                                            <  1:5402625/150‑1 (MQ=255)

                                                                             cTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCt                                                                          >  1:3727411/1‑150 (MQ=255)

                                                                              ttGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTg                                                                         >  1:5686713/1‑150 (MQ=255)

                                                                              ttGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTg                                                                         <  1:679903/150‑1 (MQ=255)

                                                                               tGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGa                                                                        <  2:5687550/150‑1 (MQ=255)

                                                                               tGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGa                                                                        <  2:336998/150‑1 (MQ=255)

                                                                                gTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCAt                                                                                                                                        >  2:5400520/1‑85 (MQ=255)

                                                                                        gTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCa                                                               >  2:7347492/1‑150 (MQ=255)

                                                                                        gTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCa                                                               >  2:6732945/1‑150 (MQ=255)

                                                                                         tGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCAc                                                              <  2:354589/150‑1 (MQ=255)

                                                                                           ggACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCa                                                            >  2:7128558/1‑150 (MQ=255)

                                                                                              ccATCGCCTCCCTGTGGGGCTCGATGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGACATGACGCCGATGTGCACGACGCACCCTATTTGAGCCCGACGGTGTTCACAGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAAcc                                                         <  1:214634/150‑1 (MQ=255)

                                                                                              ccATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAAcc                                                         <  2:1782540/150‑1 (MQ=255)

                                                                                               cATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCAt                                                                                                                                        <  2:3836488/70‑1 (MQ=255)

                                                                                                aTCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCAGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCg                                                       >  1:2603802/1‑150 (MQ=255)

                                                                                                   gccgccCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGaaa                                                    <  1:4257848/150‑1 (MQ=255)

                                                                                                    ccgccCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAg                                                   >  1:1622278/1‑150 (MQ=255)

                                                                                                      gccCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACGTTGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGTAAGAc                                                 <  2:4396538/150‑1 (MQ=255)

                                                                                                        ccTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGAccc                                               <  1:2726572/150‑1 (MQ=255)

                                                                                                          tgtgCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCcgc                                             <  1:673259/150‑1 (MQ=255)

                                                                                                           catgggCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGGTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCcgcg                                            >  1:4190469/5‑150 (MQ=255)

                                                                                                              cGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCgcca                                         <  1:6039130/150‑1 (MQ=255)

                                                                                                                  tCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCgccagcca                                     >  2:6581622/1‑150 (MQ=255)

                                                                                                                  tCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCgccagcca                                     <  2:2252365/150‑1 (MQ=255)

                                                                                                                    gctgctCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGt                                   >  2:6795483/1‑150 (MQ=255)

                                                                                                                        ctCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCg                               <  2:4609791/150‑1 (MQ=255)

                                                                                                                         tCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCgg                              >  1:6685819/1‑150 (MQ=255)

                                                                                                                          ggtcttcATGGCCAGGGCCGGCGGCCAGGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGc                             >  2:4902736/7‑150 (MQ=255)

                                                                                                                              ttcATGGCCAGGGCCGGCGGCCAGGGCCATGGAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGc                         >  2:6881670/3‑150 (MQ=255)

                                                                                                                                 aTGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGAc                      <  2:3272413/150‑1 (MQ=255)

                                                                                                                                    gCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGgc                   >  1:1684516/1‑150 (MQ=255)

                                                                                                                                     cccTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCg                               <  2:1214726/137‑1 (MQ=255)

                                                                                                                                     ccagggCCGGCGGCCAGGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGgcg                  >  2:3907547/5‑150 (MQ=255)

                                                                                                                                       cTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGcc                <  2:2958865/150‑1 (MQ=255)

                                                                                                                                        tGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCg               >  1:6543894/1‑150 (MQ=255)

                                                                                                                                         ggCCGGCGGCCAGGGCCATGGAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCGa              >  1:3356448/1‑150 (MQ=255)

                                                                                                                                             gccgGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCGACTAc          <  1:6344554/150‑1 (MQ=255)

                                                                                                                                             gccgGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCGACTAc          <  1:3757591/150‑1 (MQ=255)

                                                                                                                                              ccgGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCGACTAcc         <  2:5231593/150‑1 (MQ=255)

                                                                                                                                               cgGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCGACTACCt        >  1:819723/1‑150 (MQ=255)

                                                                                                                                                   ccTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCGACTACCTGAAg    <  1:2987242/150‑1 (MQ=255)

                                                                                                                                                    cTGGCCATGAAGACCATGTGCACGACGCACCCGGGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCGACTACCTGAAGa   >  1:1213647/1‑150 (MQ=255)

                                                                                                                                                    cTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCGACTACCTGAAGa   <  2:5644450/150‑1 (MQ=255)

ACGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCGACTACCTGAAGAC  >  NC_002947/5507997‑5508295

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