BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,287,420	G→C	5.6%	intergenic (‑92/+118)	tyrV ← / ← tyrT	tRNA‑Tyr/tRNA‑Tyr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,287,420	0	G	C	5.6%	392.2 / 9.5	117	intergenic (‑92/+118)	tyrV/tyrT	tRNA‑Tyr/tRNA‑Tyr
Reads supporting (aligned to +/- strand): ref base G (56/54); new base C (3/4); total (59/58)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 119 total aligned reads displayed.

CGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGA  >  NC_000913/1287271‑1287567

cGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCtg                                                                                                                                                     >  2:2039680/1‑150 (MQ=255)

 gAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTTCCtgt                                                                                                                                                    >  2:2409678/1‑150 (MQ=255)

   aGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTgg                                                                                                                                                  >  1:432664/1‑150 (MQ=255)

    ttCGATGACGCCAGATCTACAGCCTCCGCCCTTTGGCAGATCTTTAACCCCACCAGTCGTAATTCATATTTTGCCGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGc                                                                                                                                                 <  1:1375080/149‑1 (MQ=25)

    gTCGATGTCGGCCGATGTCCTGGCTGCCCCCTTTGGCCGCACCGGTCCCCCCCCAGGGGTAAGTCACATCTGGCCGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGc                                                                                                                                                 <  1:1366356/150‑1 (MQ=25)

       gATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAAc                                                                                                                                              >  1:484185/1‑150 (MQ=255)

        aTGACGGAACATGTACACGCTGCCCCCTTTGGCCGCTCTGTAACCCCACCAGGGGTAATTCAAATTTTGCCGTAATGATTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACg                                                                                                                                             <  2:1857378/150‑1 (MQ=37)

          gACGTCAGATTTACAGTCTGCTCCCTTTGGCCGATCTGTAAACCCACGAGGGGTAATTCATTTTTTGTTGTTTTGCTTGTGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCACAGGCAACggc                                                                                                                                           <  2:1018356/150‑3 (MQ=25)

            cgccAGATCTACAGGCTGCCCCCTTTGGCCGATCCGTAACCCCACCAGGGGTAATTCAAATTTGGACGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCACAGGCAACggccc                                                                                                                                         <  2:1851440/147‑5 (MQ=255)

            cGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTtct                                                                                                                                         >  1:2317614/1‑150 (MQ=255)

            cGGCAGATTTACAGGCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGCGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCACAGGCAACggccc                                                                                                                                         <  1:734178/150‑5 (MQ=38)

             ggCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTtctc                                                                                                                                        <  1:1375727/150‑1 (MQ=255)

             ggCAGATGTACAGTCTGCGCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTATGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTtctc                                                                                                                                        <  1:989823/150‑1 (MQ=255)

              tatgATTTCCACTCTGCTCCCTTTGGCATCTCTTTAACTTCAGCAGGGGTAATTAAATTTTTGTCGTATTGATTGTGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTtctct                                                                                                                                       <  2:1626072/147‑1 (MQ=21)

               cAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTtctctc                                                                                                                                      <  1:138088/150‑1 (MQ=255)

               cAGATCTCCAGGCTGCGCCCTTTGGCCGCTCGTGAACCCCACCAGGGGTAATTCAAATCTTGAGGTAATGAATGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTtctctc                                                                                                                                      <  1:1509093/150‑1 (MQ=38)

                aGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCg                                                                                                                                     <  1:1674078/150‑1 (MQ=255)

                 gATTTCCAGCCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCACATCTGGCCGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGc                                                                                                                                    <  1:347895/150‑1 (MQ=38)

                 gATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGGTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGc                                                                                                                                    >  2:1319615/1‑150 (MQ=255)

                 gATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCACAGGCAACgg                                                                                                                                            >  1:899126/1‑141 (MQ=255)

                  aTTTACAGGCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCt                                                                                                                                   <  1:1451520/150‑1 (MQ=255)

                        aGTCTGCTCCCTTTGGCCGCTCTTTAACCCCACCAGGGGTAATTCAAATTTTGATGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCt                                                                                                                             <  2:2313991/150‑1 (MQ=255)

                        aGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCt                                                                                                                             >  2:1166727/1‑150 (MQ=255)

                        aGTCTGCGCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCt                                                                                                                             <  1:727606/150‑1 (MQ=255)

                         gtatGCTCCCTTTGTCCGCTCCGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTAGTCGCTTCGCTCATAACCATTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCtc                                                                                                                            >  2:745416/4‑150 (MQ=38)

                          tCTGCTCCCTTTGGCAGCTCGTGAAGCCCACGAGGTGTATTTAATTTTTTGTTGTAATGATTGTGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCtcg                                                                                                                           <  2:1969648/150‑1 (MQ=25)

                          gcTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCtcg                                                                                                                           <  1:1720621/149‑1 (MQ=255)

                             gCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCtcgaat                                                                                                                        >  2:1511071/1‑150 (MQ=255)

                               tCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCtcgaatcg                                                                                                                      >  2:2444052/1‑150 (MQ=255)

                                cccTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCtcgaatcga                                                                                                                     >  2:2110762/1‑150 (MQ=255)

                                     tGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCtt                                                                                                                >  2:225301/1‑150 (MQ=255)

                                              tttaaCCCCACCAGGGGTATTTCAAATTTTGATGTTATGATTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAgg                                                                                                       <  2:509263/147‑1 (MQ=37)

                                              gggAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTAGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAgg                                                                                                       >  2:12374/1‑150 (MQ=255)

                                               ggAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGATCCTAACCCTTCGGGCCGTTGCATGTGGCAAAGTTCTCTCGATTTCGCTCGAATCGGACCTTAGTCGAAGGt                                                                                                      >  2:689341/1‑150 (MQ=38)

                                                   ccccACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTtctc                                                                                                  <  1:683291/150‑1 (MQ=255)

                                                   ccccACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTtctc                                                                                                  >  2:1256837/1‑150 (MQ=255)

                                                    cccACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCa                                                                                                 <  2:484185/150‑1 (MQ=255)

                                                     cccccaGGGGTAATTCAAATTTTGAGCTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTCGACTAAGGTtcgatt                                                                                                <  1:1970238/147‑5 (MQ=37)

                                                     ccaccaGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCAc                                                                                                >  1:2314210/1‑150 (MQ=255)

                                                      caccaGTGGTAATTCAATTTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTCGACTAAGGTtcgattc                                                                                               <  1:2193852/150‑6 (MQ=37)

                                                      caccaGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCAcc                                                                                               >  1:847241/1‑150 (MQ=255)

                                                      caccaGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCATCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAAACTTAGTCGAAGGTTCTCAcc                                                                                               >  1:1646368/1‑150 (MQ=255)

                                                       accaGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTTGTGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCTCTCGAATCGAACCTTAGTCGAAGGTTCTCAccc                                                                                              <  2:2288176/150‑1 (MQ=255)

                                                        ccaGTGGTTATTCATATCTTGCTGTAATGCTAGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTCGACTAAGGTtcgattcga                                                                                             <  1:2090424/150‑8 (MQ=21)

                                                        ccaGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCAcccc                                                                                             >  1:1823344/1‑149 (MQ=255)

                                                        ccaGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCt                                                                                             >  1:871543/1‑150 (MQ=255)

                                                           ggggTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTcc                                                                                          >  2:1797341/1‑150 (MQ=255)

                                                                aaTTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGAt                                                                                     <  1:1511071/150‑1 (MQ=255)

                                                                  ttCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGATCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGa                                                                                   >  2:63708/1‑150 (MQ=255)

                                                                    cAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGt                                                                                 >  1:197621/1‑150 (MQ=255)

                                                                       aTTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCa                                                                              <  1:2444052/150‑1 (MQ=255)

                                                                          ttgccgTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAAc                                                                           <  1:1791875/145‑1 (MQ=255)

                                                                             aGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGCGTGCAAACttt                                                                        >  2:256698/1‑150 (MQ=255)

                                                                               gTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTca                                                                      <  1:421835/150‑1 (MQ=255)

                                                                                tAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTcaa                                                                     >  2:2063127/1‑149 (MQ=255)

                                                                                tAATGCTTGAGATGGGGGGGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTGCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTcac                                                                     <  1:1458437/150‑1 (MQ=255)

                                                                                 aaTGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTTCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATAGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTcaca                                                                    >  2:1478454/1‑150 (MQ=255)

                                                                                    gCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAAtc                                                                 <  2:2317614/150‑1 (MQ=255)

                                                                                         agaTGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCAcc                                                            <  1:1166727/150‑1 (MQ=255)

                                                                                         agaTGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCAcc                                                            >  2:1806115/1‑150 (MQ=255)

                                                                                         agaTGGTGGTGGGGGAAGGATTATTCGTAGCTAAGATCCTCACCCTTCGGTACGTTTCATGTGGCAACGTTCTCTAGCTTTCGATCGAATAGAACCTTAGTCGAAGGTTCTCAACCTACCACGATGAGTGAAAACTTTCACAATCTCAcc                                                            >  2:2013354/1‑150 (MQ=25)

                                                                                          gaTGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCg                                                           <  2:1112305/150‑1 (MQ=255)

                                                                                           aTGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGa                                                          <  2:859000/150‑1 (MQ=255)

                                                                                               tggtggGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGtt                                                      >  2:836300/1‑150 (MQ=255)

                                                                                                ggtggGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTa                                                     >  2:2476900/1‑150 (MQ=255)

                                                                                                ggtggGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTa                                                     <  2:1170009/150‑1 (MQ=255)

                                                                                                ggtggGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAAACTCACCGAAGTTa                                                     >  2:628125/1‑150 (MQ=255)

                                                                                                 gtggGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTAc                                                    <  1:2110762/150‑1 (MQ=255)

                                                                                                 gtggGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTAc                                                    <  1:1244313/150‑1 (MQ=255)

                                                                                                  tggGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTAcc                                                   >  2:618900/1‑150 (MQ=255)

                                                                                                    ggggAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACcac                                                 >  1:2360346/1‑150 (MQ=255)

                                                                                                    ggggAAGGATTATTCGTCGCTTCGCTCCTCACCATTCGGGCCTTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACcac                                                 <  1:194012/150‑1 (MQ=255)

                                                                                                     gggAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCACACTTTCACAATCTCACCGAAGTTACcaca                                                >  1:2489965/1‑150 (MQ=255)

                                                                                                     gggAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACcaca                                                <  1:2049762/150‑1 (MQ=255)

                                                                                                      ggAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACAt                                               >  2:1393806/1‑150 (MQ=255)

                                                                                                      ggAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACAt                                               >  1:1376331/1‑150 (MQ=255)

                                                                                                       gAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATc                                              <  1:649804/150‑1 (MQ=255)

                                                                                                         aGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGc                                            >  1:1054621/1‑150 (MQ=255)

                                                                                                          ggATTATTCGTCGCTTCGCTCCTCACCCTACGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCt                                           >  2:1756823/1‑150 (MQ=255)

                                                                                                           gATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCtg                                          >  1:565134/1‑150 (MQ=255)

                                                                                                               attCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTgtggt                                      <  1:1334252/150‑1 (MQ=255)

                                                                                                               attCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCCTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTgtggt                                      <  2:1823344/150‑1 (MQ=255)

                                                                                                                    tcgcttcgctCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCGCCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAAtt                                 >  1:1766100/1‑150 (MQ=255)

                                                                                                                    tcgcttcgctCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAAtt                                 <  1:1620468/150‑1 (MQ=255)

                                                                                                                      gcttcgttCCTCTCCCTTCGGGCCGTTGGCTGTGGCAACGTTCTCTGGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTAt                               <  2:2294692/150‑1 (MQ=255)

                                                                                                                       cttcgctCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTAtg                              >  2:81382/1‑150 (MQ=255)

                                                                                                                         tcgctCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTAtggt                            <  1:2039680/150‑1 (MQ=255)

                                                                                                                           gctCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCCCCGAAGTTACCACATCGCTGTGGTGAATTAtggtgg                          >  1:307911/1‑150 (MQ=255)

                                                                                                                              cctcACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTAtggtggtgg                       <  2:2449408/150‑1 (MQ=255)

                                                                                                                                tcACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTgggg                     >  1:919786/1‑150 (MQ=255)

                                                                                                                                tcACCCTTAGGGCCGTTGCCTGTGGCAACGTTCTATCGCTTTCGATCGAATAGAACCTTAGTAGAAGGTTATCACACTTAACCGAAGAGTGCAAACTTTCACAATCTCACCGAAGTTACAACATAGCTGTGGTGAATTATGGTGGTgggg                     >  2:2267497/1‑150 (MQ=255)

                                                                                                                                      ttCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGa               >  1:810253/1‑150 (MQ=255)

                                                                                                                                      ttCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGa               <  1:1393806/150‑1 (MQ=255)

                                                                                                                                         gggCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTc            >  1:1896687/1‑150 (MQ=255)

                                                                                                                                          ggCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCg           <  1:2409678/150‑1 (MQ=255)

                                                                                                                                           aacGTTGCCACAGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGa          >  1:1063410/3‑150 (MQ=255)

                                                                                                                                            ccGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGaa         >  1:810717/1‑150 (MQ=255)

                                                                                                                                             cGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAAc        >  2:1672798/1‑150 (MQ=255)

                                                                                                                                                   cTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGTGGGAAGGATTCGAACCTTCGa  >  1:1116621/1‑150 (MQ=255)

                                                                                                                                                   cTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGa  <  1:81382/150‑1 (MQ=255)

CGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGA  >  NC_000913/1287271‑1287567

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