BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,554,088	+C	coding (2/624 nt)	yhgKJ ←	RNA phosphage cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,554,088	1	.	C	99.5%	1355.8 / NA	368	S133W (TCG→TGG)	yhgKJ	RNA phosphage cyclase
Reads supporting (aligned to +/- strand): ref base . (0/0); major base C (169/197); minor base A (1/1); total (170/198)
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 429 total aligned reads displayed.

CTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑G‑AGGCGGCGCCGACGGGTTATCGGTCGCCG‑CGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCA  >  NC_000913/3553993‑3554177

cTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTGCTATTTTCGCCAGCAGCGGCTCCAGCACCCGTCGGATAAAATCC‑g                                                                                             >  1:6686705/1‑96 (MQ=255)

   cGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑gcag                                                                                          <  2:5527500/96‑4 (MQ=255)

   cGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑gcag                                                                                          >  1:12158698/1‑93 (MQ=255)

      cAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAggcg                                                                                       >  2:12268977/1‑96 (MQ=255)

        gggTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCCCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAggcggc                                                                                     <  1:12472786/96‑1 (MQ=255)

         ggTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGgcggcg                                                                                    <  2:3458373/96‑1 (MQ=255)

          gTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAggcggcac                                                                                   <  1:10514141/96‑3 (MQ=255)

          gTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGgcgc                                                                                   <  1:4440159/96‑1 (MQ=255)

           tAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGcc                                                                                  >  1:8217483/1‑96 (MQ=255)

            aaaaaCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCg                                                                                 <  1:10927350/96‑1 (MQ=255)

            aaaaaCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCg                                                                                 >  2:11571085/1‑96 (MQ=255)

              aaaCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGAc                                                                               >  1:590053/1‑96 (MQ=255)

              aaaCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGAc                                                                               >  2:11329717/1‑96 (MQ=255)

                 ccGTCGCGTAACAGCGTGGCTTGCGGATGAATTCCTATTTTCGCCAGCAGCGCCCCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACggg                                                                            <  1:4165762/96‑1 (MQ=255)

                  cGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTGCAGCACCCGGCGGATAAAATTC‑GCAAGCGGCGCCGACgggc                                                                           >  1:6291140/1‑95 (MQ=255)

                  cGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGt                                                                           >  1:7939639/1‑96 (MQ=255)

                   gTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGtt                                                                          <  1:5554366/96‑1 (MQ=255)

                   gTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGtt                                                                          <  1:4711913/96‑1 (MQ=255)

                    tGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTa                                                                         >  2:7182624/1‑96 (MQ=255)

                     ggCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTAt                                                                        <  1:12745309/96‑1 (MQ=255)

                     ggCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGCTAt                                                                        >  1:7779620/1‑96 (MQ=255)

                      gCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATc                                                                       <  1:9501078/96‑1 (MQ=255)

                        gTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGGTATCgg                                                                     >  2:8929249/1‑96 (MQ=255)

                         tAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCgct                                                                    <  1:8993404/96‑3 (MQ=255)

                         tAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGt                                                                    >  2:10978904/1‑96 (MQ=255)

                          aaCAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAAGCGGCGCCGACGGGGTATTggcg                                                                   >  2:5764912/1‑94 (MQ=255)

                            cAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGtgcc                                                                 <  2:220853/96‑4 (MQ=255)

                             aGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGtgccg                                                                <  2:2756321/96‑5 (MQ=255)

                              gCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑c                                                             <  2:4868111/96‑1 (MQ=255)

                              gCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑c                                                             <  2:5364203/96‑1 (MQ=255)

                              gCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGGTATCGGTggcgc                                                               >  1:1334343/1‑94 (MQ=255)

                               cGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑cc                                                            <  2:6566968/96‑2 (MQ=255)

                               cGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATAGGt‑gccg‑cc                                                            <  1:9227724/96‑2 (MQ=255)

                               agTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑cc                                                            <  1:2085248/95‑2 (MQ=255)

                                gTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑ccg                                                           >  1:2453598/1‑94 (MQ=255)

                                  ggTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑ccgct                                                         <  2:8027267/96‑5 (MQ=255)

                                  ggTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑ccgct                                                         >  2:8366213/1‑92 (MQ=255)

                                       gCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCCCACtt                                                     >  1:4767286/1‑96 (MQ=255)

                                       gCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGGTATCGGT‑GCCGCCGCTCACtt                                                     >  1:7169867/1‑96 (MQ=255)

                                        cggatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGACGGCGGCGACGGTGTATCGGT‑GCCGCCGCTCACTTc                                                    >  1:10259785/3‑96 (MQ=255)

                                         tgatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCa                                                   >  2:14196973/1‑96 (MQ=255)

                                          gatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGTCGGATAAAATCC‑GCAGGTGGCGCCGACGGGTTATCGGGGCCTC‑CGCTCACTTCaa                                                  >  1:13164499/1‑96 (MQ=255)

                                           atgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAac                                                 >  1:9957006/1‑96 (MQ=255)

                                            tgaATTCCTACTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAACCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAaca                                                <  2:10716143/96‑1 (MQ=255)

                                              caTTCCTATTTTCGCCAGCAGCGGCTCCAGCCCCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACAcg                                              <  1:5076707/95‑1 (MQ=255)

                                              aaTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACAcg                                              <  2:7063232/96‑1 (MQ=255)

                                              aaTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACAcg                                              <  1:467008/96‑1 (MQ=255)

                                                ttCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATACAATCC‑GCAGGCGGCGCCGACGGGCTATCGG‑GGCCGCCGCTCACATCAACAcgcg                                            >  2:12747598/1‑96 (MQ=255)

                                                ttCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACAcgcg                                            <  2:4374848/96‑1 (MQ=255)

                                                   cTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAg                                         >  1:6830631/1‑96 (MQ=255)

                                                   cTATTTTCGCCAGCAGCGCCTCAAGCGCCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGCTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAg                                         <  1:706284/96‑1 (MQ=255)

                                                    tATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAgg                                        >  2:10495964/1‑96 (MQ=255)

                                                    tATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGG‑GGCCGCCGCTCACTTCAACACGCGAAgg                                        >  1:7408961/1‑96 (MQ=255)

                                                     aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGt                                       <  2:5899062/96‑1 (MQ=255)

                                                     aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGt                                       >  1:3085209/1‑96 (MQ=255)

                                                      ttgtCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGTTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTc                                      <  2:195256/93‑1 (MQ=255)

                                                         tCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGCCCAt                                   >  1:3488504/1‑96 (MQ=255)

                                                           gCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCg                                 >  2:10364858/1‑96 (MQ=255)

                                                            ccAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCgg                                <  2:1632156/96‑1 (MQ=255)

                                                             cagcagCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGc                               <  2:5129895/96‑1 (MQ=255)

                                                             cagcagCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGc                               <  1:14066065/96‑1 (MQ=255)

                                                             cagcagCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGc                               >  1:8724433/1‑96 (MQ=255)

                                                             cagcagCGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGCGGTATCTGT‑GCCGCCGCTCACTTCAACACCCGAAGGGCCATCGGc                               >  1:6560582/1‑96 (MQ=255)

                                                               gcagcGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaa                             >  1:4354928/1‑96 (MQ=255)

                                                                 ggcGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAc                           >  1:8875279/2‑96 (MQ=255)

                                                                 agcGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAc                           >  2:3746571/1‑96 (MQ=255)

                                                                 agcGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGGTATCGGT‑GCCGCCGCTCACTTCAACACGCGGAGGTCCATCGGGAAAc                           >  1:2666795/1‑96 (MQ=255)

                                                                  gcGGCTCCAGCCCCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAcc                          <  1:6487122/96‑1 (MQ=255)

                                                                  gcGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAcc                          <  1:11600048/96‑1 (MQ=255)

                                                                   cGGCTCCAGCCCCCGGCGGATAAACTCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACca                         <  1:389975/96‑1 (MQ=255)

                                                                   cGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACca                         <  1:9741168/96‑1 (MQ=255)

                                                                   cGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACca                         >  2:13462408/1‑96 (MQ=255)

                                                                   cGGCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACca                         >  1:11364755/1‑96 (MQ=255)

                                                                   cGGATCCAGCACCCGGCGGATAAAAACC‑GCAGGCGGCGCAGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACca                         <  1:1229798/96‑1 (MQ=255)

                                                                    ggCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATGGGCAAACcac                        >  1:11932144/1‑96 (MQ=255)

                                                                    ggCTCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGGTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAAGTCCATCGGCAAACcac                        >  1:1850661/1‑96 (MQ=255)

                                                                       tCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgc                     <  2:5586098/96‑1 (MQ=255)

                                                                       tCCAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGACAAACCACAgc                     <  2:2963352/96‑1 (MQ=255)

                                                                       tCCAGCACCCGGCGGATAAAATCC‑GAAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgc                     <  2:7711364/96‑1 (MQ=255)

                                                                       tCCAGCACCCGGCGGATAAAACCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgc                     <  2:1214023/96‑1 (MQ=255)

                                                                         gaGCACCCGGCGGGTAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgcgc                   >  1:11421765/2‑96 (MQ=255)

                                                                         cAGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgcgc                   <  2:13748889/96‑1 (MQ=255)

                                                                          aGCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGcgcg                  >  1:10034307/1‑96 (MQ=255)

                                                                           gCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCgg                 <  1:10364858/96‑1 (MQ=255)

                                                                           gCACCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCgg                 <  2:2527777/96‑1 (MQ=255)

                                                                             aCCCGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgc               >  2:7896349/1‑96 (MQ=255)

                                                                              cccGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGGTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAACACAGCGCGGgca              >  2:626345/1‑96 (MQ=255)

                                                                               cgggcgcATAAAATCCGG‑AGGCGGCGCCGACGGGTTGTCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcag             <  2:9564948/94‑1 (MQ=255)

                                                                               ccGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcag             <  2:9665984/96‑1 (MQ=255)

                                                                               ccGGCGGATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcag             <  1:1553999/96‑1 (MQ=255)

                                                                                 ggcggATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagca           <  2:9211370/96‑1 (MQ=255)

                                                                                 ggcggATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTAAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagca           <  2:13279620/96‑1 (MQ=255)

                                                                                  gcggATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc          <  2:11884741/96‑1 (MQ=255)

                                                                                   cggATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAcc         <  2:2342112/96‑1 (MQ=255)

                                                                                    ggATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCg        <  2:466987/96‑1 (MQ=255)

                                                                                     gATAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGAGCAGCACCGt       >  2:11507324/1‑96 (MQ=255)

                                                                                       tAAAGTCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGttt     <  1:12950587/96‑1 (MQ=255)

                                                                                       tAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGttt     <  2:7408961/96‑1 (MQ=255)

                                                                                       tAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGttt     >  2:861448/1‑96 (MQ=255)

                                                                                       tAAAATCC‑GCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGttt     >  2:1398216/1‑96 (MQ=255)

CTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCC‑G‑AGGCGGCGCCGACGGGTTATCGGTCGCCG‑CGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCA  >  NC_000913/3553993‑3554177

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 22 ≤` `ATCG/ATCG` `< 31 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`