BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,004,529	A→T	5.5%	Q184L (CAG→CTG)	PP_0866 →	transmembrane sensor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,004,529	0	A	T	5.5%	2307.2 / 81.6	717	Q184L (CAG→CTG)	PP_0866	transmembrane sensor
Reads supporting (aligned to +/- strand): ref base A (357/314); new base T (23/17); total (386/331)
Fisher's exact test for biased strand distribution p-value = 6.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.43e-01

Only 100 of 724 total aligned reads displayed.

GCCGGATGGCTCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCATATAGATG  >  NC_002947/1004380‑1004678

gCCGGATGGCTCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCa                                                                                                                                                       <  1:3910176/150‑1 (MQ=255)

   ggATGGCTCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAggg                                                                                                                                                    >  2:8770026/1‑150 (MQ=255)

       ggCTCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAgt                                                                                                                                                <  2:4313877/150‑1 (MQ=255)

       ggCTCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAgt                                                                                                                                                <  1:6954388/150‑1 (MQ=255)

       ggCTCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAgt                                                                                                                                                <  1:666581/150‑1 (MQ=255)

          tctcaacTGGAGCTCAATGGCAACACGCAGTTGCAGGTAGAGGTCAGAGCCGCCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGATCAGCGACGCCCATGACAGAAACAGGCCGATTGTGGTCCACACTGCCCAGGCCAGtgtg                                                                                                                                             >  2:6140014/7‑150 (MQ=255)

          tCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGtgtg                                                                                                                                             >  1:7739912/1‑150 (MQ=255)

          tCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGtgtg                                                                                                                                             >  1:6408061/1‑150 (MQ=255)

          tCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGtgtg                                                                                                                                             >  1:5478492/1‑150 (MQ=255)

          tCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGtgtg                                                                                                                                             >  1:2111553/1‑150 (MQ=255)

          tCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACAGCCCAGGGCCGtgtg                                                                                                                                             >  1:788745/1‑150 (MQ=255)

           cGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGa                                                                                                                                            <  2:6679342/150‑1 (MQ=255)

            gCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGCCAGTAACCGGCCGTTTGTGGTCGACACTGCCCTGGGCAgtgtcga                                                                                                                                           <  1:7479729/150‑4 (MQ=255)

            gCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCTGGGCAgtgtcga                                                                                                                                           <  2:1189498/150‑4 (MQ=255)

            gCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGAt                                                                                                                                           >  1:3450643/1‑150 (MQ=255)

            gCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGAt                                                                                                                                           <  1:5301612/150‑1 (MQ=255)

             cATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATc                                                                                                                                          <  1:1587377/150‑1 (MQ=255)

               tCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGt                                                                                                                                        >  2:2337604/1‑150 (MQ=255)

                     gCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAgg                                                                                                                                  >  2:6527670/1‑150 (MQ=255)

                      cTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGt                                                                                                                                 >  2:3450015/1‑150 (MQ=255)

                      cTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGt                                                                                                                                 <  1:6071825/150‑1 (MQ=255)

                             gCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCt                                                                                                                          >  2:7372209/1‑150 (MQ=255)

                              cAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTACCCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCtt                                                                                                                         >  2:5778116/1‑150 (MQ=255)

                              cAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCtt                                                                                                                         <  1:624491/150‑1 (MQ=255)

                              cAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCtt                                                                                                                         >  2:5323924/1‑150 (MQ=255)

                               aaCACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACCCTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCttt                                                                                                                        >  1:500425/1‑150 (MQ=255)

                               aaCACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCttt                                                                                                                        >  1:4028681/1‑150 (MQ=255)

                               aaCACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACAGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCgtt                                                                                                                        >  1:1743170/1‑150 (MQ=255)

                                    gCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGGTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGt                                                                                                                   >  1:5708527/1‑150 (MQ=255)

                                            aGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTgg                                                                                                           >  2:1462782/1‑150 (MQ=255)

                                             ggTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGa                                                                                                          <  2:5472063/150‑1 (MQ=255)

                                             ggTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGa                                                                                                          <  2:4830739/150‑1 (MQ=255)

                                              gTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGAc                                                                                                         >  1:8076435/1‑150 (MQ=255)

                                                 gAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACcca                                                                                                      >  2:1646281/1‑150 (MQ=255)

                                                   gTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACccagc                                                                                                    <  2:2677419/150‑1 (MQ=255)

                                                    ttCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACccagcc                                                                                                   >  1:6831021/1‑150 (MQ=255)

                                                      cAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACccagccag                                                                                                 <  1:7073592/150‑1 (MQ=255)

                                                       aGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCcagccagc                                                                                                >  1:319272/1‑150 (MQ=255)

                                                               ccAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCgtgt                                                                                        >  2:8534332/1‑150 (MQ=255)

                                                                   cGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCg                                                                                    >  2:4617368/1‑150 (MQ=255)

                                                                     gCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGctgg                                                                                  <  2:3832371/150‑4 (MQ=255)

                                                                       aCGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACCGTACCCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGGGGCGGTgg                                                                                <  2:3615781/150‑1 (MQ=255)

                                                                             gcCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAgg                                                                          <  1:7990613/150‑1 (MQ=255)

                                                                             gcCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAgg                                                                          <  2:227464/150‑1 (MQ=255)

                                                                                         gTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCggg                                                              >  2:4766663/1‑150 (MQ=255)

                                                                                          tGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGACCAGGGGTCGGTGCGGGt                                                             >  2:5150764/1‑150 (MQ=255)

                                                                                              gTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAg                                                         <  1:4590786/150‑1 (MQ=255)

                                                                                               tGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAgg                                                        <  1:3651090/150‑1 (MQ=255)

                                                                                               tGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAgg                                                        <  1:4265297/150‑1 (MQ=255)

                                                                                               tGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAgg                                                        <  2:7689098/150‑1 (MQ=255)

                                                                                               tGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAgg                                                        >  2:4186582/1‑150 (MQ=255)

                                                                                                 aTGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGc                                                      <  1:7571062/150‑1 (MQ=255)

                                                                                                    ttCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTAACCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAg                                                   >  2:3925467/1‑150 (MQ=255)

                                                                                                    ttCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACAGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTTACGTCCGCCAGGACCCAGCCAGCACCCCTGTGGCGGGGGAACAGGGGTCGGTGCGGGTGCAGGGCAAg                                                   >  1:455317/1‑150 (MQ=255)

                                                                                                     tCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAgg                                                  <  2:1785923/150‑1 (MQ=255)

                                                                                                      cAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAggg                                                 <  1:2444779/150‑1 (MQ=255)

                                                                                                      cAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAggg                                                 <  2:2047468/150‑1 (MQ=255)

                                                                                                       aGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGc                                                <  2:6823432/150‑1 (MQ=255)

                                                                                                       aGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGc                                                <  2:2955838/150‑1 (MQ=255)

                                                                                                        gCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCg                                               <  2:8001139/150‑1 (MQ=255)

                                                                                                         cGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGa                                              <  2:1122434/150‑1 (MQ=255)

                                                                                                         cGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGa                                              <  1:18735/150‑1 (MQ=255)

                                                                                                          gTCGCCCATGACAGTAACCGGCCGTTTTTGGTCGACACAGCCCACGGCACTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCCGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGAc                                             >  1:1018081/1‑150 (MQ=255)

                                                                                                           tCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACt                                            >  2:5023097/1‑150 (MQ=255)

                                                                                                              cccATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGc                                         <  2:2498450/150‑1 (MQ=255)

                                                                                                              cccATGACAGTAACCCGCCGTTTGTGGTCGAAACTGCCCACGGCAGTGTGATCGTCACAGGTACCCGCATTGACGTCCGCCTGGACCCAGACAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGc                                         >  2:7033647/1‑150 (MQ=255)

                                                                                                                cATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCg                                       <  2:1916297/150‑1 (MQ=255)

                                                                                                                cATGACAGTAACCGGCCGTTTGTGGTCGACACAGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCCTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCg                                       >  1:2282323/1‑150 (MQ=255)

                                                                                                                  tGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGgc                                     >  2:1329490/1‑150 (MQ=255)

                                                                                                                    aCAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGgcgc                                   >  2:3377563/1‑150 (MQ=255)

                                                                                                                      aGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCaa                                 <  1:7309509/150‑1 (MQ=255)

                                                                                                                      aGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCaa                                 <  1:881606/150‑1 (MQ=255)

                                                                                                                      aGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCGTTGACGTCCGCCGGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACACGCCGGCGCaa                                 >  1:3322157/1‑150 (MQ=255)

                                                                                                                        tAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACt                               >  2:6256737/1‑150 (MQ=255)

                                                                                                                        tAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTAACCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGGGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACt                               >  1:1415882/1‑150 (MQ=255)

                                                                                                                            cGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACg                           >  1:5830615/1‑150 (MQ=255)

                                                                                                                            cGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACg                           >  1:4346256/1‑150 (MQ=255)

                                                                                                                             ggCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACgg                          <  1:2934192/150‑1 (MQ=255)

                                                                                                                              gCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGc                         >  1:2216563/1‑150 (MQ=255)

                                                                                                                                  tttGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCggg                     <  1:2681546/150‑1 (MQ=255)

                                                                                                                                      ttgTCGACACTGCCCATGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCATCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCTCCGGCGCAACTGACGGCCGGGCAgg                 <  1:5601853/148‑1 (MQ=255)

                                                                                                                                      tgGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAgg                 >  2:8165852/1‑150 (MQ=255)

                                                                                                                                       ggTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAggg                <  2:3141180/150‑1 (MQ=255)

                                                                                                                                        gatcacacTGCCCTGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAgggg               >  2:7846575/5‑150 (MQ=255)

                                                                                                                                        gTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAgggg               >  2:8328943/1‑150 (MQ=255)

                                                                                                                                         atcacacTGCCCTGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGt              >  1:8641725/4‑150 (MQ=255)

                                                                                                                                          tcacacTGCCCTGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTc             >  2:7513344/3‑150 (MQ=255)

                                                                                                                                          tcacacTGCCCTGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTc             >  1:2821052/3‑150 (MQ=255)

                                                                                                                                           gACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCg            >  1:8338873/1‑150 (MQ=255)

                                                                                                                                            acacTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGc           <  2:4429399/150‑1 (MQ=255)

                                                                                                                                             cacTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCa          <  2:5559852/150‑1 (MQ=255)

                                                                                                                                             cacTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCa          >  1:6489017/1‑150 (MQ=255)

                                                                                                                                             cacTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCa          <  1:4876456/150‑1 (MQ=255)

                                                                                                                                                 gCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCAtata      <  2:602538/150‑1 (MQ=255)

                                                                                                                                                  cccAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCatat       >  1:4711330/1‑148 (MQ=255)

                                                                                                                                                  cccAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCATATAg     >  1:8419944/1‑150 (MQ=255)

                                                                                                                                                   ccAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCATATAGa    <  1:8167393/150‑1 (MQ=255)

                                                                                                                                                   ccAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCATATAGa    >  1:5127445/1‑150 (MQ=255)

                                                                                                                                                   ccAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCATATAGa    <  1:1176847/150‑1 (MQ=255)

                                                                                                                                                    cAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCATATAGAt   <  2:7173555/150‑1 (MQ=255)

GCCGGATGGCTCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCATATAGATG  >  NC_002947/1004380‑1004678

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