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.1%	350.4 / ‑4.6	113	Y52H (TAC→CAC)	umuD	UmuD
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (51/61); total (52/61)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 122 total aligned reads displayed.

GTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAGCGCTATTACCGCCAGC  >  NC_000913/1230053‑1230232

gTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGactcac                                                                                         <  1:2229730/93‑1 (MQ=255)

gTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGactcac                                                                                         <  2:3545726/93‑1 (MQ=255)

gTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGactcac                                                                                         <  1:2852665/93‑1 (MQ=255)

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

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

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

   cAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTc                                                                                      <  2:1682869/93‑1 (MQ=255)

   cAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTc                                                                                      <  1:703386/93‑1 (MQ=255)

     gtgtGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGt                                                                                    <  2:2242662/93‑1 (MQ=255)

     gtgtGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGt                                                                                    >  1:4034517/1‑93 (MQ=255)

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

      tgtgGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTc                                                                                   >  2:3048223/1‑93 (MQ=255)

       gtgGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCa                                                                                  <  1:3739831/93‑1 (MQ=255)

        tgGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCaa                                                                                 >  2:1057771/1‑93 (MQ=255)

        tgGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCaa                                                                                 >  2:1033006/1‑93 (MQ=255)

         ggCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCaaa                                                                                <  2:1197840/93‑1 (MQ=255)

          gCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAg                                                                               >  2:130769/1‑93 (MQ=255)

          gCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAg                                                                               >  2:1394478/1‑93 (MQ=255)

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

          gCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGAGCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAg                                                                               <  1:631178/93‑1 (MQ=255)

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

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

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

                 ttCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAgtgg                                                                        >  1:4610556/1‑93 (MQ=255)

                  tCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAgtggt                                                                       >  2:4395696/1‑93 (MQ=255)

                  tCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAgtggt                                                                       <  2:1545810/93‑1 (MQ=255)

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

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

                      cGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGAtt                                                                   >  2:911757/1‑93 (MQ=255)

                       ggCAGCAGATTACGTTGAACAGCTCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTc                                                                  >  1:1865013/1‑93 (MQ=255)

                        gcagcaGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCt                                                                 >  1:4078386/1‑93 (MQ=255)

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

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

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

                              gATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTAtgatt                                                           >  2:847758/1‑93 (MQ=255)

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

                               aTTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCGAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTAtgattg                                                          >  2:4548093/1‑93 (MQ=255)

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

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

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

                                        aaCAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGaa                                                 <  2:4554913/93‑1 (MQ=255)

                                        aaCAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGaa                                                 <  2:1865598/93‑1 (MQ=255)

                                        aaCAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGaa                                                 <  1:4118217/93‑1 (MQ=255)

                                        aaCAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGaa                                                 <  2:2844727/93‑1 (MQ=255)

                                         aCAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAAt                                                <  2:1456445/93‑1 (MQ=255)

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

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

                                              gcATCGATCTGAATCAACTGTTGATCCAGCATCCCACCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTg                                           >  1:1904556/1‑93 (MQ=255)

                                               cATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCCCTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGa                                          <  2:4407388/93‑1 (MQ=255)

                                               cATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGa                                          <  2:4110458/93‑1 (MQ=255)

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

                                                aTCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGAc                                         >  1:827299/1‑93 (MQ=255)

                                                aTCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGAc                                         >  2:1795786/1‑93 (MQ=255)

                                                aTCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGAc                                         >  2:307659/1‑93 (MQ=255)

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

                                                  cGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACgg                                       <  2:3631160/93‑1 (MQ=255)

                                                    aTCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTg                                     <  1:3153009/93‑1 (MQ=255)

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

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

                                                      cTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGAt                                   <  2:3126628/93‑1 (MQ=255)

                                                      cTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGAt                                   <  1:3065204/93‑1 (MQ=255)

                                                      cTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGAt                                   <  2:2784769/93‑1 (MQ=255)

                                                         aaTCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTa                                >  2:2048198/1‑93 (MQ=255)

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

                                                          aTCAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTAc                               <  1:2949252/93‑1 (MQ=255)

                                                            cAACTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTg                             >  2:4615267/1‑93 (MQ=255)

                                                              aCTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGAt                           >  2:389099/1‑93 (MQ=255)

                                                              aCTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGAt                           >  2:389100/1‑93 (MQ=255)

                                                               cTGTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGAtt                          >  2:3185726/1‑93 (MQ=255)

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

                                                                 gTTGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGt                        >  2:2503447/1‑93 (MQ=255)

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

                                                                  gtGATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTc                       <  1:2251595/92‑1 (MQ=255)

                                                                    gATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGa                     >  1:2427740/1‑93 (MQ=255)

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

                                                                    gATCCAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGa                     >  2:1560433/1‑93 (MQ=255)

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

                                                                     aTCCAGCATCCCAGCGCGGCTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGAt                    >  1:1593002/1‑93 (MQ=255)

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

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

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

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

                                                                       ccAGCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAg                  <  1:3727199/93‑1 (MQ=255)

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

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

                                                                          gCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAgcgc               >  2:696518/1‑93 (MQ=255)

                                                                          gCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAgcgc               <  1:1295863/93‑1 (MQ=255)

                                                                          gCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAgcgc               <  2:1051913/93‑1 (MQ=255)

                                                                          gCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAgcgc               <  2:944782/93‑1 (MQ=255)

                                                                          gCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAgcgc               >  2:3847397/1‑93 (MQ=255)

                                                                          gCATCCCAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAgcgc               <  1:466672/93‑1 (MQ=255)

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

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

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

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

                                                                               ccAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAGCGCTAttc          >  1:966687/1‑92 (MQ=255)

                                                                               ccAGCGCGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAGCGCTATTa          <  2:3398120/93‑1 (MQ=255)

                                                                                  gcgcGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAGCGCTATTACCg       >  1:3555050/1‑93 (MQ=255)

                                                                                  gcgcGACTCACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAGCGCTATTACCg       >  1:3961834/1‑93 (MQ=255)

                                                                                       actcacttCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAGCGCTATTACCgccagc  >  1:1215529/1‑93 (MQ=255)

GTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATTGTCGATAGCGCTATTACCGCCAGC  >  NC_000913/1230053‑1230232

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