BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,518,992	0	C	.	32.1%	262.6 / 118.6	106	pseudogene (272/504 nt)	insB1	IS1 transposase B,IS, phage, Tn, Transposon‑related functions, extrachromosomal, transposon related
Reads supporting (aligned to +/- strand): ref base C (35/37); new base . (19/15); total (54/52)
Fisher's exact test for biased strand distribution p-value = 5.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 110 total aligned reads displayed.

TGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACAT  >  NC_000913/4518859‑4519125

tGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCt                                                                                                                                  >  1:158060/1‑139 (MQ=255)

        cGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTc                                                                                                                          >  1:339629/1‑139 (MQ=255)

        cGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTc                                                                                                                          <  2:158060/139‑1 (MQ=255)

             gTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGt                                                                                                                     <  1:633706/139‑1 (MQ=255)

              tGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCg                                                                                                                                      >  2:838487/1‑121 (MQ=255)

              tGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCg                                                                                                                                      <  1:838487/121‑1 (MQ=255)

                                 tggctgTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                 <  2:1357753/139‑1 (MQ=255)

                                   gctgTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                               >  1:1138170/1‑139 (MQ=255)

                                             gcgTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGcc                                                                                                >  2:312016/1‑128 (MQ=255)

                                             gcgTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGcc                                                                                                <  1:312016/128‑1 (MQ=255)

                                             gcgTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAAt                                                                                     >  2:679905/1‑139 (MQ=255)

                                                 aTGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACg                                                                                 <  2:105715/139‑1 (MQ=255)

                                                   gACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGcc                                                                               >  2:229217/1‑139 (MQ=255)

                                                       ggATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGaa                                                                           <  2:1147268/139‑1 (MQ=255)

                                                           aCGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGta                                                                                                                    >  1:331235/1‑94 (MQ=255)

                                                           aCGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGta                                                                                                                    <  2:331235/94‑1 (MQ=255)

                                                            cGGAAGACGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGaaa                                                                     >  2:1039678/1‑139 (MQ=11)

                                                             ggaggaCGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGaaa                                                                     >  1:4908/1‑139 (MQ=255)

                                                                  aCGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCg                                                                                                                         <  2:1058075/82‑1 (MQ=255)

                                                                  aCGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCg                                                                                                                         >  1:1058075/1‑82 (MQ=255)

                                                                    ggTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCAc                                                              >  1:51169/1‑139 (MQ=255)

                                                                       tgtgGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCg                                                                                                                         <  1:1007370/77‑1 (MQ=255)

                                                                       tgtgGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGcc                                                                                                <  1:217181/102‑1 (MQ=255)

                                                                       tgtgGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGcc                                                                                                >  2:217181/1‑102 (MQ=255)

                                                                         ttgcgcACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGcc                                                                                                <  1:773736/97‑1 (MQ=14)

                                                                         ttgcgcACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGcc                                                                                                >  2:773736/3‑99 (MQ=14)

                                                                         tgGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTc                                                                                                                          <  2:950863/74‑1 (MQ=255)

                                                                         tgGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTc                                                                                                                          >  1:950863/1‑74 (MQ=255)

                                                                           gcgcACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTc                                                                                                                                       <  2:1133405/59‑1 (MQ=255)

                                                                           gcgcACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTc                                                                                                                                       >  1:1133405/1‑59 (MQ=255)

                                                                           gcgcACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCg                                                                                                                         <  2:783323/73‑1 (MQ=255)

                                                                           gcgcACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCg                                                                                                                         >  1:651393/1‑73 (MQ=255)

                                                                           gcgcACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCg                                                                                                                         <  2:651393/73‑1 (MQ=255)

                                                                           gcgcACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCg                                                                                                                         >  1:783323/1‑73 (MQ=255)

                                                                           gcgcACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggct                                                                                                    <  2:729738/94‑1 (MQ=255)

                                                                             gcACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGc                                                     <  1:337676/139‑1 (MQ=255)

                                                                             gcACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGc                                                     <  2:825314/139‑1 (MQ=255)

                                                                              cACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTc                                                                                                                          >  1:484830/1‑69 (MQ=255)

                                                                              cACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTc                                                                                                                          <  2:484830/69‑1 (MQ=255)

                                                                                cgtattCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGAtggct                                                                                                    >  1:1408193/5‑88 (MQ=14)

                                                                                cGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCg                                                                                                                         >  1:1053609/1‑68 (MQ=255)

                                                                                cGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCg                                                                                                                         <  2:1053609/68‑1 (MQ=255)

                                                                                cGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGtataaa                                                                                                                <  2:767813/77‑3 (MQ=255)

                                                                                cGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGtataaa                                                                                                                >  1:767813/1‑75 (MQ=255)

                                                                                cGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGtataaa                                                                                                                <  2:691699/77‑3 (MQ=255)

                                                                                cGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGtataaa                                                                                                                >  1:691699/1‑75 (MQ=255)

                                                                                cGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                 <  1:213627/92‑1 (MQ=255)

                                                                                cGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                 >  2:213627/1‑92 (MQ=255)

                                                                                cGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGc                                                                                              <  2:607593/95‑1 (MQ=255)

                                                                                cGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGc                                                                                              >  1:607593/1‑95 (MQ=255)

                                                                                cGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAg                                                  <  1:1115527/139‑1 (MQ=255)

                                                                                cGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGAtggct                                                                                                    <  2:1408193/88‑1 (MQ=17)

                                                                                     tCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTa                                             <  1:679905/139‑1 (MQ=255)

                                                                                        gTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCgatatacg                                         >  1:1334811/1‑132 (MQ=14)

                                                                                             cGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAgg                                                                                                                            >  2:222251/1‑52 (MQ=255)

                                                                                             cGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAgg                                                                                                                            <  1:222251/52‑1 (MQ=255)

                                                                                             cGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTc                                                                                                                          <  2:542752/54‑1 (MQ=255)

                                                                                             cGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTc                                                                                                                          >  1:542752/1‑54 (MQ=255)

                                                                                             cGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGtataaa                                                                                                                <  1:1207149/64‑3 (MQ=255)

                                                                                             cGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGtataaa                                                                                                                >  2:1207149/1‑62 (MQ=255)

                                                                                             cGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTc                                                                                                                          >  1:891736/1‑53 (MQ=14)

                                                                                             cGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTc                                                                                                                          <  2:891736/53‑1 (MQ=11)

                                                                                             cGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                               <  1:218729/80‑1 (MQ=18)

                                                                                             cGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCt                                                                                             <  2:1325116/82‑1 (MQ=17)

                                                                                             cGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCt                                                                                             >  1:1325116/1‑82 (MQ=17)

                                                                                                aCTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCa                                  <  2:1173206/139‑1 (MQ=255)

                                                                                                  tctGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTc                                                                                                                                       >  1:1057449/1‑36 (MQ=255)

                                                                                                  tctGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTc                                                                                                                                       <  2:1057449/36‑1 (MQ=255)

                                                                                                     ggCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAg                             >  2:203042/1‑139 (MQ=255)

                                                                                                      gcgacgcTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCg                           >  1:1138092/7‑139 (MQ=11)

                                                                                                        cacaCTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGa                          <  1:203042/139‑1 (MQ=255)

                                                                                                         acgcTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGa                              >  1:866242/4‑133 (MQ=14)

                                                                                                         acgcTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGa                              <  2:866242/130‑1 (MQ=11)

                                                                                                          cgcTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGa                              >  2:825071/3‑132 (MQ=14)

                                                                                                          cgcTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGa                              <  1:825071/130‑1 (MQ=14)

                                                                                                          cacTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACa                        <  2:51169/139‑1 (MQ=255)

                                                                                                            cTGGAGTGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGta                                                                                                                    <  1:957570/44‑1 (MQ=255)

                                                                                                            cTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                 <  2:197708/63‑1 (MQ=37)

                                                                                                            cTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                 >  1:197708/1‑63 (MQ=37)

                                                                                                            cTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCggtataa                     >  2:1069770/1‑133 (MQ=11)

                                                                                                                aGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCg                                                                                                                         >  1:476547/1‑36 (MQ=255)

                                                                                                                aGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCg                                                                                                                         <  2:476547/36‑1 (MQ=255)

                                                                                                                aGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATaatct                  >  1:286504/1‑139 (MQ=255)

                                                                                                                aGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGt                                                                                                                     <  2:699850/39‑1 (MQ=255)

                                                                                                                aGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGt                                                                                                                     >  1:699850/1‑39 (MQ=255)

                                                                                                                    tcttctGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                 >  2:1012920/1‑55 (MQ=25)

                                                                                                                    tcttctGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                 <  1:1012920/55‑1 (MQ=37)

                                                                                                                    tcttctGAGCCTACTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggct                                                                                                    >  2:251449/1‑53 (MQ=255)

                                                                                                                      ttctGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAg                                                                    <  2:659924/82‑1 (MQ=255)

                                                                                                                      ttctGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAg                                                                    >  1:659924/1‑82 (MQ=255)

                                                                                                                         tGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGAtggct                                                                                                    >  1:87451/1‑47 (MQ=255)

                                                                                                                         tGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGAgggct                                                                                                    <  2:87451/47‑1 (MQ=255)

                                                                                                                              ctgctgT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCa                                                                                   <  2:472835/59‑1 (MQ=21)

                                                                                                                              ctgctgT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCa                                                                                   >  1:472835/1‑59 (MQ=21)

                                                                                                                               tgctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGacaaca   >  2:1366836/1‑139 (MQ=255)

                                                                                                                                gctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATaatctgaat              <  2:43486/127‑1 (MQ=255)

                                                                                                                                gctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATaatctgaat              >  1:43486/1‑127 (MQ=255)

                                                                                                                                gctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACAt  >  2:340977/1‑139 (MQ=255)

                                                                                                                                 ctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggct                                                                                                    >  2:67230/1‑40 (MQ=255)

                                                                                                                                 ctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggct                                                                                                    <  1:67230/40‑1 (MQ=255)

TGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACAT  >  NC_000913/4518859‑4519125

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