BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	386,819	(G)_5→6	coding (625/828 nt)	tauC →	TauC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	386,814	1	.	G	99.5%	594.9 / ‑3.8	185	coding (620/828 nt)	tauC	TauC
Reads supporting (aligned to +/- strand): ref base . (1/0); new base G (88/96); total (89/96)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

Only 100 of 208 total aligned reads displayed.

CCAGTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCT‑GGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCGGGTGAAT  >  NC_000913/386722‑386900

ccAGTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCt                                                                                         <  1:3795260/93‑1 (MQ=255)

 cAGTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCT‑g                                                                                       >  1:5894902/1‑93 (MQ=255)

 cAGTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCT‑g                                                                                       <  1:4858602/93‑1 (MQ=255)

 cAGTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCT‑g                                                                                       >  2:5713804/1‑93 (MQ=255)

 cAGTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCT‑g                                                                                       <  1:2583074/93‑1 (MQ=255)

  aGTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCT‑gg                                                                                      >  2:2620123/1‑93 (MQ=255)

  aGTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCT‑gg                                                                                      >  2:1944919/1‑93 (MQ=255)

  aGTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCT‑gg                                                                                      <  1:3145034/93‑1 (MQ=255)

   gTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCT‑ggg                                                                                     <  2:1838688/93‑1 (MQ=255)

   gTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCT‑ggg                                                                                     >  1:2117750/1‑93 (MQ=255)

       cTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCT‑ggggggt                                                                                 <  2:2953506/93‑3 (MQ=255)

        tGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCT‑ggggggtg                                                                                <  2:1360326/93‑4 (MQ=255)

        tGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCT‑ggggggtg                                                                                >  2:5569618/1‑90 (MQ=255)

          ggTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCT‑Gggggggggg                                                                              >  1:2810078/1‑92 (MQ=255)

             gCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTg                                                                            >  1:2456474/1‑93 (MQ=255)

             gCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTg                                                                            >  2:3203060/1‑93 (MQ=255)

             gCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTg                                                                            >  1:2966851/1‑93 (MQ=255)

              ccAGCCGTGCGCTGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTTGGCTgg                                                                           <  1:1582113/93‑1 (MQ=255)

                aGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTc                                                                         <  2:603425/93‑1 (MQ=255)

                aGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTc                                                                         >  2:2344283/1‑93 (MQ=255)

                 gCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCt                                                                        <  1:3121796/93‑1 (MQ=255)

                  ccGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTa                                                                       <  1:4563765/93‑1 (MQ=255)

                   cGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTAc                                                                      <  1:6201880/93‑1 (MQ=255)

                   cGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTAc                                                                      <  1:5454948/93‑1 (MQ=255)

                    gTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACg                                                                     >  1:1260876/1‑93 (MQ=255)

                     tGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGc                                                                    >  1:2288646/1‑93 (MQ=255)

                         cAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCtggt                                                                <  2:548651/93‑1 (MQ=255)

                          aGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCtggtg                                                               >  1:4896888/1‑93 (MQ=255)

                           ggTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCtggtgg                                                              >  2:5659025/1‑93 (MQ=255)

                            gTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTggc                                                             >  2:1108645/1‑93 (MQ=255)

                            gTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGCGGGCTGGTCTACGCTGGTggc                                                             >  2:4067188/1‑93 (MQ=255)

                             tGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTggcg                                                            <  2:6277298/93‑1 (MQ=255)

                              gCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTggcgg                                                           >  2:1204816/1‑93 (MQ=255)

                              gCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTggcgg                                                           >  2:5388742/1‑93 (MQ=255)

                               cTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTggcggc                                                          >  2:3472493/1‑93 (MQ=255)

                                tgtgGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGgcggcg                                                         >  1:1675962/1‑93 (MQ=255)

                                 gtgGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGcggcgg                                                        >  2:1528275/1‑93 (MQ=255)

                                  tggttggtCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGa                                                       <  2:2830137/93‑1 (MQ=255)

                                     tttGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCt                                                    >  1:4743255/1‑93 (MQ=255)

                                       tGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGa                                                  <  2:420913/93‑1 (MQ=255)

                                        gTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGAt                                                 >  1:761718/1‑93 (MQ=255)

                                        gTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGAt                                                 <  1:462752/93‑1 (MQ=255)

                                         tCATTTTGCCCGGTGCGCTGCCGGAAATGCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGAtt                                                <  1:6058055/93‑1 (MQ=255)

                                         tCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGAtt                                                <  2:3383034/93‑1 (MQ=255)

                                         tCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGAtt                                                >  1:3664268/1‑93 (MQ=255)

                                           aTTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGc                                              <  1:3397698/93‑1 (MQ=255)

                                            ttttGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGcc                                             >  2:736580/1‑93 (MQ=255)

                                             tttGCCCGGTGCGCTGCCGGAAATCCTCGCCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCcg                                            <  1:2337436/93‑1 (MQ=255)

                                             tttGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCcg                                            <  1:949870/93‑1 (MQ=255)

                                                gCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGa                                         >  1:182240/1‑93 (MQ=255)

                                                 cccGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGAc                                        <  2:5258652/93‑1 (MQ=255)

                                                  ccGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGAcg                                       >  1:3210433/1‑93 (MQ=255)

                                                   cGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGAcgc                                      <  1:143136/93‑1 (MQ=255)

                                                    ggTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGTCGCGAcgcg                                     >  2:5048750/1‑93 (MQ=255)

                                                     gTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGAcgcgc                                    >  2:835801/1‑93 (MQ=255)

                                                       gcgcTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCgg                                  <  2:2687169/93‑1 (MQ=255)

                                                         gcTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGtt                                >  1:5880508/1‑93 (MQ=255)

                                                         gcTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGtt                                <  2:1815871/93‑1 (MQ=255)

                                                          cTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGttt                               <  1:1135789/93‑1 (MQ=255)

                                                          cTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGgtt                               >  2:1448786/1‑93 (MQ=255)

                                                           tGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTa                              <  2:4595321/93‑1 (MQ=255)

                                                           tGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTa                              <  2:3833245/93‑1 (MQ=255)

                                                            gCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAg                             >  2:1021200/1‑93 (MQ=255)

                                                              cGGAAATCGTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGa                           <  2:5981899/93‑1 (MQ=255)

                                                              cGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGa                           >  2:3693373/1‑93 (MQ=255)

                                                              cGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGa                           <  1:220241/93‑1 (MQ=255)

                                                               ggAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGACTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGAt                          >  1:2642203/1‑93 (MQ=255)

                                                                gAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGAtt                         >  1:1555953/1‑93 (MQ=255)

                                                                gAAATCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGAtt                         <  2:3661179/93‑1 (MQ=255)

                                                                 aaaTCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGAttt                        >  1:4849083/1‑93 (MQ=255)

                                                                 aaaTCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGAttt                        <  1:5022565/93‑1 (MQ=255)

                                                                  aaTCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTa                       <  2:3797618/93‑1 (MQ=255)

                                                                   aTCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTAt                      <  1:1079047/93‑1 (MQ=255)

                                                                   aTCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTAt                      <  1:4016032/93‑1 (MQ=255)

                                                                   aTCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTAt                      <  2:4397998/93‑1 (MQ=255)

                                                                   aTCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTAt                      >  1:59492/1‑93 (MQ=255)

                                                                    tCCTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTATg                     >  2:741509/1‑93 (MQ=255)

                                                                    tCCTCACCCGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTATg                     >  2:6155147/1‑93 (MQ=255)

                                                                     ccTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTATgg                    <  1:5637508/93‑1 (MQ=255)

                                                                     ccTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTATgg                    <  1:3311928/93‑1 (MQ=255)

                                                                      cTCACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTATGGt                   >  2:209268/1‑93 (MQ=255)

                                                                        cACCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTtc                 <  1:4276347/93‑1 (MQ=255)

                                                                         aCCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTtca                >  2:3520970/1‑93 (MQ=255)

                                                                         aCCGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTtca                <  1:1527606/93‑1 (MQ=255)

                                                                           cGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTtcagt              <  1:5271781/93‑1 (MQ=255)

                                                                           cGGATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTtcagt              <  2:1465336/93‑1 (MQ=255)

                                                                             gATTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTtcagtca            >  2:3067044/1‑93 (MQ=255)

                                                                              aTTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTtcagtcag           <  2:1566361/93‑1 (MQ=255)

                                                                              aTTACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTtcagtcag           >  2:4872018/1‑93 (MQ=255)

                                                                               ttACGTGTTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGCGCTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGc          >  1:5792145/1‑93 (MQ=255)

                                                                               ttACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCTACGCGCGGTTTAGGATTTATGGTTCAGTCAGc          >  1:5130663/1‑93 (MQ=255)

                                                                               ttACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGc          <  2:446729/93‑1 (MQ=255)

                                                                                tACGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCg         >  2:943366/1‑93 (MQ=255)

                                                                                  cGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCggg       <  1:3619183/93‑1 (MQ=255)

                                                                                  cGTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCggg       >  1:5534145/1‑93 (MQ=255)

                                                                                   gTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCGGGt      <  1:715143/93‑1 (MQ=255)

                                                                                   gTATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCGGGt      <  2:1729147/93‑1 (MQ=255)

                                                                                    tATTGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCGGGTg     >  2:2377437/1‑93 (MQ=255)

                                                                                      ttGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCGGGTGaa   <  1:2223673/93‑1 (MQ=255)

                                                                                       tGGTCTGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCGGGTGAAt  <  1:3019960/93‑1 (MQ=255)

CCAGTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCT‑GGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCGGGTGAAT  >  NC_000913/386722‑386900

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