BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,277,940	A→C	5.4%	noncoding (127/140 nt)	ryjA ←	small RNA RyjA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,277,940	0	A	C	5.4%	339.1 / 9.2	112	noncoding (127/140 nt)	ryjA	small RNA RyjA
Reads supporting (aligned to +/- strand): ref base A (50/56); new base C (3/3); total (53/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

Only 100 of 113 total aligned reads displayed.

TGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCACCGACG  >  NC_000913/4277797‑4278089

tGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAggg                                                                                                                                                  >  1:553689/1‑150 (MQ=255)

tGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAggg                                                                                                                                                  >  2:624633/1‑150 (MQ=255)

 ggCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGgc                                                                                                                                                 >  2:332336/1‑150 (MQ=255)

   cgcTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGgcgc                                                                                                                                               >  1:691180/1‑150 (MQ=255)

   cgcTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCACGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGATGGTACCGGCGCACGCTTGCTGGTAGATGAACAAAACTAAAGCGCCACAAGGgcgc                                                                                                                                               <  2:783823/150‑1 (MQ=255)

             cGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAg‑tttgt                                                                                                                                    >  2:1074858/1‑150 (MQ=255)

             cGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAAC‑AAACTAAAGCGCCCTTGTGGCGCTTTAg‑ttttgt                                                                                                                                   <  1:1453085/150‑4 (MQ=38)

                aCTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCCTTGTGGCGCTTTAg‑ttttgttc                                                                                                                                 <  2:1477169/150‑6 (MQ=39)

                    gACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCgg                                                                                                                             >  1:567173/1‑150 (MQ=255)

                      cGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGgtc                                                                                                                           <  1:188735/150‑1 (MQ=255)

                        gATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGCCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGgtctt                                                                                                                         >  2:680183/1‑150 (MQ=255)

                        gATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGgtctt                                                                                                                         <  2:1469600/150‑1 (MQ=255)

                        gATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGGTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGgtctt                                                                                                                         <  2:162397/150‑1 (MQ=255)

                          tGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGgtctttg                                                                                                                       <  1:589381/150‑1 (MQ=255)

                            tATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGgtctttgtc                                                                                                                     >  1:72779/1‑150 (MQ=255)

                             aTTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGgtctttgtct                                                                                                                    >  1:1230383/1‑150 (MQ=255)

                                ggTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGgtctttgtcttta                                                                                                                 >  2:1508/1‑149 (MQ=255)

                                ggTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTtc                                                                                                                 <  2:1213623/150‑1 (MQ=255)

                                ggTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTtc                                                                                                                 >  1:1059474/1‑150 (MQ=255)

                                ggTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTtc                                                                                                                 >  1:1252264/1‑150 (MQ=255)

                                    gtgGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTtctctc                                                                                                             <  2:840263/150‑1 (MQ=255)

                                     tgGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAACCTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTctctct                                                                                                            <  2:643491/150‑1 (MQ=255)

                                      ggCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTa                                                                                                           >  2:131209/1‑150 (MQ=255)

                                      ggCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTa                                                                                                           >  2:1270205/1‑150 (MQ=255)

                                        cTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATc                                                                                                         <  1:288875/150‑1 (MQ=255)

                                             tttCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCt                                                                                                    <  2:422666/150‑1 (MQ=255)

                                              ttCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTg                                                                                                   <  2:567173/150‑1 (MQ=255)

                                               tCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTgg                                                                                                  <  2:1230383/150‑1 (MQ=255)

                                                    cAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTacac                                                                                             >  2:830996/1‑150 (MQ=255)

                                                      gTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACag                                                                                           <  2:1075078/150‑1 (MQ=255)

                                                          ttGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACaggagg                                                                                       <  2:799668/150‑1 (MQ=255)

                                                          ttGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGGCTTTGTCTTTCTCTCTATCCCGCTGGTACACaggagg                                                                                       <  1:154404/150‑1 (MQ=255)

                                                            gCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGt                                                                                     >  1:1466887/1‑150 (MQ=255)

                                                             cccGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGtt                                                                                    <  1:1074858/150‑1 (MQ=255)

                                                               cGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTc                                                                                  >  1:1151930/1‑150 (MQ=255)

                                                                  tGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTcccc                                                                               >  1:1207831/1‑150 (MQ=255)

                                                                  tGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTcccc                                                                               <  1:709620/150‑1 (MQ=255)

                                                                    cGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCg                                                                             <  2:1059474/150‑1 (MQ=255)

                                                                         cccGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTc                                                                        >  1:957210/1‑150 (MQ=255)

                                                                         cccGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTc                                                                        >  1:1509522/1‑150 (MQ=255)

                                                                          ccGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCa                                                                       >  1:607199/1‑150 (MQ=255)

                                                                           cGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCaa                                                                      <  2:145557/150‑1 (MQ=255)

                                                                            gggCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAac                                                                     <  1:412727/150‑1 (MQ=255)

                                                                              gCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAacac                                                                   <  1:1337404/150‑1 (MQ=255)

                                                                              gCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAacac                                                                   <  1:802707/150‑1 (MQ=255)

                                                                                gACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAacacac                                                                 >  1:1183752/1‑150 (MQ=255)

                                                                                 aCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACAcc                                                                >  1:758973/1‑150 (MQ=255)

                                                                                   cGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTc                                                              >  2:638771/1‑150 (MQ=255)

                                                                                      ttAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGACCAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCAtt                                                           <  2:521494/150‑1 (MQ=255)

                                                                                       tAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTc                                                          <  2:358630/150‑1 (MQ=255)

                                                                                        aGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGTGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCg                                                         <  2:1509522/150‑1 (MQ=255)

                                                                                        aGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCg                                                         <  2:1300180/150‑1 (MQ=255)

                                                                                        aGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCg                                                         <  2:1294898/150‑1 (MQ=255)

                                                                                         ggAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGa                                                        >  2:1157890/1‑150 (MQ=255)

                                                                                          gagaAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTCGTACACAGGAGGGTTTCCCCCGACGTCAACACAACTCATTCGaa                                                       >  2:1253883/1‑149 (MQ=255)

                                                                                           agaagaAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAga                                                      >  2:796301/1‑149 (MQ=255)

                                                                                           agaagaAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGc                                                      >  2:1193517/1‑150 (MQ=255)

                                                                                              agaagGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCCTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACg                                                   >  2:509761/1‑150 (MQ=255)

                                                                                                  ggTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACgtggt                                               >  2:23327/1‑150 (MQ=255)

                                                                                                    tACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGgtctt                                                                                                                         <  1:525368/74‑1 (MQ=255)

                                                                                                    tACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGgtctt                                                                                                                         >  2:525368/1‑74 (MQ=255)

                                                                                                     aCCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGa                                            <  1:23327/150‑1 (MQ=255)

                                                                                                      ccGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAg                                           >  2:155489/1‑150 (MQ=255)

                                                                                                        ggCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGt                                         >  1:528809/1‑150 (MQ=255)

                                                                                                        ggCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGt                                         >  1:1476852/1‑150 (MQ=255)

                                                                                                             aggCTTGGTGGAAGATGAACAAAACTTAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCgg                                    <  1:1253883/148‑1 (MQ=255)

                                                                                                              cGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGt                                   >  2:143740/1‑150 (MQ=255)

                                                                                                              cGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGt                                   >  1:1387787/1‑150 (MQ=255)

                                                                                                                cTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTg                                 <  2:528809/150‑1 (MQ=255)

                                                                                                                  tGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGgt                               >  1:793012/1‑150 (MQ=255)

                                                                                                                  tGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGgt                               <  2:691180/150‑1 (MQ=255)

                                                                                                                   gCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTAGAGCACGTGGTCGAGGTACCGGTTGgtg                              >  1:616279/1‑150 (MQ=255)

                                                                                                                    cTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑GTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCAATCGAGCACGTGGTGGAGGTTCCGGTTGgtgt                             >  2:599535/1‑150 (MQ=255)

                                                                                                                     tGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGtt                            <  2:1315171/150‑1 (MQ=255)

                                                                                                                      tgcttaTGAAAAAACCTAACGCGCCACACGGGGGCTTTTG‑TTTGTTTTCGGGTC‑TTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGa                          <  1:815793/145‑1 (MQ=255)

                                                                                                                      tgAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTg                           <  2:1245035/149‑1 (MQ=255)

                                                                                                                      ggAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTg                           <  2:1245266/150‑1 (MQ=255)

                                                                                                                      ggAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTg                           <  2:1183752/150‑1 (MQ=255)

                                                                                                                        aaGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGAt                         <  1:32716/150‑1 (MQ=255)

                                                                                                                          gATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGc                       >  2:1255151/1‑150 (MQ=255)

                                                                                                                           aTGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCt                      <  1:155489/150‑1 (MQ=255)

                                                                                                                            tGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCtt                     <  2:1207831/150‑1 (MQ=255)

                                                                                                                             gAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCttt                    <  1:131209/150‑1 (MQ=255)

                                                                                                                             gAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGAACGTGGTGGAGGTTCCGGTTGGTGTTGATGCttt                    <  1:796301/150‑1 (MQ=255)

                                                                                                                              aaCAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCtttt                   >  2:70127/1‑149 (MQ=255)

                                                                                                                              aaCAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTa                   <  1:1193517/150‑1 (MQ=255)

                                                                                                                                acaaaCTAAAGCGCCCTTGTGGCGCTTTAGTTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTaa                  >  2:264472/3‑150 (MQ=38)

                                                                                                                                 caaaCTAAAGCGCCCTTGTGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAAtt                <  2:1145714/149‑1 (MQ=255)

                                                                                                                                 aaaaCTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAAtt                <  2:283005/150‑1 (MQ=255)

                                                                                                                                 aaaaCTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAAtt                >  1:1022578/1‑150 (MQ=255)

                                                                                                                                  aaaCTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTg               <  2:1252264/150‑1 (MQ=255)

                                                                                                                                   aaCTAAAGCGCCCTTGTGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGt              >  2:756158/1‑150 (MQ=255)

                                                                                                                                   aaCTAAAGCGCCCTTGTGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGt              >  2:886549/1‑150 (MQ=255)

                                                                                                                                   aaCTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGt              >  2:213209/1‑150 (MQ=255)

                                                                                                                                    aCTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTa             <  1:624633/150‑1 (MQ=255)

                                                                                                                                    aCTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTa             <  2:72779/150‑1 (MQ=255)

                                                                                                                                     cTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTAt            >  2:21615/1‑150 (MQ=255)

                                                                                                                                     cTAAAGCGCAACAAGTGCGCTGTAG‑ATTGTATTCAAGTCTTTGTCGTTATCTCAATCCCGCTGGTACACACGACGGTTTCCCCCGACGTAAACCCTCCTCATTCGAGCACGTGGTGGAGGGTCCGGTTGGTGTTGATTCTTTAATTGTAt            >  2:832572/1‑150 (MQ=255)

                                                                                                                                      tAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATg           <  1:544733/150‑1 (MQ=255)

                                                                                                                                               acaAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCACCGACg  <  2:1466887/150‑1 (MQ=255)

TGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAG‑TTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCACCGACG  >  NC_000913/4277797‑4278089

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