BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,736,520	A→G	Q218R (CAG→CGG)	purR →	PurR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,736,520	0	A	G	99.4%	527.1 / ‑4.8	172	Q218R (CAG→CGG)	purR	PurR
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (93/77); total (95/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

Only 100 of 185 total aligned reads displayed.

TGGAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGC  >  NC_000913/1736430‑1736608

tGGAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCggg                                                                                        <  2:1023552/93‑1 (MQ=255)

 ggAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGTTTGTGCgggg                                                                                       >  2:130897/1‑93 (MQ=255)

 ggAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCgggg                                                                                       <  2:1186968/93‑1 (MQ=255)

 ggAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCgggg                                                                                       <  1:3015281/93‑1 (MQ=255)

 ggAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCgggg                                                                                       <  1:189484/93‑1 (MQ=255)

 ggAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGGGCCGGAAAGCTGGATTGTGCgggg                                                                                       >  2:1100057/1‑93 (MQ=255)

  gAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCggggg                                                                                      >  2:5298394/1‑92 (MQ=255)

  gAACATAACACCGGCGCAGGCCGTCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGt                                                                                      <  2:6282716/93‑1 (MQ=255)

   aaCGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTg                                                                                     <  2:5544855/93‑1 (MQ=255)

   aaCGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTg                                                                                     <  2:5544857/93‑1 (MQ=255)

      gTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACt                                                                                  <  1:1762053/93‑1 (MQ=255)

       tAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACtt                                                                                 <  2:7145766/93‑1 (MQ=255)

        aaCACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACttt                                                                                >  2:3352981/1‑93 (MQ=255)

           acCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGaa                                                                             >  1:3272151/1‑93 (MQ=255)

           acCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGaa                                                                             <  2:3669792/93‑1 (MQ=255)

             cGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAAcc                                                                           <  1:4389857/93‑1 (MQ=255)

              ggCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCt                                                                          <  2:2507490/93‑1 (MQ=255)

              ggCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCt                                                                          >  1:435750/1‑93 (MQ=255)

               gcgcAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCCGGGTGACTTTGAACCTg                                                                         >  1:1457102/1‑93 (MQ=255)

                cgcAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGa                                                                        <  2:6631138/93‑1 (MQ=255)

                 gcAGGCCTCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGaa                                                                       <  2:7842837/93‑1 (MQ=255)

                 gcAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGaa                                                                       <  2:3646370/93‑1 (MQ=255)

                 gcAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGaa                                                                       >  2:2403138/1‑93 (MQ=255)

                 gcAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGaa                                                                       >  2:2520701/1‑93 (MQ=255)

                 gcAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGaa                                                                       >  1:6654840/1‑93 (MQ=255)

                 gcAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGaa                                                                       >  2:2131375/1‑93 (MQ=255)

                 gcAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGaa                                                                       <  1:2829760/93‑1 (MQ=255)

                   aGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATc                                                                     <  1:3534291/93‑1 (MQ=255)

                   aGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATc                                                                     <  1:1395707/93‑1 (MQ=255)

                    ggCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATcc                                                                    <  2:1403437/93‑1 (MQ=255)

                     gccgccTGGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCg                                                                   <  2:4297402/93‑1 (MQ=255)

                       cgccTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGt                                                                 >  2:6692976/1‑93 (MQ=255)

                       cgccTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGGGACTTTGACCCTCACTCCGGt                                                                 >  1:6140572/1‑93 (MQ=255)

                        gccTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGtt                                                                >  2:4019514/1‑93 (MQ=255)

                        gccTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGtt                                                                >  1:3271251/1‑93 (MQ=255)

                         ccTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTa                                                               >  1:306204/1‑93 (MQ=255)

                          cTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTAt                                                              >  1:5442738/1‑93 (MQ=255)

                          cTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTAt                                                              >  2:7664082/1‑93 (MQ=255)

                          cTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTAt                                                              >  1:7599259/1‑93 (MQ=255)

                           ttGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATc                                                             >  1:7368887/1‑93 (MQ=255)

                            tGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATcg                                                            >  1:3349037/1‑93 (MQ=255)

                            tGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATcg                                                            <  1:846982/93‑1 (MQ=255)

                                ggTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGcc                                                        >  2:5794908/1‑93 (MQ=255)

                                 gTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGcct                                                       >  1:1586348/1‑92 (MQ=255)

                                 gTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCa                                                       <  2:1768599/93‑1 (MQ=255)

                                 gTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCa                                                       <  1:2469821/93‑1 (MQ=255)

                                  ttttATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCAt                                                      <  2:738056/93‑1 (MQ=255)

                                  ttttATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACATGAATCCGGTTATCGCGCCAt                                                      <  1:1373169/93‑1 (MQ=255)

                                   tttATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATg                                                     >  1:2930964/1‑93 (MQ=255)

                                    ttATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATgc                                                    >  2:3298404/1‑93 (MQ=255)

                                     tATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATgca                                                   >  2:272596/1‑93 (MQ=255)

                                       tGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATgcagc                                                 <  2:3486823/93‑1 (MQ=255)

                                        gAAGGCGATGGAAGAAGCGAAGATCAAGGGGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTACCGCGCCATggagca                                                >  2:6771811/1‑93 (MQ=255)

                                         aaGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCaa                                               <  2:2195227/93‑1 (MQ=255)

                                          aGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCaaa                                              >  2:2754631/1‑93 (MQ=255)

                                          aGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCaaa                                              >  1:6899537/1‑93 (MQ=255)

                                          aGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCaaa                                              <  2:3690007/93‑1 (MQ=255)

                                          aGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCaaa                                              >  2:156439/1‑93 (MQ=255)

                                            gCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATc                                            >  1:5995927/1‑93 (MQ=255)

                                               aTGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTg                                         >  2:3150172/1‑93 (MQ=255)

                                                tGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGt                                        <  2:5019258/93‑1 (MQ=255)

                                                tGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGt                                        <  2:1263622/93‑1 (MQ=255)

                                                tGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGt                                        >  2:4972429/1‑93 (MQ=255)

                                                tGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGt                                        >  2:921009/1‑93 (MQ=255)

                                                 ggAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTc                                       >  1:1142463/1‑93 (MQ=255)

                                                  gaagaaGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCg                                      >  1:3986470/1‑93 (MQ=255)

                                                    agaagCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCa                                    >  1:5270351/1‑93 (MQ=255)

                                                    agaagCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCa                                    <  1:258441/93‑1 (MQ=255)

                                                     gaagCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAg                                   >  2:2686462/1‑93 (MQ=255)

                                                     gaagCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAg                                   <  2:549785/93‑1 (MQ=255)

                                                     gaagCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAg                                   <  1:7384897/93‑1 (MQ=255)

                                                       agCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGGGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGcc                                 >  1:6078791/1‑93 (MQ=255)

                                                          gatgatCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCa                              >  1:3719452/1‑93 (MQ=255)

                                                          gatgatCAAGGTGCCGGAAAGCTGGATTGTGCGGGGGGACTTTGAACCTGAATCCGGCTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCa                              >  2:4974090/1‑93 (MQ=255)

                                                            tgatCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATc                            >  1:1833558/1‑93 (MQ=255)

                                                            tgatCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATc                            >  1:640502/1‑93 (MQ=255)

                                                            tgatCAAGGTGCCGGAAAGCTGGATTGTGCGGGGGGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATc                            >  1:5945533/1‑93 (MQ=255)

                                                             gatCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCg                           >  1:6193558/1‑93 (MQ=255)

                                                               tCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGcc                         >  2:4131307/1‑93 (MQ=255)

                                                                cAAGGTGCCGGAAAGCTGGATTGTGCGGGGTTACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGccc                        <  1:5560328/93‑1 (MQ=255)

                                                                 aaGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCt                       >  2:4584031/1‑93 (MQ=255)

                                                                  aGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTa                      <  1:2915387/93‑1 (MQ=255)

                                                                   ggTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTAc                     <  1:2287289/93‑1 (MQ=255)

                                                                   ggTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTAc                     >  1:4696817/1‑93 (MQ=255)

                                                                   ggTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTAc                     <  2:5162776/93‑1 (MQ=255)

                                                                   ggTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTAc                     <  1:1015061/93‑1 (MQ=255)

                                                                    gTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACt                    <  1:7719675/93‑1 (MQ=255)

                                                                    gTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACt                    <  1:6778553/93‑1 (MQ=255)

                                                                        cGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCg                >  1:7007548/1‑93 (MQ=255)

                                                                         ggAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTCCCgc               >  1:5611846/1‑92 (MQ=255)

                                                                          gAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGtc              <  2:67925/93‑1 (MQ=255)

                                                                          gAAAGCTGGATTGTGCGGGGGGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGtc              >  1:649317/1‑93 (MQ=255)

                                                                             aGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGtcttc           >  1:2988568/1‑93 (MQ=255)

                                                                              gCTGGATTGTGCGGTGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGtcttct          <  2:461756/93‑1 (MQ=255)

                                                                                tGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCtgt        >  1:2802852/1‑93 (MQ=255)

                                                                                 ggATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCtgtg       >  1:5957981/1‑93 (MQ=255)

                                                                                 ggATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCtgtg       <  1:7562208/93‑1 (MQ=255)

                                                                                  gATTGTGCGGGGTGACTTTGAACTTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTgtgg      <  2:5653505/93‑1 (MQ=255)

                                                                                  gATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTgtgg      >  1:6487611/1‑93 (MQ=255)

                                                                                   aTTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAACAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTgtggt     <  1:2126912/93‑1 (MQ=255)

TGGAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGC  >  NC_000913/1736430‑1736608

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: `‑`