BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,563,585	(C)_5→6	coding (1320/1404 nt)	b2451 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,563,580	1	.	C	94.7%	285.7 / 8.9	95	coding (1325/1404 nt)	b2451	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (4/1); new base C (43/47); total (47/48)
Fisher's exact test for biased strand distribution p-value = 2.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

Only 100 of 115 total aligned reads displayed.

GCGAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGT‑C‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAA  >  NC_000913/2563488‑2563667

 cGAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGT‑c                                                                                         <  1:3055922/93‑1 (MQ=255)

  gAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGT‑c‑c                                                                                       >  1:4175523/1‑93 (MQ=255)

  gAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGT‑c‑c                                                                                       >  1:4656925/1‑93 (MQ=255)

   aaTACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGT‑c‑cc                                                                                      >  1:2921193/1‑93 (MQ=255)

    aTACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGT‑c‑ccc                                                                                     >  1:1711471/1‑93 (MQ=255)

    aTACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGT‑c‑ccc                                                                                     >  2:604379/1‑93 (MQ=255)

    aTACCGCATTCAGGAAGGAAATGCGAGCGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGT‑c‑ccc                                                                                     >  2:676613/1‑93 (MQ=255)

     tACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGT‑c‑cccc                                                                                    <  1:3778966/93‑1 (MQ=255)

     tACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGT‑c‑cccc                                                                                    >  1:1819887/1‑93 (MQ=255)

     tACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGT‑c‑cccc                                                                                    >  2:1635309/1‑93 (MQ=255)

     aaCCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGT‑c‑cccc                                                                                    <  1:168059/92‑1 (MQ=255)

      aCCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGT‑c‑ccccc                                                                                   <  1:3889415/93‑2 (MQ=255)

         gCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGT‑c‑cccccgcg                                                                                >  2:2870534/1‑89 (MQ=255)

         gCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGT‑c‑cccccgcg                                                                                <  2:619274/93‑5 (MQ=255)

         gCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGT‑c‑cccccgcg                                                                                >  1:1911227/1‑89 (MQ=255)

               aggaaggaAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACa                                                                           >  1:56429/1‑93 (MQ=255)

                ggaaggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACaa                                                                          <  2:3613325/93‑1 (MQ=255)

                ggaaggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACaa                                                                          >  1:2664900/1‑93 (MQ=255)

                ggaaggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACaa                                                                          >  2:2939936/1‑93 (MQ=255)

                 gaaggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAAt                                                                         <  1:4739628/93‑1 (MQ=255)

                 gaaggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAAt                                                                         <  2:2374809/93‑1 (MQ=255)

                 gaaggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAAt                                                                         <  2:20233/93‑1 (MQ=255)

                   aggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCa                                                                       >  1:4025974/1‑93 (MQ=255)

                   aggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCa                                                                       <  1:1714023/93‑1 (MQ=255)

                   aggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCa                                                                       >  1:679299/1‑93 (MQ=255)

                   aggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCa                                                                       <  2:12583/93‑1 (MQ=255)

                   aggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCa                                                                       <  2:2277162/93‑1 (MQ=255)

                     gaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACt                                                                     <  1:1796950/93‑1 (MQ=255)

                         tGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGt                                                                 >  1:1141865/1‑93 (MQ=255)

                          gCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTc                                                                <  2:4634226/93‑1 (MQ=255)

                              gTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATg                                                            >  2:391044/1‑93 (MQ=255)

                              gTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATg                                                            >  1:1504290/1‑93 (MQ=255)

                                gATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGAc                                                          >  1:2816842/1‑93 (MQ=255)

                                gATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGAc                                                          >  1:4442000/1‑93 (MQ=255)

                                 aTTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACt                                                         <  2:4361313/93‑1 (MQ=255)

                                 aTTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACt                                                         <  1:1484190/93‑1 (MQ=255)

                                  tttCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACtc                                                        >  1:603690/1‑92 (MQ=255)

                                   ttCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGc                                                       >  1:4644140/1‑93 (MQ=255)

                                    tcaccgtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGcc                                                      <  1:3540705/93‑1 (MQ=255)

                                    tcaccgtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGcc                                                      <  2:2164206/93‑1 (MQ=255)

                                    tcaccgtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGcc                                                      >  2:2624404/1‑93 (MQ=255)

                                         gtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACgg                                                 >  2:3879288/1‑93 (MQ=255)

                                         gtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACgg                                                 >  2:1376905/1‑93 (MQ=255)

                                         gtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACgg                                                 <  2:281227/93‑1 (MQ=255)

                                         gtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGT‑C‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACggg                                                >  2:4109513/1‑93 (MQ=255)

                                          tcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACggg                                                <  1:3554553/93‑1 (MQ=255)

                                          tcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACggg                                                <  2:2255042/93‑1 (MQ=255)

                                            cccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAg                                              <  2:2823657/92‑1 (MQ=255)

                                               ggCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAgttg                                           >  2:4301518/1‑93 (MQ=255)

                                                 cacaACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAgttgtt                                         <  2:4138152/93‑1 (MQ=255)

                                                  acaACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGttgttg                                        >  2:1946721/1‑93 (MQ=255)

                                                  acaACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGttgttg                                        >  1:2586335/1‑93 (MQ=255)

                                                   caACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTtgttgt                                       <  1:2038197/93‑1 (MQ=255)

                                                   caACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTtgttgt                                       <  1:2753001/93‑1 (MQ=255)

                                                   caACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTtgttgt                                       >  1:2689537/1‑93 (MQ=255)

                                                    aaCCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTa                                      <  2:214811/93‑1 (MQ=255)

                                                     aCCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAg                                     <  1:2852482/93‑1 (MQ=255)

                                                       cGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCt                                   <  1:3143773/93‑1 (MQ=255)

                                                        gATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCtg                                  <  1:1250298/93‑1 (MQ=255)

                                                        gATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCtg                                  >  2:168925/1‑93 (MQ=255)

                                                         aTCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCtgc                                 >  2:227935/1‑92 (MQ=255)

                                                           ccGCCAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAAACACTTCGTCAATAATTGCCAACGGGAGATGTTGTAGTAGTgg                               >  2:2973746/1‑93 (MQ=255)

                                                               cAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGgcgc                           >  2:1858344/1‑93 (MQ=255)

                                                               cAAAAAGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGgcgc                           >  1:4618155/1‑93 (MQ=255)

                                                                aaaaaGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGca                          >  1:2139268/1‑93 (MQ=255)

                                                                 aaaaGAGGCGTACCAATGTCGATATAGT‑C‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGcaac                        <  1:2326760/93‑1 (MQ=255)

                                                                  aaaGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGcaac                        >  2:1751708/1‑93 (MQ=255)

                                                                  aaaGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTATTGTAGCTGTGGGCGcaac                        <  1:2911263/93‑1 (MQ=255)

                                                                   aaGAGGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGcaaca                       >  2:2500258/1‑93 (MQ=255)

                                                                      agGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAacaaca                    <  1:2804802/93‑1 (MQ=255)

                                                                      agGCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAacaaca                    <  2:1048312/93‑1 (MQ=255)

                                                                      agGCGTACCAATGTCGATATAGTCC‑CCCCGCGCAGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAacaaca                    <  1:4577026/93‑1 (MQ=255)

                                                                       ggCGTAGCAATGTCGATATAGT‑CGCCCCGCGCGGACAATCCCTTCGTCAATGACTGCCAGCGGGAGTTGTTGTAGCTGTGGGCGCAACAACAt                   <  1:4249949/93‑1 (MQ=255)

                                                                       ggCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACAt                   <  1:4323184/93‑1 (MQ=255)

                                                                        gCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATg                  <  2:1763848/93‑1 (MQ=255)

                                                                        gCGTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATg                  <  2:1833505/93‑1 (MQ=255)

                                                                         ggTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGc                 <  1:3872002/92‑1 (MQ=255)

                                                                          gTACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGcc                >  2:3671652/1‑93 (MQ=255)

                                                                           tACCAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGccc               >  2:1732808/1‑93 (MQ=255)

                                                                            aCCAATGTCGATATAGTCC‑CCCCGCGCGTACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCa              <  1:3240403/93‑1 (MQ=255)

                                                                             ccAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCag             >  1:2863132/1‑93 (MQ=255)

                                                                             ccAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCag             <  2:601129/93‑1 (MQ=255)

                                                                             ccAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCag             <  2:4175523/93‑1 (MQ=255)

                                                                             ccAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCag             >  2:3503772/1‑93 (MQ=255)

                                                                             ccAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGAGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCag             >  2:1474492/1‑93 (MQ=255)

                                                                             ccAATGTCGATATAGT‑C‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGATGTTGTAGCTGTGGGCGCAACAACATGCCCaga            >  1:2690409/1‑93 (MQ=255)

                                                                              cAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCaga            <  2:3696256/93‑1 (MQ=255)

                                                                              cAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCaga            <  1:606031/93‑1 (MQ=255)

                                                                              cAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCaga            >  1:1867356/1‑93 (MQ=255)

                                                                              cAATGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCaga            <  1:3759448/93‑1 (MQ=255)

                                                                                aTGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGc          <  1:2596246/93‑1 (MQ=255)

                                                                                 tGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCt         <  1:2683305/93‑1 (MQ=255)

                                                                                 tGTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCt         >  1:544781/1‑93 (MQ=255)

                                                                                  gTCGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCtt        >  2:430426/1‑93 (MQ=255)

                                                                                    cGATATAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTa      <  2:4299399/93‑1 (MQ=255)

                                                                                      atatAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTAcc    >  2:1580664/1‑93 (MQ=255)

                                                                                      atatAGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTAcc    <  2:587137/93‑1 (MQ=255)

                                                                                       tataGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCa   >  1:3629229/1‑93 (MQ=255)

                                                                                       tataGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCa   >  2:3557668/1‑93 (MQ=255)

                                                                                        ataGTCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCaa  >  1:3179963/1‑93 (MQ=255)

GCGAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGT‑C‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAA  >  NC_000913/2563488‑2563667

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