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.6%	873.4 / ‑4.8	267	F124L (TTC→CTC)	nikC	NikC
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (134/132); total (134/133)
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.66e-01

Only 100 of 278 total aligned reads displayed.

GGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGTTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTATGCACGCATGGTG  >  NC_000913/3614092‑3614271

 gcAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGctct                                                                                        >  1:6083455/1‑93 (MQ=255)

 gcAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGctct                                                                                        <  2:4761447/93‑1 (MQ=255)

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

     cgcCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTAt                                                                                    >  2:3695623/1‑93 (MQ=255)

      gcCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATg                                                                                   <  2:6174967/93‑1 (MQ=255)

         gggTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGtt                                                                                <  1:1287577/93‑1 (MQ=255)

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

            ttGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGc                                                                             <  2:1179743/93‑1 (MQ=255)

             tGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCg                                                                            <  1:2091685/93‑1 (MQ=255)

             tGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCg                                                                            <  1:2332261/93‑1 (MQ=255)

              gATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGt                                                                           >  2:4197599/1‑93 (MQ=255)

              gATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGt                                                                           >  1:318890/1‑93 (MQ=255)

               aTTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTg                                                                          <  2:4434212/93‑1 (MQ=255)

                ttGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGc                                                                         >  2:803994/1‑93 (MQ=255)

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

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

                   gcgggcgCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCg                                                                      <  1:3142255/93‑1 (MQ=255)

                   gcgggcgCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCg                                                                      >  2:1026502/1‑93 (MQ=255)

                   gcgggcgCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCg                                                                      <  1:6370231/93‑1 (MQ=255)

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

                        cgcgTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCAcc                                                                 <  2:1373374/93‑1 (MQ=255)

                        cgcgTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCAcc                                                                 <  2:4133771/93‑1 (MQ=255)

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

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

                            ttGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGGACCGGGc                                                             >  1:932972/1‑93 (MQ=255)

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

                            ttGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGc                                                             >  2:703648/1‑93 (MQ=255)

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

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

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

                               aTCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTgt                                                          <  1:4524325/93‑2 (MQ=255)

                                 cAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGAcc                                                        >  1:4230127/1‑93 (MQ=255)

                                 cAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGAcc                                                        >  1:1031012/1‑93 (MQ=255)

                                 cAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGAcc                                                        >  1:6187529/1‑93 (MQ=255)

                                 cAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGAcc                                                        >  1:5059478/1‑93 (MQ=255)

                                  aGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCa                                                       >  2:297749/1‑93 (MQ=255)

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

                                   ggCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCaa                                                      <  1:3556414/93‑1 (MQ=255)

                                    gCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAAc                                                     <  1:5996123/93‑1 (MQ=255)

                                      caccaTGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGt                                                   <  2:5215405/93‑1 (MQ=255)

                                       accaTGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGGTGACGTGCTCGGCACCGGGCTGACCAACGTa                                                  >  1:3714762/1‑93 (MQ=255)

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

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

                                          aTGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAAtt                                               >  2:713191/1‑93 (MQ=255)

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

                                              gcgTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTatcg                                           >  1:3024309/1‑93 (MQ=255)

                                              gcgTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTatcg                                           <  1:109680/93‑1 (MQ=255)

                                               cgTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTatcgc                                          >  1:1164397/1‑93 (MQ=255)

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

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

                                                  cGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTatcgccat                                       <  1:1689352/93‑1 (MQ=255)

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

                                                    ccGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTatcgccatcg                                     <  1:2640649/93‑1 (MQ=255)

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

                                                     cGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGCCGTGCTCGGCACCGGGCTAACCAACGTAATTatcgccatcgc                                    <  2:4695149/93‑1 (MQ=255)

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

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

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

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

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

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

                                                        tatGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGGTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCt                                 >  2:271786/1‑93 (MQ=255)

                                                         atGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTg                                >  2:236607/1‑93 (MQ=255)

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

                                                          tGTTTATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGt                               <  2:3080603/93‑1 (MQ=255)

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

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

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

                                                            tttATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGttg                             >  1:2184049/1‑91 (MQ=255)

                                                             ttATGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGc                            <  1:140961/93‑1 (MQ=255)

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

                                                               aTGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCAc                          >  2:2683201/1‑93 (MQ=255)

                                                                tGACCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACt                         >  2:6344650/1‑93 (MQ=255)

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

                                                                  aCCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTgg                       <  1:5091859/93‑1 (MQ=255)

                                                                  aCCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTgg                       <  2:3352536/93‑1 (MQ=255)

                                                                  aCCTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTgg                       <  1:1360602/93‑1 (MQ=255)

                                                                   ccTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTggg                      <  2:5661487/93‑1 (MQ=255)

                                                                   ccTTCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATAGCCATCGCCCTGTCGCACTggg                      >  1:4888593/1‑93 (MQ=255)

                                                                      tCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGt                   <  2:1408973/93‑1 (MQ=255)

                                                                      tCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGt                   <  2:6126599/93‑1 (MQ=255)

                                                                      tCCCGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGt                   <  1:378285/93‑1 (MQ=255)

                                                                       cccGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTg                  <  2:2454143/93‑1 (MQ=255)

                                                                        ccGACCTCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTgg                 >  2:895462/1‑93 (MQ=255)

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

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

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

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

                                                                              tCGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTATGCa           >  1:333897/1‑93 (MQ=255)

                                                                               cGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTATGCAc          >  2:2163613/1‑93 (MQ=255)

                                                                               cGATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGGATGCAc          >  2:172330/1‑93 (MQ=255)

                                                                                gATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTATGCACg         >  1:2533707/1‑93 (MQ=255)

                                                                                gATTCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTATGCACg         >  1:4358801/1‑93 (MQ=255)

                                                                                  ttCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTATGCACGCa       >  2:5374207/1‑93 (MQ=255)

                                                                                  ttCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTATGCACGCa       >  1:4207232/1‑93 (MQ=255)

                                                                                   tCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTATGCACGCAt      <  1:362513/93‑1 (MQ=255)

                                                                                   tCTGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTATGCACGCAt      <  1:2061371/93‑1 (MQ=255)

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

                                                                                     tGTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACAAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTATGCACGCATgg    <  2:4227460/93‑1 (MQ=255)

                                                                                      gTCGCTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTATGCACGCATGGt   >  1:245738/1‑93 (MQ=255)

GGCAGCGCCGGGTTGATTGGCGGGCGCGTTGATCAGGCCACCATGCGCGTCGCCGATATGTTTATGACCTTCCCGACCTCGATTCTGTCGTTCTTTATGGTTGGCGTGCTCGGCACCGGGCTGACCAACGTAATTATCGCCATCGCCCTGTCGCACTGGGCGTGGTATGCACGCATGGTG  >  NC_000913/3614092‑3614271

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