BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,521,079	Δ1 bp	noncoding (29/76 nt)	valX →	tRNA‑Val

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,521,079	0	T	.	99.0%	386.3 / ‑1.7	96	noncoding (29/76 nt)	valX	tRNA‑Val
Reads supporting (aligned to +/- strand): ref base T (0/1); new base . (56/39); total (56/40)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.29e-01

CAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTT  >  NC_000913/2520943‑2521206

cAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTccc                                                                                                                               <  2:1602541/140‑1 (MQ=255)

    cgggCGAGCTCCTCCCCTAAAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTGCTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGTACC‑CCCTTACa                                                                                                                          <  2:345795/139‑1 (MQ=37)

        agagCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACcccc                                                                                                                                >  2:58953/1‑131 (MQ=255)

            cACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGGGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑gggg                                                                                                                 >  2:1541817/1‑140 (MQ=255)

              ccTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGt                                                                                                               <  1:658326/140‑1 (MQ=255)

                      aCAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGGTTAGCTCAGCTGGGAGAGCACC‑CCCCTACAAGGA‑GGGGGCcggcgggt                                                                                                       >  1:743020/1‑138 (MQ=38)

                          ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGAt                                                                                                   <  2:658764/140‑1 (MQ=255)

                            aggGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTTCGCCGGGTAGCAAGAAATTGAGAAAAGGGGGATTAGCTCCGCTGGGAGAGCACC‑CCCTTACAAGGC‑CGGGGGCGGCGGTTCGGTcc                                                                                                 >  2:465515/1‑140 (MQ=37)

                            aggGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATcc                                                                                                 >  1:540302/1‑140 (MQ=255)

                               gggTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGt                                                                                              <  1:1126777/140‑1 (MQ=255)

                                 gTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGtca                                                                                            <  2:540302/140‑1 (MQ=255)

                                 gTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCCTACAAGGA‑GGGGGTCGGCGGTTCGATCCCCtcc                                                                                            >  1:1062236/1‑139 (MQ=255)

                                  tCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGtcat                                                                                           <  2:791364/140‑1 (MQ=255)

                                  tCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGGGATTAGCTCAGCTGGGAGAGCACC‑CCCCTACAACGA‑GGGGGTCGGCGGATCTGTCCCGtccc                                                                                           >  2:395435/1‑138 (MQ=37)

                                    ggcggTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGtcatca                                                                                         <  1:813412/140‑1 (MQ=255)

                                      cggTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAATGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGGCGGCGGTTCGAACCCCTCATcacc                                                                                       >  2:7722/1‑140 (MQ=37)

                                         ttCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATcacccac                                                                                    >  2:1608281/1‑140 (MQ=255)

                                          tCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATcacccacc                                                                                   <  2:1406039/140‑1 (MQ=255)

                                          tCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATcacccacc                                                                                   <  1:7722/140‑1 (MQ=255)

                                                   tcatcaCCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCCCTTTCTCg                                                                          >  2:183465/1‑140 (MQ=255)

                                                   tcatcaCCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCg                                                                          <  1:1608281/140‑1 (MQ=255)

                                                   tcatcaCCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGAACCCGTCATCACCCACCACTTTCTCg                                                                          >  2:1197442/1‑140 (MQ=255)

                                                                 ctactTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCt                                                            >  2:799695/1‑140 (MQ=255)

                                                                       tATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTTGATCCCGTCATCCCCCCCCCCCTTCTCGCCCTTTAAATTTTTTGTaaa                                                      >  2:398072/1‑140 (MQ=21)

                                                                       tATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCACCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTaaa                                                      >  2:514862/1‑140 (MQ=255)

                                                                        aTGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCATCTAAATTTTTTGTaaaa                                                     >  2:288553/1‑140 (MQ=255)

                                                                        aTGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTaaaa                                                     >  2:767217/1‑140 (MQ=255)

                                                                        aTGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTaaaa                                                     >  1:947836/1‑140 (MQ=255)

                                                                        aTGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTaaaa                                                     >  2:1217265/1‑140 (MQ=255)

                                                                        aTGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTaaaa                                                     >  1:1712725/1‑140 (MQ=255)

                                                                        aTGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCATTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTaaaa                                                     >  2:1515479/1‑140 (MQ=255)

                                                                            aGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAAtgt                                                 >  1:324250/1‑140 (MQ=255)

                                                                            aGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAAtgt                                                 >  2:1471054/1‑140 (MQ=255)

                                                                            aGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGCCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAAtgt                                                 >  2:501589/1‑140 (MQ=255)

                                                                             ggccccgccgTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAAtgtg                                                <  2:1242168/136‑1 (MQ=255)

                                                                             gTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAAtgtg                                                <  1:799695/140‑1 (MQ=255)

                                                                             gTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAAtgtg                                                <  1:767217/140‑1 (MQ=255)

                                                                             gTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAAtgtg                                                <  2:1345470/140‑1 (MQ=255)

                                                                             gTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAAtgtg                                                <  1:1217265/140‑1 (MQ=255)

                                                                             gTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAAtgtg                                                <  1:121785/140‑1 (MQ=255)

                                                                             gTCTCCGCAGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAAtgtg                                                <  1:1629965/140‑1 (MQ=255)

                                                                              tctcCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCCCTTTCTCGCCAGCTAAATTTCTTGTaaaaaattgg                                               >  2:690681/1‑135 (MQ=255)

                                                                              tctcCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGa                                               >  1:641793/1‑140 (MQ=255)

                                                                              tctcCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGa                                               >  2:1470539/1‑140 (MQ=255)

                                                                              tctcCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGa                                               >  1:1457438/1‑140 (MQ=255)

                                                                              tctcCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGa                                               >  1:959449/1‑140 (MQ=255)

                                                                              tctcCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGa                                               >  1:956284/1‑140 (MQ=255)

                                                                              tctcCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGa                                               >  1:263518/1‑140 (MQ=255)

                                                                              tctcCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGa                                               >  1:596615/1‑140 (MQ=255)

                                                                                tcCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAg                                             >  1:850940/1‑140 (MQ=255)

                                                                                  cgccgTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTa                                           >  1:446607/1‑140 (MQ=255)

                                                                                  cgccgTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCCTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTa                                           >  1:1558595/1‑140 (MQ=255)

                                                                                  cgccgTGTAGCAAGAAATTGAGAAGGGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTa                                           <  2:1637371/140‑1 (MQ=255)

                                                                                   gccgTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTTTTGTAAAAATGTGAAGTAc                                          >  2:459209/1‑140 (MQ=255)

                                                                                   gccgTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTAc                                          >  1:869470/1‑140 (MQ=255)

                                                                                   gccgTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTAc                                          >  1:1228104/1‑140 (MQ=255)

                                                                                   gccgTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTAc                                          >  2:220019/1‑140 (MQ=255)

                                                                                   gccgTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGGCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCCGCTAAATTTCTTGTAAAAATGTGaacttc                                          >  2:100170/1‑136 (MQ=21)

                                                                                     cgTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCg                                        >  1:344449/1‑140 (MQ=255)

                                                                                     cgTATAGCAATAAACCGAGCAGGGGGTGTTTAGCGCAGCTGGGAGGGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCg                                        <  2:1525026/140‑1 (MQ=21)

                                                                                      gtgtAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGa                                       <  1:1533248/140‑1 (MQ=255)

                                                                                       tgtAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGaa                                      >  1:1651230/1‑140 (MQ=255)

                                                                                       tgtAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGaa                                      <  1:690681/140‑1 (MQ=255)

                                                                                         tAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGt                                    <  1:465515/140‑1 (MQ=255)

                                                                                            cAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTggg                                 >  1:1495974/1‑140 (MQ=255)

                                                                                               gAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGa                              <  2:344449/140‑1 (MQ=255)

                                                                                                 aaTTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTTTTGTAAAAATGTGAATTACACAAATGGGGGAtt                            >  2:387656/1‑140 (MQ=18)

                                                                                                 aaTTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGGGAtt                            >  1:1287802/1‑140 (MQ=255)

                                                                                                     gagaAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCt                        <  1:1197442/140‑1 (MQ=255)

                                                                                                     gagaAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCt                        <  1:183465/140‑1 (MQ=255)

                                                                                                      agaAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTc                       >  2:711130/1‑140 (MQ=255)

                                                                                                      agaAGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGGCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTTTCAAAGAGGGGGATTAggtt                       >  2:1599661/1‑137 (MQ=18)

                                                                                                         aGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGc                    <  1:220019/140‑1 (MQ=255)

                                                                                                         aGTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGGCGGCGGTTCGGTCCCCTCCACCCCCCCCACTTTCTCACCAGCTAAATTTCTTGTAAAAATGGGAAGTACCGCAGTGGGGGATTGACTCAGc                    >  2:919446/1‑140 (MQ=18)

                                                                                                          gTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTTTTGTAAAAAGGTGAAGTTCCGAAAAGTGGGATTAGGTCCGCt                   >  2:796906/1‑140 (MQ=18)

                                                                                                          gTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCt                   <  1:1541817/140‑1 (MQ=255)

                                                                                                          gTGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCt                   <  2:596615/140‑1 (MQ=255)

                                                                                                           tGGGTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTg                  >  1:488492/1‑140 (MQ=255)

                                                                                                             ggTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTTTTGTAAAAATGTGAAGTACCGAAAAGGGGGATTAGGTCAGCTggg                >  2:1293386/1‑140 (MQ=21)

                                                                                                             ggTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTggg                >  1:607541/1‑140 (MQ=255)

                                                                                                             ggTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTggg                >  1:691401/1‑140 (MQ=255)

                                                                                                             ggTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTggg                >  1:1715249/1‑140 (MQ=255)

                                                                                                             ggTGATTAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTggg                >  1:720106/1‑140 (MQ=255)

                                                                                                                  ttAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGc           <  2:1715249/140‑1 (MQ=255)

                                                                                                                  ttAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGc           <  1:951051/140‑1 (MQ=255)

                                                                                                                  ttAGCCCAGCTGGGGAAGC‑CCCCCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCCCTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGc           <  2:923224/140‑1 (MQ=37)

                                                                                                                   tAGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCa          <  1:395435/140‑1 (MQ=255)

                                                                                                                    aGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTTCCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCAc         <  2:641793/140‑1 (MQ=255)

                                                                                                                    aGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCAc         <  2:1712725/140‑1 (MQ=255)

                                                                                                                    aGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCAc         <  1:1471054/140‑1 (MQ=255)

                                                                                                                    aGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCAc         <  2:1495974/140‑1 (MQ=255)

                                                                                                                    aGCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCAc         <  2:1228104/140‑1 (MQ=255)

                                                                                                                     gCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCAcc        <  2:1457438/140‑1 (MQ=255)

                                                                                                                     gCTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCAcc        <  1:514862/140‑1 (MQ=255)

                                                                                                                      cTCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCt       <  1:711130/140‑1 (MQ=255)

                                                                                                                      cTCAGCTGGGAGAGCACC‑CCATTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCt       <  1:1515479/140‑1 (MQ=255)

                                                                                                                       tCAGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCCTCCTCACCCACCCCTTTCTCGCCGGCTAAATTTTTTTTAAAAATGTGAACTACCGAAGAGGGGGATTAGCTCAGCTGGGAGGGCACCTc      >  2:1475984/1‑140 (MQ=21)

                                                                                                                         aGCTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTATCGAAAAGGGGGATTAGGTCAGCAGGGAGGGCACCTccc    >  2:389345/1‑140 (MQ=25)

                                                                                                                           cTGGGAGAGCACC‑CCCTTACAAGGAGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTTTTGTAAAAATGTGAAGTATCGAAATGGGGGATTAGCTCATCTGGGAGAGCACCTCCCt   >  2:1220020/1‑140 (MQ=25)

                                                                                                                           cTGGGAGAGCACC‑CCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCCGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAATGGGGGATTAGCTCAGCTGGGAGAGCACCTCCCtt  >  1:1756881/1‑140 (MQ=37)

CAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGA‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTT  >  NC_000913/2520943‑2521206

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 15 ≤` `ATCG/ATCG` `< 22 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`