BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,811,545	G→A	A303A (GCG→GCA)	emrB →	EmrB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,811,545	0	G	A	89.9%	707.9 / 43.6	247	A303A (GCG→GCA)	emrB	EmrB
Reads supporting (aligned to +/- strand): ref base G (9/16); new base A (101/121); total (110/137)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.79e-01

Only 100 of 259 total aligned reads displayed.

GCTGCTTGTGTATCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCGACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCGGCCGCTTCGCGCATAAACTGGATAT  >  NC_000913/2811454‑2811634

gctgctTGTGTATCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCGa                                                                                          >  2:1427146/1‑93 (MQ=255)

 ctgctTGTGTATCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAAc                                                                                         <  1:3264857/93‑1 (MQ=255)

   gctTGTGTATCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCt                                                                                       >  2:5005730/1‑93 (MQ=255)

   gctTGTGTATCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCt                                                                                       >  2:5447793/1‑93 (MQ=255)

     ttGTGTATCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTgg                                                                                     >  2:1484426/1‑93 (MQ=255)

     ttGTGTATCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTgg                                                                                     <  1:4381830/93‑1 (MQ=255)

      tgtgTATCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTggg                                                                                    >  2:9609075/1‑93 (MQ=255)

      tgtgTATCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTggg                                                                                    >  2:836572/1‑93 (MQ=255)

      tgtgTATCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTggg                                                                                    >  2:8343107/1‑93 (MQ=255)

       gtgtATCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGc                                                                                   <  2:6628610/93‑1 (MQ=255)

          tATCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAgg                                                                                <  1:6373799/93‑1 (MQ=255)

           aTCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGt                                                                               >  2:4822179/1‑93 (MQ=255)

           aTCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGt                                                                               >  2:7862235/1‑93 (MQ=255)

           aTCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGt                                                                               <  2:8683026/93‑1 (MQ=255)

            tCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGtt                                                                              <  1:7904373/93‑1 (MQ=255)

            tCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGtt                                                                              >  2:2986105/1‑93 (MQ=255)

             cAGCCTGGCGTACATGCTCTACTTCGGCGCTATTGTTCTGCTGCTGCAGTGGCTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGttt                                                                             <  2:4445708/93‑1 (MQ=255)

                 cTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGcc                                                                         >  1:8549750/1‑93 (MQ=255)

                 cTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGcc                                                                         <  1:4256295/93‑1 (MQ=255)

                 cTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGcc                                                                         <  1:9117256/93‑1 (MQ=255)

                  tCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCct                                                                        <  2:6112651/93‑1 (MQ=255)

                  tCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCct                                                                        >  1:2876672/1‑93 (MQ=255)

                   cgcgTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCctc                                                                       >  2:9316742/1‑93 (MQ=255)

                    gcgTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCctct                                                                      >  2:9282594/1‑93 (MQ=255)

                    gcctataTGCTCTACTTCGGCACTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAGCCTGGGCAGGTTTGGCctct                                                                      <  2:8041490/90‑1 (MQ=255)

                     cgTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTg                                                                     <  1:6061354/93‑1 (MQ=255)

                      gTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTgc                                                                    <  1:8006682/93‑1 (MQ=255)

                      gTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTgc                                                                    <  1:5662674/93‑1 (MQ=255)

                       tataTGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTgcg                                                                   >  2:2828647/1‑93 (MQ=255)

                        atatGCTCTACTTCGGCGCTATGGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTgcgc                                                                  <  2:6872130/93‑1 (MQ=255)

                        atatGCTCTACTTCGGCGCTAGTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTgcgc                                                                  <  2:2207736/93‑1 (MQ=255)

                         tatGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGcc                                                                 >  1:710800/1‑93 (MQ=255)

                         tatGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGcc                                                                 >  2:6902237/1‑93 (MQ=255)

                          atGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCGGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCg                                                                >  1:1424289/1‑93 (MQ=255)

                             gtctGCTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTa                                                             <  1:1640777/92‑1 (MQ=255)

                             ctctACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGtg                                                             >  2:3057762/1‑92 (MQ=255)

                             ctctACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGtg                                                             >  2:7455756/1‑92 (MQ=255)

                             ctctACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTa                                                             <  2:5242961/93‑1 (MQ=255)

                               ctACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAgg                                                           <  1:9653029/93‑1 (MQ=255)

                                     cGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGattat                                                     >  1:9173048/1‑93 (MQ=255)

                                          cTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGCTTGGCCTCTGCGCCGGTAGGGATTATTCCgg                                                >  2:1073689/1‑93 (MQ=255)

                                          cTATTGTTCTGCCGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCgg                                                >  2:7086029/1‑93 (MQ=255)

                                           tATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGTGCCGGTAGGGATTATTCCGGt                                               <  1:7597002/93‑1 (MQ=255)

                                             ttGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCGACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGa                                             <  1:8433700/93‑1 (MQ=255)

                                               gTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATc                                           >  2:8936902/1‑93 (MQ=255)

                                                ttCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATcc                                          >  1:372990/1‑93 (MQ=255)

                                                ttCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATcc                                          >  1:4085262/1‑93 (MQ=255)

                                                ttCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATcc                                          >  2:7685558/1‑93 (MQ=255)

                                                  ctgctgCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCGACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTg                                        >  1:6171602/1‑93 (MQ=255)

                                                  ctgctgCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTg                                        <  1:5848868/93‑1 (MQ=255)

                                                   tgctgcCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGGCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGt                                       <  1:4821746/93‑1 (MQ=255)

                                                   tgctgcCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGt                                       <  1:9398785/93‑1 (MQ=255)

                                                   tgctgcCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGt                                       <  2:4780253/93‑1 (MQ=255)

                                                    gctgcCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTc                                      >  2:9611922/1‑93 (MQ=255)

                                                      tgcCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGc                                    >  2:4375331/1‑93 (MQ=255)

                                                      tgcCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGc                                    >  1:4110466/1‑93 (MQ=255)

                                                      tgcCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGc                                    >  2:4539549/1‑93 (MQ=255)

                                                      tgcCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCCCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGc                                    <  1:3076773/93‑1 (MQ=255)

                                                       gcCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGcc                                   <  2:567825/93‑1 (MQ=255)

                                                       gcCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGcc                                   >  1:1990194/1‑93 (MQ=255)

                                                        ccGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCg                                  >  1:5537407/1‑93 (MQ=255)

                                                          gCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCggt                                >  2:7146932/1‑91 (MQ=255)

                                                          gCAGTTGTTGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGAt                                <  2:6347712/93‑1 (MQ=255)

                                                             gttgttGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTAt                             >  1:35704/1‑93 (MQ=255)

                                                             gttgttGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTAt                             >  1:6963789/1‑93 (MQ=255)

                                                             gttgttGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTAt                             >  2:6683401/1‑93 (MQ=255)

                                                             gttgttGCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTAt                             >  2:3885641/1‑93 (MQ=255)

                                                              ttgttgCAGGAGGTCTACGGTTACACGGCGACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATc                            >  1:8782436/1‑93 (MQ=255)

                                                              ttgttgCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATc                            <  1:89895/93‑1 (MQ=255)

                                                                gttgCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCgg                          >  1:4681992/1‑93 (MQ=255)

                                                                gttgCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCgg                          <  2:1992183/93‑1 (MQ=255)

                                                                gttgCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCgg                          >  1:3152548/1‑93 (MQ=255)

                                                                 ttgCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCGGc                         <  2:3392941/93‑1 (MQ=255)

                                                                 gtgCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCGGc                         <  1:5283875/92‑1 (MQ=255)

                                                                  tgCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCGGcc                        <  2:6458690/93‑1 (MQ=255)

                                                                   gCAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCGGCCg                       >  1:8710668/1‑93 (MQ=255)

                                                                    cAGGAGGTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCGGCCGc                      <  2:7051604/93‑1 (MQ=255)

                                                                      ggaggTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCGGCCGCtt                    >  1:2493801/1‑93 (MQ=255)

                                                                      ggaggTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCGGCCGCtt                    >  1:1010808/1‑93 (MQ=255)

                                                                       gaggTCTACGGTTACACGGCGACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCGGCCGCTTc                   >  2:9782214/1‑93 (MQ=255)

                                                                       gaggTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCGGCCGCTTc                   <  1:6395115/93‑1 (MQ=255)

                                                                       gaggTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCGGCCGCTTc                   <  2:3510029/93‑1 (MQ=255)

                                                                         ggTCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCGGCCGCTTcgc                 <  2:7195592/93‑1 (MQ=255)

                                                                           tCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCGGCCGCTTCgcgc               <  2:6203071/93‑1 (MQ=255)

                                                                           tCTACGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCGGCCGCTTCgcgc               >  2:1017441/1‑93 (MQ=255)

                                                                              aCGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCGGCCGCTTCGCGCATa            >  1:7298682/1‑93 (MQ=255)

                                                                              aCGGTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCGGCCGCTTCGCGCATa            <  1:3424020/93‑1 (MQ=255)

                                                                                ggatACACGGCGACCTGGGAAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCGGCCGCTTCGCGCATaaa          <  2:9206920/90‑1 (MQ=255)

                                                                                ggTTACACGGCGACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCGGCCGCTTCGCGCATaaa          <  1:3491347/93‑1 (MQ=255)

                                                                                ggTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCGGCCGCTTCGCGCATaaa          <  2:4194351/93‑1 (MQ=255)

                                                                                ggTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCGGCCGCTTCGCGCATaaa          <  2:2532793/93‑1 (MQ=255)

                                                                                 gTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCGGCCGCTTCGCGCATAAAc         >  1:2258484/1‑93 (MQ=255)

                                                                                 gTTACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCGGCCGCTTCGCGCATAAAc         <  1:2035984/93‑1 (MQ=255)

                                                                                  ttACACGGCGACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCGGCCGCTTCGCGCATAAACt        <  1:4591456/93‑1 (MQ=255)

                                                                                  ttACACGGCGACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCGGCCGCTTCGCGCATAAACt        <  2:6312604/93‑1 (MQ=255)

                                                                                   tACACGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCGGCCGCTTCGCGCATAAACTg       <  2:5791775/93‑1 (MQ=255)

                                                                                    acacGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCGGCCGCTTCGCGCATAAACTgg      <  1:1750596/93‑1 (MQ=255)

                                                                                      acGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCGGCCGCTTCGCGCATAAACTGGat    <  2:1811144/93‑1 (MQ=255)

                                                                                       cGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCGGCCGCTTCGCGCATAAACTGGata   <  1:5249082/93‑1 (MQ=255)

                                                                                       cGGCAACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCGGCCGCTTCGCGCATAAACTGGata   <  1:2804899/93‑1 (MQ=255)

GCTGCTTGTGTATCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCGACCTGGGCAGGTTTGGCCTCTGCGCCGGTAGGGATTATTCCGGTGATCCTGTCGCCGATTATCGGCCGCTTCGCGCATAAACTGGATAT  >  NC_000913/2811454‑2811634

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 20 ≤` `ATCG/ATCG` `< 28 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`