BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	780,841:1	+T	7.0%	intergenic (+1/‑2)	valT → / → lysW	tRNA‑Val/tRNA‑Lys

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	780,841	1	.	T	7.0%	688.9 / 43.4	183	intergenic (+1/‑2)	valT/lysW	tRNA‑Val/tRNA‑Lys
Reads supporting (aligned to +/- strand): ref base . (70/100); new base T (4/6); total (77/106)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

Only 100 of 193 total aligned reads displayed.

GGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCG  >  NC_000913/780692‑780989

gggCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCcaccac                                                                                                                                                       <  2:3581013/150‑1 (MQ=255)

 ggCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCaccac‑c                                                                                                                                                     >  2:3363014/1‑150 (MQ=255)

  gCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑Cg                                                                                                                                                    <  2:2385115/150‑1 (MQ=255)

   cGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑Cgg                                                                                                                                                   >  2:2036906/1‑150 (MQ=255)

      gaCGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTc                                                                                                                                                <  2:1847622/150‑1 (MQ=255)

      gaCGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTc                                                                                                                                                >  1:2700381/1‑150 (MQ=255)

      gaCGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGACTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTc                                                                                                                                                >  1:522711/1‑150 (MQ=255)

        cGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGt                                                                                                                                              >  1:3649131/1‑150 (MQ=255)

        cGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGCTCga                                                                                                                                              >  1:2927552/1‑149 (MQ=255)

          aaGTCGAGTCATCCTGCACGACCCACCACTAACATCGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTa                                                                                                                                            <  1:1933965/150‑1 (MQ=255)

          aaGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTa                                                                                                                                            <  1:3434695/150‑1 (MQ=255)

          aaGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTa                                                                                                                                            <  2:1046039/150‑1 (MQ=255)

          aaGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTa                                                                                                                                            <  2:72880/150‑1 (MQ=255)

          aaGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTa                                                                                                                                            <  2:502256/150‑1 (MQ=255)

           aGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAg                                                                                                                                           <  2:1693062/150‑1 (MQ=255)

            gTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CCGGTGCTTAgg                                                                                                                                          >  2:285041/1‑149 (MQ=255)

             tCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CCGGTCGTTAggt                                                                                                                                         >  1:622650/1‑148 (MQ=255)

                aGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAg                                                                                                                                      >  1:4281294/1‑150 (MQ=255)

                  tCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGtt                                                                                                                                    <  2:3991619/150‑1 (MQ=255)

                  tCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGtt                                                                                                                                    <  1:156391/150‑1 (MQ=255)

                   cATCCTGCACGACCCACCACTAACATCGTTAGTTGTAGTATCCAGCGTAGGATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGGCGGCGGTGCGATCCCGTCATCACCCACCCC‑CGGGGCGTTAGCGCAgggg                                                                                                                                   >  1:851906/1‑147 (MQ=18)

                     tCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGTGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGATCAGTTGGt                                                                                                                                 >  2:410295/1‑150 (MQ=255)

                     tCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCCCCCACCAC‑CGGGTCGATAGCTCAGTTGGt                                                                                                                                 >  2:4213433/1‑150 (MQ=255)

                      ccTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACCAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTa                                                                                                                                <  1:4716182/150‑1 (MQ=255)

                      ccTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACATGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACAAC‑CGGGTCGTTAGCTCAGTTGGTa                                                                                                                                <  2:3633058/150‑1 (MQ=255)

                      ccTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTa                                                                                                                                <  1:4370292/150‑1 (MQ=255)

                      ccTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTa                                                                                                                                <  1:3151334/150‑1 (MQ=255)

                         gCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTaga                                                                                                                              >  2:3925239/1‑150 (MQ=255)

                          cACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTgg                                                                                                                                  <  2:1729232/145‑1 (MQ=255)

                          cACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTagag                                                                                                                             >  1:1689708/1‑150 (MQ=255)

                           aCGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCa                                                                                                                           <  2:2278747/150‑1 (MQ=255)

                           aCGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCa                                                                                                                           <  2:3371528/150‑1 (MQ=255)

                            cGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTATCTAAGTTGGTAGATCAg                                                                                                                          >  2:4045430/1‑150 (MQ=255)

                            cGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAg                                                                                                                          <  2:3494815/150‑1 (MQ=255)

                              aCCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGtt                                                                                                                        <  2:3818650/150‑1 (MQ=255)

                                 caccacTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTATCTCAGTTGGTAGAGCAGTTGAc                                                                                                                     >  1:1055295/1‑150 (MQ=255)

                                  accacTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACt                                                                                                                    <  1:2002012/150‑1 (MQ=255)

                                     acTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACtttt                                                                                                                 >  2:930324/1‑150 (MQ=255)

                                      cTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTa                                                                                                                <  1:1875560/150‑1 (MQ=255)

                                      cTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTa                                                                                                                <  1:4758817/150‑1 (MQ=255)

                                       tAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTTCAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCTGTTGACTTTTaa                                                                                                               <  1:3499/150‑1 (MQ=255)

                                       tAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACGCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTaa                                                                                                               <  1:679772/150‑1 (MQ=255)

                                       tAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTaa                                                                                                               <  1:4567657/150‑1 (MQ=255)

                                       tAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTaa                                                                                                               <  1:4215289/150‑1 (MQ=255)

                                         acctagtttgtTGTAGTATACAGAGTATTTTGGGGTGATTATCTCTGCTGTGAGAGCACCTGCCTGACCTGGAGGGGGTCGTCGGTTCGATCCCGTCATCACGAAACAG‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATc                                                                                                             <  1:3322541/147‑1 (MQ=14)

                                         aCATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATc                                                                                                             <  1:1155118/150‑1 (MQ=255)

                                           aTAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAg                                                                                                                                           <  2:639459/119‑1 (MQ=19)

                                           aTAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCaa                                                                                                           <  1:2281674/150‑1 (MQ=255)

                                            tagttagtTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAAt                                                                                                          <  1:151832/150‑1 (MQ=255)

                                             agttagttGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAAt                                                                                                          <  1:1822154/150‑1 (MQ=17)

                                               ttagttGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTgg                                                                                                       <  2:1303197/150‑1 (MQ=39)

                                               ttagttGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCACGTGGTCGAGCAGTTGACATTTAATCAAATgg                                                                                                       >  2:1924227/1‑150 (MQ=14)

                                                      tagtaTCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAgg                                                                                                >  1:3491859/1‑150 (MQ=35)

                                                        gtaTCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGtt                                                                                              <  1:2144219/150‑1 (MQ=35)

                                                         taTCCATAGTAGTATCGTGTGATTAGCTCCGCTGTGATAGCACCTCCCTTACCTGGAGGGGGTCGTCGGTTCGATCCCGTCATAACGCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTc                                                                                             <  2:303711/150‑1 (MQ=11)

                                                         taTCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTc                                                                                             <  2:2362345/150‑1 (MQ=35)

                                                                 gtagttTCGGGTGATTAGCTCAGCTGGGAGATCACCTCCCTTACTAGGAGGGGGCCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTg                                                                                     <  1:1260558/150‑1 (MQ=255)

                                                                       tCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGAc                                                                               >  2:4340437/1‑150 (MQ=33)

                                                                         gggTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGAccc                                                                             >  1:1046262/1‑150 (MQ=33)

                                                                         gggTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGAccc                                                                             >  1:1046106/1‑150 (MQ=33)

                                                                         gggTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGGTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGAccc                                                                             >  2:3822699/1‑150 (MQ=18)

                                                                            tGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACccacc                                                                          >  1:2882224/1‑150 (MQ=33)

                                                                             gATTAGCTCATCTGGGAGAGATCCTCCCTTACCTGGAGTGGGTCGGCGGTTCGATCCCGTCCTAACCCACCAG‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACccacca                                                                         <  2:3268481/150‑1 (MQ=11)

                                                                             gATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCaccac‑cta                                                                                                                                                   >  1:1654477/1‑74 (MQ=14)

                                                                              aTTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAg                                                                        >  1:2028795/1‑150 (MQ=33)

                                                                              aTTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAg                                                                        >  1:2770/1‑150 (MQ=33)

                                                                               ttAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGt                                                                       >  2:1799109/1‑150 (MQ=33)

                                                                                          gggAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGa                                                            >  2:4272885/1‑150 (MQ=35)

                                                                                           ggAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGaa                                                           <  2:4376038/150‑1 (MQ=35)

                                                                                           ggAGAGCAACTCACTTACTAGGAGGGGTTCGGCGGTTCGATCCCGTCATCACCCACCAG‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGaa                                                           <  1:3924853/150‑1 (MQ=11)

                                                                                             agagCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGa                                                         <  1:4633625/150‑1 (MQ=38)

                                                                                             agagCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGa                                                         <  2:4058395/150‑1 (MQ=38)

                                                                                                 cACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGa                                                     >  2:1776301/1‑150 (MQ=40)

                                                                                                  aCCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCAGCACGACCCACCAGTTTTAACATCGAAGAAAGAt                                                    >  1:2480668/1‑150 (MQ=255)

                                                                                                  aCCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTCTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGAt                                                    >  2:3153339/1‑150 (MQ=255)

                                                                                                   ccTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGGAGAGCCGTTGACTTTTAATCAATTGGTCGCCGGTTCGAATACTGCACGACCCCCCAGTTTTAACAACGAAGACAGagg                                                   >  1:999222/1‑148 (MQ=14)

                                                                                                        cTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAg                                              <  2:4153035/150‑1 (MQ=40)

                                                                                                         ttACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCCCGTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCCGGGTCCAAGCCAGCACCACCCACCCGGTTTAACATCCAAGACCGATGGGAAGc                                             >  1:4335221/1‑150 (MQ=11)

                                                                                                         ttACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGGAGAGCAGTTGACTTTTAATAAATTGGTCGCAGGTTCGAATACTGCACGCCCCACCAGGTTTAACATCGAAGACAGATGTTAAGc                                             >  1:3871327/1‑150 (MQ=14)

                                                                                                              aggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTa                                        <  1:4010120/150‑1 (MQ=255)

                                                                                                              aggaggGGGTCGGCGGTTCGATCCCGTCATCACCCAACAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTa                                        <  1:2318642/150‑1 (MQ=255)

                                                                                                               ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAg                                       <  2:4710254/150‑1 (MQ=255)

                                                                                                               ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAg                                       <  2:1055295/150‑1 (MQ=255)

                                                                                                                 aggGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTATGa                                     <  1:4045430/150‑1 (MQ=255)

                                                                                                                 aggGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGa                                     <  2:1196696/150‑1 (MQ=255)

                                                                                                                     ggTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAAc                                 <  2:3649131/150‑1 (MQ=255)

                                                                                                                          gcggTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGc                            <  1:4790633/150‑1 (MQ=255)

                                                                                                                            ggTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGCAGACAGATGTTAAGCGTGTAGGATAACGTTGCGt                          <  1:99235/150‑1 (MQ=255)

                                                                                                                            ggTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGt                          <  2:4484938/150‑1 (MQ=255)

                                                                                                                            ggTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGt                          <  1:3518634/150‑1 (MQ=255)

                                                                                                                            ggTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGt                          <  1:372464/150‑1 (MQ=255)

                                                                                                                                  aTCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGCTCGAATCCTGCACGACCCACCAGTTTTaa                                                                  <  1:3216336/104‑1 (MQ=25)

                                                                                                                                     ccGTCATCACCCACCAC‑GGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACgg                 <  1:2876932/150‑1 (MQ=255)

                                                                                                                                         catcaCCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTACTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCg             <  2:3491859/150‑1 (MQ=255)

                                                                                                                                         catcaCCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCg             <  2:664248/150‑1 (MQ=255)

                                                                                                                                              ccccccAC‑CGTGGCGTTAGCTCAGTTTTTAGTGCTGTTGTCTTTTAATCTATTGGTGTCAGGTTCGATTCCTGCACGACCCACCAGTTTTACCATCGAAGACAGATGTTAAGCGTGTAGGTTAACGTTGCGTCAGCAACGGCCCGTAggg        <  2:2350428/150‑1 (MQ=14)

                                                                                                                                                 accac‑cGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGTCAGATGTTAAGCGTGAAGGATAACGTTGCGACAGCAACGGCCCGTAGGgcga     >  2:111142/1‑150 (MQ=16)

                                                                                                                                                 accac‑cGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGgcgg     >  2:2672219/1‑149 (MQ=255)

                                                                                                                                                  ccac‑cGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGgcgag    <  2:2770/150‑1 (MQ=255)

                                                                                                                                                    ac‑cGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTATTAAAATCGTACACAGATGTTAAGAGTGTAGGATAACGTAGCGGAAGCAAAGGCCCGTAGGgcgagcg  >  2:1571824/1‑150 (MQ=11)

GGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCG  >  NC_000913/780692‑780989

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