BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,458,752	T→C	R86R (CGT→CGC)	b1396 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,458,752	0	T	C	98.5%	387.4 / ‑2.4	130	R86R (CGT→CGC)	b1396	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/2); new base C (67/61); total (67/63)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.49e-01

Only 100 of 139 total aligned reads displayed.

ATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGTCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGAAATTGTTAA  >  NC_000913/1458662‑1458841

aTTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGccc                                                                                         >  2:1301386/1‑93 (MQ=255)

  tttCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAg                                                                                       <  1:2290202/93‑1 (MQ=255)

      ccTGCCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCtt                                                                                   <  1:2984125/92‑1 (MQ=255)

       cTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCttt                                                                                  >  2:4537491/1‑93 (MQ=255)

       cTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCttt                                                                                  <  2:2646745/93‑1 (MQ=255)

         ggCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGc                                                                                <  2:2625070/93‑1 (MQ=255)

            tGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCgg                                                                             >  1:303473/1‑93 (MQ=255)

             gATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGga                                                                            <  1:2827218/93‑1 (MQ=255)

             gATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGga                                                                            >  1:852707/1‑93 (MQ=255)

              aTACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGgag                                                                           >  1:1769867/1‑93 (MQ=255)

               tACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGgaga                                                                          >  2:3760237/1‑93 (MQ=255)

                 cTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGaca                                                                        <  1:4318508/93‑1 (MQ=255)

                  tGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGacac                                                                       <  1:444226/93‑1 (MQ=255)

                  tGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGacac                                                                       >  1:1669270/1‑93 (MQ=255)

                   gCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACAcc                                                                      <  1:1155528/93‑1 (MQ=255)

                   gCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACAcc                                                                      >  2:3067389/1‑93 (MQ=255)

                    ccTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCt                                                                     >  1:637984/1‑93 (MQ=255)

                    ccTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCt                                                                     >  1:2778649/1‑93 (MQ=255)

                     cTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCtt                                                                    >  1:193693/1‑93 (MQ=255)

                      tttGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTa                                                                   >  2:3459686/1‑93 (MQ=255)

                         gCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACt                                                                <  1:1841267/93‑1 (MQ=255)

                         cccTCCGCCTGCAATAGTCAGGGGCTGGCACCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACt                                                                <  1:2653/92‑1 (MQ=255)

                          ccTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTg                                                               <  2:492983/93‑1 (MQ=255)

                            tACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCt                                                             >  2:3659145/1‑93 (MQ=255)

                             aCGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTa                                                            <  2:3138954/93‑1 (MQ=255)

                              cGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTAc                                                           <  1:502220/93‑1 (MQ=255)

                               gCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACt                                                          >  2:2148662/1‑93 (MQ=255)

                               gCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACt                                                          >  1:4253163/1‑93 (MQ=255)

                                ccTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTg                                                         >  2:1575668/1‑93 (MQ=255)

                                 cTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTgc                                                        >  2:1928780/1‑93 (MQ=255)

                                 cTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTgc                                                        >  1:309978/1‑93 (MQ=255)

                                   gCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTgcgc                                                      >  1:97960/1‑93 (MQ=255)

                                     aaTAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAg                                                    <  1:4239995/93‑1 (MQ=255)

                                      aTAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAgg                                                   >  2:698694/1‑93 (MQ=255)

                                        agtcAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGGGCAGGTa                                                 >  1:3044904/4‑93 (MQ=255)

                                          tcAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACg                                               >  2:220260/2‑93 (MQ=255)

                                            aGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGtc                                             <  2:1948705/93‑1 (MQ=255)

                                            aGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGtc                                             <  1:271620/93‑1 (MQ=255)

                                             ggggCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGtca                                            >  1:531704/1‑93 (MQ=255)

                                               ggCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGTCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGtcatc                                          <  2:1544878/93‑1 (MQ=255)

                                               ggCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGtcatc                                          >  2:3339928/1‑93 (MQ=255)

                                               ggCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGtcatc                                          <  1:890085/93‑1 (MQ=255)

                                                gCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGtcatca                                         <  2:3628992/93‑1 (MQ=255)

                                                gCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGtcatca                                         <  2:1867878/93‑1 (MQ=255)

                                                 cTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAg                                        <  1:1560663/93‑1 (MQ=255)

                                                 cTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAg                                        <  1:887212/93‑1 (MQ=255)

                                                  tGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAgg                                       >  1:204277/1‑93 (MQ=255)

                                                   ggCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAggg                                      >  1:4018678/1‑93 (MQ=255)

                                                   ggCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAggg                                      <  2:4539447/93‑1 (MQ=255)

                                                      aGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGgcaa                                   >  2:1462533/1‑93 (MQ=255)

                                                         cGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGgcaagca                                >  2:2437051/1‑93 (MQ=255)

                                                          gTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGgcaagcaa                               <  2:2531142/93‑1 (MQ=255)

                                                           tCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGgcaagcaga                              >  2:1999171/1‑91 (MQ=255)

                                                           tCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCaaa                              >  2:111266/1‑93 (MQ=255)

                                                           tCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCaaa                              >  2:3028761/1‑93 (MQ=255)

                                                           tCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCaaa                              <  1:2651583/93‑1 (MQ=255)

                                                            cGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAAc                             <  2:4651888/93‑1 (MQ=255)

                                                            cGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAAc                             <  2:1282223/93‑1 (MQ=255)

                                                            cGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAAc                             <  2:3444729/93‑1 (MQ=255)

                                                             gCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAAcc                            >  2:101348/1‑93 (MQ=255)

                                                             ccTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAAcc                            <  1:2103956/92‑1 (MQ=255)

                                                              cTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCg                           >  1:1178528/1‑93 (MQ=255)

                                                              cTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCg                           >  2:3474408/1‑93 (MQ=255)

                                                               ttCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCgg                          >  2:640875/1‑93 (MQ=255)

                                                               ttCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCgg                          <  2:1887889/93‑1 (MQ=255)

                                                               ttCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCgg                          >  1:2487247/1‑93 (MQ=255)

                                                               ttCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCgg                          >  1:1540566/1‑93 (MQ=255)

                                                                tCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGgt                         <  1:3785201/93‑1 (MQ=255)

                                                                  tgcctgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGgtgt                       <  1:705119/93‑1 (MQ=255)

                                                                    cctgcACGATTGATTTTTTGCGCCCAGGCTTTGCCTGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGttt                     <  2:4026255/93‑1 (MQ=255)

                                                                    cctgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGttt                     >  1:2570227/1‑93 (MQ=255)

                                                                     ctgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTa                    >  1:1480951/1‑93 (MQ=255)

                                                                     ctgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTa                    <  2:577542/93‑1 (MQ=255)

                                                                     ctgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTa                    <  1:1365314/93‑1 (MQ=255)

                                                                      tgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTAc                   <  1:1681294/93‑1 (MQ=255)

                                                                       gcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACg                  >  2:4122795/1‑93 (MQ=255)

                                                                        cACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGa                 <  1:3285743/93‑1 (MQ=255)

                                                                         aCGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCCGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGAc                >  2:2645982/1‑93 (MQ=255)

                                                                         aCGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGAc                >  1:2421996/1‑93 (MQ=255)

                                                                         aCGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGAc                >  1:2574135/1‑93 (MQ=255)

                                                                          cGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACa               >  2:4049906/1‑93 (MQ=255)

                                                                          cGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACa               >  2:3303283/1‑93 (MQ=255)

                                                                           gattgattTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACAt              <  1:4671731/93‑1 (MQ=255)

                                                                           gattgattTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACAt              <  2:1195916/93‑1 (MQ=255)

                                                                            attgattTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATc             <  1:1742384/93‑1 (MQ=255)

                                                                            attgattTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATc             >  2:3834514/1‑93 (MQ=255)

                                                                               gattTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGaa          >  2:2816355/1‑93 (MQ=255)

                                                                                attTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGaaa         <  2:4635514/93‑1 (MQ=255)

                                                                                attTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGaaa         <  2:587835/93‑1 (MQ=255)

                                                                                 ttttttGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGAAAt        <  1:3788330/93‑1 (MQ=255)

                                                                                 ttttttGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGAAAt        >  1:1454474/1‑93 (MQ=255)

                                                                                 ttttttGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGAAAt        <  2:3312612/93‑1 (MQ=255)

                                                                                 ttttttGCGCCCAGGCTTTGCCGGAGACACCTTAACGGCTACGGCGGAGGTGCGTCATCAGGGCAAGCAAAGCGGTGTTTACGACATCGAAAt        <  2:3046932/93‑1 (MQ=255)

                                                                                  tttttGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGAAAtt       <  2:2696865/93‑1 (MQ=255)

                                                                                    tttGCGCCCGGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGCCAAGCAAACCGGTGTTTACGACATCGAAATTGt     <  2:1228033/93‑1 (MQ=255)

                                                                                     ttGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGAAATTGtt    <  2:3572201/93‑1 (MQ=255)

                                                                                     ttGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGAAATTGtt    >  1:4519996/1‑93 (MQ=255)

                                                                                      tGCGCCCAGGCTTTGCCGGAGACTCCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGAAATTGTTa   <  2:3103153/93‑1 (MQ=255)

                                                                                      tGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGAAATTGTTa   <  1:1011601/93‑1 (MQ=255)

                                                                                       gcgcCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGAAATTGTTaa  >  2:2869947/1‑93 (MQ=255)

ATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGTCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGAAATTGTTAA  >  NC_000913/1458662‑1458841

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