BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,563,580	T→C	D442G (GAC→GGC)	b2451 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,563,580	0	T	C	92.6%	512.8 / 24.2	176	D442G (GAC→GGC)	b2451	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (7/6); new base C (88/75); total (95/81)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.69e-01

Only 100 of 192 total aligned reads displayed.

GAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑T‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAAAG  >  NC_000913/2563490‑2563669

gAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑c‑cc                                                                                         >  2:1698558/1‑93 (MQ=255)

 aaTACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑c‑ccc                                                                                        >  2:2557067/1‑93 (MQ=255)

  aTACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑c‑ccgc                                                                                       >  2:264577/1‑91 (MQ=255)

  aTACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑cccc                                                                                       >  1:5765460/1‑93 (MQ=255)

  aTACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑cccc                                                                                       >  1:3765657/1‑93 (MQ=255)

  aTACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑cccc                                                                                       >  1:2791711/1‑93 (MQ=255)

   tACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑c‑ccccc                                                                                      >  2:764064/1‑93 (MQ=255)

   tACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑c‑ccccc                                                                                      >  2:52208/1‑93 (MQ=255)

   tACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGAAATAG‑c‑ccccc                                                                                      >  2:1630172/1‑93 (MQ=255)

     ccGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCcgc                                                                                    >  1:2470063/1‑93 (MQ=255)

      cGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCcgcg                                                                                   >  1:4902208/1‑93 (MQ=255)

        cATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCcgcgcg                                                                                 >  1:5360204/1‑93 (MQ=255)

        cATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCcgcgcg                                                                                 >  1:4796599/1‑93 (MQ=255)

          ttCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGa                                                                               >  2:1400727/1‑93 (MQ=255)

          ttCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGa                                                                               >  2:1853851/1‑93 (MQ=255)

           tCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGAc                                                                              >  1:6053224/1‑93 (MQ=255)

            cAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACa                                                                             >  1:4210406/1‑93 (MQ=255)

             aggaaggaAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACaa                                                                            >  1:5106291/1‑93 (MQ=255)

              ggaaggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAAt                                                                           >  2:4775264/1‑93 (MQ=255)

                aaggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCa                                                                         >  1:4948040/1‑93 (MQ=255)

                  ggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACt                                                                       >  1:2001895/1‑93 (MQ=255)

                   gaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACtt                                                                      >  2:2343482/1‑93 (MQ=255)

                   gaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACtt                                                                      <  2:3576055/93‑1 (MQ=255)

                     aaTGCGAGTTATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCg                                                                    <  2:5132987/93‑1 (MQ=255)

                     aaTGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCg                                                                    <  2:5339215/93‑1 (MQ=255)

                     aaTGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCg                                                                    <  1:3406968/93‑1 (MQ=255)

                      aTGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGt                                                                   >  2:5708325/1‑93 (MQ=255)

                      aTGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAACAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGt                                                                   >  1:713796/1‑93 (MQ=255)

                       tGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTc                                                                  <  1:5476698/93‑1 (MQ=255)

                        gCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCa                                                                 <  1:5616070/93‑1 (MQ=255)

                        gCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCa                                                                 <  2:3932874/93‑1 (MQ=255)

                        gCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCa                                                                 >  1:3491304/1‑93 (MQ=255)

                         cGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCaa                                                                <  2:5347830/93‑1 (MQ=255)

                          gAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCCCCCGCGCGGACAATCACTTCGTCaa                                                                <  1:2846397/93‑1 (MQ=255)

                          gAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAAt                                                               >  2:1909198/1‑93 (MQ=255)

                               aTTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTg                                                          >  1:4551786/1‑93 (MQ=255)

                                 ttCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGcc                                                        <  1:56758/93‑1 (MQ=255)

                                  tcaccgtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCa                                                       >  2:344328/1‑93 (MQ=255)

                                    accgtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑T‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAAc                                                     >  2:2194351/1‑93 (MQ=255)

                                     ccgtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACg                                                    >  2:109975/1‑93 (MQ=255)

                                       gtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACggg                                                  <  2:2660731/93‑1 (MQ=255)

                                       ctccccgGCAACACCAGTCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCCCGTACAATCACTTCGTCAATGACTGCCAACggg                                                  <  1:1466391/89‑1 (MQ=255)

                                         caccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAg                                                >  2:5951552/1‑93 (MQ=255)

                                          cccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAgt                                               <  2:3849983/92‑1 (MQ=255)

                                           ccgGAACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCGCTTCGTCAATGACTGCCAACGGGAgtt                                              <  2:1413403/93‑1 (MQ=255)

                                             ggCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAgttgt                                            <  2:4103841/93‑1 (MQ=255)

                                              gCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAgttgtt                                           <  1:4308537/93‑1 (MQ=255)

                                               cacaACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGttgttg                                          <  1:38256/93‑1 (MQ=255)

                                                acaACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTtgttgt                                         <  2:4250798/93‑1 (MQ=255)

                                                acaACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTtgttgt                                         >  2:5748131/1‑93 (MQ=255)

                                                acaACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTtgttgt                                         >  1:4785770/1‑93 (MQ=255)

                                                 caACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTa                                        >  2:5845510/1‑93 (MQ=255)

                                                 caACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTa                                        <  1:1055657/93‑1 (MQ=255)

                                                  aaCCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAg                                       >  2:828112/1‑93 (MQ=255)

                                                  aaCCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAg                                       >  1:457342/1‑93 (MQ=255)

                                                   aCCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGc                                      <  2:4763081/93‑1 (MQ=255)

                                                   aCCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGc                                      <  1:129661/93‑1 (MQ=255)

                                                    ccGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCt                                     <  1:4734598/93‑1 (MQ=255)

                                                     cGATCCGCCAAAAGGGGGCGTACCGATGTCGATCTAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCtg                                    <  2:5485157/93‑1 (MQ=255)

                                                     cGATCCGCCAAAAAGAGGCGTACGAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCtg                                    <  1:3509136/93‑1 (MQ=255)

                                                      gATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCtgt                                   <  2:3243144/93‑1 (MQ=255)

                                                       aTCCACCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCtgtg                                  <  1:4953651/93‑1 (MQ=255)

                                                        tCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTgg                                 >  1:4607891/1‑93 (MQ=255)

                                                        tCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTgg                                 <  1:5406372/93‑1 (MQ=255)

                                                        tCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTgg                                 >  1:5064782/1‑93 (MQ=255)

                                                         ccGCCAAAAAGAGGCGTCCCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCCCTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTggg                                <  2:4925980/93‑1 (MQ=255)

                                                         ccGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTggg                                <  2:278775/93‑1 (MQ=255)

                                                         ccGCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTggg                                <  1:1747704/93‑1 (MQ=255)

                                                           gCCAAAAAGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGgcg                              >  2:4653129/1‑93 (MQ=255)

                                                              aaaaaGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGGGGACAATCACTTCGTCAATGACTGCCAACGGGCGTTGTTGTAGCTGTGGgcgcca                           >  1:3255589/1‑91 (MQ=255)

                                                              aaaaaGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGcaa                           >  2:1081044/1‑93 (MQ=255)

                                                               aaaaGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGcaac                          >  2:2630211/1‑93 (MQ=255)

                                                               aaaaGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGcaac                          >  2:5923925/1‑93 (MQ=255)

                                                               aaaaGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGcaac                          >  1:1290580/1‑93 (MQ=255)

                                                               aaaaGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGcaac                          <  1:4156790/93‑1 (MQ=255)

                                                               aaaaGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGcaac                          <  1:145008/93‑1 (MQ=255)

                                                                aaaGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGcaaca                         >  2:1763911/1‑93 (MQ=255)

                                                                 aaGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGcaacga                        >  2:4563168/1‑91 (MQ=255)

                                                                 aaGAGGCGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGcaacaa                        <  2:3710960/93‑1 (MQ=255)

                                                                       cGTACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGcc                  <  2:2742735/93‑1 (MQ=255)

                                                                         tACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCa                <  2:1972830/93‑1 (MQ=255)

                                                                         tACCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCa                <  2:4818824/93‑1 (MQ=255)

                                                                          aCCAATGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCag               <  2:894500/93‑1 (MQ=255)

                                                                             aaTGTCGATATAGCC‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCagag             <  2:5189169/93‑1 (MQ=255)

                                                                             aaTGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCACCAACATGCCCAGAGc            >  2:1301626/1‑93 (MQ=255)

                                                                             aaTGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGc            >  2:454045/1‑93 (MQ=255)

                                                                               tGTCGATATAG‑T‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCtt          >  1:787204/1‑93 (MQ=255)

                                                                               tGTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCtt          >  1:1960814/1‑93 (MQ=255)

                                                                                gTCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCttt         >  1:1041218/1‑93 (MQ=255)

                                                                                 tCGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTa        <  2:671566/93‑1 (MQ=255)

                                                                                  cGATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTAc       >  2:5684670/1‑93 (MQ=255)

                                                                                   gATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTAcc      <  1:2792178/93‑1 (MQ=255)

                                                                                   gATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTAcc      >  1:1791637/1‑93 (MQ=255)

                                                                                   gATATAG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTAcc      <  1:3727323/93‑1 (MQ=255)

                                                                                      ataG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCaaa   <  1:2037604/93‑1 (MQ=255)

                                                                                      ataG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCaaa   <  2:3030979/93‑1 (MQ=255)

                                                                                       taG‑T‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAAAg  >  2:426849/1‑93 (MQ=255)

                                                                                       taG‑T‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAAAg  <  1:5590466/93‑1 (MQ=255)

                                                                                       taG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACTAAAg  >  2:1527599/1‑93 (MQ=255)

                                                                                       taG‑C‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAAAg  <  1:5313600/93‑1 (MQ=255)

GAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑T‑CCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAAAG  >  NC_000913/2563490‑2563669

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