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	99.5%	604.3 / ‑5.6	191	D442G (GAC→GGC)	b2451	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/1); major base C (100/89); minor base A (1/0); total (101/90)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 228 total aligned reads displayed.

GCGAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑T‑C‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAA  >  NC_000913/2563488‑2563667

gCGAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑t                                                                                           >  1:3626260/1‑93 (MQ=255)

  taaTACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑T‑c‑c                                                                                       <  1:2815025/92‑1 (MQ=255)

  gAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑T‑c‑c                                                                                       >  2:5750727/1‑93 (MQ=255)

   aaTACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑T‑c‑cc                                                                                      >  2:2511567/1‑93 (MQ=255)

   aaTACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑T‑c‑cc                                                                                      <  2:3058826/93‑1 (MQ=255)

    aTACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑T‑c‑ccc                                                                                     >  2:74567/1‑93 (MQ=255)

      aCCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑T‑c‑ccccc                                                                                   >  2:5021583/1‑92 (MQ=255)

      aCCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑T‑c‑ccccc                                                                                   <  2:5141288/93‑2 (MQ=255)

      aCCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑T‑c‑ccccc                                                                                   <  1:2633037/93‑2 (MQ=255)

       ccGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑T‑c‑cccccg                                                                                  <  2:3029765/93‑3 (MQ=255)

       ccGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑T‑c‑cccccg                                                                                  <  1:1031829/93‑3 (MQ=255)

        cGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑T‑c‑cccccgc                                                                                 <  1:4721237/93‑4 (MQ=255)

         gCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑T‑c‑cccccgcg                                                                                >  1:3119190/1‑89 (MQ=255)

         gCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑T‑c‑cccccgcg                                                                                >  1:4424398/1‑89 (MQ=255)

               aggaaggaAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACa                                                                           >  2:4497955/1‑93 (MQ=255)

               aggaaggaAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑T‑CTCCCCGCGCGGACa                                                                           >  2:1635730/1‑93 (MQ=255)

                 gaaggaaATGCGGGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCCTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAAt                                                                         >  2:1743414/1‑93 (MQ=255)

                 gaaggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAAt                                                                         <  2:708071/93‑1 (MQ=255)

                 gaaggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAAt                                                                         >  1:1997225/1‑93 (MQ=255)

                 gaaggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAAt                                                                         >  1:3201823/1‑93 (MQ=255)

                  aaggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATc                                                                        >  1:6109050/1‑93 (MQ=255)

                   aggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCa                                                                       >  2:6043182/1‑93 (MQ=255)

                   aggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCa                                                                       <  1:4076397/93‑1 (MQ=255)

                    ggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCAc                                                                      >  2:20657/1‑93 (MQ=255)

                    ggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCAc                                                                      <  1:929671/93‑1 (MQ=255)

                     gaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACt                                                                     <  1:831855/93‑1 (MQ=255)

                      gaaTGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAGAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACtt                                                                    <  2:197981/92‑1 (MQ=255)

                       aaTGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTc                                                                   <  2:1694383/93‑1 (MQ=255)

                       aaTGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTc                                                                   <  1:4818580/93‑1 (MQ=255)

                       aaTGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTc                                                                   <  2:2892475/93‑1 (MQ=255)

                        aTGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCg                                                                  <  2:1317415/93‑1 (MQ=255)

                        aTGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCg                                                                  >  1:2483571/1‑93 (MQ=255)

                         tGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGt                                                                 >  2:5941947/1‑93 (MQ=255)

                           cGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCa                                                               >  2:1108673/1‑93 (MQ=255)

                            gAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCaa                                                              <  1:3050465/93‑1 (MQ=255)

                            gAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCaa                                                              <  1:5000065/93‑1 (MQ=255)

                             aGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAAt                                                             >  1:3401930/1‑93 (MQ=255)

                             aGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGGCAAt                                                             >  1:5994714/1‑93 (MQ=255)

                             aGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAAACAACTCGTCCAt                                                             >  1:3501444/1‑93 (MQ=255)

                               tGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGa                                                           >  1:2508974/1‑93 (MQ=255)

                               tGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGa                                                           <  2:608020/93‑1 (MQ=255)

                                 aTTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACt                                                         <  2:3559579/93‑1 (MQ=255)

                                 aTTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACt                                                         >  1:5882658/1‑93 (MQ=255)

                                  tttCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTg                                                        >  1:5825035/1‑93 (MQ=255)

                                  tttCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTg                                                        >  1:2433607/1‑93 (MQ=255)

                                   ttCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGc                                                       <  2:4631474/93‑1 (MQ=255)

                                    tcaccgtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGcc                                                      >  2:319107/1‑93 (MQ=255)

                                    tcaccgtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGcc                                                      >  2:2992057/1‑93 (MQ=255)

                                      accgtcaccgGCACGACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGCC‑C‑CCCCGCGCGGCCAATCACTTCGTCAATGACTGCCaa                                                    <  2:4076781/93‑1 (MQ=255)

                                      accgtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCaa                                                    <  2:4023310/93‑1 (MQ=255)

                                       ccgtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAAc                                                   <  2:5133861/93‑1 (MQ=255)

                                       ccgtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAAc                                                   >  1:1193615/1‑93 (MQ=255)

                                        cgtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACg                                                  <  2:1940395/93‑1 (MQ=255)

                                        cgtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACg                                                  >  1:3743354/1‑93 (MQ=255)

                                         gtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCTATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACgg                                                 <  1:2781094/93‑1 (MQ=255)

                                           caccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGa                                               <  1:385066/93‑1 (MQ=255)

                                             ccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAgt                                             >  2:3416505/1‑93 (MQ=255)

                                             ccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCCCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAgt                                             >  2:3431534/1‑93 (MQ=255)

                                                   caACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTtgttgt                                       <  1:5473901/93‑1 (MQ=255)

                                                    aaCCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTtgttgtg                                      >  2:1198288/1‑92 (MQ=255)

                                                    aaCCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTa                                      <  1:1161156/93‑1 (MQ=255)

                                                    aaCCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTa                                      <  1:1751113/93‑1 (MQ=255)

                                                     aCCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAg                                     >  2:3663456/1‑93 (MQ=255)

                                                     aCCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAg                                     >  1:1671782/1‑93 (MQ=255)

                                                     aCCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAg                                     >  2:3618612/1‑93 (MQ=255)

                                                       cGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCt                                   >  2:2730300/1‑93 (MQ=255)

                                                          tCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCtgtg                                <  2:1772871/93‑1 (MQ=255)

                                                          tCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCtgtg                                <  2:3477342/93‑1 (MQ=255)

                                                           ccGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTgg                               <  2:3831530/93‑1 (MQ=255)

                                                           ccGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTgg                               <  2:1249927/93‑1 (MQ=255)

                                                            cGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTggg                              <  1:1587919/93‑1 (MQ=255)

                                                               cAAAAAGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGgcgc                           >  1:5589678/1‑93 (MQ=255)

                                                                aaaaaGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGca                          >  1:4888285/1‑93 (MQ=255)

                                                                aaaaaGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGca                          >  1:1299842/1‑93 (MQ=255)

                                                                aaaaaGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGca                          >  2:4776964/1‑93 (MQ=255)

                                                                aaaaaGAGGCGTACCAATGTAGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGca                          >  2:6119166/1‑93 (MQ=255)

                                                                 aaaaGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGcaa                         >  2:665553/1‑93 (MQ=255)

                                                                  aaaGAGGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGcaca                        >  1:1908977/1‑91 (MQ=255)

                                                                     gagGCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCaacaac                     >  2:1819206/1‑93 (MQ=255)

                                                                       ggCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACAt                   >  2:4308754/1‑93 (MQ=255)

                                                                        gCGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATg                  <  1:2064135/93‑1 (MQ=255)

                                                                         cGTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGc                 >  2:4826862/1‑93 (MQ=255)

                                                                          gTACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGcc                >  1:1574894/1‑93 (MQ=255)

                                                                           tACCAATGTCGATATAG‑TCC‑CCCCGCGCTGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGccc               <  2:5948179/93‑1 (MQ=255)

                                                                           tACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGccc               >  1:529829/1‑93 (MQ=255)

                                                                           tACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGccc               >  2:2291931/1‑93 (MQ=255)

                                                                           tACCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGATGTGGGCGCAACAACATGccc               >  2:1610845/1‑93 (MQ=255)

                                                                            tccAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCa              <  1:2118178/92‑1 (MQ=255)

                                                                            aCCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCa              >  1:4231566/1‑93 (MQ=255)

                                                                            aCCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCa              <  1:4718917/93‑1 (MQ=255)

                                                                            aCCAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCa              <  1:6247556/93‑1 (MQ=255)

                                                                              cAATGTCGATATAG‑TCC‑CCCCGCGCGGCCAATCACTTCGTCAATTGCTGCCCACGGGAGTTGTTCTAGCTGGGGCCGCAACATGATGCCCaaa            >  2:5092604/1‑91 (MQ=255)

                                                                              cAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCaga            <  1:4471778/93‑1 (MQ=255)

                                                                              cAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCaga            <  1:4178193/93‑1 (MQ=255)

                                                                              cAATGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCaga            <  1:5422432/93‑1 (MQ=255)

                                                                                aTGTCGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGc          <  1:5695128/93‑1 (MQ=255)

                                                                                    cGATATAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTa      <  1:4533707/93‑1 (MQ=255)

                                                                                      atatAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTAcc    <  1:1518794/93‑1 (MQ=255)

                                                                                      atatAG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTAcc    <  1:21523/93‑1 (MQ=255)

                                                                                        ataG‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCaa  >  2:5864712/1‑93 (MQ=255)

GCGAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑T‑C‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAA  >  NC_000913/2563488‑2563667

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