BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	780,841:1	+T	5.4%	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	5.4%	483.9 / 26.2	130	intergenic (+1/‑2)	valT/lysW	tRNA‑Val/tRNA‑Lys
Reads supporting (aligned to +/- strand): ref base . (38/85); new base T (3/4); total (41/89)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.79e-01

Only 100 of 134 total aligned reads displayed.

GGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGC  >  NC_000913/780692‑780988

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

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

    gagaCGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCCGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑Cggg                                                                                                                                                 >  1:1017080/1‑150 (MQ=255)

     agaCGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGt                                                                                                                                                <  2:1795038/150‑1 (MQ=255)

     agaCGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGt                                                                                                                                                >  2:2591206/1‑150 (MQ=255)

       aCGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGCGTCg                                                                                                                                              >  1:2131206/1‑150 (MQ=255)

        cGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGGTCCCGTCATCACACACCAC‑CGGGTCGt                                                                                                                                             >  1:329363/1‑150 (MQ=255)

         gAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCggt                                                                                                                                            >  1:2367507/1‑148 (MQ=255)

          aaGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTCTCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTa                                                                                                                                           <  2:1634289/150‑1 (MQ=255)

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

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

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

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

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

          aaGTCGAGTCATCCTGCACGACCCACCACGACCATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTa                                                                                                                                           <  2:1639352/150‑1 (MQ=255)

            gTCGAGTCATCCTGCTCGTCCCCCCTCTCAGCTAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAg                                                                                                                                          <  2:2819260/150‑1 (MQ=14)

               gAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCa                                                                                                                                      <  2:2815091/150‑1 (MQ=255)

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

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

                    aTCCTGCACGACCCACCACTAACATCGTTAGTTGTAGTATACAGCGTAGTATCGGGTGATTAGCTAAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCACGTCATCACCCACCAC‑CGGGTCGTTAGCTCAggggg                                                                                                                                 >  1:1078580/1‑146 (MQ=18)

                    aTCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTg                                                                                                                                  >  1:12293/1‑150 (MQ=255)

                    aTCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGGAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTg                                                                                                                                  >  2:140639/1‑150 (MQ=255)

                     tCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGt                                                                                                                                <  2:1225955/150‑1 (MQ=255)

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

                        tGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTaga                                                                                                                             <  2:436120/150‑1 (MQ=255)

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

                             gACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CCGGTCCTTAGCGCAGGTGGTATAGccat                                                                                                                        >  1:1097371/1‑147 (MQ=17)

                                ccaccaCTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGa                                                                                                                     >  1:2157198/1‑150 (MQ=255)

                                 caccacTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGAc                                                                                                                    >  2:1072964/1‑150 (MQ=255)

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

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

                                  accacTAACATAGTTAGTTGTAGTATCCAGCGCAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGAc                                                                                                                    <  1:140639/150‑1 (MQ=255)

                                     acTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACATCCCAAAAGCGGAGGGGGTCGACGGATCGATCCCGACGTCGAACAACAC‑CGGGTGCTAAGGACAGTTGGTAGAGCGGTTgggtggg                                                                                                                >  2:1301047/1‑144 (MQ=9)

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

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

                                       tAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTaa                                                                                                              >  1:962743/1‑150 (MQ=255)

                                       caaCCTAGTTAGTTGTAGTATCCAGCGTCGTCCAGGGTGATTCGCTCAGCTGGGAGAGCACCTCGCTTACTAGGAGGGGGTCGGCGTTTCGATCCGGTAATCACGCACCAG‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTaa                                                                                                              <  2:2296560/149‑1 (MQ=14)

                                         aCATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATc                                                                                                            >  1:2563051/1‑150 (MQ=255)

                                         aCATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATc                                                                                                            >  1:1491014/1‑150 (MQ=255)

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

                                            tagttagtTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAAt                                                                                                         >  1:1952781/1‑150 (MQ=255)

                                            tagttagtTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTGGGTAGAGCAGTTGACTTTTAATCAAt                                                                                                         >  1:1877902/1‑150 (MQ=255)

                                             agttagttGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAAtt                                                                                                        <  1:2176033/150‑1 (MQ=40)

                                             agttagttGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAAtt                                                                                                        <  1:1672246/150‑1 (MQ=40)

                                               ttagttGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACCAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTgg                                                                                                      <  2:2999908/150‑1 (MQ=255)

                                               ttagttGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTg                                                                                                       >  2:2749938/1‑150 (MQ=11)

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

                                                  gttGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCg                                                                                                   <  1:1779793/150‑1 (MQ=37)

                                                               gCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATcc                                                                                      >  1:753856/1‑150 (MQ=35)

                                                               gCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATcc                                                                                      >  1:8086/1‑150 (MQ=35)

                                                               gCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATcc                                                                                      <  1:1744973/150‑1 (MQ=35)

                                                               gCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATcc                                                                                      >  1:204288/1‑150 (MQ=35)

                                                                   agtaTCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCa                                                                                  <  1:2295508/150‑1 (MQ=33)

                                                                   agtaTCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCa                                                                                  <  2:454986/150‑1 (MQ=33)

                                                                      aTCGGGTGTTTAGCTCAGCTGGGAGAGCACCTCCCTTACATGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGa                                                                               <  2:1658369/150‑1 (MQ=18)

                                                                      aTCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGa                                                                               <  2:2217153/150‑1 (MQ=33)

                                                                          ggTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACcca                                                                           >  2:1044686/1‑150 (MQ=33)

                                                                           gTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACccac                                                                          >  1:184874/1‑150 (MQ=33)

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

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

                                                                                   cTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTa                                                                  <  2:1833896/150‑1 (MQ=33)

                                                                                   cTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTa                                                                  <  2:2559046/150‑1 (MQ=33)

                                                                                        cTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATc                                                             <  1:2591206/150‑1 (MQ=33)

                                                                                        cTGGGAGAGCACCTCCCTGACACGGAGGGGGGCGGCGGTTCGATCCCGTCATCACCCACCTC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATc                                                             <  1:790952/150‑1 (MQ=255)

                                                                                           ggAGAGCACCTCCCTTACACGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTTGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGaa                                                          <  1:1732413/150‑1 (MQ=25)

                                                                                           ggAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCCC‑CGGGTCGTTAGCTCTGTTGGTAGAGCAGTTGACTTTTACTCAATTGGTCGCTGGTTCGATTCCTGCACGACCCACCAGTTTTAACATCGaa                                                          <  2:2799494/150‑1 (MQ=11)

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

                                                                                            gagaGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAg                                                         >  2:1215296/1‑150 (MQ=37)

                                                                                              gagaACCTCCATTCCATGGAGGTGGTCGTGGGTTCGATCCCTTACTAACGAACCAG‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGAc                                                       <  2:335396/150‑1 (MQ=11)

                                                                                              gagCACCGCCCTTACCAGCAGGTGGTCGGTGCTTAGTTCCCGTCCGCACCAAACAG‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGAc                                                       <  2:335423/150‑1 (MQ=11)

                                                                                               agCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGTCTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACa                                                      <  1:2350123/150‑1 (MQ=255)

                                                                                                 cACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGa                                                    <  2:939098/150‑1 (MQ=40)

                                                                                                 cACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGa                                                    <  2:2001217/150‑1 (MQ=40)

                                                                                                  aCCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGAt                                                   >  2:1315850/1‑150 (MQ=40)

                                                                                                  aCCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGAt                                                   <  1:3044314/150‑1 (MQ=40)

                                                                                                  aCCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGAt                                                   <  1:454701/150‑1 (MQ=40)

                                                                                                     tCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGGTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCCCCACCCACCAGGTTTAACATCGAAGACAGagggt                                                >  2:1845052/1‑146 (MQ=14)

                                                                                                     tCCCTTACAAGGAGGGGGTCGGCGGTTCGATACCGAAATCACCCACCAC‑CGGGACGTTAGCTCAGTAGGTAGAGCAGATCACTTTTAAACAAATCGTAGAAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGAGGtt                                                >  2:1740567/1‑150 (MQ=14)

                                                                                                       ccTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGGTGGTTGAGCAGTTGACTTTTTATGAACTGGTCGCAGGTTAGAATCCAGCACGACCCACCAGGTTTAACCTCCAAGCCCGGTGttca                                              >  1:460204/1‑148 (MQ=11)

                                                                                                        cTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAg                                             >  2:1301270/1‑150 (MQ=40)

                                                                                                           aCAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCgt                                          <  1:287130/150‑1 (MQ=40)

                                                                                                               ggaggGGGTCGGCGGTTCGATCCCGTCCTCACCCACCAA‑CGGGTCGTTAGCTCCGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAg                                      <  1:726194/150‑1 (MQ=255)

                                                                                                                gaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCGGAGTTGGTAGTGCACTTGACTTTTTTTCAATTGGTCGCAGGTTCGAATACTGCACGACCCCCCAGTTTTAAAATGGAAGACAGATGTTTAGCGGGGAgg                                     >  1:941976/1‑150 (MQ=14)

                                                                                                                  ggggggCGTCGGTCCGACCCGGTCCGCACGCACCAG‑GGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTGAACATCGAAGACAGATGTTAAGCGTGTAGGAt                                   <  1:2377339/150‑1 (MQ=14)

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

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

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

                                                                                                                                 gATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCa                    >  2:2440658/1‑150 (MQ=255)

                                                                                                                                 gATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCa                    >  2:2282791/1‑150 (MQ=255)

                                                                                                                                  aTCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCaa                   <  2:9571/150‑1 (MQ=255)

                                                                                                                                   tCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAAc                  <  1:1215296/150‑1 (MQ=255)

                                                                                                                                     ccGTCATCACCCTCCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTAGCCTTGTAAGCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACgg                <  2:2419149/150‑1 (MQ=14)

                                                                                                                                     ccGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACgg                <  2:2419011/150‑1 (MQ=255)

                                                                                                                                        tcatcaCCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGccc             <  2:2563051/150‑1 (MQ=255)

                                                                                                                                        tcatcaCCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGccc             <  2:1491014/150‑1 (MQ=255)

                                                                                                                                         cagcacccaccAC‑CGGGTCGTTAGCTCAGGTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCg            <  2:2146849/147‑1 (MQ=255)

                                                                                                                                           tcaCCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTa          <  2:2679521/150‑1 (MQ=255)

                                                                                                                                             acccaccAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAgg        >  2:1571582/1‑150 (MQ=255)

                                                                                                                                              cccaccAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATAAATTGTTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAggg       <  2:1008021/150‑1 (MQ=255)

                                                                                                                                                caccac‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGgcg     >  2:2685145/1‑150 (MQ=255)

GGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGC  >  NC_000913/780692‑780988

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