BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	781,238	T→A	8.3%	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	8.3%	1436.5 / 47.7	479	intergenic (+16/‑131)	lysY/lysZ	tRNA‑Lys/tRNA‑Lys
Reads supporting (aligned to +/- strand): ref base T (230/193); new base A (37/11); total (272/207)
Fisher's exact test for biased strand distribution p-value = 3.16e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 487 total aligned reads displayed.

AGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGG  >  NC_000913/781090‑781387

agCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGATTTAACATCAAACTc                                                                                                                                                       >  2:3195523/1‑150 (MQ=255)

  cACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACACAg                                                                                                                                                     >  1:2549747/1‑150 (MQ=37)

   aCCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGa                                                                                                                                                    >  1:5161412/1‑150 (MQ=255)

      tCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCCAACTCAGATGt                                                                                                                                                 >  2:7217488/1‑150 (MQ=255)

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

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

      tCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATACTGCACGACCCACCAGTTTTAAAATCAAACTCAGATGt                                                                                                                                                 >  1:1729864/1‑150 (MQ=255)

       cccTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATAGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGaagat                                                                                                                                                >  2:9659650/1‑146 (MQ=255)

        ccTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTa                                                                                                                                               >  1:9377509/1‑150 (MQ=255)

          ttACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTTGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAg                                                                                                                                             <  2:9045542/150‑1 (MQ=255)

                ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGaa                                                                                                                                       <  2:7793725/150‑1 (MQ=255)

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

                ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGaa                                                                                                                                       <  2:3007400/150‑1 (MQ=255)

                ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGAATTAACATCAAACTCAGATGTTAAGCGTGaa                                                                                                                                       >  1:6584079/1‑150 (MQ=255)

                 gaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTATCTCAGTTGGTATAGCAGTTGTATTTTAATGAATTGGTCGCAGGTTCGAATACTGCACGACCCACCATTTTTAAAATCAAAATATGATGTTTATAGTGAAg                                                                                                                                      >  1:4206247/1‑150 (MQ=11)

                  aggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAAACAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCACTTTTAACATCAAACACAGATGTTAAGCGTGAAgg                                                                                                                                     >  1:7074037/1‑150 (MQ=255)

                   gggggTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAAAACAGATGTTAAGCGTGAAGGa                                                                                                                                    >  1:3167873/1‑150 (MQ=255)

                     ggctCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCTTTAGCTCAGCTGGTAGAGCAGTTGACTTATAATCTATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATa                                                                                                                                  <  2:6900187/147‑1 (MQ=14)

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

                      ggTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCACAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGGTTTAACATCAAACACAGATGTTAAGCGTGAAGGATaa                                                                                                                                 >  1:4847996/1‑150 (MQ=255)

                        tCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACg                                                                                                                               >  1:7823205/1‑150 (MQ=255)

                        tCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCACTTTAAACATCCAACACCGAAGTGAACCGTGACGGAAAACg                                                                                                                               >  2:1422859/1‑150 (MQ=11)

                          ggcggTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGAACCAGGTTCGCAGCGGGAGCGACACACCGGGTTTAACATCGTACACCGATGTTCAGCGGGAAGGATAACGtt                                                                                                                             >  2:4885601/1‑150 (MQ=11)

                           gcggTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTg                                                                                                                            >  2:7559554/1‑150 (MQ=255)

                               ttCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACCACCCACCAGTTTTAACATCAAACACAGATGTTAAGCGTGAAGGATAACGTTGCGcc                                                                                                                        >  1:5273619/1‑150 (MQ=255)

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

                                 cGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAg                                                                                                                      <  1:5554993/150‑1 (MQ=255)

                                        gTCATCACCCACCACTCGGGTCGTTAGGTAAGTTGGTATAGCAGTTTACTTTTTATAAATTGGTGGGAGGTTCGAATACTGCACGACCAAACAGTTTTAACATCAAAATAATATTTTTTGGGTGTAGGATAAAGTTTGGGCAGAAACGGc                                                                                                               >  1:9491512/1‑150 (MQ=255)

                                                    cacTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTCAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGGAGGgcgag                                                                                                   >  2:6443412/1‑150 (MQ=255)

                                                        cGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGagcgaa                                                                                               <  2:1078961/150‑1 (MQ=255)

                                                         gggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGagcgaag                                                                                              >  2:5771756/1‑150 (MQ=255)

                                                         gggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGagcgaag                                                                                              <  2:8778588/150‑1 (MQ=255)

                                                            tCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGagcgaag‑cga                                                                                          >  1:5768284/1‑150 (MQ=255)

                                                              gTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGtt                                                                                                                                                >  2:9568420/1‑95 (MQ=40)

                                                              gTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAg                                                                                                                                      >  1:1318190/1‑105 (MQ=255)

                                                                   cTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATcc                                                                                   >  1:9859296/1‑150 (MQ=255)

                                                                    tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCTGTAGATAACAGATAAGGATAACGTAGCGACAGCAACGGCCCGTAGGGCGAGAGAAG‑AGATACATCCt                                                                                  >  2:847144/1‑150 (MQ=14)

                                                                    tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCt                                                                                  >  2:3267961/1‑150 (MQ=255)

                                                                    tCAGTTGGTAGAGCAGTTGACTTTTAATCAAATGGTCGCAGGTTATAATCCTGCACCACCCACAAGTTTTAACATAAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCCAAC‑CGAGACATccg                                                                                  >  1:9003564/1‑149 (MQ=14)

                                                                      aGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGc                                                                                <  2:1816387/150‑1 (MQ=255)

                                                                        ttGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGTTGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGTCGAGTCATCCTGCa                                                                               <  2:9618860/150‑1 (MQ=255)

                                                                        ttGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACAAAGATGTTAAGCGTGAAGCATAACGTTGCGACAGCAACGGCCCGAAGGGCGAGCGAAG‑CGAGTCATACTGCAc                                                                              >  1:4625877/1‑150 (MQ=11)

                                                                            tAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAg                                                                                                                                             <  1:8403046/84‑1 (MQ=255)

                                                                            tAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGAcc                                                                          <  2:786324/150‑1 (MQ=255)

                                                                               agCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACccac                                                                       >  1:1644344/1‑150 (MQ=255)

                                                                                    ttGACTTTTAATCAATTGGTCGCTGGTTCGATCACTTCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGTAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATc                                                                  <  2:1624445/150‑1 (MQ=14)

                                                                                    ttGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTATCTTCATTGTCTGATGTTATGGGTGATGGTTAACGTTGCGCCTGCTAGGGCCGGTAGGGCGAGCGAAG‑CGAGTCTTCCTGCACGACCCACCAAtg                                                                  >  2:7273661/1‑149 (MQ=11)

                                                                                    ttGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAAATCAGAAGATAAGCGTGAAGGATAACGAAGCGCCACCAACGGCCCGAAGGGCGAGCGACG‑CGAGTCATCCTGCACGACCCACCAATc                                                                  >  2:6270747/1‑150 (MQ=14)

                                                                                    atGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑TGAGTCATCCTGCACGACCCACCAATc                                                                  >  2:2501593/2‑150 (MQ=255)

                                                                                       aCTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTa                                                               <  2:9282789/150‑1 (MQ=255)

                                                                                         ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTaaa                                                             >  1:10043726/1‑150 (MQ=255)

                                                                                          tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAg                                                            >  1:2823913/1‑150 (MQ=255)

                                                                                          tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAg                                                            >  2:3625056/1‑150 (MQ=255)

                                                                                          tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTATAACATCAAACTCAGATGTTAAACGTGAAGGATAAAGTTCAGCCAGCAACGGCACGAAGTGCGAGCGAAG‑CGAGACATCATGCACGACCCACCAATCAAAAAg                                                            >  1:9578721/1‑150 (MQ=14)

                                                                                             aaTCTATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTATGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGAtt                                                         <  2:3933583/150‑1 (MQ=255)

                                                                                             aaTCAATTTGTCTCAGGTTCGAATCCTGCAGGTCCCTGCAGTTTTTACATCCACCGCTGTTGTTATGCGTGATGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGAtt                                                         <  1:1580395/150‑1 (MQ=14)

                                                                                              aTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTg                                                        <  1:5404834/150‑1 (MQ=255)

                                                                                               tCATTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTgg                                                       <  1:6504109/150‑1 (MQ=255)

                                                                                               tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGGGAAGCGTGAAGGATAACGTGGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAgggtgg                                                       >  1:8121191/1‑145 (MQ=255)

                                                                                                ccaggggTCGCAGGTTCGAATGCTGCACGACCCACCAGGTTTAACATCAATCTCCGCTGTTTAGCGTGAAGGATATCGTTGCGCCTGAAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGc                                                      <  2:5291747/145‑1 (MQ=17)

                                                                                                cAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGc                                                      <  2:8124370/150‑1 (MQ=255)

                                                                                                 aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGcc                                                     >  1:3119051/1‑150 (MQ=255)

                                                                                                 aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACACAGAAGTTAACCGTGAAGGATAACGTTGCTCAAGCACCGGCCCGAAGGGCGGGCGAAG‑AGAGTCATACAGCACCACCCCCCAATCATAAAGGTTGGcc                                                     >  1:3628278/1‑150 (MQ=14)

                                                                                                   ttGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGcccc                                                   <  2:9377509/150‑1 (MQ=255)

                                                                                                     ggTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGa                                                 <  1:8885916/150‑1 (MQ=255)

                                                                                                       tCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGt                                               >  2:1230648/1‑150 (MQ=255)

                                                                                                       tCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGt                                               >  1:6242834/1‑150 (MQ=255)

                                                                                                       tCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGt                                               >  1:9106111/1‑150 (MQ=255)

                                                                                                       tCGCAGGTTCGAATCCTGCACGACCCACCAGTATTAACATCAAACTCAGATGTTAAGCGTTAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGt                                               >  1:5965880/1‑150 (MQ=255)

                                                                                                             gTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTGAAGATTGGCCCCGAGTAAAAAt                                         >  1:2692742/1‑150 (MQ=255)

                                                                                                             gTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAAt                                         <  1:7449038/150‑1 (MQ=255)

                                                                                                              ttCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATc                                        <  1:473020/150‑1 (MQ=255)

                                                                                                              ttCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATc                                        >  1:6473173/1‑150 (MQ=255)

                                                                                                                  aaTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGCCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCATAAAGATTGGCCCCGAGTAAAAATCTTTc                                    >  2:5606991/1‑150 (MQ=255)

                                                                                                                    tCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTCAAGATTGGCCCCGAGTAAAAATCTTTCAg                                  <  2:2692742/150‑1 (MQ=255)

                                                                                                                    tCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAg                                  <  2:5149982/150‑1 (MQ=255)

                                                                                                                     ccTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAgg                                 >  1:7907385/1‑150 (MQ=255)

                                                                                                                       tGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCTGGTa                               >  1:5873375/1‑150 (MQ=255)

                                                                                                                        gCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTaa                              <  2:9896090/150‑1 (MQ=255)

                                                                                                                            gACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACAcc                          >  1:3380431/1‑150 (MQ=255)

                                                                                                                             aCCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACAccc                         >  1:6913602/1‑150 (MQ=255)

                                                                                                                               ccaccaGTTTTAACATCTAACACAGATGTTAAGCGTGTAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCGCCGAGTAAAAAGCTTTGAGGTTCCAccgct                       >  2:4732989/1‑147 (MQ=17)

                                                                                                                                caccaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTa                      >  1:3997208/1‑150 (MQ=255)

                                                                                                                                 accaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTAt                     <  1:1111045/150‑1 (MQ=255)

                                                                                                                                 accaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTAt                     >  1:3834937/1‑150 (MQ=255)

                                                                                                                                 accaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTAt                     >  1:9712778/1‑150 (MQ=255)

                                                                                                                                  ccaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATg                    <  2:2865607/150‑1 (MQ=255)

                                                                                                                                   caGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATgg                   >  1:7603537/1‑150 (MQ=255)

                                                                                                                                     gTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGt                 <  2:4226545/150‑1 (MQ=255)

                                                                                                                                        ttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGt              <  2:2312956/150‑1 (MQ=255)

                                                                                                                                        ccaaCCTCCAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGt              <  2:8286899/148‑1 (MQ=255)

                                                                                                                                          aaCATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTa            <  2:4637107/150‑1 (MQ=255)

                                                                                                                                          aaCATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGGATGGGTCGTTa            >  2:7971237/1‑150 (MQ=255)

                                                                                                                                           aCATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCTTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAg           <  2:8854274/150‑1 (MQ=255)

                                                                                                                                           aCATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAg           >  2:4434911/1‑150 (MQ=255)

                                                                                                                                             aTCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCt         <  1:1995660/150‑1 (MQ=255)

                                                                                                                                               cAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCa       <  2:6093707/150‑1 (MQ=255)

                                                                                                                                                 aaCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGt     <  2:7443651/150‑1 (MQ=255)

                                                                                                                                                    tCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTgg  >  1:9221685/1‑150 (MQ=255)

                                                                                                                                                    tCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAG‑CGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTgg  >  1:4980067/1‑150 (MQ=255)

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