BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,851,205	G→C	7.3%	G116R (GGG→CGG) ‡	ansA →	L‑asparaginase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,851,205	0	G	C	7.3%	412.8 / 17.7	124	G116R (GGG→CGG) ‡	ansA	L‑asparaginase 1
Reads supporting (aligned to +/- strand): ref base G (58/57); new base C (3/6); total (61/63)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

Only 100 of 125 total aligned reads displayed.

AAGACTGGCAGCATATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCC  >  NC_000913/1851056‑1851352

aaGACTGGCAGCATATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAg                                                                                                                                                     <  2:1509628/150‑1 (MQ=255)

    cTGGCAGCATATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTc                                                                                                                                                 >  2:748515/1‑150 (MQ=255)

    cTGGCACCATATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATTGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTc                                                                                                                                                 <  2:441124/150‑1 (MQ=255)

      ggCAGCATATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTcac                                                                                                                                               <  1:315276/150‑1 (MQ=255)

       gcagcaTATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTcaca                                                                                                                                              >  2:1599607/1‑150 (MQ=255)

         agcaTATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACCCTGTCACaat                                                                                                                                            <  1:1410743/150‑2 (MQ=255)

         agcaTATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACAGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACaaa                                                                                                                                            >  1:925802/1‑150 (MQ=255)

          gcaTATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACCCTGTCACaatg                                                                                                                                           <  1:404270/150‑3 (MQ=255)

           caTATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATc                                                                                                                                          >  2:1523167/1‑150 (MQ=255)

             tatTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACCCTGTCACaatgacc                                                                                                                                        <  1:669737/150‑6 (MQ=255)

             tatTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACCCTGTCACaatgacc                                                                                                                                        <  1:1392081/150‑6 (MQ=255)

             tatTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGGCATTGTGACCCTGTCACaatgacc                                                                                                                                        <  2:405125/150‑6 (MQ=255)

              atTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACaaa                                                                                                                                            >  2:45362/1‑145 (MQ=255)

              atTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCg                                                                                                                                       >  2:791909/1‑150 (MQ=255)

               ttGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgc                                                                                                                                      >  1:541792/1‑150 (MQ=255)

                tGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTc                                                                                                                                                 <  2:896241/138‑1 (MQ=255)

                tGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTc                                                                                                                                                 >  1:896241/1‑138 (MQ=255)

                tGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgct                                                                                                                                     >  2:928848/1‑150 (MQ=255)

                  cTGAAGATATTAAAGCGCACTCTGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgctgc                                                                                                                                   <  1:1620912/150‑2 (MQ=255)

                  cTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgctgg                                                                                                                                   >  1:432904/1‑150 (MQ=255)

                   tGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgctggc                                                                                                                                  >  2:940097/1‑150 (MQ=255)

                    gAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgctggct                                                                                                                                 >  2:1608291/1‑150 (MQ=255)

                    gAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgctggct                                                                                                                                 >  1:515980/1‑150 (MQ=255)

                         tatTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgctggctg                                                                                                                                >  1:661616/1‑146 (MQ=255)

                         tatTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgctggctg                                                                                                                                <  2:661616/146‑1 (MQ=255)

                         tatTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGtt                                                                                                                            >  2:183455/1‑150 (MQ=255)

                          atTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTa                                                                                                                           >  2:126807/1‑150 (MQ=255)

                           ttAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCAGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACCCTGTCACaatgacc                                                                                                                                        <  1:92291/136‑6 (MQ=255)

                           ttAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCAGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACCCTGTCACaatgacc                                                                                                                                        >  2:92291/1‑131 (MQ=255)

                            tAAAGCGCACTATGACGACTATGATGGTTTTGTCATTATGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACg                                                                                                                         >  2:1053619/1‑150 (MQ=255)

                             aaaGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGc                                                                                                                        >  1:221487/1‑150 (MQ=255)

                              aaGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGct                                                                                                                       <  1:1667623/150‑1 (MQ=255)

                                  gcACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGa                                                                                                                   <  1:734228/150‑1 (MQ=255)

                                    aCTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACg                                                                                                                 <  1:1066261/150‑1 (MQ=255)

                                       aTGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGAc                                                                                                              <  2:1003285/150‑1 (MQ=255)

                                         gacgacTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACaa                                                                                                            >  2:1427724/1‑150 (MQ=255)

                                              cTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTaa                                                                                                       >  1:1163544/1‑150 (MQ=255)

                                                atgatgGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAAta                                                                                                     >  1:42857/1‑149 (MQ=255)

                                                atgatgGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATc                                                                                                     >  1:959438/1‑150 (MQ=255)

                                                atgatgGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATc                                                                                                     <  1:1013799/150‑1 (MQ=255)

                                                   atgGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATctgc                                                                                                  >  2:119352/1‑150 (MQ=255)

                                                   atgGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATctgc                                                                                                  <  1:1082354/150‑1 (MQ=255)

                                                   atgGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATctgc                                                                                                  >  1:244011/1‑150 (MQ=255)

                                                    tgGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATctgct                                                                                                 >  2:110647/1‑150 (MQ=255)

                                                     ggTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATctgctg                                                                                                >  1:1261238/1‑150 (MQ=255)

                                                      gTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGa                                                                                               >  2:700685/1‑150 (MQ=255)

                                                           gTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCg                                                                                          >  2:1561429/1‑150 (MQ=255)

                                                           gTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCg                                                                                          >  2:505035/1‑150 (MQ=255)

                                                            tCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGt                                                                                         >  2:958597/1‑150 (MQ=255)

                                                            tCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGt                                                                                         <  1:156335/150‑1 (MQ=255)

                                                            tCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGt                                                                                         >  1:1435833/1‑150 (MQ=255)

                                                             cATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGtt                                                                                        >  2:934210/1‑150 (MQ=255)

                                                              aTTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTg                                                                                       >  2:456487/1‑150 (MQ=255)

                                                               ttCTGCACGGCACCGACACGATGGCGTATACCGCCTCTTCGCTGTCGTTCATGCTCGAGAATCTCTGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGt                                                                                      <  1:37218/150‑1 (MQ=255)

                                                                tCTGCACGGGACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTa                                                                                     <  2:291063/150‑1 (MQ=255)

                                                                   gCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGt                                                                                  >  2:1146322/1‑150 (MQ=255)

                                                                         cACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTgcggc                                                                            <  1:1453378/150‑1 (MQ=255)

                                                                            cGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGaa                                                                         <  1:742275/150‑1 (MQ=255)

                                                                            cGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGaa                                                                         <  1:1200252/150‑1 (MQ=255)

                                                                             gACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTgcggcg                                                                           >  1:1447403/1‑147 (MQ=255)

                                                                              acacGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAAtt                                                                       <  2:1056497/150‑1 (MQ=255)

                                                                               cacGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTa                                                                      >  2:558875/1‑150 (MQ=255)

                                                                                 cGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATc                                                                    >  1:1662969/1‑150 (MQ=255)

                                                                                  gATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATcc                                                                   <  1:928848/150‑1 (MQ=255)

                                                                                     ggCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGAt                                                                >  2:1624591/1‑150 (MQ=255)

                                                                                           taCCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGa                                                          <  1:1146322/150‑1 (MQ=255)

                                                                                           taCCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGa                                                          >  2:1504692/1‑150 (MQ=255)

                                                                                            aCCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGaa                                                         <  2:541792/150‑1 (MQ=255)

                                                                                             ccGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTTACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAg                                                        <  1:603222/150‑1 (MQ=255)

                                                                                             ccGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAg                                                        <  2:1156494/150‑1 (MQ=255)

                                                                                             ccGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAg                                                        <  2:1645210/150‑1 (MQ=255)

                                                                                                  tctGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACg                                                   >  1:177288/1‑150 (MQ=255)

                                                                                                    tGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGct                                                 <  1:126807/150‑1 (MQ=255)

                                                                                                     gcgcTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGctc                                                >  1:583359/1‑150 (MQ=255)

                                                                                                     gcgcTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGctc                                                >  1:1566158/1‑150 (MQ=255)

                                                                                                           tCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTc                                          <  1:1605877/150‑1 (MQ=255)

                                                                                                            cGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCa                                         >  1:1343333/1‑150 (MQ=255)

                                                                                                             gTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCaa                                        <  2:1163544/150‑1 (MQ=255)

                                                                                                             gTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCaa                                        <  1:940097/150‑1 (MQ=255)

                                                                                                                  tGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAAcc                                   <  1:1029596/150‑1 (MQ=255)

                                                                                                                  tGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAAcc                                   <  2:515980/150‑1 (MQ=255)

                                                                                                                      cGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGAtt                               >  2:532764/1‑150 (MQ=255)

                                                                                                                       gagaATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTg                              <  1:542980/150‑1 (MQ=255)

                                                                                                                         gaATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTTCTGTATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTa                            <  2:921304/150‑1 (MQ=255)

                                                                                                                             ctcGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcgc                        <  2:244011/150‑1 (MQ=255)

                                                                                                                               cGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCgg                      <  1:558875/150‑1 (MQ=255)

                                                                                                                                   aaaCCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAAc                  >  1:1236220/1‑150 (MQ=255)

                                                                                                                                    aaCCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAAcc                 <  1:1445579/150‑1 (MQ=255)

                                                                                                                                     aCCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCg                >  2:444151/1‑150 (MQ=255)

                                                                                                                                      ccGGTCATTGTGACAGGGTCTCAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAGATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGc               <  1:139323/150‑1 (MQ=255)

                                                                                                                                       cGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGAACAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCa              >  1:492183/1‑150 (MQ=255)

                                                                                                                                        gggatttGTGACCCTGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCAc             >  2:1059120/6‑150 (MQ=255)

                                                                                                                                        ggTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCAc             <  2:1261238/150‑1 (MQ=255)

                                                                                                                                         gTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACt            <  1:958597/150‑1 (MQ=255)

                                                                                                                                         gTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACt            <  1:58634/150‑1 (MQ=255)

                                                                                                                                          gatttGTGACCCTGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTa           >  1:1235618/4‑150 (MQ=255)

                                                                                                                                           cATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTAc          <  1:1565968/150‑1 (MQ=255)

                                                                                                                                             ttGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCa        >  2:1252289/1‑150 (MQ=255)

                                                                                                                                               gtgACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCaaa      >  1:1634589/1‑150 (MQ=255)

                                                                                                                                                tgACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAACGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAg     >  1:1012545/1‑150 (MQ=255)

AAGACTGGCAGCATATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCC  >  NC_000913/1851056‑1851352

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`