BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	781,254	A→T	8.1%	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	8.1%	392.4 / 13.2	125	intergenic (+32/‑115)	lysY/lysZ	tRNA‑Lys/tRNA‑Lys
Reads supporting (aligned to +/- strand): ref base A (51/61); new base T (4/8); total (55/70)
Fisher's exact test for biased strand distribution p-value = 5.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 131 total aligned reads displayed.

GGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTT  >  NC_000913/781106‑781403

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

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

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

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

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

    ggggTCGGCGGTTCGATCCCGTCATCACCCACCACTCGTGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGAt                                                                                                                                                  <  2:750305/150‑1 (MQ=255)

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

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

      ggTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATaa                                                                                                                                                >  2:2151269/1‑150 (MQ=255)

      ggTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAAATCAGATGTTAAGCGTGAAGGaaaa                                                                                                                                                >  2:2151174/1‑150 (MQ=255)

       gTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAAc                                                                                                                                               >  1:1602921/1‑150 (MQ=255)

        tCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACg                                                                                                                                              <  1:889369/150‑1 (MQ=255)

        tCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCGCAGTTGGTAGAGCAGTTGACTTTTAAACAATTTGTCGCAGGTTCTAAAACTGCACGACCCACCACTTTTAACAACAAAAACCGAGGTTAAGCGGGAAGGATAACg                                                                                                                                              >  1:2244738/1‑150 (MQ=11)

          ggcggTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGtt                                                                                                                                            >  2:485165/1‑150 (MQ=255)

           gcggTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTg                                                                                                                                           <  1:71552/150‑1 (MQ=255)

           gcggTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAAATCAGATGTTAAGCGTGAAGGATAACGTTg                                                                                                                                           >  1:472011/1‑150 (MQ=255)

             ggTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACCTTgcg                                                                                                                                         >  1:827663/1‑150 (MQ=255)

             ggTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAAAATCAAAATAAGTAGTTTATCGGGAAGGATAACGTTgcg                                                                                                                                         >  1:159337/1‑150 (MQ=14)

                 cGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATCGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAg                                                                                                                                     <  2:2169163/150‑1 (MQ=255)

                   aTCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCTAACTCAGATGTTAAGCGGGAAGGATAACGTTGCGCCAGCa                                                                                                                                   >  2:1424835/1‑150 (MQ=255)

                   aTCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAATCTCAGATGTTAAGCGTGACGGATAACGTTGCGCCAGCa                                                                                                                                   >  2:29402/1‑150 (MQ=255)

                    tCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCCAAATCAGATGTTAAGCGGGAAGGATAACGTTGCGCCAGCaa                                                                                                                                  >  1:2271928/1‑150 (MQ=255)

                    tCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAAAGTTGCGCCAGCaa                                                                                                                                  >  1:656129/1‑150 (MQ=255)

                        gTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCGGGGTCGAATCCTGCCCGACCCACCAGTTTTAACATCAAACACAGATGTTTAGCGTGAAGGATAACGTTGCGCCAGCAACGGc                                                                                                                              >  1:1084319/1‑150 (MQ=14)

                            tcaCCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCgg                                                                                                                          >  1:1384141/1‑149 (MQ=255)

                             caccaaccACTCGGGTCGTTAGCTAAGTTGGTATAGCTTTTTAATTTTAATCAATTGGTCGCTGTTTCGATTCCTGCACGACCCACCAGTTTTAACATCAAAATCAGATTTTAAGCGTGAAGTATAACGTTGCGCCAGCAACGGCCCGTa                                                                                                                         <  2:614685/150‑1 (MQ=14)

                                       tCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTAAACATCAAAAGAAGATGATAACCTTGATGGATAACGTTCCGACAGCAACGGCCCGTAGGgcgtgcga                                                                                                               >  2:1772043/1‑150 (MQ=11)

                                            tCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGagcgaagcga                                                                                                          <  2:1643270/150‑1 (MQ=255)

                                               ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTc                                                                                                       >  1:2220668/1‑150 (MQ=255)

                                               ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTc                                                                                                       >  1:2220704/1‑150 (MQ=255)

                                                tAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCa                                                                                                      >  1:939211/1‑150 (MQ=255)

                                                tAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCa                                                                                                      >  1:2240220/1‑150 (MQ=255)

                                                tAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCa                                                                                                      >  1:939304/1‑150 (MQ=255)

                                                tAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCa                                                                                                      >  1:2399735/1‑150 (MQ=255)

                                                tAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCa                                                                                                      >  1:380410/1‑150 (MQ=255)

                                                tAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAACGATAACGATGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCa                                                                                                      >  1:1426706/1‑150 (MQ=255)

                                                              gagCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGATGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCAGCACGACcca                                                                                        >  2:1729873/1‑150 (MQ=255)

                                                              gagCAGTTGACTTTTAATCAATTGGTCGCAGGTGCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCTTGAAGGATAACGTTGCGCAAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCATGCACGACcca                                                                                        >  2:1095783/1‑150 (MQ=255)

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

                                                                                cAATTGGCCCCAGGTTCGAATCCTGCACGACCCACCACTATTAACAACAAACTCAGATGTTAAGCTTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTTAAGATTggt                                                                      <  1:1753326/150‑2 (MQ=17)

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

                                                                                     ggTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCTAACACAGATGTTAAGCGTGTAGGATAACGTTGCGACAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGa                                                                 <  2:129641/150‑1 (MQ=255)

                                                                                       tCGCTGGTTCGAATCCTGCACGACCCACCAGTTTTAACATAAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGt                                                               <  1:1106035/150‑1 (MQ=255)

                                                                                       tCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCATAAAGAATGGCCACGAGt                                                               >  1:409147/1‑150 (MQ=255)

                                                                                               tCGAATCCTGCACGACCCACCAGTTTTAACATCTAACTCAGATGTTAAGCGTGTAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTaaa                                                            <  1:1544498/145‑1 (MQ=255)

                                                                                               tCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTaaa                                                            >  2:1544498/1‑145 (MQ=255)

                                                                                               tCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTaaa                                                            >  1:694831/1‑145 (MQ=255)

                                                                                                cGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGATAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATACTTCACGACCCACCAATCTTAAAGATTGGCCCCGAGAAAAAATCtt                                                      >  1:333761/1‑150 (MQ=255)

                                                                                                  aaTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCCATCTTCAAGATTGGCCCCGAGTAAAAATcggtc                                                    >  2:1631080/1‑146 (MQ=255)

                                                                                                  aaTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGTCCCCGAGTAAAAATCGTTc                                                    >  2:1757882/1‑150 (MQ=255)

                                                                                                  aaTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGGGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCCATCTTAAAGATTGGCCCCGAGTAAAAATcgggc                                                    >  2:2006313/1‑146 (MQ=255)

                                                                                                  aaTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGAAGTTAAGCGTGAGGGATAACGTTGCGACAGCAACGGCCCGTAGGGCGAGAGAAGCGAGTCATCCTGCGCGACCCACCCATCTTAAAGATTGAGGCAGAGTaaaaagcgtgc                                                    >  2:1538297/1‑144 (MQ=17)

                                                                                                   aTCCTGCACGACCCACCAGTTTTAACATCCACCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCTGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCa                                                   <  2:704801/150‑1 (MQ=255)

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

                                                                                                    tCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTGGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTCAAAATCTGTCAg                                                  >  2:2259029/1‑150 (MQ=255)

                                                                                                     ccTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTaaaaa                                                          <  1:1423840/141‑1 (MQ=255)

                                                                                                     ccTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTaaaaa                                                          >  2:1423840/1‑141 (MQ=255)

                                                                                                       tGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCt                                                       <  2:406744/142‑1 (MQ=14)

                                                                                                       tGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaaaagcg                                                                        <  2:1788339/125‑6 (MQ=255)

                                                                                                       tGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTa                                               <  2:1297361/150‑1 (MQ=255)

                                                                                                       tGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTa                                               <  2:1748779/150‑1 (MQ=255)

                                                                                                        gCACGACCCTCCAGCTTTAACATGAAAATCAGATGTTTTGCGTTTAGTTTATTTTTGCGGCTGCATCTTCCGTTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTaa                                              <  1:138383/150‑1 (MQ=16)

                                                                                                        gCACGACCCCCCAGTTTTAACAGCTAACTCAGATGTTTAGCGTGAAGGATAACGTTGCGCCAGCCACGGCCGGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTaa                                              <  2:2081252/150‑1 (MQ=21)

                                                                                                        gCACGACCCACAAGTTTTAACATCCAACACAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTaa                                              <  2:1290375/150‑1 (MQ=255)

                                                                                                                caccaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTa                                      <  1:105620/150‑1 (MQ=255)

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

                                                                                                                 accaGTTTTAACACCATACTCAGATGTTAAGCGTGTAGGTTAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTAt                                     <  1:1824822/150‑1 (MQ=255)

                                                                                                                 aacaGTTTTAACATCTTACTCAGATGTTAAGCGTGTAGGATATCGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTAt                                     <  1:1247197/148‑1 (MQ=255)

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

                                                                                                                      ttttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTc                                >  2:325131/1‑150 (MQ=255)

                                                                                                                        ttAACATCCAAGTCAGATGTTAAGCGTGTAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGt                              <  2:1981542/150‑1 (MQ=255)

                                                                                                                        ttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCATTCAGGTAACAACCGTATGGGTCGt                              >  2:1531023/1‑150 (MQ=255)

                                                                                                                             aTCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTATCAGGTAACACCCGTATGGGTCGTTAGCt                         >  1:1660747/1‑150 (MQ=255)

                                                                                                                              tcgaagacAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTc                        >  1:1149454/8‑150 (MQ=255)

                                                                                                                                 aaCTCAGATGTTAAGCGTGAAGGATACCGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCATAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGt                     <  1:2402577/150‑1 (MQ=255)

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

                                                                                                                                    tCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTgg                  >  1:2190314/1‑150 (MQ=255)

                                                                                                                                      tgctgatAAGCTTGCAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTa                <  1:300858/144‑1 (MQ=255)

                                                                                                                                      aGATTTTAAGCTTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTa                <  1:432443/150‑1 (MQ=255)

                                                                                                                                      aGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTa                <  2:1149454/150‑1 (MQ=255)

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

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

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

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

                                                                                                                                        aTGTTAAGCGTGAAGGATCACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTaga              <  1:751193/150‑1 (MQ=255)

                                                                                                                                        aTGTTAAGCGTGAAGGATAACGTTGCGCCCGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTaga              <  2:637265/150‑1 (MQ=255)

                                                                                                                                        aTGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTaga              <  2:1482999/150‑1 (MQ=255)

                                                                                                                                          gTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGc            <  1:2050708/150‑1 (MQ=255)

                                                                                                                                          gTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGc            <  1:202808/150‑1 (MQ=255)

                                                                                                                                           ttAAGCGTGAAGTATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCa           <  2:1426706/150‑1 (MQ=255)

                                                                                                                                           ttAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCa           <  2:939304/150‑1 (MQ=255)

                                                                                                                                           ttAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCa           <  2:2240220/150‑1 (MQ=255)

                                                                                                                                           ttAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCa           <  2:939211/150‑1 (MQ=255)

                                                                                                                                           ttAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCa           <  2:2399735/150‑1 (MQ=255)

                                                                                                                                             aaGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGt         <  2:436096/150‑1 (MQ=255)

                                                                                                                                             aaGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGt         <  2:1768030/150‑1 (MQ=255)

                                                                                                                                              aGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGtt        <  2:1084319/150‑1 (MQ=255)

                                                                                                                                                cGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGa      <  2:2220668/150‑1 (MQ=255)

                                                                                                                                                cGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGa      <  2:2220826/150‑1 (MQ=255)

                                                                                                                                                    aaGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGTTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACttt  <  1:818415/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: `‑`