BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	781,254	A→T	6.2%	intergenic (+32/‑115)	lysY → / → lysZ	tRNA‑Lys/tRNA‑Lys

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	781,254	0	A	T	6.2%	1207.2 / 32.9	375	intergenic (+32/‑115)	lysY/lysZ	tRNA‑Lys/tRNA‑Lys
Reads supporting (aligned to +/- strand): ref base A (193/155); new base T (13/14); total (206/169)
Fisher's exact test for biased strand distribution p-value = 5.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 380 total aligned reads displayed.

GGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTT  >  NC_000913/781106‑781403

ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTGAGGTGGTAGTGCAGGTGACTTTTAATCAATTGGTCGCGGGTTCGAATCCTGCACGTCCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGaa                                                                                                                                                      >  1:7793279/1‑150 (MQ=14)

    ggggTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGAt                                                                                                                                                  <  2:7258750/150‑1 (MQ=255)

        tCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATAAATTGGTAGCAGGTTCGAATACTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATaaag                                                                                                                                              >  1:5718257/1‑148 (MQ=255)

          ggcggTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATAAATTGGTCGCAGGTTCGAATCCTGCTCGACCCACCAGGTTTAACTTCAAACTCTGTTGTTAAGCGTGAAGGATACCGtt                                                                                                                                            >  2:2521366/1‑150 (MQ=14)

           gcggTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACAACAAACTCAGATGTTAAGCGTGAAGGATAACGTTg                                                                                                                                           >  1:8371989/1‑150 (MQ=255)

                  gATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGc                                                                                                                                    >  2:582880/1‑150 (MQ=255)

                  gATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGc                                                                                                                                    >  1:6315618/1‑150 (MQ=255)

                  gATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAAATCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGc                                                                                                                                    >  1:4714011/1‑150 (MQ=255)

                   aTCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAAAATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCa                                                                                                                                   >  1:1864969/1‑150 (MQ=255)

                    tCCCGTCATCACCCACGACTCGGGTCGTTATCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCaa                                                                                                                                  <  1:6414508/150‑1 (MQ=255)

                     cccGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATAAATTGGTCGCCGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAAc                                                                                                                                 >  1:1499809/1‑150 (MQ=255)

                      ccGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACg                                                                                                                                >  2:1965132/1‑150 (MQ=255)

                       cGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTTGTAGAGCAGTTGACTTTTAATAAATTTGTAGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAAAATAAAACTCAGATGTTAATCGTGAAGGATAAAGTTGCGCCAGCAACgg                                                                                                                               >  1:3371341/1‑150 (MQ=14)

                        gTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATGAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTTAAATCTATCTGAGATGTTAATCGTGAAGGATAACGTTTCGCCAGCAACGGc                                                                                                                              >  1:3214478/1‑150 (MQ=14)

                               cccaccACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGGAgg                                                                                                                       >  1:4855324/1‑150 (MQ=255)

                                ccaccaCTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAggg                                                                                                                      >  2:1779993/1‑150 (MQ=255)

                                ccaccaCTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAAATGGTCGCAGGTTCGAATACTGCACGACCCACCAGTTTTAAAAACAAAAACAGATGTTAAGAGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAggg                                                                                                                      >  1:6954066/1‑150 (MQ=14)

                                 caccacTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGgc                                                                                                                     <  2:8416050/150‑1 (MQ=255)

                                            tCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAg                                                                                                                                                     >  2:4558144/1‑107 (MQ=255)

                                            tCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGagcgaagcga                                                                                                          >  1:4571337/1‑150 (MQ=255)

                                             cGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAg                                                                                                         <  2:143297/150‑1 (MQ=255)

                                               ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACACAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTc                                                                                                       >  1:1632039/1‑150 (MQ=37)

                                                      aGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCCAACTCAGATGTTAAGCGTGTAGGATAACGTTGCGACAGCAACGGAGCGGATGGCGCTCGAGGCGGGTCATCCTga                                                                                                >  2:6414508/1‑149 (MQ=11)

                                                       gTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCa                                                                                               <  2:4386013/150‑1 (MQ=255)

                                                       gTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCa                                                                                               <  2:7612228/150‑1 (MQ=255)

                                                       gTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGACAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCa                                                                                               >  2:4916390/1‑150 (MQ=255)

                                                          ggTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGa                                                                                            <  2:7081697/150‑1 (MQ=255)

                                                            tAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATACCGTTGCGCCCGCAACGGCCCGTGGGGCGAGCGGAGCGAGTCATCCGGCACGGcc                                                                                          >  2:275657/1‑150 (MQ=14)

                                                            tAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGAcc                                                                                          <  2:3507428/150‑1 (MQ=255)

                                                              gagcttttGACTTTTAAAAATTTGGGCGCTGTTTCGAACCCTCCACCACCCACCATTTTTAAAAACAAACTCATATTTTATGGGTTAAGTATAAGGTGGCGCCAGAAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACcca                                                                                        <  1:3960952/144‑1 (MQ=9)

                                                                    ttGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATc                                                                                  >  1:2552465/1‑150 (MQ=255)

                                                                     tGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCt                                                                                 >  1:4048516/1‑150 (MQ=255)

                                                                         ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTaaa                                                                             >  1:1560265/1‑150 (MQ=255)

                                                                         ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTaaa                                                                             >  1:5608320/1‑150 (MQ=255)

                                                                          tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAATCTCTGATGTTAAGCGTGATGTATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGATTCATCCTGCACGACCCACCAATCTTaaat                                                                            >  2:114027/1‑149 (MQ=17)

                                                                          tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAACTTAAAg                                                                            >  1:4229686/1‑150 (MQ=255)

                                                                              aTCAATTGGTCGCAGGTTCGAATCCTGAACGACCCACCACTTTAAACATCAAAATCACATCTTAAGCGACAACTATAACTTAGAGACAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATACTGCACGACCCACCAATCTTAAAGATTg                                                                        >  1:1011927/1‑150 (MQ=14)

                                                                                 aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGTATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGCCATCCAGCACGACCCACCAATCTGAAAGATGGGcc                                                                     >  2:6164886/1‑150 (MQ=17)

                                                                                 aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTTAAGATTggtg                                                                     <  2:5408714/150‑3 (MQ=255)

                                                                                 aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGcc                                                                     >  1:6544602/1‑150 (MQ=255)

                                                                                  tttGGTCGCTGGTTCTAATCCTGCACTACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGccc                                                                    <  1:3928125/149‑1 (MQ=255)

                                                                                  aTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGccc                                                                    >  1:7540316/1‑150 (MQ=255)

                                                                                   ttGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTAAACATCAAACTCAGAAGTTAAGCGTGAAGGATAACGTTGCGCCAGCAAAGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGcccc                                                                   >  2:3062198/1‑150 (MQ=255)

                                                                                       tCGCAGGTTCGAATCCTGCACGACCCACCAGGGTAAACATCAAACTCAGATGTTAACCGGGACGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGACGCGAGTCATCCTGCACGACCCACCAATCATAAAGAGTGGCCCCGAGt                                                               >  1:611072/1‑150 (MQ=17)

                                                                                        cGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCTAAGTCAGATGTTAAGCGTGTAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTaa                                                                              <  2:2613094/134‑1 (MQ=255)

                                                                                         gCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGTGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTaa                                                             >  2:1627824/1‑150 (MQ=255)

                                                                                         gCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTaa                                                             >  2:2837677/1‑150 (MQ=255)

                                                                                         gCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGAAGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCCATCTTAAAGATTGGCCCCGAGaaa                                                             >  1:5204244/1‑150 (MQ=255)

                                                                                          cAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAACACAGATGTTAAGCGTGTAGGATAACGTTGCGACAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTaaa                                                            >  2:3959040/1‑150 (MQ=255)

                                                                                          cAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTaaa                                                            <  2:3113817/150‑1 (MQ=255)

                                                                                               tCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCt                                                       >  2:8247665/1‑150 (MQ=255)

                                                                                                cGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGTATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCtt                                                      <  1:2953182/150‑1 (MQ=255)

                                                                                                cGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCtt                                                      <  1:1200675/150‑1 (MQ=255)

                                                                                                cGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCtt                                                      <  1:4309392/150‑1 (MQ=255)

                                                                                                  aaTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCCACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTc                                                    >  2:6057377/1‑150 (MQ=255)

                                                                                                  aaTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGCGTCATCCTGCACGACACACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTc                                                    >  2:4860853/1‑150 (MQ=255)

                                                                                                  aaTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTCTTAAAGATTGGCCCCGAGTAAAAATCTTTc                                                    >  2:7452535/1‑150 (MQ=255)

                                                                                                  aaTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCATGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTc                                                    >  1:3767027/1‑150 (MQ=255)

                                                                                                  aaTCCTGCACGACCCACCAGTTTAAACGTCGAAAACGGTAGTTAAGGGGTAGGGATAAAGATCCGACAGCAGCCGCCGGTGGGGCGAGGGAAGCGAGACTTCATGCACTACCCGCAAATCTTAAAGATTGGCCCCGGGTAAGTATCTGTc                                                    >  2:905938/1‑150 (MQ=9)

                                                                                                   aTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTct                                                   >  1:5125100/1‑149 (MQ=255)

                                                                                                    tCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTGGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAg                                                  <  2:4781376/150‑1 (MQ=255)

                                                                                                     ccTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAgg                                                 <  2:22900/150‑1 (MQ=255)

                                                                                                      cTGCACGACCCACCAGTTTTAACATCTAACACAGATGTTAAGCGTGTAGGATAACGTTGCGACAGCAACGGCCCGGAGGGCGAGCGAAGCGAGTCCTCATGCACGACCCGCCAATCTTAAAGGTTGGCGCCGAGTAAAAATATTTCAGGt                                                >  2:807033/1‑150 (MQ=17)

                                                                                                        gCACGACCCACCAGTTTTAACTTCAACCTCTGATGTTATGCGTGATGTATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTaa                                              <  1:3464516/150‑1 (MQ=21)

                                                                                                        gCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTaa                                              <  2:3371341/150‑1 (MQ=255)

                                                                                                               ccaccaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTAGCCCCGAGTAAAAAACGTTCAGGTACCACCCgg                                       >  2:1644267/1‑149 (MQ=255)

                                                                                                                caccaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTa                                      >  1:1960249/1‑150 (MQ=255)

                                                                                                                 tccaGTTTTAACTTCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTAt                                     <  1:7860917/149‑1 (MQ=255)

                                                                                                                 accaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTAt                                     <  1:7040693/150‑1 (MQ=255)

                                                                                                                 accaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTAt                                     >  1:6894573/1‑150 (MQ=255)

                                                                                                                 accaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTAt                                     <  1:1779993/150‑1 (MQ=255)

                                                                                                                  ccaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCTGGTAACTCCCGTATg                                    >  1:923505/1‑150 (MQ=255)

                                                                                                                   caGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATgg                                   >  1:5335919/1‑150 (MQ=255)

                                                                                                                   caGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATgg                                   <  1:6755527/150‑1 (MQ=255)

                                                                                                                   caGTTTTAACACCAAACACAGATGTTAAGCGTGTAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATgg                                   <  2:2938773/150‑1 (MQ=255)

                                                                                                                    aGTTTTAACATCCAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATggg                                  <  1:4784506/150‑1 (MQ=255)

                                                                                                                    aGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATggg                                  >  2:4332925/1‑150 (MQ=255)

                                                                                                                     gTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGt                                 <  2:6788371/150‑1 (MQ=255)

                                                                                                                      ttttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGGATGGGTc                                >  2:623129/1‑150 (MQ=255)

                                                                                                                        ttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGt                              <  1:735648/150‑1 (MQ=255)

                                                                                                                         tAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGtt                             >  2:3632656/1‑150 (MQ=255)

                                                                                                                             aTCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCt                         >  2:2234166/1‑150 (MQ=255)

                                                                                                                             aTCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACCACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCCTCTGGGTCGTTAGCt                         >  2:6667603/1‑150 (MQ=255)

                                                                                                                              tCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTc                        >  1:3100737/1‑150 (MQ=255)

                                                                                                                               cAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCa                       <  2:1315681/150‑1 (MQ=255)

                                                                                                                               cAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCa                       <  1:5325287/150‑1 (MQ=255)

                                                                                                                                 aaCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGt                     >  2:3185207/1‑150 (MQ=255)

                                                                                                                                  aCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGGATGGGTCGTTAGCTCAGtg                    >  2:4125616/1‑149 (MQ=255)

                                                                                                                                     cAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGt                 <  1:267059/150‑1 (MQ=255)

                                                                                                                                     cAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGt                 <  1:126811/150‑1 (MQ=255)

                                                                                                                                      aGATGTTAAGCGTGAAGGATACCGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTa                <  1:8172570/150‑1 (MQ=255)

                                                                                                                                      aGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTa                <  1:7994663/150‑1 (MQ=255)

                                                                                                                                      aGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTa                >  1:8066742/1‑150 (MQ=255)

                                                                                                                                       gATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGATTCATCCTGAACGACCCACCAATCTTAAAGATAGTACCCGACTAAAAATCTTTCAGGTAACACCCTTATGGGTCGTTAGCTAAGAAGGTag               >  2:5606054/1‑150 (MQ=21)

                                                                                                                                       gATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTag               >  1:5865203/1‑150 (MQ=255)

                                                                                                                                       gATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTag               >  1:4944523/1‑150 (MQ=255)

                                                                                                                                       gATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTag               >  2:3749845/1‑150 (MQ=255)

                                                                                                                                                 gTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGAc     <  2:2052166/150‑1 (MQ=255)

                                                                                                                                                   gAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGAATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACtt   >  1:2295776/1‑150 (MQ=255)

                                                                                                                                                    aaGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACttt  <  2:8000385/150‑1 (MQ=255)

GGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTT  >  NC_000913/781106‑781403

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