BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,554,113	Δ1 bp	coding (372/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,113	0	C	.	98.3%	510.2 / 1.6	119	coding (372/396 nt)	yhgKJ	RNA phosphage cyclase
Reads supporting (aligned to +/- strand): ref base C (0/0); major base . (48/68); minor base G (0/2); total (48/71)
Fisher's exact test for biased strand distribution p-value = 5.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 131 total aligned reads displayed.

TGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCG‑AGGCGGCGCCGACGGGTTATCGGTCGCCG‑CGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGG  >  NC_000913/3554026‑3554198

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

 ggTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑cc                                                                                 <  2:4439807/93‑2 (MQ=255)

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

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

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

         gatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACtt                                                                          <  1:2766718/93‑1 (MQ=255)

         gatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACtt                                                                          <  1:2766717/93‑1 (MQ=255)

         gatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCAACCGGCGGATAAAATCCGCAGGAGGCGCCGACGGGGTATCGGTGCCGG‑CGCGCACtt                                                                          >  2:458800/1‑93 (MQ=255)

          atgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGCGTATCGGTGCCGG‑CGCTCACTTc                                                                         >  2:3435127/1‑93 (MQ=255)

            gaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCaa                                                                       <  1:3828318/93‑1 (MQ=255)

            gaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCaa                                                                       >  1:3119911/1‑93 (MQ=255)

              aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATTCGCAGGCGGCGCCGACGGGTTATCGGTGGCGC‑CGCTCACTTCAaca                                                                     >  2:4091944/1‑93 (MQ=255)

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

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

              aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAaca                                                                     >  2:3283508/1‑93 (MQ=255)

                 ccTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCCCCCTTCACCAcgc                                                                  >  2:2450606/1‑93 (MQ=255)

                  cTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACAcgcg                                                                 <  1:4229768/93‑1 (MQ=255)

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

                   tATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGT‑GCCGCCGCCCACTTCAACACGCGa                                                                >  1:1242478/1‑93 (MQ=255)

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

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

                     ttttCGCCAGCAGCGGCTCCAGCACCCGGTGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAg                                                              >  2:3591060/1‑93 (MQ=255)

                     ttttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAg                                                              <  1:2846258/93‑1 (MQ=255)

                      tttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAgg                                                             <  1:2282234/93‑1 (MQ=255)

                      tttCGCCAGCAGCCGCTCCAGCCCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAgg                                                             <  1:2282232/93‑1 (MQ=255)

                       ttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGt                                                            >  1:3262255/1‑93 (MQ=255)

                        tCGCCAGCAGCGGCTCCAGCCCCCGGCGGATAAAAGCCGCAGGCGGCGCCGACGGGTGATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTc                                                           <  2:423878/93‑1 (MQ=255)

                        tCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTc                                                           <  2:2585618/93‑1 (MQ=255)

                         cGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTcc                                                          <  1:995531/93‑1 (MQ=255)

                         cGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTcc                                                          <  2:1297010/93‑1 (MQ=255)

                         cGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTcc                                                          <  1:2228602/93‑1 (MQ=255)

                          gCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCa                                                         <  2:1442330/93‑1 (MQ=255)

                           ccAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCAt                                                        <  1:161103/93‑1 (MQ=255)

                           ccAGCAGCGGCTCCAGCACCCGGAGCACAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCAt                                                        <  2:2725390/93‑1 (MQ=255)

                           ccAGCAGCGGCTCCAGCACCCCGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGTGCCAt                                                        >  1:1005292/1‑93 (MQ=255)

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

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

                              gcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCgg                                                     <  1:2613394/93‑1 (MQ=255)

                              gcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCgg                                                     >  2:211734/1‑93 (MQ=255)

                              gcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCgg                                                     <  1:4345312/93‑1 (MQ=255)

                               cagcGGCTCCAGCCCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGc                                                    <  1:419116/93‑1 (MQ=255)

                               cagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGc                                                    >  1:1814034/1‑93 (MQ=255)

                               cagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGc                                                    <  2:1326308/93‑1 (MQ=255)

                                agcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCAATCGGCa                                                   >  2:138280/1‑93 (MQ=255)

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

                                 gcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaa                                                  <  1:2525006/93‑1 (MQ=255)

                                  cGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaaa                                                 <  1:3014300/93‑1 (MQ=255)

                                    gCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAcc                                               <  2:4336504/93‑1 (MQ=255)

                                    gCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAcc                                               <  2:3824765/93‑1 (MQ=255)

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

                                       ccAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcaca                                            >  1:1388452/1‑93 (MQ=255)

                                       ccAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcaca                                            >  2:2256637/1‑93 (MQ=255)

                                        cAGCCCCCGGCGGATAACATCCGCAGGCGGCGCCGACGGGTTATCGG‑AGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAg                                           <  2:4216350/93‑1 (MQ=255)

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

                                           ctcccGGCAGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgcgc                                        <  1:3698984/91‑1 (MQ=255)

                                            aCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGcgcg                                       >  1:3513648/1‑93 (MQ=255)

                                             cccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAACCCACAGCGCgg                                      <  2:4676865/93‑1 (MQ=255)

                                               cggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgc                                    <  1:4363881/93‑1 (MQ=255)

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

                                                 gcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcag                                  <  2:2660219/93‑1 (MQ=255)

                                                   ggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACCGCGCGGgcagca                                >  2:3127015/1‑93 (MQ=255)

                                                   ggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagcc                                >  2:1433318/1‑92 (MQ=255)

                                                    gATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCCC‑CGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc                               <  1:1685437/93‑1 (MQ=255)

                                                      tAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCg                             >  1:837997/1‑93 (MQ=255)

                                                       aaaaTCCGCAGGCGGCGCCGCCGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGt                            <  2:3583466/93‑1 (MQ=255)

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

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

                                                        aaaTCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGtt                           <  2:4350820/93‑1 (MQ=255)

                                                        aaaTCCGCAGGCGGCGCCGACGGATTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGtt                           <  2:4350821/93‑1 (MQ=255)

                                                            ccgcaGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgca                       <  1:4156678/89‑1 (MQ=255)

                                                               cagtcggcgCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgcagca                    <  2:1284819/89‑1 (MQ=255)

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

                                                                aGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCAc                   <  2:3439832/93‑1 (MQ=255)

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

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

                                                                    ggcgCCGACGGGTTATCGGTGCCTC‑CGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGc               <  1:4039897/93‑1 (MQ=255)

                                                                    ggcgCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGc               <  1:4039898/93‑1 (MQ=255)

                                                                    ggcgCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGc               <  1:2769477/93‑1 (MQ=255)

                                                                    ggcgCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGc               <  1:1651159/93‑1 (MQ=255)

                                                                     gcgcCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCg              <  1:735614/93‑1 (MQ=255)

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

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

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

                                                                        ccGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTAc           >  2:3852614/1‑93 (MQ=255)

                                                                        ccGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTAc           >  2:2272133/1‑93 (MQ=255)

                                                                           aCGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAAc        >  1:3346831/1‑93 (MQ=255)

                                                                           aCGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAAc        >  1:3660836/1‑93 (MQ=255)

                                                                           aCGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAAc        >  2:2142212/1‑93 (MQ=255)

                                                                           aCGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAAc        >  2:3917258/1‑93 (MQ=255)

                                                                            cGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACt       <  2:1120925/93‑1 (MQ=255)

                                                                            cGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACt       <  2:1157397/93‑1 (MQ=255)

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

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

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

                                                                               gTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTcc    >  1:1427244/1‑93 (MQ=255)

                                                                               gTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTcc    >  1:1921360/1‑93 (MQ=255)

                                                                               gTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTcc    >  2:653365/1‑93 (MQ=255)

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

                                                                                 tATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCgg  <  1:3406815/93‑1 (MQ=255)

                                                                                 tATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCgg  <  2:456872/93‑1 (MQ=255)

TGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCG‑AGGCGGCGCCGACGGGTTATCGGTCGCCG‑CGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGG  >  NC_000913/3554026‑3554198

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