BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,949,344	0	G	A	62.4%	372.7 / inf	371	Q174Q (CAG→CAA)	yecD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (80/59); new base A (118/113); total (199/172)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.40e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 400 total aligned reads displayed.

ATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTA  >  NC_000913/1949252‑1949436

  cGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAgcaaca                                                                                            <  2:4260211/93‑1 (MQ=255)

   ggTGTTGAATCCATCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGGGATTTCCGAAGACGACTGTAGTGCCGCCAGCGCCGagcagccc                                                                                           >  2:4526401/1‑91 (MQ=255)

    gtgtTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAAcaca                                                                                          >  2:3427798/1‑93 (MQ=255)

     tgtTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACaa                                                                                         <  1:1902251/93‑1 (MQ=255)

      gtTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAAt                                                                                        <  1:3453157/93‑1 (MQ=255)

      gtTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAAt                                                                                        <  1:1476669/93‑1 (MQ=255)

      gtTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAAt                                                                                        <  1:3272941/93‑1 (MQ=255)

      gtTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAAt                                                                                        <  2:2119234/93‑1 (MQ=255)

        tGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATaa                                                                                      <  2:3483295/93‑1 (MQ=255)

         gAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAAc                                                                                     <  2:3008837/93‑1 (MQ=255)

          aaTCCACGGCCCGCAATGGCTGGGAACTTGGCTTTAAGCTGGTGATTGCCGAAGACGCCTGTGGTGCCGCCAGCGCCGAGCAACACAATAACa                                                                                    <  1:3843926/93‑1 (MQ=255)

          aaTCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACa                                                                                    >  2:2851318/1‑93 (MQ=255)

          aaTCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACa                                                                                    >  2:3295517/1‑93 (MQ=255)

            tCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAgc                                                                                  >  2:4171343/1‑92 (MQ=255)

            accACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGt                                                                                  <  1:1238307/92‑1 (MQ=255)

             ccACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAgca                                                                                 <  1:447188/93‑3 (MQ=255)

              cACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCAt                                                                                >  2:4099707/1‑93 (MQ=255)

                ccgcccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCAtta                                                                              <  1:845134/93‑3 (MQ=255)

                  gcccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTaat                                                                            >  1:650312/1‑92 (MQ=255)

                   cccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTaatc                                                                           <  2:2176244/93‑3 (MQ=255)

                      gcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCata                                                                        <  1:4262366/93‑1 (MQ=255)

                      gcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCata                                                                        <  1:3595403/93‑1 (MQ=255)

                       cAATGCCTGGGAACTTGGCTTTACTCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCatat                                                                       <  2:2584203/93‑1 (MQ=255)

                        aaTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATc                                                                      >  2:3185064/1‑93 (MQ=255)

                        aaTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATc                                                                      >  2:2201110/1‑93 (MQ=255)

                        aaTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATc                                                                      >  1:2665561/1‑93 (MQ=255)

                        aaTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATc                                                                      >  1:3500325/1‑93 (MQ=255)

                        aaTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATc                                                                      <  2:2689854/93‑1 (MQ=255)

                        aaTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATc                                                                      >  1:3897677/1‑93 (MQ=255)

                        aaTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATc                                                                      >  1:4022155/1‑93 (MQ=255)

                         aTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCt                                                                     >  1:4458242/1‑93 (MQ=255)

                         aTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCt                                                                     >  1:978618/1‑93 (MQ=255)

                         aTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCt                                                                     <  1:3698583/93‑1 (MQ=255)

                          tGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTa                                                                    >  1:3340445/1‑93 (MQ=255)

                          tGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTa                                                                    >  1:1422608/1‑93 (MQ=255)

                           gCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTAc                                                                   >  2:4737758/1‑93 (MQ=255)

                             cTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTAccc                                                                 >  2:39251/1‑93 (MQ=255)

                              tGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCcg                                                                >  1:1134948/1‑93 (MQ=255)

                              tGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCcg                                                                >  1:1709012/1‑93 (MQ=255)

                               gggAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCcgc                                                               >  1:277218/1‑93 (MQ=255)

                                ggAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCACCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCcgcc                                                              >  1:3153893/1‑92 (MQ=255)

                                 gAACTTGGCTTTAATCTGGTGACTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCgcgc                                                             <  2:2041755/93‑1 (MQ=255)

                                   aCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCAt                                                           <  2:3750/93‑1 (MQ=255)

                                   aCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCAt                                                           <  1:3765198/93‑1 (MQ=255)

                                   aCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCAt                                                           >  1:59985/1‑93 (MQ=255)

                                   aCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCAt                                                           >  1:2267307/1‑93 (MQ=255)

                                    cTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATc                                                          <  2:561507/93‑1 (MQ=255)

                                     ttGGCTTTAATCTGGTGATTGCCGATGACGCCTGTAGTGCCGTCAGCGCCGAGCAACATAATAACAGTATTAATCATATCTACTCGCGCATCg                                                         >  2:4220956/1‑93 (MQ=255)

                                      tGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGc                                                        >  2:4605605/1‑93 (MQ=255)

                                      tGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGc                                                        >  2:2819714/1‑93 (MQ=255)

                                       ggCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGcc                                                       <  2:3473944/93‑1 (MQ=255)

                                         cTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCg                                                     <  2:3814608/93‑1 (MQ=255)

                                          tttAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGAATTAATCATATCTACCCGCGCATCGCCCgg                                                    >  1:3562782/1‑92 (MQ=255)

                                            tAATCTGGTGATTGCCGAGGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCgtgt                                                  <  2:1918699/93‑1 (MQ=255)

                                             aaTCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGtgtg                                                 <  1:521864/93‑1 (MQ=255)

                                                 tGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTa                                             <  1:304736/93‑1 (MQ=255)

                                                  ggTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAg                                            >  1:1533347/1‑93 (MQ=255)

                                                    tGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCg                                          <  1:2129892/93‑1 (MQ=255)

                                                     gATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGt                                         >  1:2504703/1‑93 (MQ=255)

                                                     gATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGt                                         <  1:161679/93‑1 (MQ=255)

                                                      aTTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTg                                        >  2:3151594/1‑93 (MQ=255)

                                                      aTTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTg                                        <  1:2598936/93‑1 (MQ=255)

                                                      aTTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTg                                        >  1:2202450/1‑93 (MQ=255)

                                                       ttGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTgg                                       >  2:4486667/1‑93 (MQ=255)

                                                         gCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGaa                                     <  1:2030198/93‑1 (MQ=255)

                                                           cGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAaga                                   <  2:2479921/93‑1 (MQ=255)

                                                            gAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAagag                                  >  2:89670/1‑93 (MQ=255)

                                                            gAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAagag                                  >  2:1028913/1‑93 (MQ=255)

                                                             aaGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAgaga                                 <  2:2713525/93‑1 (MQ=255)

                                                             aaGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAgaga                                 >  2:1272975/1‑93 (MQ=255)

                                                               gACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATc                               <  2:2071238/93‑1 (MQ=255)

                                                                 cGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCt                             >  1:777304/1‑93 (MQ=255)

                                                                  gCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTCCCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCtt                            <  2:4715764/93‑2 (MQ=255)

                                                                    cTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATccgtga                          >  1:909796/1‑89 (MQ=255)

                                                                    cTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCttca                          >  2:2344656/1‑90 (MQ=255)

                                                                    cTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCttca                          >  2:3962901/1‑90 (MQ=255)

                                                                     tGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAAc                         >  2:3661089/1‑93 (MQ=255)

                                                                       gaGTGCCGCCAGCGTAGAGCAGCACAATAACAGTATCAACCATATCTCCCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAAcgc                       <  2:72727/92‑1 (MQ=255)

                                                                         gTGGCGCCAGCGCCGTGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGAGTAGAGTGGTAGAGATCCTTCACGCGt                     >  2:1538386/1‑93 (MQ=255)

                                                                           gccgccAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTa                   <  1:2843486/93‑1 (MQ=255)

                                                                            ccgccAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTAt                  <  2:1867003/93‑1 (MQ=255)

                                                                               ccAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGAt               >  2:3036642/1‑93 (MQ=255)

                                                                               ccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGAt               <  1:3981634/93‑1 (MQ=255)

                                                                                cAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCTCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGAtt              >  1:4373560/1‑93 (MQ=255)

                                                                                cAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGAtt              >  1:2340561/1‑93 (MQ=255)

                                                                                cAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGAtt              <  1:1768072/93‑1 (MQ=255)

                                                                                  gcgcCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTa            >  1:1082692/1‑93 (MQ=255)

                                                                                   cgcCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTAc           >  1:2008401/1‑93 (MQ=255)

                                                                                   cgcCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTAc           <  2:80694/93‑1 (MQ=255)

                                                                                    gcCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACa          >  2:1456629/1‑93 (MQ=255)

                                                                                    gcCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACa          >  1:195139/1‑93 (MQ=255)

                                                                                    gcCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACa          <  1:4164913/93‑1 (MQ=255)

                                                                                    gcCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACa          <  1:4719406/93‑1 (MQ=255)

                                                                                     ccGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACAt         <  1:1871282/93‑1 (MQ=255)

                                                                                      cGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATc        <  2:441023/93‑1 (MQ=255)

                                                                                        agcagcACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCgg      <  1:460286/93‑1 (MQ=255)

                                                                                        agcagcACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCgg      <  2:178613/93‑1 (MQ=255)

                                                                                        agcaacACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCgg      <  1:2411963/93‑1 (MQ=255)

                                                                                         gcaacaCAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGt     <  2:4359667/93‑1 (MQ=255)

                                                                                           agcaCAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCt   >  2:2923095/1‑93 (MQ=255)

ATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTA  >  NC_000913/1949252‑1949436

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