BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,287,124	C→A	6.7%	noncoding (113/171 nt) R18L (CGC→CTC)	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,124	0	C	A	6.7%	639.8 / 9.6	233	noncoding (113/171 nt) R18L (CGC→CTC)	rttR tpr	rtT sRNA, processed from tyrT transcript protamine‑like protein
Reads supporting (aligned to +/- strand): ref base C (124/79); new base A (21/2); total (152/81)
Fisher's exact test for biased strand distribution p-value = 4.78e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 239 total aligned reads displayed.

GAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTC  >  NC_000913/1286976‑1287271

gAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGACAGAATATTTTATTGCGGATTCGCTTGAGAGTTCAGGGACTAATGAAATTGATCGTGCAGGGTGAAGGATTACTCAGCGCTGCGAGATTCTCACTACGGGTCGTTGCCTGCg                                                                                                                                                    >  1:2342783/1‑150 (MQ=14)

 aaCCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCgg                                                                                                                                                   <  1:6921344/150‑1 (MQ=255)

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

        gTCGACGCTTCTCCTCCGGCCCCGCATGGGCAGAATATTGGAGTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGct                                                                                                                                            <  2:7745827/150‑1 (MQ=14)

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

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

         tCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGctc                                                                                                                                           >  1:5150678/1‑150 (MQ=255)

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

           gAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGTGAGTTCAGGGACTTTTGAATGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGctctc                                                                                                                                         >  1:6158151/1‑150 (MQ=255)

           gAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCACCGctctc                                                                                                                                         >  1:5712033/1‑150 (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)

                    cATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGAGCTTCGCCCTTCGGGTCGTTGCCTGAGGCAAAGCTCTCTCGCTGgcg                                                                                                                                >  2:7099251/1‑150 (MQ=255)

                     aTCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGATAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAGTCAGAGCTTAGAGCATCTCCCATCGGGTAGATTCCTGTGGCAAAGGACTCTCGCttaccc                                                                                                                               >  2:430408/1‑145 (MQ=14)

                     aTCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGgcgc                                                                                                                               >  2:1732415/1‑150 (MQ=255)

                         ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGa                                                                                                                           <  1:9076721/150‑1 (MQ=255)

                          tCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCCACCCGCGCACGCGGGcgcgcgcg                                                                                                                          >  1:803652/1‑148 (MQ=14)

                           ccccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGATTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGt                                                                                                                         <  1:2419859/150‑1 (MQ=255)

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

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

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

                                  tGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCt                                                                                                                  >  1:2148590/1‑150 (MQ=255)

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

                                     gCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCTCTTCTCCCTTCGGGTCGTTGCCTGCGGCAACGATCTCTCGCTGTCGCTCAACGCGAACCATgc                                                                                                               >  2:3606602/1‑149 (MQ=14)

                                     gCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCCCATCGCCCTTCGGGTCGTTGCCGGCCCCAAAGCTCTCTCGCTCGCCCGCCACTCGAACCTTgg                                                                                                               >  2:6362844/1‑150 (MQ=14)

                                     gCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACACAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTgg                                                                                                               >  2:1420472/1‑150 (MQ=255)

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

                                        gAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCTCGCTGGCGCTCGAGTCGAACCTTGGTCg                                                                                                            >  2:5700578/1‑150 (MQ=255)

                                         aaTATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCGCGAGTCGAACCGTGGTCGa                                                                                                           >  2:1875581/1‑150 (MQ=255)

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

                                                    gCGGATTCGCTTGAGAGTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATc                                                                                                <  2:2019399/150‑1 (MQ=255)

                                                    gCGGATTCGATTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATc                                                                                                <  2:1855315/150‑1 (MQ=255)

                                                     cgtaTTCGATTGAGATTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATcc                                                                                               <  2:2155834/147‑1 (MQ=255)

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

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

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

                                                         ttCGCTTGAGAGTTCAGGGACTTTCGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGAGGCAACGCTCACTCGCTGGCGCTCGAGTCGAACCATTGTCGAAGCTTCTCAAccaaac                                                                                           >  1:7603849/1‑146 (MQ=14)

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

                                                             cTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAATCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGATATCTCGCAGTCGCTCGAGTCGAAACTTGGTCGAAGCTTCTCAACCATaccccc                                                                                       >  1:1449013/1‑148 (MQ=14)

                                                             cTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAATCAGCGCTGCGCGCTTCGCCCTTCGGGACGTTGCCAGCGGCAACGCTCTCACGCTGTCGCTCGAGTCGAACCTTGGTAGAAGCTTCTCATCCTTCCCCGc                                                                                       >  1:5923017/1‑150 (MQ=14)

                                                              ttGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCCACGCTCACACGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTATCAACCATCCCCCCa                                                                                      >  1:1891515/1‑150 (MQ=14)

                                                              ttGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCa                                                                                      >  2:1434336/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                                             cTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTa                                                                       >  2:6826048/1‑150 (MQ=38)

                                                                              ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTAGAATCTTCTCATACTTCCCCGCAAGGGCAGAATATTTaa                                                                      >  2:3864790/1‑150 (MQ=255)

                                                                              ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAAAGCACTCTCGCGGGCGCGCGAGTCGTACCAGGGTCGAATCATCTCAAGCATCACCGCAAGGGCCGAATAgtttt                                                                      >  1:1555171/1‑148 (MQ=11)

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

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

                                                                                      gTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGGGGCAACCCTCTCTCGCAGGCGCGCGAATCGAACCATGGTCGAAGCTTCTCAAACATCCCCCCAAGGTCCGAATATGTAAATTCgggt                                                              >  1:6577678/1‑148 (MQ=11)

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

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

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

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

                                                                                                     gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCa                                               <  1:4413030/150‑1 (MQ=255)

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

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

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

                                                                                                            ttcTCAGCGCGGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACtt                                        <  1:78239/148‑1 (MQ=255)

                                                                                                              cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCCCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACtttt                                      <  1:641242/150‑1 (MQ=255)

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

                                                                                                               tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTg                                     >  2:2853778/1‑150 (MQ=255)

                                                                                                                 aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCTGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGAGTCGTTGGGAAGTTCAGGGACTTTTGaa                                   >  2:2448177/1‑150 (MQ=255)

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

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

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

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

                                                                                                                        cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtg                            <  2:6428492/150‑1 (MQ=255)

                                                                                                                           cggtCGCCCGTCGGCTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtg                         <  1:4958282/147‑1 (MQ=255)

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

                                                                                                                                gCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTggg                    <  1:5030192/150‑1 (MQ=255)

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

                                                                                                                                       agagcGTTGCCGCAGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGa                              <  2:2618153/129‑1 (MQ=255)

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

                                                                                                                                          tCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTc          <  1:1420472/150‑1 (MQ=255)

                                                                                                                                          ccgcggCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGCACCAGCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAAGTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTc          <  2:1265975/145‑1 (MQ=14)

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

                                                                                                                                           cGTTGCCGCAGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATACTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCg         >  1:5831416/1‑150 (MQ=255)

                                                                                                                                             ttGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGaa       >  2:6868792/1‑150 (MQ=255)

                                                                                                                                             ttGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACttt                                       <  2:3486924/118‑1 (MQ=255)

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

                                                                                                                                                ccTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCCGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGTAAGTGATGGTGGTGGGGGAAGGATTCGAACCt    >  1:5766860/1‑150 (MQ=255)

                                                                                                                                                ccTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCt    >  2:1907813/1‑150 (MQ=255)

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

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

                                                                                                                                                  tGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTc  >  2:1667984/1‑150 (MQ=255)

                                                                                                                                                  tGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCATCTAAAACATACCCCCAAGGGCAGAATATTTAAATTCGGGTTAGTTGGGAAGTTAATGGACTTTTGAAAGTGATTGTTGTGGGGGGAGGATTAGaaacaac  >  1:7247849/1‑145 (MQ=12)

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

                                                                                                                                                    cGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGt                                                  <  1:2866983/100‑1 (MQ=255)

                                                                                                                                                    cGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGt                                                  >  2:2971193/1‑100 (MQ=255)

GAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTC  >  NC_000913/1286976‑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: `‑`