BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	934,731	A→G	Q762R (CAG→CGG)	ftsK →	cell division protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	934,731	0	A	G	99.3%	473.4 / NA	143	Q762R (CAG→CGG)	ftsK	cell division protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (101/41); total (102/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 187 total aligned reads displayed.

TTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATC  >  NC_000913/934641‑934820

ttACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑gccgcggc                                                                                         <  1:6289967/93‑1 (MQ=255)

ttACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑gccgcggc                                                                                         <  1:1688493/93‑1 (MQ=255)

 tACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑gccgcggcc                                                                                        >  1:6346737/1‑93 (MQ=255)

  ccGCCAATGGTTGACCCTGTCAAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑gccgcggccg                                                                                       <  2:4264298/92‑1 (MQ=18)

  aCGCCAATTGTTGACCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑gccgcggccg                                                                                       <  1:417429/93‑1 (MQ=255)

    gCCAATTGTTTAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑gccgcgtccgca                                                                                     >  1:621743/1‑93 (MQ=255)

    gCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑gccgcggccgca                                                                                     <  2:237967/93‑1 (MQ=255)

    cccAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑gccgcggccgca                                                                                     <  2:3232335/92‑1 (MQ=255)

     cgaaTTGTTGAACCTGTACAGCAGCCGCAACAACCGGTGGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑Gccgcggccgcag                                                                                    <  2:4589937/91‑1 (MQ=255)

     ccAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑Gccgcggccgcag                                                                                    >  2:807161/1‑93 (MQ=255)

     ccAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑Gccgcggccgcag                                                                                    >  2:6263934/1‑93 (MQ=255)

     caatttgttgAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑Gccgcggccgcag                                                                                    <  1:5497648/89‑1 (MQ=255)

      cAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGt                                                                                   >  1:4286335/1‑93 (MQ=255)

       atttgttgAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTa                                                                                  <  1:658944/91‑1 (MQ=255)

        aTTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTAt                                                                                 >  2:4219619/1‑93 (MQ=255)

        aTTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTAt                                                                                 <  2:3855127/93‑1 (MQ=255)

         ttgttgACCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATc                                                                                <  1:771592/93‑1 (MQ=255)

          tgtggACCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATca                                                                               <  1:6002238/93‑1 (MQ=18)

             tgAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATcagca                                                                            >  1:1489028/1‑93 (MQ=255)

                acttGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGcc                                                                         <  1:4899864/90‑1 (MQ=255)

                aCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGcc                                                                         >  2:272555/1‑93 (MQ=255)

                 ccTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCCGCCGCca                                                                        >  1:1959225/1‑93 (MQ=255)

                 ccTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCca                                                                        <  1:3464418/93‑1 (MQ=255)

                 ccTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCca                                                                        >  1:1350591/1‑93 (MQ=255)

                   tGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCcaca                                                                      >  2:225563/1‑93 (MQ=255)

                   tGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCcaca                                                                      >  1:4612060/1‑93 (MQ=255)

                    gTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCacaa                                                                     >  2:2774989/1‑93 (MQ=255)

                    gTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCAGCCGCAGTATCAGCAGCCacaa                                                                     >  1:200493/1‑93 (MQ=255)

                     tACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCacaac                                                                    <  2:3180544/93‑1 (MQ=255)

                      acacCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCacaaca                                                                   <  1:3615323/93‑1 (MQ=255)

                      aCAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCacaaca                                                                   >  2:6064268/1‑93 (MQ=255)

                        agcagcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGc                                                                 <  1:4072035/93‑1 (MQ=255)

                          cagcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCg                                                               >  2:2924352/1‑93 (MQ=255)

                           agcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCgg                                                              >  2:4554574/1‑93 (MQ=255)

                           agcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCgg                                                              >  1:1726899/1‑93 (MQ=255)

                            gcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGt                                                             >  1:2299591/1‑93 (MQ=255)

                            gcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCACAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGt                                                             <  2:5718283/93‑1 (MQ=255)

                             ccGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGtt                                                            <  1:1266634/93‑1 (MQ=255)

                              cGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTg                                                           >  1:533767/1‑93 (MQ=255)

                              cGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTg                                                           >  2:846362/1‑93 (MQ=255)

                              cGCAACAACCGGTTGCACCGCAGCAGCAATACCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTg                                                           <  2:3484964/93‑1 (MQ=255)

                              cGCAACAACCGGTTGCACCGCAGCAGCAAGATAAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTg                                                           <  2:2652989/93‑1 (MQ=17)

                               gCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTgc                                                          >  1:2743988/1‑93 (MQ=255)

                                caacaaCCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTgcg                                                         >  1:2976030/1‑93 (MQ=255)

                                 aacaccCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTgcgc                                                        <  2:5541030/93‑1 (MQ=255)

                                  acaacCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGcgcc                                                       >  1:2255347/1‑93 (MQ=255)

                                  acaacCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGcgcc                                                       >  1:5004503/1‑93 (MQ=255)

                                    aacCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGcgccgc                                                     >  1:577950/1‑93 (MQ=255)

                                    aacCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGcgccgc                                                     >  1:1463051/1‑93 (MQ=255)

                                     acCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCa                                                    >  1:912988/1‑93 (MQ=255)

                                        ggTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGcc                                                 >  2:1428571/1‑93 (MQ=255)

                                        ggTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGcc                                                 >  2:4863430/1‑93 (MQ=255)

                                        ggTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGcc                                                 >  2:4444649/1‑93 (MQ=255)

                                        ggTTCCACCCCAACACAAATATCAGCAGCCGCAACAACCAGTTCG‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGcc                                                 <  2:138001/93‑1 (MQ=18)

                                         gTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAACCGGTTCCGCCGCAGCca                                                >  1:3844063/1‑93 (MQ=255)

                                            gCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACaa                                             <  2:1489028/93‑1 (MQ=255)

                                              cccGCAGCAGAAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAata                                           <  1:2382195/92‑1 (MQ=255)

                                               ccGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAatat                                          >  2:2122416/1‑93 (MQ=255)

                                                 gcggcagcaATATCAGCAGCCGCAACAACCAGTCCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATca                                        <  2:2615368/90‑1 (MQ=255)

                                                 gcagaagaaATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATca                                        <  2:471722/93‑1 (MQ=255)

                                                  cagcagcaATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATcag                                       >  2:685867/1‑93 (MQ=255)

                                                   agcagcaAAATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATcagc                                      <  2:42895/93‑1 (MQ=255)

                                                    gcagcaATATCAGCAGCCGCAACACCCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATTTcagca                                     <  1:5362050/93‑1 (MQ=255)

                                                    gcagcaATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATcagca                                     >  1:3738842/1‑93 (MQ=255)

                                                    cgagcaAAACCACCACCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATcagca                                     <  2:2239847/91‑1 (MQ=21)

                                                     cagcaATATCAGCGGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCAACAATATcagcag                                    <  1:6475051/93‑1 (MQ=255)

                                                      agcaATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCagcagc                                   >  2:1855801/1‑93 (MQ=255)

                                                       gcaATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGcc                                  >  2:945139/1‑93 (MQ=255)

                                                       gcaATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGcc                                  >  1:1177533/1‑93 (MQ=255)

                                                       gcaATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTCCGCCGCAGCCACAATATCAGCAGcc                                  >  2:3880160/1‑93 (MQ=255)

                                                        aaaTATCGGCAGCCCAAAAACCCATTCCC‑GCCGCGTCCGCAGTACCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCg                                 <  2:3859137/92‑1 (MQ=21)

                                                          atatCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGca                               >  2:5829102/1‑93 (MQ=255)

                                                           tatCAGCAGCCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGcaa                              <  2:1621273/93‑1 (MQ=255)

                                                              cggcagcCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGcaacaa                           <  2:5857021/91‑1 (MQ=255)

                                                              cagcagCCGAAACACCAGGTCCC‑GCCGCGCCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGcaacaa                           <  1:1628711/93‑1 (MQ=18)

                                                               agcagcCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCaacaac                          >  2:2005153/1‑93 (MQ=255)

                                                               agcagcCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCaaaaac                          >  1:5136197/1‑93 (MQ=255)

                                                               agcagcCGCAACAACCAGCTCC‑GCCGCGGCCGCAGTATCAGCAGCCACGACAGCCGGCTGCGCCGCAGCCACAATATCAGCAGCCGCaacaac                          >  2:269220/1‑93 (MQ=21)

                                                                gcagcCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAAcc                         >  1:1536825/1‑93 (MQ=255)

                                                                gcagcCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAAcc                         >  2:5668088/1‑93 (MQ=255)

                                                                gcagcCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGACACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAAcc                         >  1:4146367/1‑93 (MQ=255)

                                                                 cagcCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCg                        >  2:2795603/1‑93 (MQ=255)

                                                                 caaccGAAACACCCATTTCC‑GCCCCGCCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCg                        <  1:4541975/90‑1 (MQ=18)

                                                                   gcCGCAACAACCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGt                      >  1:3295317/1‑93 (MQ=255)

                                                                    ccGCAAAAACC‑GTTCCGGCGGCGCCCGCAGTACCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGtt                     <  2:1870371/93‑1 (MQ=18)

                                                                      ccaacaaCCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTgc                   <  2:2525487/92‑1 (MQ=255)

                                                                       caacaaCCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTgcg                  >  1:1114276/1‑93 (MQ=255)

                                                                         acaacCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGcc                >  1:286453/1‑93 (MQ=255)

                                                                         acaacCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGcc                >  2:3593603/1‑93 (MQ=255)

                                                                          caacCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCca               >  2:1243951/1‑93 (MQ=255)

                                                                           aacCAGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCcac              >  2:1000708/1‑93 (MQ=255)

                                                                            acCAGTCCC‑GCGGCGGCCTCAGTATGACCAGCCACAACAGCCGGTTGCTCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCcaca             <  1:2830665/93‑1 (MQ=18)

                                                                               aGTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCAcagca          >  1:2691665/1‑93 (MQ=255)

                                                                                gTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCAcagcag         <  1:5448855/93‑1 (MQ=255)

                                                                                gTTCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCAcagcag         >  1:1651799/1‑93 (MQ=255)

                                                                                  tCC‑GCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAgcagca       <  1:3424009/93‑1 (MQ=255)

                                                                                   cc‑gccgCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACAGGTTGCGCCACAGCAGCaa      <  1:950380/93‑1 (MQ=255)

                                                                                    c‑gccgcGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAat     >  2:3971664/1‑93 (MQ=255)

                                                                                      gccgcggccgcaGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAata    >  1:3358587/1‑93 (MQ=255)

                                                                                      gccgcggccgaaGTATCACCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAata    <  2:5607898/93‑1 (MQ=255)

TTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCC‑GCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATC  >  NC_000913/934641‑934820

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