BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,588,897	A→G	Y354H (TAT→CAT)	ugpB ←	UgpB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,588,897	0	A	G	97.1%	525.8 / 2.9	170	Y354H (TAT→CAT)	ugpB	UgpB
Reads supporting (aligned to +/- strand): ref base A (2/3); new base G (93/72); total (95/75)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.78e-01

Only 100 of 176 total aligned reads displayed.

ACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTCGCGAGGAAATCGAGGAACTTCGCCA  >  NC_000913/3588807‑3588986

acggcaacggcGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCg                                                                                         >  1:535978/1‑93 (MQ=255)

  ggcaacggcGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCgct                                                                                       <  1:3788136/93‑1 (MQ=255)

    caacggcGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCgctgc                                                                                     >  2:2349314/1‑93 (MQ=255)

     aacggcGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCgctgct                                                                                    >  1:1926653/1‑93 (MQ=255)

      acggcGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCtt                                                                                   >  1:3904830/1‑93 (MQ=255)

      acggcGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCtt                                                                                   <  1:997818/93‑1 (MQ=255)

       cggcggCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCttt                                                                                  <  1:1643507/93‑1 (MQ=255)

        ggcggcTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTg                                                                                 <  2:3197348/93‑1 (MQ=255)

        ggcgccTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTg                                                                                 <  2:2481761/93‑1 (MQ=255)

           ggcTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTg                                                                              >  2:662349/1‑93 (MQ=255)

           ggcTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTg                                                                              >  1:2079220/1‑93 (MQ=255)

           ggcTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTg                                                                              <  2:2379155/93‑1 (MQ=255)

           ggcTTATTCAGCATCTGACGCGTCGCGGTATCCCCCCCTGGGTTTTTCTCTTAAAAGCCCTCCTCACGGGCCAGGCCAGGGCCTGCTTTgtta                                                                              >  1:1799506/1‑90 (MQ=255)

             cTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGAt                                                                            <  1:396483/93‑1 (MQ=255)

              ttATTCAGCATCTGACGCGTCGCTGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGAtt                                                                           <  2:2785903/93‑1 (MQ=255)

               tATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTg                                                                          <  1:4244491/93‑1 (MQ=255)

               tATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTg                                                                          <  2:551119/93‑1 (MQ=255)

                    aGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGa                                                                     >  2:2223742/1‑93 (MQ=255)

                    aGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGa                                                                     <  1:2213661/93‑1 (MQ=255)

                     gCATCTGACGCGTCGCGGTAGCCGCCCCTGGGTTTTTCTCATAAAAGCCCTACTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGAt                                                                    <  2:4508847/93‑1 (MQ=255)

                      cATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATa                                                                   >  1:1058746/1‑93 (MQ=255)

                        tCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAAc                                                                 <  2:3840047/93‑1 (MQ=255)

                         cTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAAcc                                                                <  2:4666621/93‑1 (MQ=255)

                         cTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAAcc                                                                <  2:3818254/93‑1 (MQ=255)

                         cTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAAcc                                                                >  1:696032/1‑93 (MQ=255)

                          tGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCg                                                               <  2:4497178/93‑1 (MQ=255)

                          tGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCg                                                               <  1:2283719/93‑1 (MQ=255)

                           ggcgcgTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCgg                                                              <  2:3892138/91‑1 (MQ=255)

                           gACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCgg                                                              >  1:3734398/1‑93 (MQ=255)

                           gACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCgg                                                              >  1:3734264/1‑93 (MQ=255)

                            aCGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGt                                                             <  2:2855581/93‑1 (MQ=255)

                             cgcgTCGCGGTATCCGCCCCTGGGTTTTGCTCATAAAAGTCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGtt                                                            <  1:3292665/93‑1 (MQ=255)

                              gcgTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGttt                                                           >  2:3932099/1‑93 (MQ=255)

                               cgTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGtttt                                                          <  2:1539416/93‑1 (MQ=255)

                                gTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAATCGGTTTTc                                                         >  1:2365107/1‑93 (MQ=255)

                                gTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTc                                                         <  2:928771/93‑1 (MQ=255)

                                gTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTc                                                         >  2:2468767/1‑93 (MQ=255)

                                 tCGCGGTATCCGCCCCTGGGGTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCt                                                        >  1:3497122/1‑93 (MQ=255)

                                  cgcgGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTg                                                       >  2:4045308/1‑93 (MQ=255)

                                   gcgGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGa                                                      <  2:378179/93‑1 (MQ=255)

                                   gcgGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGa                                                      >  1:485093/1‑93 (MQ=255)

                                   gcgGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGa                                                      >  1:3475743/1‑93 (MQ=255)

                                    cgGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGAt                                                     <  1:3665107/93‑1 (MQ=255)

                                    cgGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGAt                                                     >  1:3792633/1‑93 (MQ=255)

                                    cgGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGAt                                                     <  2:3676459/93‑1 (MQ=255)

                                      gTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGc                                                   <  1:1350390/93‑1 (MQ=255)

                                       tATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGcc                                                  <  1:3475769/93‑1 (MQ=255)

                                        aTCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCa                                                 >  1:3203934/1‑93 (MQ=255)

                                         tCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCAc                                                >  1:791394/1‑93 (MQ=255)

                                          ccGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACt                                               >  2:326661/1‑93 (MQ=255)

                                            gCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCg                                             <  2:750178/93‑1 (MQ=255)

                                            gCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCg                                             >  1:4675323/1‑93 (MQ=255)

                                             ccccTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCgg                                            >  1:1890748/1‑93 (MQ=255)

                                              cccTGGGTTTTTCTCATAAAAGGCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGCGATTGGCAGATAACCGGTTTTCTGATGCAACTCGGc                                           >  2:4229973/1‑93 (MQ=255)

                                               ccTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCa                                          <  2:1713932/93‑1 (MQ=255)

                                                cTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAg                                         <  1:4266418/93‑1 (MQ=255)

                                                cTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAg                                         >  1:118753/1‑93 (MQ=255)

                                                 tGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGc                                        <  2:2152282/93‑1 (MQ=255)

                                                  gggTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCg                                       >  2:1132547/1‑93 (MQ=255)

                                                  gggTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCg                                       <  1:1535674/93‑1 (MQ=255)

                                                   ggTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGt                                      <  1:662349/93‑1 (MQ=255)

                                                   ggTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGt                                      >  2:3116322/1‑93 (MQ=255)

                                                    gTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGtt                                     >  1:2380284/1‑93 (MQ=255)

                                                     tttttCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGttt                                    >  2:2637286/1‑93 (MQ=255)

                                                     tttttCTCATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGttt                                    <  1:2649608/93‑1 (MQ=255)

                                                         tctcATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTg                                <  2:2723436/93‑1 (MQ=255)

                                                         tctcATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGCTTTCTGATGCCACTCGGCAGCGTTTTCTg                                >  1:3938468/1‑93 (MQ=255)

                                                           tcATAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGc                              >  1:2196142/1‑93 (MQ=255)

                                                              tAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTc                           >  1:3090393/1‑93 (MQ=255)

                                                              tAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTc                           <  1:3901857/93‑1 (MQ=255)

                                                              tAAAAGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTc                           >  1:1953468/1‑93 (MQ=255)

                                                                aaaGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTcgc                         <  1:2819268/93‑1 (MQ=255)

                                                                aaaGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTcgc                         <  2:1072145/93‑1 (MQ=255)

                                                                 aaGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTcgcg                        >  2:3547257/1‑93 (MQ=255)

                                                                  aGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTCGCGa                       >  2:1483872/1‑93 (MQ=255)

                                                                  aGCCCTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTCGCGa                       <  1:4549992/93‑1 (MQ=255)

                                                                    cccTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTCGCGAgg                     <  2:818580/93‑1 (MQ=255)

                                                                     ccTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTCGCGAGGa                    <  2:1320261/93‑1 (MQ=255)

                                                                      cTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGTTTAGCGAGGaa                   >  1:1476623/1‑93 (MQ=255)

                                                                      cTGCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTCGCGAGGaa                   <  2:562899/93‑1 (MQ=255)

                                                                       tGCTCACGGGTCAGGTCATGCGCTGCTTTGTTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTCGCGAGGaaa                  >  1:1013672/1‑93 (MQ=255)

                                                                        gCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTCGCGAGGAAAt                 >  1:3532305/1‑93 (MQ=255)

                                                                        gCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTCGCGAGGAAAt                 <  2:3269146/93‑1 (MQ=255)

                                                                        gCTCACGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTCGCGAGGAAAt                 >  2:3187633/1‑93 (MQ=255)

                                                                         cTCACGGGTCAGGTCATGCACTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTCGCGAGGAAATc                <  2:4214527/93‑1 (MQ=255)

                                                                            aCGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTCGCGAGGAAATCGAg             <  1:1260124/93‑1 (MQ=255)

                                                                             cGGGTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTCGCGAGGAAATCGAgg            <  1:2414290/93‑1 (MQ=255)

                                                                              gggTCAGGTCATGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCAGATGCCACTCGGCAGCGTTTTCTGGCTTCGCGAGGAAATCGAGGa           >  1:2064143/1‑93 (MQ=255)

                                                                               ggtcaggtcaTGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTCGCGAGGAAATCGAGGaa          >  2:4172292/1‑93 (MQ=255)

                                                                               ggtcaggtcaTGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTCGCGAGGAAATCGAGGaa          >  1:3399015/1‑93 (MQ=255)

                                                                                gtcaggtcaTGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTCGCGAGGAAATCGAGGAAc         >  2:3956005/1‑93 (MQ=255)

                                                                                gtcaggtcaTGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTCGCGAGGAAATCGAGGAAc         <  1:3436946/93‑1 (MQ=255)

                                                                                   aggtcaTGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTCGCGAGGAAATCGAGGAACTTc      >  1:2966777/1‑93 (MQ=255)

                                                                                    ggtcaTGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTCGCGAGCAAATCGAGGAACTTCg     <  2:4506409/93‑1 (MQ=255)

                                                                                     gtcaTGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTCGCGAGGAAATCGAGGAACTTCGc    >  1:2492869/1‑93 (MQ=255)

                                                                                     gtcaTGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTCGCGAGGAAATCGAGGAACTTCGc    >  2:4056673/1‑93 (MQ=255)

                                                                                     gtcaTGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTCGCGAGGAAATCGAGGAACTTCGc    >  1:2911096/1‑93 (MQ=255)

                                                                                      tcaTGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTCGCGAGGAAATCGAGGAACTTCGcc   >  1:1118026/1‑93 (MQ=255)

                                                                                      gcaTGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTCGCGAGGAAATCGAGGAACTTCGcc   <  1:1441448/92‑1 (MQ=255)

                                                                                       caTGCGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTCGCGAGGAAATCGAGGAACTTCGCca  >  2:1498058/1‑93 (MQ=255)

ACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCTGGCTTCGCGAGGAAATCGAGGAACTTCGCCA  >  NC_000913/3588807‑3588986

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