BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,287,122	T→G	11.3%	noncoding (115/171 nt) R19R (AGG→CGG)	rttR ← tpr ←	rtT sRNA, processed from tyrT transcript protamine‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,287,122	0	T	G	11.3%	455.0 / 26.2	229	noncoding (115/171 nt) R19R (AGG→CGG)	rttR tpr	rtT sRNA, processed from tyrT transcript protamine‑like protein
Reads supporting (aligned to +/- strand): ref base T (110/76); new base G (28/2); total (150/79)
Fisher's exact test for biased strand distribution p-value = 1.48e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 235 total aligned reads displayed.

GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTC  >  NC_000913/1286973‑1287271

 ccGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTg                                                                                                                                                       >  1:7904659/1‑150 (MQ=255)

    aaCCTTAGTCGAAGCTTCTCATCCGTCCCCGCATGGGCAGATTATTTGATTGGGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCgg                                                                                                                                                    <  2:3911735/150‑1 (MQ=255)

       cTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGGGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCGTCGGGTCGTTGCCGGCGGCac                                                                                                                                                 >  2:2162299/1‑149 (MQ=255)

        ttggTCGAAGCTTCTCCTCCGTCCCCGCCTGGGCAGAATATTTGAGTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAc                                                                                                                                                <  2:767046/147‑1 (MQ=11)

        ttAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCCCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAc                                                                                                                                                <  1:9000184/150‑1 (MQ=255)

           gTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGct                                                                                                                                             >  2:1327226/1‑150 (MQ=255)

            tCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGctc                                                                                                                                            >  1:5135203/1‑150 (MQ=255)

            tCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGcgc                                                                                                                                            >  1:8429405/1‑148 (MQ=255)

             cGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGACTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGctct                                                                                                                                           <  1:1204629/150‑1 (MQ=255)

              gAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGctctc                                                                                                                                          >  2:8829727/1‑150 (MQ=255)

              gAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGctctc                                                                                                                                          >  2:3705986/1‑150 (MQ=255)

                   gtctcAGCCTTCCCCGCTTGGGGAGTCTATTTGTTTGCGGATTCGCTTGAGCGTTGCGGGAGTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCt                                                                                                                                     <  2:4792804/149‑1 (MQ=11)

                    tctcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCATCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCCCAGCGCTCTCGcgg                                                                                                                                    >  2:3272245/1‑148 (MQ=14)

                    tctcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTg                                                                                                                                    >  1:7688072/1‑150 (MQ=255)

                    tctcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCTg                                                                                                                                    >  1:6443519/1‑150 (MQ=255)

                        aTCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGGCTGCGGCAACGCTCTCTCGCtttagc                                                                                                                                >  2:3565096/1‑145 (MQ=255)

                        aTCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTGGCCGGCGGCAACGCGCACTCGCTGgcgc                                                                                                                                >  1:2386492/1‑150 (MQ=255)

                        aTCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTGGAGAGTTCAGGGACTTTTGAGATAGGTGGTGGTGGGGGAAGGGATACTCAGCGCTGCGCGCTTCGCCCGTCGGGTCGAGGCCTGAGGCAACGCTCTCTCGGTGgcgc                                                                                                                                >  2:641242/1‑150 (MQ=14)

                            ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTGGCCTGCGGCCACCCTCTCACGCGGGCGCGCGa                                                                                                                            >  1:2421398/1‑150 (MQ=14)

                            ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCGGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCACAAACGCTGGCGCTCGa                                                                                                                            >  1:1710671/1‑150 (MQ=14)

                            ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGATTGATAGTTCAGGTACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTAAGAGCTGAGCGCTTCGCCCTTCGGGTCGTTGGCTGGGGCAACGCAATCTAGCTTGGGCTAGa                                                                                                                            >  1:407634/1‑150 (MQ=14)

                               cccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTc                                                                                                                         >  2:4516071/1‑150 (MQ=255)

                                ccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCATAGGCAACGCTCACTCGCTTGCGCTCGAGtag                                                                                                                        >  1:4870309/1‑148 (MQ=14)

                                  gCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGaa                                                                                                                      <  1:3454332/150‑1 (MQ=255)

                                   cATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAAc                                                                                                                     <  1:3135718/150‑1 (MQ=255)

                                     tGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTGGCCGGCGGCAACCCGCTCACGCTGGCGCTCGAGTCGAACCt                                                                                                                   >  1:4650132/1‑150 (MQ=14)

                                        gCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTAGTTGCCTGCGGCAACGCTATCTCGCTGGAGCTCGAGTCGAAACATgg                                                                                                                >  2:1420487/1‑150 (MQ=14)

                                         cAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCCGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGt                                                                                                               <  1:2260113/150‑1 (MQ=255)

                                           gAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCATCGCCCTTCGGTTCGATTCCTCACGCAAAGCACTCTCGCTGGCCCTAGAGTCGAACCTTGGTCg                                                                                                             >  2:1038486/1‑150 (MQ=14)

                                            aaTATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCTGGCGCGCGCGGCGAACCTTGGTcca                                                                                                            >  2:1875580/1‑148 (MQ=255)

                                                tttGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCTGGCGCTCGTCTCGAACCGTGGTAGAAGCt                                                                                                        >  2:5492153/1‑150 (MQ=14)

                                                tttGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGACCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCt                                                                                                        >  2:3063353/1‑150 (MQ=255)

                                                     ttGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGAGGCAACGATCTCTCGCTGTCGCTCGAATCGAACCTTTGTCGAAtcttctaa                                                                                                   >  2:3608033/1‑148 (MQ=14)

                                                      tGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTACGGGTCGTTGCCTGAGGCAAAGCTCTCTCGCTGGCCCACGAGTCGAACCTTGGTCGAAGCTTcacaa                                                                                                  >  2:2230615/1‑149 (MQ=14)

                                                        cGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTTGCGCTCGAGTCGAACCTTGGTCGAAGCATCTCATcc                                                                                                >  2:5071224/1‑150 (MQ=255)

                                                        cGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATcc                                                                                                >  2:6883461/1‑150 (MQ=255)

                                                          gATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAACCtt                                                                                              >  2:5296866/1‑150 (MQ=255)

                                                            ttCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTcc                                                                                            <  2:6050265/150‑1 (MQ=255)

                                                            ttCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTAGAATCTTATCaaacatcc                                                                                            >  1:7603900/1‑143 (MQ=14)

                                                               gCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCGCGCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTATCATACGTccccc                                                                                         >  1:8180784/1‑149 (MQ=14)

                                                                 ttGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACCCGCTCTCGCTGGCGCTCGCGTCGAACCTTGGTAGAAGCTTCTCATACATCCCCGCa                                                                                       >  1:7255798/1‑150 (MQ=14)

                                                                 ttGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCAACACAACCTTCCCCGcc                                                                                       >  1:6187380/1‑149 (MQ=14)

                                                                   gagaGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAAGATCTCTCGCTTTCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAACCTTACCCgaagg                                                                                     >  2:4965258/1‑146 (MQ=14)

                                                                       gTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTTGCGCTCGAGTCGAACCTGGGTCGAAGCTTCTCATCCATCACCGCA‑GGGGCa                                                                                >  2:4448695/1‑150 (MQ=14)

                                                                       gTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCa                                                                                >  2:8335584/1‑150 (MQ=255)

                                                                       gTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCa                                                                                <  2:9217429/150‑1 (MQ=255)

                                                                       gTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAACCTTCCCCGCT‑TGGGcc                                                                                >  2:6932973/1‑149 (MQ=255)

                                                                        ttCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGACGTTGCCTGCGGCAACCCTCTCACGCGGGCCCACCAGGAGAACCCTGGTCGAAGCAGCGCAACCTTCCCCGCC‑TGGGccg                                                                               >  2:1893735/1‑148 (MQ=14)

                                                                          cAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACCCTCTCAAGCTTGCGCTCCAGTCGAACCGTGGTCGAAGATTCAAATCCGACCCCGCA‑gggggcgcc                                                                             >  1:3576834/1‑145 (MQ=14)

                                                                           aGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCAGAat                                                                            >  2:8131963/1‑150 (MQ=40)

                                                                           aGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCATCGCCCGTCGGGTCGTTGCGTGAGGCAAAGCTCTAAAGGTGGCGCGCGACTCGAGCCGTGGTCGAAGCTTCTCATCCTTCCGCGCC‑Tggggcgggg                                                                            >  2:7867689/1‑144 (MQ=14)

                                                                             ggACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGAGGCAAAGCTATCTCGCTTGCGCTCGAGTCGAAACATGGTCGAATCTTCTCATACTTCCCCGCA‑TGGGCAGAatat                                                                          >  2:2487158/1‑150 (MQ=14)

                                                                             ggACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACCCTCGCTCGCGGGCGCGCGAGTCGAACCAGGGTCGCAGCCTCGCATACATCCCCCCC‑GGGGCCGCatat                                                                          >  1:1037386/1‑150 (MQ=14)

                                                                               aCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCAGAATAttt                                                                        >  2:1763773/1‑150 (MQ=38)

                                                                                 ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCATCGCTGCGCTCTTCGCCCTTCGGGACGTGGACTTCGACAACGCGATCTCGCTGCCGCTAGAGTAGAACCTTGGTCGACGCTTCACAAACTTCCCCGCA‑TGGGCAGAATAGTTaa                                                                      >  2:6843759/1‑150 (MQ=11)

                                                                                 ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCATCGCCCTTCGGGACGTTGCCTGAGGCAACGCTCTCTCGCTGGCGCTCGAATCGAACCGTGGTCGAAGCATATCCTAATACCCCCCT‑GGGGCAAAATAGTTaa                                                                      >  2:5845878/1‑150 (MQ=11)

                                                                                  tttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAATCGAACCTTGGTCGAAGCTTCTCAACCTTCCCCCCATTGGGCaaaaattaaaa                                                                      >  2:195212/1‑143 (MQ=11)

                                                                                    tGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTACGCCCTTCGGGTCGTTTCCTGCGGCAACGCTCTATCGCAGGCGCGCGAGTCGAACCAAGGTCGAAGCTTCTCAACCATCCCCCCC‑TGGGCCGAAAAATTAAATg                                                                   >  1:8000192/1‑150 (MQ=11)

                                                                                     gAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCAGAATATTTAATTGc                                                                  <  1:3167044/150‑1 (MQ=39)

                                                                                          tGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGGGGCAAAGCGCTCTCGCTGGCGCCCGCGTCGAACCTTTGTCGAAGCTTCTCAACCATCCCCCCA‑CGGGGCGAAAATTTAAATTCGGaat                                                             >  1:2672410/1‑148 (MQ=11)

                                                                                          cgATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCCAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCAGAATATTTAATTGCGGAtt                                                             >  2:5726925/2‑150 (MQ=18)

                                                                                            aTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCCGCCACCCACACACGCAGGCCCACCAGTCGAACCATGGTCGAAACTTCACAACCATCCCCCCC‑GGGGCCGAAAAATTAATTGCGGAGTCg                                                           >  1:7561592/1‑150 (MQ=2)

                                                                                             tggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCAGAATATTTAATTGCGGATTCGt                                                          <  1:5726925/150‑1 (MQ=37)

                                                                                                tggtgtGGGATTGATTTCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCAGAATATTTAATTGCGGATTCGTTgg                                                       <  1:7034629/150‑1 (MQ=255)

                                                                                                tggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCAGAATATTTAATTGCGGATTCGTTgg                                                       <  1:6826048/150‑1 (MQ=38)

                                                                                                tggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCAGAATATTTAATTGCGGATTCGTTgg                                                       >  2:5690665/1‑150 (MQ=38)

                                                                                                       ggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCTCTCTCGATGGCGCTCGAGGCGCACCTTGGGCGACGCTTCTCATCCGTCCCCGCA‑TGGGCCGGATAGTTCATTGCGGGTTCGGTggggcggtc                                                >  2:4328849/1‑144 (MQ=11)

                                                                                                         aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAg                                              <  1:4493808/150‑1 (MQ=255)

                                                                                                          aGGATTTCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAgg                                             <  1:3713622/150‑1 (MQ=255)

                                                                                                          aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAgg                                             <  1:948256/150‑1 (MQ=255)

                                                                                                               tACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACtt                                        <  2:7258917/150‑1 (MQ=255)

                                                                                                                  tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGGCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCCGAATAGTTAATTGCGGATTCGGTGGGAAGTTCAGGGACGTTTg                                     >  2:6731167/1‑150 (MQ=14)

                                                                                                                   cAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGa                                    <  1:6079605/150‑1 (MQ=255)

                                                                                                                    aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA‑GGGGCCGAACATTGAAGTGCGGATTCGGGGGGGAGGTCCGGGACGTTTGaa                                   >  2:9115873/1‑150 (MQ=14)

                                                                                                                    aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAACCATCCCCGCA‑TGGGCAGAAGATTGAAATGCGGGTTCGGGGGGAAGTTCAGGGACTGTTGaa                                   >  2:6161588/1‑150 (MQ=14)

                                                                                                                    aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCGTGGTCGCAGCTTCTCATCCGTCCCCGCC‑TGGGGAGAATAGTTTATTGCGGGTTCGGTGGGGAGTTCAGGGACGGTTGat                                   >  2:78239/1‑149 (MQ=14)

                                                                                                                     gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGaaa                                  <  2:3993171/150‑1 (MQ=255)

                                                                                                                     gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCTCACTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTACCCCCA‑GGGGCAGAATATTTTATTGCGGGTTCGTTGGGCAGTTCAGGGACGTTTGtaa                                  >  2:1171574/1‑150 (MQ=14)

                                                                                                                      cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAg                                 >  2:6600878/1‑150 (MQ=255)

                                                                                                                         tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTccc                                                                                           <  1:2190004/90‑1 (MQ=33)

                                                                                                                                ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtggt                       <  1:5921249/150‑1 (MQ=255)

                                                                                                                                  cGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtggtgg                     <  2:2148590/150‑1 (MQ=255)

                                                                                                                                   gCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTggg                    <  2:2909825/150‑1 (MQ=255)

                                                                                                                                       ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAACCTTCCCCGCA‑TGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                >  2:422062/1‑150 (MQ=255)

                                                                                                                                         gagagcGTTGCCGCAGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtggt                       <  1:6627584/136‑1 (MQ=255)

                                                                                                                                            gTCGTTGCCTGCGGCACCGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAtt           <  1:2093254/150‑1 (MQ=255)

                                                                                                                                            gTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAACCTTCCCCGCA‑GGGGCACAAAATTTACATGCGGACTCGTGGGGACGGTCAGGGCCGTTGGAAAGGGCTGGTGGGGGGGGCAGGAta           >  2:7366746/1‑149 (MQ=17)

                                                                                                                                             tCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTc          <  1:1420487/150‑1 (MQ=255)

                                                                                                                                             ccgccgCCCGCGCCAACGCCCTCTCGCTGTCGCTCGAGTCGCACCTCCGTCGAATCTTCTCAGCCTTCCCCTCT‑TGGGCTGAACTACTAACTTCCGATCCGTTGGTCACTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTc          <  2:1266001/145‑1 (MQ=9)

                                                                                                                                               gTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGa        >  2:162943/1‑150 (MQ=255)

                                                                                                                                                ttGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGt                                                                                                                          <  1:1548851/36‑1 (MQ=40)

                                                                                                                                                ttGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCAGAATATTAAATTGCGGATGCGTTGTGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGGAGGATTAGaa       >  2:2718442/1‑150 (MQ=17)

                                                                                                                                                 tGCCTGCGGCAACGCTCTCTCGCTGGCGCTCTAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtg                         <  1:4958323/131‑1 (MQ=255)

                                                                                                                                                 tGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAACCGTCCCCCCA‑TGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGGAGGATTAGAAc      >  1:129630/1‑150 (MQ=17)

                                                                                                                                                  gCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAAcc     <  1:5690665/150‑1 (MQ=255)

                                                                                                                                                   ccTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTACCCGCA‑TGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGACGGATTCGAACCt    >  1:5767456/1‑150 (MQ=255)

                                                                                                                                                   ccTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCATCTCAAACAACCCCCCA‑AGGGCAGAAAATTTAATTGCGGTTTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGGAGGATTAGTAcca    >  1:7683045/1‑149 (MQ=17)

                                                                                                                                                    cTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCtt   <  2:803652/150‑1 (MQ=255)

                                                                                                                                                     tGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCAGAATATTTAATTGCGGATTCGTGGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTc  >  2:1684540/1‑150 (MQ=255)

                                                                                                                                                     tGCGGCAACGCTCTCTCGCTGCCGCTCGAGCCGCACCTTCGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCAGAATATTTCATTGCGGATTCGCTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTc  <  2:6257276/150‑1 (MQ=14)

GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA‑TGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTC  >  NC_000913/1286973‑1287271

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