BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	339,498	T→C	V37A (GTG→GCG)	yahI →	putative kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	339,498	0	T	C	99.6%	871.2 / ‑4.9	276	V37A (GTG→GCG)	yahI	putative kinase
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (147/127); total (148/128)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.40e-01

Only 100 of 286 total aligned reads displayed.

TTGCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGTGCTGGAAATGCTGGCTTCCGATTACGACATTGTGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGCCGTGCGGAGATT  >  NC_000913/339408‑339586

ttGCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGgcgc                                                                                        >  2:7730401/1‑93 (MQ=255)

  gCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTg                                                                                      >  1:10006843/1‑93 (MQ=255)

   ccATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTgg                                                                                     >  2:5424028/1‑93 (MQ=255)

    cATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGa                                                                                    >  1:754891/1‑93 (MQ=255)

    cATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGa                                                                                    <  2:6349858/93‑1 (MQ=255)

     aTTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGaa                                                                                   >  1:8192639/1‑93 (MQ=255)

      ttGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGaaa                                                                                  <  1:7885534/93‑1 (MQ=255)

      ttGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGaaa                                                                                  >  1:5953267/1‑93 (MQ=255)

       tggtggCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAAt                                                                                 <  1:3864931/93‑1 (MQ=255)

          tggCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCt                                                                              <  1:2095720/93‑1 (MQ=255)

          tggAAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGTACTGCt                                                                              >  1:9649129/1‑93 (MQ=255)

           ggCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTg                                                                             >  1:9305022/1‑93 (MQ=255)

           ggCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTg                                                                             >  1:7083237/1‑93 (MQ=255)

           ggCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATAGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTg                                                                             <  2:8945590/93‑1 (MQ=255)

            gCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTgg                                                                            <  1:1382090/93‑1 (MQ=255)

             cAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGGTGGc                                                                           >  1:969696/1‑93 (MQ=255)

             cAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGc                                                                           >  2:7845481/1‑93 (MQ=255)

             cAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGc                                                                           <  2:1911451/93‑1 (MQ=255)

              aaCAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCt                                                                          >  1:9108891/1‑93 (MQ=255)

              aaCAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCt                                                                          >  1:4290738/1‑93 (MQ=255)

              aaCAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCt                                                                          >  1:1890794/1‑93 (MQ=255)

               aCAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCtt                                                                         >  2:2487097/1‑93 (MQ=255)

                 aGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTcc                                                                       <  1:5042489/93‑1 (MQ=255)

                 aGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTcc                                                                       >  2:1219964/1‑93 (MQ=255)

                 aGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTcc                                                                       <  2:8108082/93‑1 (MQ=255)

                  gCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCg                                                                      <  1:6150885/93‑1 (MQ=255)

                  gCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCg                                                                      >  2:6712544/1‑93 (MQ=255)

                   cATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGa                                                                     >  1:568590/1‑93 (MQ=255)

                    aTTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGAt                                                                    >  1:4542459/1‑93 (MQ=255)

                     ttATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGAtt                                                                   >  1:3644474/1‑93 (MQ=255)

                     ttATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGAtt                                                                   <  2:7897532/93‑1 (MQ=255)

                      tATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGAttc                                                                  >  2:8661136/1‑92 (MQ=255)

                      tATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTa                                                                  >  1:6017840/1‑93 (MQ=255)

                      tATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTa                                                                  <  2:3508167/93‑1 (MQ=255)

                      tATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTa                                                                  >  2:4879232/1‑93 (MQ=255)

                         cAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGa                                                               <  1:5324972/93‑1 (MQ=255)

                          aaaGATAACGCCAGCGAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGAc                                                              >  2:9919813/1‑93 (MQ=255)

                           aaGATACCGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACa                                                             <  2:7115748/93‑1 (MQ=255)

                            aGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACAt                                                            <  2:818725/93‑1 (MQ=255)

                              aTAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATtg                                                          >  1:6602467/1‑93 (MQ=255)

                               tAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATtgc                                                         <  2:2645253/93‑2 (MQ=255)

                               tAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATtgc                                                         >  1:3486463/1‑92 (MQ=255)

                               tAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATtgc                                                         >  1:2638648/1‑92 (MQ=255)

                               caaCGCCAGACAGTCGATTGAGCATCAGGAGGAGGCGGGGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATtgc                                                         <  1:8112094/92‑2 (MQ=255)

                                aaCGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATtgcg                                                        <  2:4089375/93‑3 (MQ=255)

                                  ctccagccagTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCt                                                      <  2:3849266/91‑1 (MQ=255)

                                  cGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCt                                                      <  1:2934849/93‑1 (MQ=255)

                                  cGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCt                                                      <  1:6242125/93‑1 (MQ=255)

                                   gccagccaGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTg                                                     >  2:1612553/1‑93 (MQ=255)

                                    ccagccagTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGa                                                    >  1:7764810/1‑93 (MQ=255)

                                    ccagccagTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGa                                                    <  2:5575575/93‑1 (MQ=255)

                                    ccagccagTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGa                                                    >  2:360435/1‑93 (MQ=255)

                                     cagccagTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGAc                                                   >  2:5774782/1‑93 (MQ=255)

                                      agccagTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGAcc                                                  <  1:7136411/93‑1 (MQ=255)

                                      agccagTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGAcc                                                  <  2:9614039/93‑1 (MQ=255)

                                        ccagTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCa                                                <  2:4522364/93‑1 (MQ=255)

                                         cagTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCAc                                               <  2:4687770/93‑1 (MQ=255)

                                         cagTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCAc                                               <  1:1161889/93‑1 (MQ=255)

                                           gTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGCTTACGACATTGCGCGGACCCCCgc                                             >  1:551959/1‑92 (MQ=255)

                                           gTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACgg                                             <  2:7321259/93‑1 (MQ=255)

                                           gTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACgg                                             >  2:6796039/1‑93 (MQ=255)

                                           gTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACgg                                             >  2:9228389/1‑93 (MQ=255)

                                           gTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACgg                                             >  1:1505167/1‑93 (MQ=255)

                                           gTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATCACGACATTGCGCTGACCCACgg                                             >  1:7005181/1‑93 (MQ=255)

                                             cGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCa                                           >  2:7002118/1‑93 (MQ=255)

                                              gATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCaa                                          <  2:6546125/93‑1 (MQ=255)

                                                ttGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACg                                        <  1:6335406/93‑1 (MQ=255)

                                                ttGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACg                                        >  2:4792948/1‑93 (MQ=255)

                                                 tGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACgg                                       <  2:5262856/93‑1 (MQ=255)

                                                 tGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACgg                                       <  1:3157448/93‑1 (MQ=255)

                                                 tGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACgg                                       >  2:6663526/1‑93 (MQ=255)

                                                   aGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGc                                     <  1:8151133/93‑1 (MQ=255)

                                                      aTCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGc                                  <  1:6379272/93‑1 (MQ=255)

                                                       tCAGGCGGAGGCGGTGAAAGCCTTCGCCGATACGGGGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGcg                                 >  1:1816865/1‑92 (MQ=255)

                                                       tCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCa                                 <  1:4035935/93‑1 (MQ=255)

                                                         aggcggaggcggGCAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAgg                               <  1:2872770/93‑1 (MQ=255)

                                                           gcggaggcggTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTc                             >  2:4359531/1‑93 (MQ=255)

                                                            gggaggcggTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCg                            <  2:2926984/92‑1 (MQ=255)

                                                            cggaggcggTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCg                            <  1:4641006/93‑1 (MQ=255)

                                                                ggcggTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGAGCTGACCCACGGCAACGGGCCGCAGGTCGGGCt                        >  1:9543955/1‑93 (MQ=255)

                                                                   ggTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGATTCCGATTACGACATTGCGATGACCCACGGCAACGGGCAGCAGGGCGGGCTGGa                     >  2:2233440/1‑93 (MQ=255)

                                                                    gTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGAt                    >  2:3379065/1‑93 (MQ=255)

                                                                       aaaGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTa                 <  1:9960121/93‑1 (MQ=255)

                                                                        aaGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTAc                >  1:2370474/1‑93 (MQ=255)

                                                                         aGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACg               >  1:5192316/1‑93 (MQ=255)

                                                                          gCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGc              <  2:9694994/93‑1 (MQ=255)

                                                                           ccGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGcc             <  1:2487366/93‑1 (MQ=255)

                                                                             gTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGCCGt           <  2:2296249/93‑1 (MQ=255)

                                                                             gTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGCCGt           <  1:8345204/93‑1 (MQ=255)

                                                                             gTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGCCGt           >  2:1986675/1‑93 (MQ=255)

                                                                              tCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGCCGTg          >  2:6967642/1‑93 (MQ=255)

                                                                               cGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGCCGTGc         <  1:8933346/93‑1 (MQ=255)

                                                                                gCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGCCGTGCg        <  1:2218260/93‑1 (MQ=255)

                                                                                gCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGCCGTGCg        >  1:7291038/1‑93 (MQ=255)

                                                                                 ccGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGCCGTGCgg       <  2:1739939/93‑1 (MQ=255)

                                                                                  cGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGCCGTGCGga      >  2:2831642/1‑93 (MQ=255)

                                                                                   gATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGCCGTGCGgag     >  1:8530302/1‑93 (MQ=255)

                                                                                    aTACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGCCGTGCGgaga    >  2:6057312/1‑93 (MQ=255)

                                                                                     tACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGCCGTGCGGAGAt   >  1:733937/1‑93 (MQ=255)

                                                                                      aCGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGCCGTGCGGAGAtt  >  2:6084718/1‑93 (MQ=255)

TTGCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGTGCTGGAAATGCTGGCTTCCGATTACGACATTGTGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGCCGTGCGGAGATT  >  NC_000913/339408‑339586

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