BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	780,841:1	+T	6.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	6.0%	852.6 / 53.7	232	intergenic (+1/‑2)	valT/lysW	tRNA‑Val/tRNA‑Lys
Reads supporting (aligned to +/- strand): ref base . (78/140); new base T (6/7); total (85/147)
Fisher's exact test for biased strand distribution p-value = 5.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

Only 100 of 237 total aligned reads displayed.

GGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCG  >  NC_000913/780693‑780989

   gagaCGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑Cggg                                                                                                                                                   <  1:2524727/150‑1 (MQ=255)

   gagaCGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑Cggg                                                                                                                                                   <  1:3655372/150‑1 (MQ=255)

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

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

        gAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGt                                                                                                                                               <  2:1999249/150‑1 (MQ=255)

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

         aaGTCGAGACATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTa                                                                                                                                             <  2:3166921/150‑1 (MQ=255)

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

           gTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAg                                                                                                                                            <  1:1013463/150‑1 (MQ=255)

              gAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTTGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTa                                                                                                                                             <  2:1144565/146‑1 (MQ=255)

              gAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTa                                                                                                                                             >  1:1144711/1‑146 (MQ=255)

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

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

                gTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCAAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAgg                                                                                                                                      >  1:328245/1‑149 (MQ=255)

                 tCATCCTGCACGACCCACCACTAAGATAGTTCGTTGGAGTATCCAGCGTAGTAGCGGGTGGTTACCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGGTCGAGCCCGTCATCACCCCCCAC‑CGGGTCGTTAGCTcgggg                                                                                                                                     >  2:5065194/1‑146 (MQ=18)

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

                   aTCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCTTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACTAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTgg                                                                                                                                   <  1:1221539/150‑1 (MQ=255)

                   aTCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTgg                                                                                                                                   <  2:3386609/150‑1 (MQ=255)

                    tCCTGCACGACCCACCACTAACCTAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAggtggt                                                                                                                                  <  2:3894004/150‑1 (MQ=255)

                    tCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGAGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGt                                                                                                                                  >  1:695176/1‑150 (MQ=255)

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

                      cTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTa                                                                                                                                 >  1:2589478/1‑150 (MQ=255)

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

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

                       tGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACACGGGACGTTAGAGCAGtggtaat                                                                                                                                >  2:3434771/1‑144 (MQ=17)

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

                         cACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGc                                                                                                                             >  1:3980006/1‑150 (MQ=255)

                         cACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGCATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGc                                                                                                                             >  1:4278430/1‑150 (MQ=255)

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

                            gACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGt                                                                                                                          <  1:2054305/150‑1 (MQ=255)

                                 accacTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑AGGGTCGTTAGCTCAGGTGGTAGAGCAGTTGaat                                                                                                                     >  1:2498075/1‑148 (MQ=255)

                                   cacTAACTAAGTTATTTGTTGTATCCAGAGTTTTATCGGGTGATTAGCTCAGATGGGTGACCTCATCGCTTACTATGAGGGGGTCGGCGGTTCGATCGGGTCAACACGCACCAG‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACttt                                                                                                                   <  2:3940956/150‑1 (MQ=9)

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

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

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

                                       aaCATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAAt                                                                                                               <  1:95634/150‑1 (MQ=255)

                                        acctagttagtTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCcaccactc                                                                                                                                                      >  2:3260505/4‑109 (MQ=255)

                                          aTAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTTGGAGAGCACCTCCCTTACAAGGAGGGGGTCGTCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCaa                                                                                                            <  2:847531/150‑1 (MQ=255)

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

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

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

                                            cgcgagttGGAGTTTGCCGCGTTGTTTCGGGTGATGAGCTCAGCTGGGAGTTCTACGCGCGTCCGAGGAGGGGGTCGGCGGTTCGATCCGGTCCTCACGAACCAG‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAAtt                                                                                                          <  2:2486770/146‑1 (MQ=2)

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

                                             gttagttGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTCGACTTTTAATCAATTg                                                                                                         <  2:2431324/150‑1 (MQ=255)

                                               tagttGTAGTATCCAGCGTAGTATCGGGTGATTTGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGt                                                                                                       <  1:2236389/150‑1 (MQ=38)

                                                agttGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTc                                                                                                      <  2:13674/150‑1 (MQ=38)

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

                                                  ttGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGc                                                                                                    >  1:2655829/1‑150 (MQ=35)

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

                                                    gtagtaTCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACCAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAg                                                                                                  <  2:4704751/150‑1 (MQ=25)

                                                        taTCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTc                                                                                             <  2:553277/150‑1 (MQ=35)

                                                        taTCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTc                                                                                             <  1:67598/150‑1 (MQ=35)

                                                          tCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGa                                                                                           >  2:1467697/1‑150 (MQ=35)

                                                             aGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATc                                                                                        >  1:1249061/1‑150 (MQ=35)

                                                               cGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCt                                                                                      >  2:3715393/1‑150 (MQ=35)

                                                                  agtaTCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCa                                                                                   <  2:1158274/150‑1 (MQ=33)

                                                                  agtaTCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCa                                                                                   <  2:4454281/150‑1 (MQ=33)

                                                                  agtaTCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCa                                                                                   <  1:4516282/150‑1 (MQ=33)

                                                                  agtaTCGGATGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGATTCGAA‑CCTGcg                                                                                   <  2:4150869/150‑2 (MQ=11)

                                                                    taTCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CCGGGCGTTAGCGCAGGTGGGAGAGCAGTTGCCGTTTAAACCATTGGGCGCAGG‑TGCGAATCCGGCACg                                                                                 >  2:2623158/1‑150 (MQ=11)

                                                                       cGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGAcc                                                                              <  2:2915692/150‑1 (MQ=33)

                                                                        gggTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGAccc                                                                             >  1:4762453/1‑150 (MQ=33)

                                                                         ggTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACcca                                                                            >  2:1900047/1‑150 (MQ=33)

                                                                           tGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACccacc                                                                          >  2:2568766/1‑150 (MQ=33)

                                                                           tGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACccacc                                                                          >  2:126195/1‑150 (MQ=33)

                                                                               tAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACCCACCAGtt                                                                      <  2:3020599/150‑1 (MQ=33)

                                                                                  cTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACCCACCAGTTTTa                                                                   >  2:1391327/1‑150 (MQ=33)

                                                                                   tCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACCCACCAGTTTTaa                                                                  <  2:67416/150‑1 (MQ=33)

                                                                                      gcagggAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACCCACCAGTATTAACAt                                                               <  2:3767542/147‑1 (MQ=18)

                                                                                          ggcgagCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTAATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACCCACCAGTTTTAACATCGaa                                                           <  1:3780143/147‑1 (MQ=25)

                                                                                            agagCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGa                                                         <  2:3980006/150‑1 (MQ=38)

                                                                                            agaaCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGa                                                         <  1:1681650/150‑1 (MQ=38)

                                                                                                 aCCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACCCACCAGTTTTAAAATCGAAGACAGAt                                                    >  1:664054/1‑150 (MQ=255)

                                                                                                 aCCTCCCTAACAACGAGGGGGTCGGCGGTTCGATCCAGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGCTAGAGCACTTGACTTATAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGAt                                                    <  2:253840/150‑1 (MQ=14)

                                                                                                  ccTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATg                                                   >  1:1603579/1‑150 (MQ=40)

                                                                                                  ccTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATg                                                   <  1:2205625/150‑1 (MQ=40)

                                                                                                  ccTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAAATTGTCGCAGG‑TTCGAATACGGCACGACCCACCAGTTTTAACATCGAAGACAGATg                                                   >  1:462022/1‑150 (MQ=255)

                                                                                                  ccTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTATTAATCAATTTGTCGCAGG‑TTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGTCAGATg                                                   >  1:3039132/1‑150 (MQ=255)

                                                                                                      ccTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTaa                                               >  2:4506347/1‑150 (MQ=40)

                                                                                                       cTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTGGAGCAGTTGACGTTTAATCCATTGGTCGCAGG‑GTCGAATCCAGCCCGACCCACCAGTTTTAACATCCAACACAGATGGTAAg                                              >  1:555552/1‑150 (MQ=11)

                                                                                                          aCAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCgt                                           <  2:634685/150‑1 (MQ=40)

                                                                                                          aCAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCgt                                           <  1:864677/150‑1 (MQ=40)

                                                                                                            aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑Cgg                                                                                                                                                    >  2:4965731/1‑44 (MQ=38)

                                                                                                             aggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTa                                        <  2:4323348/150‑1 (MQ=255)

                                                                                                             aggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTa                                        <  2:1785622/150‑1 (MQ=255)

                                                                                                             aggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTa                                        <  2:3200811/150‑1 (MQ=255)

                                                                                                              ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCTCCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAg                                       <  2:2026862/150‑1 (MQ=255)

                                                                                                              ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAg                                       <  2:91016/150‑1 (MQ=255)

                                                                                                                         gcggTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGc                            <  1:4326714/150‑1 (MQ=255)

                                                                                                                         gcggTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATAAATTTGTCGCAGG‑TTCGAATCCTGCACCACCCACCAGATTTAACATCGAAGACAGATGTTAATAGTGTATGATAACGTTGc                            >  1:2246215/1‑150 (MQ=14)

                                                                                                                            gTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTc                         <  1:1276820/150‑1 (MQ=255)

                                                                                                                               cGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCGCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACCCACCAGTTTTAACATAGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGc                      >  1:3160570/1‑150 (MQ=255)

                                                                                                                                 aTCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCaa                    >  2:2107109/1‑150 (MQ=255)

                                                                                                                                 aTCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGGTGGTAGAGCAGTTGACTTGTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACCCCCCCGGGTTAACATCGAAGACCGAGGTGAAGCGGGTAGGATAACGGTGCGgcagcac                    >  2:3597412/1‑149 (MQ=14)

                                                                                                                                  tCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAAc                   <  1:3219404/150‑1 (MQ=255)

                                                                                                                                   cccGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACg                  <  1:2972932/150‑1 (MQ=255)

                                                                                                                                    ccGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑tt                                                                                              >  2:3455303/1‑73 (MQ=38)

                                                                                                                                    ccGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑tt                                                                                              <  1:3455303/73‑1 (MQ=38)

                                                                                                                                          ccacccaccAG‑GGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTa           <  2:2265943/149‑1 (MQ=255)

                                                                                                                                               caccac‑CGGGTCGTTAGCTCAGTTGGTAGAGCATTTGACATTTAATCAATTGGTCGCAGG‑TTCGAAACCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGgcg      <  2:4278430/150‑1 (MQ=255)

GGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGG‑TTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCG  >  NC_000913/780693‑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: `‑`