BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,366,054	A→G	H198R (CAC→CGC)	yhcG →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,366,054	0	A	G	98.8%	543.1 / ‑2.8	172	H198R (CAC→CGC)	yhcG	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (2/0); new base G (88/82); total (90/82)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.41e-01

Only 100 of 184 total aligned reads displayed.

GTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACACGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTGGAATTGAAAGATGAATATTCAG  >  NC_000913/3365965‑3366143

  cTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTc                                                                                      <  1:181679/93‑1 (MQ=255)

   tGTTCGCCAGCTAGAGCGGCAAATTGCGACCCAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTcc                                                                                     >  2:769761/1‑93 (MQ=255)

   tGTTCGCCAGCTAGAGCGGCAAATTGCGACCCAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTcc                                                                                     >  2:1877796/1‑93 (MQ=255)

   tGTTCGCCAGCTAGAGCGGCAAATTGCGACCCAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTcc                                                                                     <  1:2165313/93‑1 (MQ=255)

    gTTCGCCAGCTAGAGCGGCAAATTGCGACCCAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCt                                                                                    <  2:3365789/93‑1 (MQ=255)

     ttCGCCAGCTAGAGCGGCAAATTGCGACCCAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTg                                                                                   <  1:4038891/93‑1 (MQ=255)

      tCGCCAGCTAGAGCGGCAAATTGCGACCCAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGc                                                                                  <  1:3791823/93‑1 (MQ=255)

      tCGCCAGCTAGAGCGGCAAATTGCGACCCAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGc                                                                                  <  2:26083/93‑1 (MQ=255)

        gCCAGCTAGAGCGGCAAATTGCGACCCAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCg                                                                                >  1:91758/1‑93 (MQ=255)

        gCCAGCTAGAGCGGCAAATTGCGACCCAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCg                                                                                >  2:3504074/1‑93 (MQ=255)

         ccAGCTAGAGCGGCAAATTGCGACCCAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCga                                                                               >  1:3849585/1‑93 (MQ=255)

           aGCTAGAGCGGCAAATTGCGACCCAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCgaga                                                                             >  1:2537189/1‑93 (MQ=255)

            gCTAGAGCGGCAAATTGCGACCCAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGAc                                                                            <  2:3843589/93‑1 (MQ=255)

              tAGAGCGGCAAATTGCGACCCAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGc                                                                          >  2:3103489/1‑93 (MQ=255)

               agagCGGCAAATTGCGACCCAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCa                                                                         <  2:3752054/93‑1 (MQ=255)

                    ggCAAATTGCGACCCAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATAttc                                                                    >  1:1602274/1‑93 (MQ=255)

                      cAAATTGCGACCCAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATAttctt                                                                  <  2:1313157/93‑1 (MQ=255)

                      cAAATTGCGACCCAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATAttctt                                                                  <  1:3520892/93‑1 (MQ=255)

                       aaaTTGCGACCCAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATAttcttc                                                                 >  1:1295932/1‑93 (MQ=255)

                         aTTGCGACCCAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAATGCGCTCCTGCCGAGACGCATATTCTTCCg                                                               >  1:1990962/1‑93 (MQ=255)

                         aTTGCGACCCAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCg                                                               >  1:4397534/1‑93 (MQ=255)

                           tGCGACCCAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCa                                                             <  2:4324584/93‑1 (MQ=255)

                              gACCCAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACa                                                          >  1:1735598/1‑93 (MQ=255)

                              gACCCAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACa                                                          <  1:1020509/93‑1 (MQ=255)

                               aCCCAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAg                                                         >  2:3568155/1‑93 (MQ=255)

                               aCCCAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAg                                                         >  2:2208305/1‑93 (MQ=255)

                               aCCCAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAg                                                         <  2:874552/93‑1 (MQ=255)

                                 ccAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGc                                                       >  2:1111123/1‑93 (MQ=255)

                                  cAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCg                                                      >  2:3042326/1‑93 (MQ=255)

                                  cAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCg                                                      <  1:554660/93‑1 (MQ=255)

                                  cAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACACGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCg                                                      >  1:1575318/1‑93 (MQ=255)

                                   aGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGa                                                     >  2:2650174/1‑93 (MQ=255)

                                    gTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGAt                                                    >  1:1028080/1‑93 (MQ=255)

                                    gTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGAt                                                    <  2:1826062/93‑1 (MQ=255)

                                    gTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGAt                                                    >  1:2939318/1‑93 (MQ=255)

                                    gTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGAt                                                    <  1:2448741/93‑1 (MQ=255)

                                    gTTTTATAAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGAt                                                    <  1:2862331/93‑1 (MQ=255)

                                     ttttATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATa                                                   <  1:3312700/93‑1 (MQ=255)

                                     ttttATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATa                                                   <  2:3030656/93‑1 (MQ=255)

                                      tttATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATAc                                                  <  1:4303427/93‑1 (MQ=255)

                                       ttATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATAcg                                                 <  2:1720261/93‑1 (MQ=255)

                                        tATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATAcgc                                                <  2:2412540/93‑1 (MQ=255)

                                        tATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATAcgc                                                >  2:504128/1‑93 (MQ=255)

                                           gAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGAt                                             >  2:2852036/1‑93 (MQ=255)

                                             gCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATcc                                           <  2:2626881/93‑1 (MQ=255)

                                             gCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCACATCAGGCGATACGCGATcc                                           >  2:4209616/1‑93 (MQ=255)

                                              cGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATccc                                          <  2:630173/93‑1 (MQ=255)

                                              cGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATccc                                          <  1:4548927/93‑1 (MQ=255)

                                                gACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCtt                                        <  1:3232312/93‑1 (MQ=255)

                                                gACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCtt                                        >  2:3416590/1‑93 (MQ=255)

                                                gACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCtt                                        <  1:1176929/93‑1 (MQ=255)

                                                gACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCtt                                        >  1:3583665/1‑93 (MQ=255)

                                                gACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCtt                                        <  2:2057695/93‑1 (MQ=255)

                                                 acacTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCttt                                       >  2:714999/1‑93 (MQ=255)

                                                  cacTACTGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTtg                                      <  1:3853677/93‑1 (MQ=255)

                                                   actactGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTtgt                                     <  2:4059816/93‑1 (MQ=255)

                                                     tactGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGc                                   >  2:237926/1‑93 (MQ=255)

                                                      actGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCt                                  <  1:2086272/93‑1 (MQ=255)

                                                      actGTCACATGACAAATCAGCCATGCTGCAACAACACGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCt                                  >  2:1587218/1‑93 (MQ=255)

                                                       ctGTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTc                                 >  1:1444985/1‑93 (MQ=255)

                                                         gTCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGa                               >  2:4590401/1‑93 (MQ=255)

                                                          tCACATGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGaa                              >  1:4513056/1‑93 (MQ=255)

                                                          tCACATGAAAACTCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGaa                              <  2:2541027/93‑1 (MQ=255)

                                                            cgaTGACAAACCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAAtt                            <  2:2352815/91‑1 (MQ=255)

                                                            acaTGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAAtt                            <  2:1974220/93‑1 (MQ=255)

                                                            acaTGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAAtt                            >  2:1259818/1‑93 (MQ=255)

                                                             caTGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAAttt                           <  2:4536816/93‑1 (MQ=255)

                                                              aTGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTc                          >  2:815688/1‑93 (MQ=255)

                                                               tGACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCt                         >  1:992475/1‑93 (MQ=255)

                                                               agACAAATCAGCCATCGCGCAACAACCCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCt                         <  1:387800/92‑1 (MQ=255)

                                                                gACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTg                        <  1:2987378/93‑1 (MQ=255)

                                                                gACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTg                        >  2:314854/1‑93 (MQ=255)

                                                                gACAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTg                        >  1:4737945/1‑93 (MQ=255)

                                                                 aCAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTgg                       <  2:2502112/93‑1 (MQ=255)

                                                                 aCAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTgg                       >  2:2366229/1‑93 (MQ=255)

                                                                  cAAATCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTGGa                      >  2:3500585/1‑93 (MQ=255)

                                                                   aaaTCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTGGaa                     >  2:3175094/1‑93 (MQ=255)

                                                                    aaTCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTGGAAt                    >  1:42301/1‑93 (MQ=255)

                                                                     aTCAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTGGAAtt                   <  2:3787679/93‑1 (MQ=255)

                                                                       cAGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTGGAATTGa                 <  1:3128268/93‑1 (MQ=255)

                                                                        aGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTGGAATTGaa                >  2:1284411/1‑93 (MQ=255)

                                                                        aGCCATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTGGAATTGaa                >  2:383752/1‑93 (MQ=255)

                                                                          ccATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTGGAATTGAAAg              >  1:4140631/1‑93 (MQ=255)

                                                                           cATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTGGAATTGAAAGa             >  1:4714665/1‑93 (MQ=255)

                                                                           cATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTGGAATTGAAAGa             <  1:1829993/93‑1 (MQ=255)

                                                                           cATGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTGGAATTGAAAGa             >  1:865432/1‑93 (MQ=255)

                                                                            ctgctgcAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTGGAATTGAAAGAt            <  2:1803738/92‑1 (MQ=255)

                                                                            aTGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTGGAATTGAAAGAt            >  1:768101/1‑93 (MQ=255)

                                                                            aTGCTGCAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTGGAAGTGAAAGAt            >  2:3649052/1‑93 (MQ=255)

                                                                              gctgcAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTGGAATTGAAAGATGa          >  1:1631103/1‑93 (MQ=255)

                                                                              gctgcAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTGGAATTGAAAGATGa          <  2:3192329/93‑1 (MQ=255)

                                                                              gctgcAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTGGAATTGAAAGATGa          <  1:950749/93‑1 (MQ=255)

                                                                                tgcAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTGGAATTGAAAGATGAat        <  2:1514146/93‑1 (MQ=255)

                                                                                tgcAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTGGAATTAAAAGATGAat        >  1:1074742/1‑93 (MQ=255)

                                                                                 gcAACAACGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTGGAATTGAAAGATGAata       >  1:4203265/1‑93 (MQ=255)

                                                                                   aacaacGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTGGAATTGAAAGATGAATAtt     >  2:3253540/1‑93 (MQ=255)

                                                                                   aacaacGCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTGGAATTGAAAGATGAATAtt     <  2:1346835/93‑1 (MQ=255)

                                                                                    acaacgCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTGGAATTGAAAGATGAATATTc    <  2:1181460/93‑1 (MQ=255)

                                                                                     caacgCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTGGAATTGAAAGATGAATATTCa   <  2:1124761/93‑1 (MQ=255)

                                                                                      aacgCGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTGGAATTGAAAGATGAATATTCAg  <  2:3471685/93‑1 (MQ=255)

GTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCAGTTTTATGAGCGGACACTACTGTCACATGACAAATCAGCCATGCTGCAACAACACGCTCCTGCCGAGACGCATATTCTTCCGCAACAGGCGATACGCGATCCCTTTGTGCTCGAATTTCTGGAATTGAAAGATGAATATTCAG  >  NC_000913/3365965‑3366143

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: `‑`