BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	781,238	T→A	5.6%	intergenic (+16/‑131)	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,238	0	T	A	5.6%	867.7 / 14.2	274	intergenic (+16/‑131)	lysY/lysZ	tRNA‑Lys/tRNA‑Lys
Reads supporting (aligned to +/- strand): ref base T (143/106); new base A (15/6); total (162/112)
Fisher's exact test for biased strand distribution p-value = 2.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 278 total aligned reads displayed.

AGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGG  >  NC_000913/781090‑781387

agCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAAAATCAAACta                                                                                                                                                      >  2:2167925/1‑149 (MQ=255)

      tCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAATCACAGATGt                                                                                                                                                >  1:3436141/1‑150 (MQ=255)

      tCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGt                                                                                                                                                >  2:1845090/1‑150 (MQ=255)

           tACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGc                                                                                                                                           >  1:1880234/1‑150 (MQ=255)

            aCAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACACAGATGTTAAGCg                                                                                                                                          >  2:544354/1‑150 (MQ=37)

              aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTg                                                                                                                                        >  2:4066039/1‑150 (MQ=255)

              aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTg                                                                                                                                        >  1:4922316/1‑150 (MQ=255)

              aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGGTGGTAGAGCAGTTGACTTTTAATCAAATGGTCGCAGGTTCGAATCCTGCACGACCCACCCGTTTTAACATCAAACACAGATGTTAAGCggg                                                                                                                                        >  1:4297223/1‑148 (MQ=255)

                ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGaa                                                                                                                                      >  2:4938433/1‑150 (MQ=255)

                 gaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAg                                                                                                                                     <  1:5032587/150‑1 (MQ=255)

                    ggggTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAAAATCAAACTCAGATGTTAAGCGTGAAGGAt                                                                                                                                  >  1:1834117/1‑150 (MQ=255)

                     gggTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATAAATTGGTCGCAGGTTCGAAACCTGCACGACCCACCACTTTTAACAACAAAAACAGATGTTAAGCGTGAAGGATa                                                                                                                                 >  1:1672634/1‑150 (MQ=11)

                      gttCGGCGTTTCTACCCCTTAATAACCAACCA‑ACGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAAc                                                                                                                               <  2:4385157/148‑1 (MQ=11)

                      ggTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATaa                                                                                                                                <  2:5488632/150‑1 (MQ=255)

                         cggcggTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAAATAAGATGTTAAGCGTGAAGGATAACGt                                                                                                                             >  1:1249383/1‑150 (MQ=255)

                             ggTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGACGGATAACGTTgcg                                                                                                                         >  1:2605484/1‑150 (MQ=255)

                             ggTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTgcg                                                                                                                         >  2:5168011/1‑150 (MQ=255)

                                tCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCa                                                                                                                      <  1:2450559/150‑1 (MQ=255)

                                 cGAACCCGCCACCACCCACCACTCGGGTCGTTAGCTCAGTTGGTTGAGCATTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAg                                                                                                                     <  2:21161/150‑1 (MQ=14)

                                  gATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACGCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGc                                                                                                                    >  1:3092423/1‑150 (MQ=255)

                                   aTCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATACTGCACGACCCACCAGTTTTAAAATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCa                                                                                                                   >  1:2908678/1‑150 (MQ=255)

                                   aTCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATACTGCACGACCCACCAGTTTTAAAATCAAACACAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCa                                                                                                                   >  1:2553867/1‑150 (MQ=255)

                                    tCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATAAAACTCAGATGTTAAGCGTGTAGGATAAAGTTGCGCCAGCaa                                                                                                                  >  1:5747255/1‑150 (MQ=255)

                                    tCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCAGCACCAACCACCACTTTTAACATCAAACAAAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGcca                                                                                                                  >  1:2692416/1‑148 (MQ=14)

                                      ccGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCATTTTTAACTTCAAACTCTGTTTTTAAGCGTGTAGTTTAACGTTGCTCCAGAAACg                                                                                                                >  2:3033653/1‑150 (MQ=14)

                                       cGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACgg                                                                                                               >  2:583370/1‑150 (MQ=255)

                                          catcaCCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAAATGGTCGCAGGTTCGAATACGGCACGACCCACCACTTTTAACGTCAAACGCAGATGTTAAGCGTGAAGGATAACGTTGGGCCAGCCACGGccc                                                                                                            >  1:2903462/1‑150 (MQ=14)

                                              acccaccACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAg                                                                                                        >  2:2958356/1‑150 (MQ=255)

                                                ccaccaCTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAACCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGCATTAACATCAAACTCAGATGGGAAGCGAGAAGGGTAACGTTGCGCCGGAAACGCCACGGAggg                                                                                                      >  2:2500857/1‑150 (MQ=14)

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

                                                     acTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCTATCTCAGATGTTAAGCGGGTAGGATAACGTTGCGCCAGCGACGGCCCGTAGGgcgagc                                                                                                 >  2:3612807/1‑150 (MQ=11)

                                                               ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTc                                                                                       <  2:3185556/150‑1 (MQ=255)

                                                                  gCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCg                                                                                                                                          <  2:1180246/96‑1 (MQ=40)

                                                                  gCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCg                                                                                                                                          >  1:1180246/1‑96 (MQ=40)

                                                                    tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCATGATTAAAGTCAAACTCAGAATATAACCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGATCGAAGCGAGTCATCCt                                                                                  >  1:2990765/1‑150 (MQ=14)

                                                                    tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGagcgaagcga                                                                                          >  1:4580644/1‑142 (MQ=255)

                                                                    tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGagcgaagcga                                                                                          <  2:4580644/142‑1 (MQ=255)

                                                                    tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCt                                                                                  <  2:1227492/150‑1 (MQ=255)

                                                                    tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTATAACATCAAACTCAGAAGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCt                                                                                  >  1:4666871/1‑150 (MQ=255)

                                                                      aGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGc                                                                                <  2:2338673/150‑1 (MQ=255)

                                                                          ggTAGAGCAGTTGACTTTTAATCAATTGGTCGGAGGTTCGAATCCTGCACGACCCACCAGTTTTTACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGa                                                                            >  1:4455504/1‑150 (MQ=255)

                                                                            tAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCAACCAGCACGAcc                                                                          >  2:2291945/1‑150 (MQ=255)

                                                                                 cAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGAAGATAACCGTGAAGGATAACGTTGCGCCAGCAACGCCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacca                                                                     >  2:1211659/1‑150 (MQ=255)

                                                                                   gTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAt                                                                   <  2:1434716/150‑1 (MQ=255)

                                                                                    ttGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATc                                                                  <  2:5091807/150‑1 (MQ=255)

                                                                                      gACTTTTAATCATTTGGTCGCAGTTTCTAATCCTGCACGACCCACCTGTTTTAACATCAAACTCATTTGTTAAGCGTGATGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCtt                                                                <  1:3988309/150‑1 (MQ=14)

                                                                                      gACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACACAGATGATAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCta                                                                >  2:4802337/1‑149 (MQ=255)

                                                                                         ttttAATCAATTGGTCGCTGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTATGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTaaa                                                             <  1:5768687/150‑1 (MQ=255)

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

                                                                                         ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTaaa                                                             <  2:378450/150‑1 (MQ=255)

                                                                                         ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTaaa                                                             <  2:1131643/150‑1 (MQ=255)

                                                                                          tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAACACAGATGTTAAGCGTGTAGGATAACGTTGCGACAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAg                                                            >  1:3818768/1‑150 (MQ=255)

                                                                                            tAATCAATTGGTCGCTGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGAt                                                          <  1:3338005/150‑1 (MQ=255)

                                                                                               tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATAAAGCACGACCCAACAAAATAAAAGAtttg                                                       >  1:865307/1‑148 (MQ=14)

                                                                                                cAATCGGTCGCAGGTTCGCATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGGAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTTAAGATTggt                                                      <  2:3959084/150‑2 (MQ=255)

                                                                                                 aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGcc                                                     <  2:944834/150‑1 (MQ=255)

                                                                                                       tCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCATACTCAGATGCTAAGGTTTAATGAGAACGTTGCGCAAGCAAAGGCCCGTATGGCGAGCGAGGCGAGTCATCCTGCACGACCCACCAAACTTAAATAATCGCGCCGAGt                                               >  2:3084185/1‑150 (MQ=14)

                                                                                                       tCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGGAGCGAGTCATCCGGCACGACCCACCGATCTTAACGAGTGGCCCCGAGt                                               >  2:4702501/1‑150 (MQ=17)

                                                                                                             gTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAAt                                         >  1:5586034/1‑150 (MQ=255)

                                                                                                              ttCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTAAAGCGAGAAGGATAACGTTCCGCAAGCAACGCCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGTTTGGCCCCGAGTAAACATc                                        >  1:4107877/1‑150 (MQ=18)

                                                                                                               tCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCt                                       <  2:3388672/150‑1 (MQ=255)

                                                                                                                  aaTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAAc                                                                                                                 >  1:593640/1‑73 (MQ=255)

                                                                                                                  aaTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTc                                    >  1:1541203/1‑150 (MQ=255)

                                                                                                                  aaTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTc                                    >  1:5800714/1‑150 (MQ=255)

                                                                                                                   aTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGTAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCa                                   <  2:3818768/150‑1 (MQ=255)

                                                                                                                   aTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCa                                   <  2:1641090/150‑1 (MQ=255)

                                                                                                                    tCCTTCACGACCAACCAGTTTTAACATCAAACTCTGATGTTAAGCGTGAAGTATAAGTTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAg                                  <  1:583370/150‑1 (MQ=21)

                                                                                                                    tCCTGCTCGACCCACCAGTTTTAACATCAAACTCATATTTTTCGCGTGATTGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAg                                  <  1:2674851/150‑1 (MQ=21)

                                                                                                                    tCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAg                                  >  2:1544743/1‑150 (MQ=255)

                                                                                                                      cTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGt                                >  1:3383513/1‑150 (MQ=255)

                                                                                                                       tGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTa                               >  1:1352036/1‑150 (MQ=255)

                                                                                                                       tGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTa                               >  1:2194816/1‑150 (MQ=255)

                                                                                                                           cGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATATCGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAacac                           <  1:1463257/150‑1 (MQ=255)

                                                                                                                             aCCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACAccc                         <  2:3696984/150‑1 (MQ=255)

                                                                                                                                 accaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGTCCCCGAGTAAAAATCTTTCAGGTAACACCCGTAt                     >  2:2737120/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                                                                                                       tttAACATCCAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCg               <  1:4283998/150‑1 (MQ=255)

                                                                                                                                        ttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACTACCGAATGGGTCGt              >  2:5328701/1‑150 (MQ=255)

                                                                                                                                        tgaaCATCAACCTCAGCTGATAAGCGTGCAGGATAACGTTGCGCCAGCAACGGCCCTTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCCATCTTACAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGt              <  1:966464/148‑1 (MQ=25)

                                                                                                                                          aacaacAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTa            <  1:1051349/150‑1 (MQ=255)

                                                                                                                                          aaCATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTa            <  2:734090/150‑1 (MQ=18)

                                                                                                                                           acaacaAACACAGATGATAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAg           <  2:5612104/150‑1 (MQ=255)

                                                                                                                                           aCATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAg           <  1:644962/150‑1 (MQ=255)

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

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

                                                                                                                                               ccaaCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCa       <  1:4727407/148‑1 (MQ=255)

                                                                                                                                               cAAACTCAGATGTTAAGCGTTAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCa       <  2:5959307/150‑1 (MQ=255)

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

                                                                                                                                                 aaCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGt     <  2:1327913/150‑1 (MQ=255)

                                                                                                                                                 aaCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGt     <  1:2739879/150‑1 (MQ=255)

                                                                                                                                                   cTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTg   <  2:2762074/150‑1 (MQ=255)

                                                                                                                                                    tCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACTCCCGTATGGGTCGTTAGCTCAGTTgg  >  2:2016349/1‑150 (MQ=255)

AGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGG  >  NC_000913/781090‑781387

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