BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,974,350	G→A	R555R (CGC→CGT)	aas ←	acyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,974,350	0	G	A	99.6%	892.0 / NA	224	R555R (CGC→CGT)	aas	acyltransferase
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (108/115); total (108/116)
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 224 total aligned reads displayed.

ATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCTACCGTACCGGGTTTCGCCGCCATCGGTAC  >  NC_000913/2974252‑2974445

aTTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACg                                                                                               >  1:3849769/1‑101 (MQ=255)

    cgccgcGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCg                                                                                           >  1:3849773/1‑101 (MQ=255)

    cgccgcGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCg                                                                                           <  2:3849430/101‑1 (MQ=255)

    cgccgcGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCg                                                                                           >  2:3849772/1‑101 (MQ=255)

       cgcgAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGccc                                                                                        <  1:3849540/101‑1 (MQ=255)

        gcgAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCct                                                                                       <  2:3849246/101‑1 (MQ=255)

         cgAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCctc                                                                                      >  1:3849778/1‑101 (MQ=255)

           aaCATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCtt                                                                                    >  2:3849779/1‑101 (MQ=255)

            aCATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTc                                                                                   <  2:3849438/101‑1 (MQ=255)

            aCATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTc                                                                                   <  1:3849321/101‑1 (MQ=255)

            aCATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTc                                                                                   <  1:3849460/101‑1 (MQ=255)

              aTTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTACGCCCTCTTCGa                                                                                 >  1:3849782/1‑101 (MQ=255)

               ttCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGAt                                                                                <  1:3849509/101‑1 (MQ=255)

                  tcGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATAcc                                                                             >  1:3849784/1‑101 (MQ=255)

                    ggcggTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAg                                                                           <  1:3849425/101‑1 (MQ=255)

                    ggcggTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAg                                                                           >  1:3849786/1‑101 (MQ=255)

                      cggTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCg                                                                                           >  2:3849787/1‑83 (MQ=255)

                        gTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGAc                                                                       <  1:3849255/101‑1 (MQ=255)

                        gTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGAc                                                                       <  1:3849489/101‑1 (MQ=255)

                         tGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGAcc                                                                      <  1:3849478/101‑1 (MQ=255)

                          gggCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCg                                                                     <  2:3849447/101‑1 (MQ=255)

                          gggCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCg                                                                     >  2:3849791/1‑101 (MQ=255)

                            gCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGac                                                                   <  2:3849401/101‑1 (MQ=255)

                             cACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGaca                                                                  >  1:3849792/1‑101 (MQ=255)

                              aCTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGacaa                                                                 >  2:3849794/1‑101 (MQ=255)

                               cTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGacaac                                                                <  1:3849348/101‑1 (MQ=255)

                               cTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGacaac                                                                <  1:3849376/101‑1 (MQ=255)

                                ttCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGacaaca                                                               <  2:3849475/101‑1 (MQ=255)

                                ttCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGacaaca                                                               <  2:3849530/101‑1 (MQ=255)

                                ttCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGacaaca                                                               >  1:3849798/1‑101 (MQ=255)

                                ttCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGacaaca                                                               >  1:3849800/1‑101 (MQ=255)

                                 tCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAg                                                              <  1:3849549/101‑1 (MQ=255)

                                 tCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAg                                                              <  1:3849471/101‑1 (MQ=255)

                                  ccAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAgg                                                             >  2:3849801/1‑101 (MQ=255)

                                   cAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGgc                                                            <  2:3849461/101‑1 (MQ=255)

                                   cAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGgc                                                            <  1:3849357/101‑1 (MQ=255)

                                      tACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGgcgcg                                                         >  1:3849803/1‑101 (MQ=255)

                                       acacCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGgcgcgc                                                        >  1:3849806/1‑101 (MQ=255)

                                        cacCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCa                                                       >  2:3849809/1‑101 (MQ=255)

                                        cacCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCa                                                       >  1:3849808/1‑101 (MQ=255)

                                           cTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATcc                                                    >  1:3849813/1‑101 (MQ=255)

                                           cTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATcc                                                    >  2:3849812/1‑101 (MQ=255)

                                            tGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGcc                                                                                         >  2:3849816/1‑63 (MQ=255)

                                            tGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCa                                                   <  2:3849413/101‑1 (MQ=255)

                                             ggCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCAt                                                  <  2:3849610/101‑1 (MQ=255)

                                              gCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCAtt                                                 >  1:3849818/1‑101 (MQ=255)

                                               cTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTc                                                <  2:3849513/101‑1 (MQ=255)

                                               cTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTc                                                >  2:3849821/1‑101 (MQ=255)

                                                ttCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTcc                                               >  2:3849823/1‑101 (MQ=255)

                                                  ctcCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTg                                             >  1:3849825/1‑101 (MQ=255)

                                                   tcCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTgg                                            <  2:3849355/101‑1 (MQ=255)

                                                   tcCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTgg                                            >  2:3849827/1‑101 (MQ=255)

                                                    ccACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCAt                                                  >  2:3849831/1‑94 (MQ=255)

                                                    ccACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGt                                           >  1:3849829/1‑101 (MQ=255)

                                                     cACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTa                                          <  1:3849445/101‑1 (MQ=255)

                                                     cACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTa                                          <  2:3849522/101‑1 (MQ=255)

                                                     cACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTa                                          <  1:3849420/101‑1 (MQ=255)

                                                     cACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTa                                          >  1:3849834/1‑101 (MQ=255)

                                                     cACCCGCAGATAGCCGTTCATTATGTTCGGCCCGTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTa                                          <  1:3849523/101‑1 (MQ=255)

                                                        ccGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGaa                                       <  2:3849275/101‑1 (MQ=255)

                                                         cGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAAt                                      >  1:3849837/1‑101 (MQ=255)

                                                          gCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATa                                     <  1:3849442/101‑1 (MQ=255)

                                                           cAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATAc                                    >  2:3849838/1‑101 (MQ=255)

                                                            aGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACg                                   >  2:3849840/1‑101 (MQ=255)

                                                            aGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACg                                   <  1:3849623/101‑1 (MQ=255)

                                                              aTAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGcc                                 <  1:3849370/101‑1 (MQ=255)

                                                               tAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGccc                                >  1:3849843/1‑101 (MQ=255)

                                                               tAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGccc                                <  2:3849584/101‑1 (MQ=255)

                                                                 gCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCta                              <  2:3849488/101‑1 (MQ=255)

                                                                 gCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCta                              <  1:3849381/101‑1 (MQ=255)

                                                                   cGTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCtacc                            >  2:3849845/1‑101 (MQ=255)

                                                                    gTTCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCtaccg                           <  2:3849363/101‑1 (MQ=255)

                                                                     ttCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCtaccgt                          >  2:3849848/1‑101 (MQ=255)

                                                                     ttCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCtaccgt                          <  1:3849573/101‑1 (MQ=255)

                                                                      tCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCtaccgta                         <  1:3849390/101‑1 (MQ=255)

                                                                      tCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCtaccgta                         <  2:3849493/101‑1 (MQ=255)

                                                                      tCATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCtaccgta                         <  2:3849562/101‑1 (MQ=255)

                                                                       cATTATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCtaccgtac                        >  1:3849849/1‑101 (MQ=255)

                                                                         ttATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCtaccgtaccg                      >  2:3849853/1‑101 (MQ=255)

                                                                          tATGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCTACCGTACCgg                     >  1:3849854/1‑101 (MQ=255)

                                                                            tGTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCTACCGTACCGGGt                   >  2:3849856/1‑101 (MQ=255)

                                                                             gTTCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCTACCGTACCGGGtt                  <  1:3849483/101‑1 (MQ=255)

                                                                              ttCGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCTACCGTACCGGGttt                 <  1:3849309/101‑1 (MQ=255)

                                                                                cGGCCCTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCTACCGTACCGGGTTTcg               >  1:3849858/1‑101 (MQ=255)

                                                                                   cccTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCTACCGTACCGGGTTTcgccg            >  1:3849863/1‑101 (MQ=255)

                                                                                   cccTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCTACCGTACCGGGTTTcgccg            >  2:3849862/1‑101 (MQ=255)

                                                                                   cccTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCTACCGTACCGGGTTTcgccg            <  2:3849374/101‑1 (MQ=255)

                                                                                   cccTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCTACCGTACCGGGTTTcgccg            <  1:3849599/101‑1 (MQ=255)

                                                                                    ccTTTCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCTACCGTACCGGGTTTcgccgc           <  2:3849620/101‑1 (MQ=255)

                                                                                      tttCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCTACCGTACCGGGTTTCGCCGCCa         >  1:3849866/1‑101 (MQ=255)

                                                                                        tCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCTACCGTACCGGGTTTCGCCGCCATc       <  2:3849322/101‑1 (MQ=255)

                                                                                        tCAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCTACCGTACCGGGTTTCGCCGCCATc       >  2:3849867/1‑101 (MQ=255)

                                                                                         cAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCTACCGTACCGGGTTTCTCCGCCATCg      <  2:3849545/101‑1 (MQ=255)

                                                                                         cAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCTACCGTACCGGGTTTCGCCGCCATCg      <  2:3849278/101‑1 (MQ=255)

                                                                                         cAGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCTACCGTACCGGGTTTCGCCGCCATCg      >  1:3849869/1‑101 (MQ=255)

                                                                                          aGTTGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCTACCGTACCGGGTTTCGCCGCCATCgg     <  1:3849496/101‑1 (MQ=255)

                                                                                            ttGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCTACCGTACCGGGTTTCGCCGCCATCGGTa   >  1:3849872/1‑101 (MQ=255)

                                                                                            ttGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCTACCGTACCGGGTTTCGCCGCCATCGGTa   <  1:3849410/101‑1 (MQ=255)

                                                                                             tGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCTACCGTACCGGGTTTCGCCGCCATCGGTAc  >  1:3849875/1‑101 (MQ=255)

                                                                                             tGCAGACGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCTACCGTACCGGGTTTCGCCGCCATCGGTAc  <  1:3849556/101‑1 (MQ=255)

ATTTCGCCGCGAACATTCTCGGCGGTGGGCACTTCCAGTACACCTGGCTTCTCCACCCGCAGATAGCCGTTCATTATGTTCGGCCCTTTCAGTTGCAGGCGTCCGCCCTCTTCGATACCAGGGACCGACAACAGGCGCGCATCCATTCCTGGTAGAATACGCCCTACCGTACCGGGTTTCGCCGCCATCGGTAC  >  NC_000913/2974252‑2974445

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