BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,614,182	T→C	F124L (TTC→CTC)	nikC →	NikC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,614,182	0	T	C	99.3%	985.4 / ‑5.0	305	F124L (TTC→CTC)	nikC	NikC
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (141/162); total (142/163)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 314 total aligned reads displayed.

GGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGTTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTATGCACGCATGGTG  >  NC_000913/3614092‑3614271

ggcAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTc                                                                                         <  2:8219152/93‑1 (MQ=255)

ggcAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTc                                                                                         >  2:5642674/1‑93 (MQ=255)

 gcAGCGCCGGGTTGTTTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGctct                                                                                        <  2:7745052/93‑1 (MQ=255)

   aGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCttt                                                                                      <  2:8384185/93‑1 (MQ=255)

   aGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCttt                                                                                      <  2:1826812/93‑1 (MQ=255)

    gcgcCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGCTTCTGTCGCTCTTTa                                                                                     >  1:1102299/1‑93 (MQ=255)

    gcgcCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTa                                                                                     >  2:9660007/1‑93 (MQ=255)

    gcgcCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTa                                                                                     <  1:2574906/93‑1 (MQ=255)

    gcgcCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTa                                                                                     <  1:5861378/93‑1 (MQ=255)

    gcgcCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTa                                                                                     <  2:4088409/93‑1 (MQ=255)

       ccGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATgg                                                                                  >  2:1496130/1‑93 (MQ=255)

        cGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGt                                                                                 <  1:2716144/93‑1 (MQ=255)

        cGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGt                                                                                 >  2:8075294/1‑93 (MQ=255)

         gggTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGtt                                                                                >  1:4755271/1‑93 (MQ=255)

          ggTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTAttgttg                                                                               >  1:8430532/1‑93 (MQ=255)

              gATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATTGGTGGCgg                                                                           >  2:9166404/1‑92 (MQ=255)

               aTTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGTTTGGCggg                                                                          >  1:2363179/1‑91 (MQ=255)

                ttGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGc                                                                         <  2:2146062/93‑1 (MQ=255)

                 tGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCt                                                                        <  2:4777204/93‑1 (MQ=255)

                 tGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCt                                                                        >  2:7562322/1‑93 (MQ=255)

                 tGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCt                                                                        <  2:3982615/93‑1 (MQ=255)

                  ggcgggcgCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTc                                                                       <  2:3460140/93‑1 (MQ=255)

                  ggcgggcgCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTc                                                                       >  2:1384547/1‑93 (MQ=255)

                   gcgggcgCGTTTATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCg                                                                      <  1:8236379/93‑1 (MQ=255)

                    cgggcgCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCgg                                                                     >  2:890092/1‑93 (MQ=255)

                      ggcgCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCa                                                                   >  1:4383791/1‑93 (MQ=255)

                       gcgcGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCAc                                                                  <  1:71746/93‑1 (MQ=255)

                         gcgTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCg                                                                <  1:7780949/93‑1 (MQ=255)

                           gTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCggg                                                              <  2:6899832/93‑1 (MQ=255)

                           gTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCggg                                                              >  2:5915976/1‑93 (MQ=255)

                           gTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCggg                                                              <  1:4035212/93‑1 (MQ=255)

                           gTTGATCAGGCCACCATGCGCGTCGCAGATATGTTTATTACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCggg                                                              <  2:2626966/93‑1 (MQ=255)

                            ttGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGc                                                             <  1:5673602/93‑1 (MQ=255)

                              gATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTg                                                           >  1:7308781/1‑93 (MQ=255)

                                   ggCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCaa                                                      >  2:608302/1‑93 (MQ=255)

                                   ggCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCaa                                                      >  1:647962/1‑93 (MQ=255)

                                     ccaccaTGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACg                                                    >  1:7473784/1‑93 (MQ=255)

                                       accaTGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTa                                                  >  2:72861/1‑93 (MQ=255)

                                        ccaTGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTaa                                                 <  2:9397536/93‑1 (MQ=255)

                                         caTGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAAt                                                >  2:5975782/1‑93 (MQ=255)

                                          aTGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAAtt                                               <  2:3886432/93‑1 (MQ=255)

                                           tGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTa                                              >  1:7040872/1‑93 (MQ=255)

                                            gcgcGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTat                                             <  2:5796098/93‑1 (MQ=255)

                                             cgcgTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTatc                                            <  2:643252/93‑1 (MQ=255)

                                               cgTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGTGCTGACCAACGTAATTatcgc                                          >  1:6802994/1‑93 (MQ=255)

                                                gTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTatcgcc                                         >  2:8324296/1‑93 (MQ=255)

                                                gTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGCTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTatagca                                         >  1:1847396/1‑92 (MQ=255)

                                                 tCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTatcgcca                                        >  1:5315993/1‑93 (MQ=255)

                                                 tCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTatcgcca                                        <  2:7787926/93‑1 (MQ=255)

                                                 tCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTatcgcca                                        >  2:8994413/1‑93 (MQ=255)

                                                 tCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTatcgcca                                        >  2:3071681/1‑93 (MQ=255)

                                                  cGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCATGCTCGGCACCGGGCTGACCAACGTAATTatcgccat                                       <  2:5521846/93‑1 (MQ=255)

                                                   gCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTatcgccatc                                      >  2:7048891/1‑93 (MQ=255)

                                                   gCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTatcgccatc                                      <  2:3402455/93‑1 (MQ=255)

                                                     cGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTatcgccatcgc                                    >  2:8767661/1‑93 (MQ=255)

                                                      gATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTatcgccatcgcc                                   >  2:9256911/1‑93 (MQ=255)

                                                      catatGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTatcgccatcgcc                                   <  1:3612480/92‑1 (MQ=255)

                                                       atatGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGccc                                  >  1:6202766/1‑93 (MQ=255)

                                                       atatGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGccc                                  <  1:6006088/93‑1 (MQ=255)

                                                       atatGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGccc                                  <  1:470330/93‑1 (MQ=255)

                                                       atatGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGACATCGccc                                  >  2:4620145/1‑93 (MQ=255)

                                                        tatGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCt                                 >  1:1068804/1‑93 (MQ=255)

                                                        tatGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCt                                 >  2:9133070/1‑93 (MQ=255)

                                                          tGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGt                               >  1:1540335/1‑93 (MQ=255)

                                                           ttttATGACCCTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGTCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTc                              <  1:4953947/92‑1 (MQ=255)

                                                           gTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTc                              >  2:7149884/1‑93 (MQ=255)

                                                           gTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTc                              >  2:3603264/1‑93 (MQ=255)

                                                           gTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTc                              >  2:3262690/1‑93 (MQ=255)

                                                           gTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTc                              >  1:7913026/1‑93 (MQ=255)

                                                           gTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTc                              <  2:7667604/93‑1 (MQ=255)

                                                            tttATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCg                             >  1:3871785/1‑93 (MQ=255)

                                                              tATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCa                           <  2:2291260/93‑1 (MQ=255)

                                                               aTGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGTCCTGTCGCAc                          <  2:5868473/93‑1 (MQ=255)

                                                               aTGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCAc                          <  1:3453498/93‑1 (MQ=255)

                                                               aTGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCAc                          <  2:177005/93‑1 (MQ=255)

                                                                tGACCTTCCCGACCTCGATTCTGTCGTTCTTTAGGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACt                         <  1:2908749/93‑1 (MQ=255)

                                                                 gACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTg                        >  2:1841305/1‑93 (MQ=255)

                                                                   ccTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCCCTggg                      >  2:4348221/1‑93 (MQ=255)

                                                                    cTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGc                     <  1:1264654/93‑1 (MQ=255)

                                                                      tCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGt                   >  1:550516/1‑93 (MQ=255)

                                                                         cGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGt                >  2:2556482/1‑93 (MQ=255)

                                                                          gACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTa               <  1:7239388/93‑1 (MQ=255)

                                                                          gACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTa               >  1:8740758/1‑93 (MQ=255)

                                                                           gccTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTAt              <  2:5216777/92‑1 (MQ=255)

                                                                           aCCTCGATTCTGTCGTTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTAt              >  1:3203739/1‑93 (MQ=255)

                                                                           aCCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTAt              >  1:152332/1‑93 (MQ=255)

                                                                            ccTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTATg             >  2:4738299/1‑93 (MQ=255)

                                                                             cTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTATGc            >  2:5052255/1‑93 (MQ=255)

                                                                             cTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTATGc            <  2:5899722/93‑1 (MQ=255)

                                                                              tCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTATGCa           <  1:4859838/93‑1 (MQ=255)

                                                                              tCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTATGCa           <  1:4486673/93‑1 (MQ=255)

                                                                                  ttCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTATGCACGCa       <  1:5116260/93‑1 (MQ=255)

                                                                                    cTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTATGCACGCATg     >  2:1189837/1‑93 (MQ=255)

                                                                                    cTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTATGCACGCATg     <  1:1770883/93‑1 (MQ=255)

                                                                                    cTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTATGCACGCATg     <  2:2884493/93‑1 (MQ=255)

                                                                                     tGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTATGCACGCATgg    <  1:9256912/93‑1 (MQ=255)

                                                                                     tGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTATGCACGCATgg    <  1:6231332/93‑1 (MQ=255)

                                                                                      gTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTATGCACGCATGGt   <  1:203882/93‑1 (MQ=255)

                                                                                      gTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGACCTGTCGCACTGGGCGTGGTATGCACGCATGGt   >  1:9731319/1‑93 (MQ=255)

                                                                                       tCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTATGCACGCATGGTg  >  1:9452471/1‑93 (MQ=255)

GGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGTTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTATGCACGCATGGTG  >  NC_000913/3614092‑3614271

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