BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	780,841:1	+T	8.1%	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	8.1%	655.3 / 56.9	184	intergenic (+1/‑2)	valT/lysW	tRNA‑Val/tRNA‑Lys
Reads supporting (aligned to +/- strand): ref base . (58/111); new base T (6/7); total (66/118)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.59e-01

Only 100 of 191 total aligned reads displayed.

GGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCG  >  NC_000913/780692‑780989

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

   cGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑Cgg                                                                                                                                                   <  1:3699013/150‑1 (MQ=255)

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

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

       aCGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTc                                                                                                                                                <  2:2434920/150‑1 (MQ=255)

       aCGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCg                                                                                                                                               >  2:3508510/1‑150 (MQ=255)

       aCGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCg                                                                                                                                               >  2:1493208/1‑150 (MQ=255)

       aCGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATAAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCg                                                                                                                                               >  1:3299026/1‑150 (MQ=255)

        cGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGt                                                                                                                                              <  2:3733293/150‑1 (MQ=255)

         gAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGtt                                                                                                                                             >  1:2441680/1‑150 (MQ=255)

         gAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGtt                                                                                                                                             <  1:3493238/150‑1 (MQ=255)

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

             tCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCt                                                                                                                                         <  2:159815/150‑1 (MQ=255)

             tCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCt                                                                                                                                         <  2:608252/150‑1 (MQ=255)

                aGTCATCCTGCACGACCCACCACTAACAGAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGATCAg                                                                                                                                      >  1:1802469/1‑150 (MQ=255)

                 gTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTc                                                                                                                                                >  2:1451101/1‑139 (MQ=255)

                     tCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGt                                                                                                                                 >  2:800963/1‑150 (MQ=255)

                     tCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGGCGTTAGCGCAggcggg                                                                                                                                 >  2:1319691/1‑145 (MQ=255)

                      ccTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTa                                                                                                                                >  1:3677807/1‑150 (MQ=255)

                      ccTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTa                                                                                                                                >  2:2540885/1‑150 (MQ=255)

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

                       cTGCCCGCCCCCCCACGCACCTAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATGAGCTCAGCTGGGAGAGCACCTCCCGTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTag                                                                                                                               <  2:3038244/150‑1 (MQ=14)

                       cTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTag                                                                                                                               >  2:3888541/1‑150 (MQ=255)

                       cTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTag                                                                                                                               <  1:168625/150‑1 (MQ=255)

                        tGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTaga                                                                                                                              >  2:1018268/1‑150 (MQ=255)

                                ccaccaCTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTg                                                                                                                       >  2:1365618/1‑150 (MQ=255)

                                ccaccaCTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGa                                                                                                                      <  2:2422001/150‑1 (MQ=255)

                                ccaccaCTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGa                                                                                                                      <  2:1300586/150‑1 (MQ=255)

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

                                    cacTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACttt                                                                                                                  <  2:648203/150‑1 (MQ=255)

                                     acTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACACGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACttt                                                                                                                  >  2:221639/1‑150 (MQ=255)

                                      ctcaCATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTa                                                                                                                <  1:156459/147‑1 (MQ=255)

                                      cTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACtttt                                                                                                                 <  1:1365618/150‑1 (MQ=255)

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

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

                                       gaaCATAGTTAGTTGTAGTATCCAGCGGAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTaa                                                                                                               <  2:309580/149‑1 (MQ=255)

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

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

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

                                              gttagttGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTa                                                                                                                                            >  2:3394089/1‑115 (MQ=9)

                                               gtagttGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGGCGGTACCTCAGTTGGGAGGGCAGGTGACTTTTAATCAATTgg                                                                                                       >  1:2833877/2‑150 (MQ=11)

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

                                                    tGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCa                                                                                                  <  2:240662/150‑1 (MQ=35)

                                                    tGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCa                                                                                                  <  1:963006/150‑1 (MQ=35)

                                                    tGTAGTATCCAGCGTAGTATCGGGTGATGAGCTCAGCTGGGAGAGCACCTCCCTTACACGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCa                                                                                                  <  2:2553469/150‑1 (MQ=25)

                                                        gtaTCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGGTCAGTTGGTAGTGCAGTTGACTTTTATTGAATTGGTCGCAGGtt                                                                                              >  2:3272324/1‑150 (MQ=255)

                                                        gtaTCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGtt                                                                                              >  2:760709/1‑150 (MQ=35)

                                                           tCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGGTGACTTTTAATCAAATGGTCGCAGGTTAGa                                                                                           >  2:634455/1‑150 (MQ=255)

                                                            cccgCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGaa                                                                                          <  2:3213206/147‑1 (MQ=25)

                                                            ccAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCTCCTCCCTGACTAGGAGGGGGTCGGCGGTTCGATCCCGTCCTCACGCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGaa                                                                                          <  2:1588649/150‑1 (MQ=11)

                                                              aGCGTATTATCGGGTGATTATCTCAGCATTGAGATCTCCTCCCTTACCAGGAGTGGGTAGGCGTTTCGATCCCGTCATCACGCACCAG‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATc                                                                                        <  1:1347661/150‑1 (MQ=11)

                                                              aGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATc                                                                                        >  1:1805587/1‑150 (MQ=35)

                                                                 gtagtaTCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACCAGGAGTGGGTCGGCGGTTCGATCCCGTCATCACCCACCAG‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTg                                                                                     <  2:3669418/150‑1 (MQ=255)

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

                                                                         gggTGATTAGCTCAGCTGTGAGATCTCCTCCCTTACATTGAGTGGGTCGTCGGATCGATCCCCTCATCACGCAACAG‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGAccc                                                                             <  2:3025440/150‑1 (MQ=11)

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

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

                                                                                 aGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTCATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGttt                                                                     <  2:2835240/150‑1 (MQ=18)

                                                                                 aGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGttt                                                                     <  1:2486862/150‑1 (MQ=33)

                                                                                 aGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGttt                                                                     <  1:800963/150‑1 (MQ=33)

                                                                                  gCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCCGGGTCGAATCCTGCACGACCCACCAGtttt                                                                    >  1:9943/1‑150 (MQ=18)

                                                                                  gCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGtttt                                                                    >  1:3339438/1‑150 (MQ=33)

                                                                                  gCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGtt                                                                      <  1:2178274/148‑1 (MQ=33)

                                                                                  gCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGtt                                                                      >  2:2178274/1‑148 (MQ=33)

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

                                                                                    tCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTaa                                                                  >  2:2395373/1‑150 (MQ=33)

                                                                                          gggAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGa                                                            <  1:3555173/150‑1 (MQ=35)

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

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

                                                                                             agagCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑GGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGa                                                         <  2:2722438/150‑1 (MQ=38)

                                                                                               atcaacTCCCTTACATGGAGGGGGTCGGCGGTTCGATCCCGTCATCACGCACCAG‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACa                                                       <  2:1723785/145‑1 (MQ=255)

                                                                                               agcagcTCGCGTAGATGGAGGGGGGTGTCGGTTGGTTCGCGGGTTCACCCTCCTC‑TGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACa                                                       <  1:3782959/150‑1 (MQ=11)

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

                                                                                                     tCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGtt                                                 >  1:1548690/1‑150 (MQ=40)

                                                                                                           aCAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACATTTAATCAAATGGTCGCAGGATCGAATACAGCACGACCCACCAGGTTTAAAATAGAGGACAGGGGTTATGCgg                                           >  1:1398661/1‑149 (MQ=11)

                                                                                                              cggaggGGGTCGTCTGTTCTATCCCGTCTTCACCCACCAC‑CGTGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTa                                        <  1:3943460/149‑1 (MQ=14)

                                                                                                              aggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGTGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTa                                        <  2:2362989/150‑1 (MQ=255)

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

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

                                                                                                               ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAg                                       >  2:2944933/1‑150 (MQ=255)

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

                                                                                                               ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAg                                       >  2:2083670/1‑150 (MQ=255)

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

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

                                                                                                                gaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTTACTTTTAATAAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAAAATCGAAGACAGATGTTAAGCGTGTAgg                                      >  1:891328/1‑150 (MQ=255)

                                                                                                                gaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGGTCGAATACTGCACGACCCACCCCGTTTAACATAGAAGACAGATGTTAAGCGTGTAgg                                      >  1:2873657/1‑150 (MQ=14)

                                                                                                                   ggggTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATa                                   <  1:1018268/150‑1 (MQ=255)

                                                                                                                   ggggTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATa                                   <  1:1018238/150‑1 (MQ=255)

                                                                                                                    gggTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATaa                                  <  1:3757805/150‑1 (MQ=255)

                                                                                                                      gTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACg                                <  2:3118875/150‑1 (MQ=255)

                                                                                                                       tCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATAAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAATCGTGTAGGATAACGt                               >  1:2564203/1‑150 (MQ=255)

                                                                                                                       tCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCGCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGt                               >  1:563724/1‑150 (MQ=255)

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

                                                                                                                              ctCGATCCCGTCATCACCAACCAC‑GGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCa                        <  1:782083/149‑1 (MQ=255)

                                                                                                                               tCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAg                       <  2:563426/150‑1 (MQ=255)

                                                                                                                               tCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGCTAAGCGTGTAGGATAACGTTGCGTCAg                       <  2:3736870/150‑1 (MQ=255)

                                                                                                                                 gATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTTGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATACTGCACGACCCACCAGGTTTAAAATAGAAGACAGATGTTAATCGTGTAGGATAAAGTATCGTCAGCa                     >  1:2340396/1‑150 (MQ=14)

                                                                                                                                  aTCCCGTCAACACCCACCAG‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCaa                    <  1:2082609/150‑1 (MQ=255)

                                                                                                                                   tCCCGTCATCACCCACCCC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAAc                   <  2:2567024/150‑1 (MQ=255)

                                                                                                                                                    ac‑cGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGgcgagcg  >  2:1816069/1‑150 (MQ=255)

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