BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,554,118	+C	coding (367/396 nt)	yhgKJ ←	RNA phosphage cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,554,117	1	.	C	98.8%	582.6 / ‑1.2	174	coding (368/396 nt)	yhgKJ	RNA phosphage cyclase
Reads supporting (aligned to +/- strand): ref base . (0/2); new base C (95/76); total (96/78)
Fisher's exact test for biased strand distribution p-value = 2.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 186 total aligned reads displayed.

GGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCG‑AGGCGGCGCCGACGGGTTATCGGTCGCCG‑CGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGGC  >  NC_000913/3554027‑3554199

ggTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑cc                                                                                  >  1:2011230/1‑92 (MQ=255)

  tttGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑ccgc                                                                                <  2:535334/93‑4 (MQ=255)

   ttGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑ccgct                                                                               <  1:1096755/93‑5 (MQ=255)

     gCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCa                                                                              >  1:2641241/1‑93 (MQ=255)

       tgatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACt                                                                            >  2:3942883/1‑93 (MQ=255)

       tgatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACt                                                                            <  2:1900543/93‑1 (MQ=255)

       tgatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACt                                                                            <  1:4020427/93‑1 (MQ=255)

       tgatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACt                                                                            <  2:2460881/93‑1 (MQ=255)

        gatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACtt                                                                           >  1:2759401/1‑93 (MQ=255)

         atgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTc                                                                          <  2:4993785/93‑1 (MQ=255)

            aaTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAac                                                                       >  2:6135702/1‑93 (MQ=255)

            aaTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAac                                                                       >  1:3305822/1‑93 (MQ=255)

             aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAaca                                                                      <  1:6414941/93‑1 (MQ=255)

             aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAaca                                                                      <  2:5818638/93‑1 (MQ=255)

              ttCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAacac                                                                     >  1:5463871/1‑93 (MQ=255)

                ccTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACAcgc                                                                   >  2:4603073/1‑93 (MQ=255)

                 cTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACTGGTTATCCGT‑GCCGCCGCTCACTTCAGCAcgcg                                                                  >  2:3957892/1‑93 (MQ=255)

                 cTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACAcgcg                                                                  >  1:1351929/1‑93 (MQ=255)

                 cTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACAcgcg                                                                  >  1:2797909/1‑93 (MQ=255)

                  tATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGa                                                                 >  1:5922204/1‑93 (MQ=255)

                   aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCTC‑CGCTCACTTCAACACGCGaa                                                                <  1:5631147/93‑1 (MQ=255)

                   aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGaa                                                                <  2:2286386/93‑1 (MQ=255)

                    ttttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAg                                                               >  1:489165/1‑93 (MQ=255)

                      ttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGt                                                             <  1:5496370/93‑1 (MQ=255)

                      ttCCCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGt                                                             <  1:5097560/93‑1 (MQ=255)

                        cgcgAGCAGCGGCTCCAGCACCCGGCGGATAAAGTCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTcc                                                           <  2:1154584/93‑1 (MQ=255)

                        cGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTcc                                                           >  1:3394840/1‑93 (MQ=255)

                        cGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTcc                                                           >  1:5799561/1‑93 (MQ=255)

                        cGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTcc                                                           >  2:1033357/1‑93 (MQ=255)

                         gCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCa                                                          >  1:1951941/1‑93 (MQ=255)

                          ccAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCAt                                                         >  2:1957928/1‑93 (MQ=255)

                           cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATc                                                        <  1:930756/93‑1 (MQ=255)

                           cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATc                                                        <  1:3520429/93‑1 (MQ=255)

                            agcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCg                                                       >  2:4802773/1‑93 (MQ=255)

                            agcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCg                                                       >  2:1645717/1‑93 (MQ=255)

                            agcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCg                                                       >  1:1501884/1‑93 (MQ=255)

                                gcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaa                                                   >  2:6012811/1‑93 (MQ=255)

                                gcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaa                                                   >  2:1881167/1‑93 (MQ=255)

                                gcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGAGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGGCCATCGGCaa                                                   >  1:4295204/1‑93 (MQ=255)

                                  ggCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAc                                                 <  2:2131388/93‑1 (MQ=255)

                                  gccTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAc                                                 <  2:6161572/91‑1 (MQ=255)

                                   gCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAcc                                                >  2:2822095/1‑93 (MQ=255)

                                   gCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAcc                                                >  2:275652/1‑93 (MQ=255)

                                   gCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTAACGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAcc                                                >  1:3205282/1‑93 (MQ=255)

                                    cTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACca                                               <  2:4678407/93‑1 (MQ=255)

                                     tCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcac                                              <  1:2941674/93‑1 (MQ=255)

                                       cAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAg                                            >  1:1318593/1‑93 (MQ=255)

                                       cAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGAT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAg                                            <  1:5226284/93‑1 (MQ=255)

                                        atcACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgc                                           <  1:1990190/91‑1 (MQ=255)

                                        aGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgc                                           >  1:3513141/1‑93 (MQ=255)

                                         gCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgcg                                          >  2:3821725/1‑93 (MQ=255)

                                         gCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgcg                                          <  2:636153/93‑1 (MQ=255)

                                         gCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgcg                                          >  1:3263347/1‑93 (MQ=255)

                                          cACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgcgc                                         <  2:6393934/93‑1 (MQ=255)

                                           ccccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGcgcg                                        <  2:457501/92‑1 (MQ=255)

                                            cccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCgg                                       <  1:5838846/93‑1 (MQ=255)

                                             ccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCggg                                      >  2:4019838/1‑93 (MQ=255)

                                              cggcggATAAAACCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgc                                     <  1:3910415/93‑1 (MQ=255)

                                               ggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgca                                    >  1:5249566/1‑93 (MQ=255)

                                               ggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgca                                    >  1:4344198/1‑93 (MQ=255)

                                               ggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgca                                    <  1:2042344/93‑1 (MQ=255)

                                               ggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCAATTCAACACGCGAAGGTCCATCGGCAAACCACACCGCGGgca                                    >  2:273109/1‑93 (MQ=255)

                                                  ggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagca                                 <  2:2794805/93‑1 (MQ=255)

                                                  ggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagca                                 >  2:318246/1‑93 (MQ=255)

                                                  ggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagca                                 <  2:4691646/93‑1 (MQ=255)

                                                  ggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagca                                 >  1:5268788/1‑93 (MQ=255)

                                                   gATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc                                >  2:845697/1‑93 (MQ=255)

                                                    aTAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAcc                               <  2:5155673/93‑1 (MQ=255)

                                                      aaaaTCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGt                             <  2:21657/93‑1 (MQ=255)

                                                      aaaaTCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGt                             <  2:172815/93‑1 (MQ=255)

                                                      aaaaTCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGt                             <  1:530792/93‑1 (MQ=255)

                                                      aaaaTCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGt                             >  1:5166347/1‑93 (MQ=255)

                                                      aaaaTCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGt                             >  1:4737417/1‑93 (MQ=255)

                                                       aaaTCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGtt                            >  2:3540471/1‑93 (MQ=255)

                                                             gcaGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgcagc                      <  2:5700899/91‑1 (MQ=255)

                                                             gcaGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgcagc                      >  1:4215655/3‑93 (MQ=255)

                                                               aGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCAc                    >  2:5889447/1‑93 (MQ=255)

                                                                ggcggcGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCAcc                   >  1:3804402/1‑93 (MQ=255)

                                                                ggcggcGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCAcc                   >  2:1439326/1‑93 (MQ=255)

                                                                 gcggcgCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCa                  <  1:2332030/93‑1 (MQ=255)

                                                                    gcgcCGGCGGGTTATCGGTACCGC‑CGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCg               <  1:2280773/93‑1 (MQ=255)

                                                                    gcgcCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCg               >  2:5542215/1‑93 (MQ=255)

                                                                    gcgcCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCg               >  2:5200592/1‑93 (MQ=255)

                                                                      gcCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTa             >  1:1186980/1‑93 (MQ=255)

                                                                       ccGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTAc            <  1:1276266/93‑1 (MQ=255)

                                                                        cGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACa           >  2:6310377/1‑93 (MQ=255)

                                                                         gACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACaa          >  1:2505385/1‑93 (MQ=255)

                                                                         gACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACaa          >  2:2864187/1‑93 (MQ=255)

                                                                          aCGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAAc         <  1:4035478/93‑1 (MQ=255)

                                                                          aCGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAAc         <  2:252045/93‑1 (MQ=255)

                                                                           cGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACt        >  1:611779/1‑93 (MQ=255)

                                                                            gggTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACtt       <  1:3363922/93‑1 (MQ=255)

                                                                             ggTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTc      >  2:1709975/1‑93 (MQ=255)

                                                                             ggTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTc      <  2:3442954/93‑1 (MQ=255)

                                                                             ggTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTc      >  2:1972373/1‑93 (MQ=255)

                                                                              gTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTcc     <  1:4590381/93‑1 (MQ=255)

                                                                               ttATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCg    <  2:5454575/93‑1 (MQ=255)

                                                                               ttATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCg    >  1:3992145/1‑93 (MQ=255)

                                                                                tATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCgg   >  1:177774/1‑93 (MQ=255)

                                                                                 aTCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGgc  >  2:259483/1‑93 (MQ=255)

GGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCG‑AGGCGGCGCCGACGGGTTATCGGTCGCCG‑CGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGGC  >  NC_000913/3554027‑3554199

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 19 ≤` `ATCG/ATCG` `< 26 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`