BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	780,841	1	.	T	7.6%	423.8 / 37.1	118	intergenic (+1/‑2)	valT/lysW	tRNA‑Val/tRNA‑Lys
Reads supporting (aligned to +/- strand): ref base . (41/68); new base T (5/4); total (46/72)
Fisher's exact test for biased strand distribution p-value = 3.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-01

Only 100 of 130 total aligned reads displayed.

GGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAG  >  NC_000913/780693‑780987

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

 gCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAAGATAGTTAGTTTGAGTATCCAGAGCAGTATCGGGTGATAACCTCAGCTGGGAGAGCACCACCCTTACAAGGAGGGTCTCGTCGGTTCGATCCCGACAACACAcaccacac                                                                                                                                                   >  2:2306696/1‑148 (MQ=18)

 gCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATCGGTAGGTGGAGTATCCAGCGTAGGATCGGGGGATTAGCGCAGCTGGGAGAGCACCTCCCTGCCACGGAGCGGGGCGGCCGTTCGATCCCGGCATCACCCacccc‑cc                                                                                                                                                  >  2:1487188/1‑149 (MQ=18)

   tagaCGAAGTCGAGTCATCCTGCACGACACACCTCTTAGATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCcaccactcgg                                                                                                                                                 <  2:1890810/149‑5 (MQ=17)

   gagaCGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCACCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑Cggg                                                                                                                                                >  2:1046605/1‑150 (MQ=255)

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

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

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

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

         caGTCGAGTCATCCTCCACGGCCCCGCACTTTCCTAGTTTGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGtt                                                                                                                                           <  2:1869382/149‑1 (MQ=14)

         aaGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGtt                                                                                                                                           <  2:841053/150‑1 (MQ=255)

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

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

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

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

              gAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCCGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCa                                                                                                                                     <  1:1594385/150‑1 (MQ=255)

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

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

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

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

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

                  cATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGtt                                                                                                                                  >  2:307807/1‑150 (MQ=255)

                    tCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCCGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTGCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGt                                                                                                                               <  1:720002/150‑1 (MQ=255)

                    tCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCCGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTGCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGt                                                                                                                               <  1:521081/150‑1 (MQ=255)

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

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

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

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

                    tCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTGCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGt                                                                                                                               <  1:1224246/150‑1 (MQ=255)

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

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

                            gACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGGTGGTAGAGCAGt                                                                                                                       >  2:452392/1‑150 (MQ=255)

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

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

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

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

                                 accacTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACt                                                                                                                  >  1:145627/1‑150 (MQ=255)

                                    acTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACtttt                                                                                                               <  2:130070/150‑1 (MQ=255)

                                     cTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACCAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGCCTTTTa                                                                                                              <  2:2296728/150‑1 (MQ=255)

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

                                      tAACATAGTTTGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCa                                                                                                                         <  1:1826572/139‑1 (MQ=255)

                                      tAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCa                                                                                                                         >  2:1826572/1‑139 (MQ=255)

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

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

                                      tAACATACTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTaa                                                                                                             <  2:1088866/150‑1 (MQ=255)

                                        aCATAGTTAGTTGTAGTATCCAGCGAAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGTCGGTACGATCCCGTCATCACCCACCAC‑CGGGTCGATAGCTCAGTTGGTAGAGCAGTTGACTTTTAATc                                                                                                           >  1:2141961/1‑150 (MQ=255)

                                         cATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATc                                                                                                           >  1:840179/1‑150 (MQ=21)

                                          aTAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGATTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCaa                                                                                                         <  2:2289394/150‑1 (MQ=255)

                                           tagttagtTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCaa                                                                                                         >  2:2414117/1‑150 (MQ=17)

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

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

                                            agttagttCTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAAtt                                                                                                       <  2:954596/150‑1 (MQ=255)

                                               tagttGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGt                                                                                                    <  2:2411043/150‑1 (MQ=38)

                                                  ttGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTTGGTCAGTTGGTTGAGCAGTTGAGTTTTATTGAATTGGTGGc                                                                                                 >  1:494354/1‑150 (MQ=11)

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

                                                      agtaTCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCCGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGt                                                                                             <  1:1197121/150‑1 (MQ=25)

                                                          tCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGa                                                                                         >  2:916189/1‑150 (MQ=35)

                                                          tCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGa                                                                                         >  2:564196/1‑150 (MQ=35)

                                                            cAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTg                                                                                                      >  1:2217545/1‑135 (MQ=37)

                                                            cAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTg                                                                                                      <  2:2217545/135‑1 (MQ=37)

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

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

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

                                                                    taTCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACg                                                                               >  1:1521913/1‑150 (MQ=33)

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

                                                                            gATTAGCTCCGCTGGGAGAGCACCTCCCTTACAAGGTGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCCGTTGGTAGAGCAGTTGACTTTTAATCATTTGGTCGCAGGTTCGAATCCTGCACGACccacca                                                                       <  1:543905/150‑1 (MQ=255)

                                                                            gATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCCGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACccacca                                                                       <  1:543918/150‑1 (MQ=18)

                                                                              ttAGCTCAGCTGTGAGAGCACCTCCCTTACAAGGAGGGGGTCTGCGGTTCGATCCCGTCTACACCCAACAG‑CGGGACGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGt                                                                     <  2:2541909/150‑1 (MQ=11)

                                                                                   tCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCCGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTaa                                                                <  1:2416438/150‑1 (MQ=18)

                                                                                      gCTGGGAGAGCACCTCCCTTACTAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCAACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACAt                                                             <  2:1360599/150‑1 (MQ=18)

                                                                                                 aCCTCGCTTACTAGGAGTGGGTCGTCGGTTCGATCCCTTCCTCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGAt                                                  <  2:2399920/150‑1 (MQ=14)

                                                                                                 aCCTCCGTTTCCTGGAGGGGGTCGGTGGTCAGATCCGTTCTTCACGCACAAG‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGAt                                                  <  1:1439123/150‑1 (MQ=14)

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

                                                                                                 aCCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTAAGGTGGTTGGGCAGTTTAAGTTTAATAAAATTGTACCAGGTTACGATACGGCACGGACCACCAGTTTTAAAATCCaaacactag                                                  >  2:6216/1‑143 (MQ=2)

                                                                                                 aCCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCGCCGTTGGTAGAGCAGTTGACTTTTAATCAAATTGTCGCAGGGTCGGATCCGGCACGACCCACCAGGTTTAACAACCAAcacaaag                                                  >  2:470338/1‑147 (MQ=11)

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

                                                                                                    tCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGGTTTAACATCGAAGACAGATGtt                                               >  2:388821/1‑150 (MQ=255)

                                                                                                      ccTTACAAGGAGGGGGTCGGCGGTGCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTa                                              <  1:324884/150‑1 (MQ=12)

                                                                                                          aCAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTAAGTTGGTAGAGCAGTTGACTTTTAATAAATTGGTCGCAGGTTCGAATCCAGCACGAACCACCAGTTTTAAAATCGAAGACAGATGTTAAGCgt                                         >  1:2490761/1‑150 (MQ=14)

                                                                                                            aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCgtgt                                       >  2:773911/1‑150 (MQ=255)

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

                                                                                                               gaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATAAAATGGTCGCAGGTTAGAATACAGCACGACCCACCACGTTTTACATAGCAGACAGATGTTAAGCGTGTAgg                                    >  1:895417/1‑150 (MQ=14)

                                                                                                                      tCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACGTTTAATAAAATGGTCGCAGGTTAGAATACAGCACGACCCACCACTTTTAAAAACGAAGAAAAATGTTAAACGTGTATGATAAAGt                             >  1:2189879/1‑150 (MQ=14)

                                                                                                                        ggcggTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTg                           <  2:675860/150‑1 (MQ=255)

                                                                                                                             ttCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCa                      <  1:1866458/150‑1 (MQ=255)

                                                                                                                             ttCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCa                      <  1:564266/150‑1 (MQ=255)

                                                                                                                             ttCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGGTGGTAGAGCAGTTTTCTTTTTATCAAATGGGCGCAGGTTCGAATCCTGCCCGACCCACCAGGTTTAAAATAGAAGACAGATGTTAAGCGTGTAGGCTAACGTTGCGTCa                      >  1:615677/1‑150 (MQ=14)

                                                                                                                             atCGATCCGGTCCCCACGCACCAG‑GGGGTCGTGAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCa                      <  1:61552/149‑1 (MQ=14)

                                                                                                                                gATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGGTGGTAGAGCAGTTGACTTTTAATCAAATGGTCGCAGGTTCGAATACGGCACGACCCACCAGGTTTAAAATAGAAGACAGATGGTAATCGTGTAGGATAACGTTGCGTCAGCa                   >  1:2158542/1‑150 (MQ=14)

                                                                                                                                 aTCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGTCAGATGTTTAGCGGGTGGGATAACGTTGCGTCAGCaa                  >  2:743945/1‑150 (MQ=255)

                                                                                                                                 aTCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCaa                  >  2:758010/1‑150 (MQ=255)

                                                                                                                                          tcaCCCACCAG‑GGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTa         <  2:1917687/150‑1 (MQ=255)

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

                                                                                                                                              ccaccaC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGAACCACCAGTTTTAACATCGAAGACAGATGTTAATCGGGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGgc     >  2:5753/1‑150 (MQ=255)

                                                                                                                                                accac‑cGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGgcga   <  2:2189879/150‑1 (MQ=255)

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

                                                                                                                                                 cacc‑cGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTa                                                                 <  1:2494185/84‑1 (MQ=1)

                                                                                                                                                  ccc‑cGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCg                                                                                          <  1:329213/59‑1 (MQ=14)

                                                                                                                                                  cac‑cGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCg                                                                                          <  1:2142634/61‑1 (MQ=35)

                                                                                                                                                  cac‑cGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCg                                                                                          >  2:329213/1‑61 (MQ=35)

GGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAG  >  NC_000913/780693‑780987

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