BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,230,143	T→C	Y52H (TAC→CAC)	umuD →	UmuD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,230,143	0	T	C	99.7%	1037.3 / ‑5.2	320	Y52H (TAC→CAC)	umuD	UmuD
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (145/174); total (146/174)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.90e-01

Only 100 of 341 total aligned reads displayed.

GTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAGCGCTATTACCGCCAGC  >  NC_000913/1230053‑1230232

gTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGactcac                                                                                         >  2:6289648/1‑93 (MQ=255)

  tCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGactcactt                                                                                       <  1:7624999/93‑1 (MQ=255)

      tgtgGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTc                                                                                   <  2:1735874/93‑1 (MQ=255)

          gCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAg                                                                               <  2:8245209/93‑1 (MQ=255)

           cTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGc                                                                              <  2:4507466/93‑1 (MQ=255)

           cTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGc                                                                              >  2:265324/1‑93 (MQ=255)

            tttCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCa                                                                             <  2:35031/93‑1 (MQ=255)

             ttCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCaa                                                                            <  2:3206732/93‑1 (MQ=255)

             ttCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCaa                                                                            <  1:7523008/93‑1 (MQ=255)

             ttCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCaa                                                                            <  1:6169895/93‑1 (MQ=255)

             ttCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCaa                                                                            >  1:474747/1‑93 (MQ=255)

             ttCCGTCACCGGCAGCATAGTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCTCGACTCACTTCGTCAAATCaa                                                                            <  1:1136439/93‑1 (MQ=255)

              tCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAg                                                                           <  1:2675549/93‑1 (MQ=255)

              tCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAg                                                                           >  2:2109566/1‑93 (MQ=255)

              tCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAg                                                                           <  2:8643248/93‑1 (MQ=255)

               ccTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAgt                                                                          <  1:7530357/93‑1 (MQ=255)

                cTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAgtg                                                                         <  2:1924686/93‑1 (MQ=255)

                  tCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAgtggg                                                                       >  1:4175035/1‑92 (MQ=255)

                    aCCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGa                                                                     <  1:6370082/93‑1 (MQ=255)

                     ccGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGAt                                                                    <  2:4732673/93‑1 (MQ=255)

                     ccGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGAt                                                                    >  1:9114825/1‑93 (MQ=255)

                      cGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGAtt                                                                   <  2:4261850/93‑1 (MQ=255)

                         cagcagATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTa                                                                <  1:4243734/93‑1 (MQ=255)

                          agcagATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTAt                                                               <  1:4440656/93‑1 (MQ=255)

                           gcagATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTAtg                                                              >  1:3313567/1‑93 (MQ=255)

                           gcagATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTAtg                                                              >  2:153152/1‑93 (MQ=255)

                            cagATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTAtga                                                             <  1:9906757/93‑1 (MQ=255)

                              gATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTAtgatt                                                           >  1:2518022/1‑93 (MQ=255)

                              gATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTAtgatt                                                           <  2:120557/93‑1 (MQ=255)

                               aTTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTAtgattg                                                          <  1:9056955/93‑1 (MQ=255)

                               aTTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTAtgattg                                                          >  1:5582092/1‑93 (MQ=255)

                               aTTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTAtgattg                                                          <  1:9525196/93‑1 (MQ=255)

                                ttACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTAtgattga                                                         <  2:4895912/93‑1 (MQ=255)

                                ttACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTAtgattga                                                         <  1:3685768/93‑1 (MQ=255)

                                ttACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTAtgattga                                                         >  1:3188245/1‑93 (MQ=255)

                                ttACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTAtgattga                                                         >  2:7742248/1‑93 (MQ=255)

                                 tACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTAtgattgat                                                        >  1:9341859/1‑93 (MQ=255)

                                  aCGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGAtg                                                       <  1:1027751/93‑1 (MQ=255)

                                   cGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGAtgg                                                      <  1:185676/93‑1 (MQ=255)

                                    gTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGAtggt                                                     <  2:9609527/93‑1 (MQ=255)

                                    gTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGAtggg                                                     >  1:3612437/1‑92 (MQ=255)

                                     ttcaacaacGCAACTATCTGAATCAACTGTTGATCCACCATCCCAGCGCCACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGCtggtg                                                    <  1:7161176/90‑1 (MQ=255)

                                     ttGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGAtggtg                                                    >  1:140090/1‑93 (MQ=255)

                                     ttGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGAtggtg                                                    >  2:4053478/1‑93 (MQ=255)

                                      tGAACAGCGCATCGATCTGAACCAATGGTTGATCCAGCATTCCTTCGCGACGCTCTGCTGCAAAGCAAGTGGTGATTCTATGAGTGAtggggg                                                   >  2:8973964/1‑93 (MQ=255)

                                        aaCAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGat                                                 >  2:8298827/1‑92 (MQ=255)

                                        aaCAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGaa                                                 >  2:5751082/1‑93 (MQ=255)

                                           aGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTa                                              >  2:2856227/1‑93 (MQ=255)

                                            gcgcATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAg                                             >  2:5442484/1‑93 (MQ=255)

                                            gcgcATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAg                                             <  2:3011660/93‑1 (MQ=255)

                                            gcgcATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAg                                             >  2:6578723/1‑93 (MQ=255)

                                              gcATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTg                                           <  2:1764386/93‑1 (MQ=255)

                                              gcATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTg                                           >  1:4799253/1‑93 (MQ=255)

                                                aTCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGAc                                         <  1:6947903/93‑1 (MQ=255)

                                                 tCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACg                                        >  1:9199270/1‑93 (MQ=255)

                                                 tCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACg                                        >  2:5031660/1‑93 (MQ=255)

                                                    aTCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTg                                     >  2:9355588/1‑93 (MQ=255)

                                                     tCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGa                                    >  2:4568067/1‑93 (MQ=255)

                                                     tCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGa                                    >  1:8205321/1‑93 (MQ=255)

                                                      cTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGAt                                   >  2:384215/1‑93 (MQ=255)

                                                      cTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTACTGGTGGAATTAGTGACGGTGAt                                   >  2:6938983/1‑93 (MQ=255)

                                                       tGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGAtt                                  <  2:7037819/93‑1 (MQ=255)

                                                         aaTCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTa                                <  2:5313265/93‑1 (MQ=255)

                                                         aaTCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTa                                <  2:2323824/93‑1 (MQ=255)

                                                           tCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACt                              <  1:8071766/93‑1 (MQ=255)

                                                           tCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACt                              >  1:6677380/1‑93 (MQ=255)

                                                              aCTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGAt                           <  1:654702/93‑1 (MQ=255)

                                                               cTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGCATTAGTGACGGTGATTTACTGAtt                          >  1:5281367/1‑93 (MQ=255)

                                                               cTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGAtt                          <  2:2545511/93‑1 (MQ=255)

                                                                tGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTg                         >  2:9869449/1‑93 (MQ=255)

                                                                tGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGGGATTCTATGATTAATGGTGGAATTAGTGACGGTGATTTACTGATTg                         >  1:9699450/1‑93 (MQ=255)

                                                                 gTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGt                        <  1:6034966/93‑1 (MQ=255)

                                                                  ttGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTc                       <  2:3763120/93‑1 (MQ=255)

                                                                    taTCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGa                     <  2:7612482/92‑1 (MQ=255)

                                                                    gATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGa                     <  1:5530071/93‑1 (MQ=255)

                                                                     aTCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGAt                    <  2:5854123/93‑1 (MQ=255)

                                                                     aTCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGAt                    <  2:6645515/93‑1 (MQ=255)

                                                                      tCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATa                   <  2:30814/93‑1 (MQ=255)

                                                                      tCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATa                   >  1:7243511/1‑93 (MQ=255)

                                                                       ccAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGCGATTTACTGATTGTCGCTAg                  >  2:857899/1‑93 (MQ=255)

                                                                        cAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAgc                 <  2:8682729/93‑1 (MQ=255)

                                                                         aGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAgcg                >  1:3780717/1‑93 (MQ=255)

                                                                         aGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGAGATTTACTGATTGTCGATAgcg                >  1:5900372/1‑93 (MQ=255)

                                                                            aTCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAGCGCTa             >  2:6336580/1‑93 (MQ=255)

                                                                            aTCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAGCGCTa             >  2:8156581/1‑93 (MQ=255)

                                                                             tCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAGCGCTAt            >  1:8585178/1‑93 (MQ=255)

                                                                             tCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAGCGCTAt            >  2:9085004/1‑93 (MQ=255)

                                                                             tCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAGCGCTAt            <  2:7332732/93‑1 (MQ=255)

                                                                             tCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAGCGCTAt            <  1:700635/93‑1 (MQ=255)

                                                                              cccAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAGCGCTAtt           >  2:6060037/1‑93 (MQ=255)

                                                                               ccAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAGCGCTATTa          >  2:8383184/1‑93 (MQ=255)

                                                                               ccAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAGCGCTATTa          >  2:3667222/1‑93 (MQ=255)

                                                                               ccAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAGCGCTATTa          >  1:2297779/1‑93 (MQ=255)

                                                                                cAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAGCGCTATTAc         <  2:2192873/93‑1 (MQ=255)

                                                                                cAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAGCGCTATTAc         <  1:4946314/93‑1 (MQ=255)

                                                                                cAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAGCGCTATTAc         <  1:1732332/93‑1 (MQ=255)

                                                                                  gcgcGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAGCGCTATTACCg       <  1:8461053/93‑1 (MQ=255)

                                                                                    gcgACGCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAGCGCTATTACCgcc     >  2:5820514/1‑93 (MQ=255)

                                                                                      gACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAGCGCTATTACCgccag   <  1:7818694/93‑1 (MQ=255)

                                                                                       actcacttCGTCAAAGCAAGTGGTGATTCTATGATTTATGGTGGAATTAGTGACGGTGATTTACTGATTGTTGATAGCGCTATTACCgccagc  >  1:8963558/1‑93 (MQ=255)

GTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAGCGCTATTACCGCCAGC  >  NC_000913/1230053‑1230232

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